Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • honey
    replied
    error in running Repitools

    I tried to run the package as mentioned in fig5 of http://nar.oxfordjournals.org/conten...ar.gkr416.full
    but is getting error-



    I have changed upstream downstream distance, frequency as well as smoothing and no smoothing.
    Thanks for your help

    Error in pam(all, n.clusters, cluster.only = TRUE, do.swap = FALSE) :
    NA/NaN/Inf in foreign function call (arg 4)
    In addition: Warning messages:
    1: In min(x) : no non-missing arguments to min; returning Inf
    2: In max(x) : no non-missing arguments to max; returning -Inf

    Leave a comment:


  • honey
    replied
    K means on Chipseq

    frozenlyse,
    Thanks

    Let me clarify what I have I have two bed files with following structure following the Vignette.
    Chr start stop strand Reads Per Million (RPM) Transcript ID Gene Symbol normal: reads _cancer: reads
    chr4 56953396 56946396 - 5.04223 NM_181806 AASDH 35 106
    chr4 462999 455999 - 4.20186 NR_002451 ABCA11P 26 96
    chr7 150560251 150553251 - 4.97758 NM_007189 ABCF2 43 115
    chr21 42513151 42589572 + 3.55542 NM_207174 ABCG1 4 57
    chr4 8216460 8209460 - 3.6847 NM_001130088 ABLIM2 12 51
    chr3 130076023 130083023 + 4.46043 NR_033426 ACAD9 34 96
    chr3 133518991 133569357 + 4.2665 NM_001134194 ACPP 44 123

    One is for cancer and other is for normal. The coordinates are different for cancer and normal cells (above is snapshot of cancer bed file only). The question is if I have to do K means to show that normal and cancer show different type of clusters. Should I be taking the normal and cancer in one data frame (like Vignette). Then clustering will be based on raw reads(which perhaps may not be ideal) or RPM ( I get only for one sample- normal or cancer; further it may be one gene may be enriched in cancer but not in normal and vice versa). I looked for data(Hi1amples) from 4.3 but could not fine find structure of H1samples.

    I understand it appear that I am asking something which has to readymade, but it is kind of frustrating that I have results in seqminer but no good fig and I have not got any feedback from there.

    Thanks
    Last edited by honey; 06-20-2012, 04:09 PM.

    Leave a comment:


  • frozenlyse
    replied
    This can be done using Repitools from bioconductor - see section 5.3 of the vignette for an example.

    Leave a comment:


  • honey
    replied
    k- means clustering of chipseq peak data wrt to TSS

    Part of the problem is I dont have wiggle file and has only bed file with peak score. Moreover, Cistrome will not allow me to get list of candidate peaks in each cluster.

    Leave a comment:


  • kadircaner
    replied
    Try heatmap module under cistrome.org/ap. You can provide your own bed file along with wig files that you want to display. Picture quality is better than seqminer...

    Leave a comment:


  • honey
    started a topic k- means clustering of chipseq peak data wrt to TSS

    k- means clustering of chipseq peak data wrt to TSS

    This question has been discussed number of times in this forum- I have a paculiar situation. I have chips seq data for normal and cancer cells I have used a tool to call peaks and now has bed files which shsiblow enrichment in normal as well as cancer cells. I want to perform K means clustering using TSS as refrence of these peaks. I have used Seqminer to perform clustering I do have results which specifically show that different clusters for normal and cancer cells. The problem is I cannot get a good quality picture. I have tried all differenrt (everything) possible combinations. but picyure come out be very fuzzy. I used 5KB+/- t0 500bp+/- various other adjustments, of no use.

    I want to generate fig5 of http://nar.oxfordjournals.org/conten...ar.gkr416.full
    or Fig 2E of http://www.ncbi.nlm.nih.gov/pubmed/18992931
    Fig1A of http://www.jbc.org/content/early/201...bc.M111.266254
    I am also open to use any other tool which can allow me to use bed files for doing K means clustering. I want to stick with my peak calling alogrithm. I undertsand MACS in Galaxy can give such clustering but I have to use MACs peak calling to use such that tool.

    Thanks for your help.
    Last edited by honey; 06-18-2012, 10:34 PM.

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
30 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 09:45 AM
0 responses
202 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 08:54 AM
0 responses
212 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-02-2024, 03:00 PM
0 responses
194 views
0 likes
Last Post seqadmin  
Working...
X