I am curious if Seqmonk can give me coverage of specific genes/ transcripts from DNA seq data as is given by coverageBed in bedtools http://bedtools.readthedocs.org/en/latest/
I want to calculate coverage per gene/ chromosomes
Thanks
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by Aspadia View PostI want to use the 'relative quantitation' option, I create probes in an IP and input sample from ChIP then apply the running window probe generator and would then like to use the relative quantitation. However, after the probe generation it does not show me the option for relative quantitation nor when I go to Data --> quantitate existing probes....
Can you tell me how I can get to this relative quantitation?
Leave a comment:
-
Relative quantitation
Hi Simon,
Thanks for your previous quick reply, I got seqmonk running now
This brings a next problem. I want to use the 'relative quantitation' option, I create probes in an IP and input sample from ChIP then apply the running window probe generator and would then like to use the relative quantitation. However, after the probe generation it does not show me the option for relative quantitation nor when I go to Data --> quantitate existing probes....
Can you tell me how I can get to this relative quantitation?
Thanks a lot!
Leave a comment:
-
Originally posted by markf View PostDear Simon,
I have a question about seqmonk: Is it possible to create my own probesets & quantifications outside of seqmonk, and then load them? For example from GFF, and assign the GFF score as the quantification score for that probe.
What you would then do is to import the data for that ChIP and then use the standard quantitation tools to assign quantitated values to the peaks you imported. What you can't do is to import pre-quantitated data - we have thought about this a lot since we've had a fair few requests for it, but it causes so many problems in other parts of the data model that we've decided against this. There are other programs around which deal solely with pre-quantitated data (IGV etc) and the point about SeqMonk is that it has access to the raw data so you can do things which those programs can't do.
Originally posted by markf View PostPS - keep up the good work - seqmonk is great!
Leave a comment:
-
Data import
Dear Simon,
I have a question about seqmonk: Is it possible to create my own probesets & quantifications outside of seqmonk, and then load them? For example from GFF, and assign the GFF score as the quantification score for that probe.
I have my own set of peak calls (from a 4C experiment, with my own preprocessing, and using R3Cseq) - and I was hoping to use seqmonk for some sanity checking, visualization and report generation.
Thanks.
Mark
PS - keep up the good work - seqmonk is great!
Leave a comment:
-
If java is not found on your command line the most likely cause we've seen is that you have installed the 32 bit version of java on a 64 bit machine. Oracle only show the 32 bit version of java on the front page of java.com, to get the 64 bit version you need to click on "See all java downloads" and then get the 64 bit version.
You can install 32 and 64 bit side by side so adding the 64 bit version won't mess up your browser, and you'll need the 64 bit version to be able to use more than 2GB memory.
Hope this helps
Leave a comment:
-
Seqmonk and Java
Simon,
The Java and seqmonk problem is now creating a real headache. I am also in same boat. Perhaps at some point in time Seqmonk and Java dependency need to be thoroughly investigated. I agree it may be because of Java updates but we want Seqmonk to work.
Thanks
Leave a comment:
-
Hi Simon,
Seqmonk sounds really awesome but I do not manage to run it
It is downloaded and when I try to run it it says it cannot find java. When I type java -version in cmd it says 'java' is not recognized as an internal or external command, operable program or batch file.
In previous reply you mentioned this:
Java isn't installed properly, or the java binary isn't in your path. Open a command prompt and type 'java -version' if you get an error saying this isn't a recognised command then this is the problem.
So I guess the second part is my problem but how can I solve this? I just installed the latest version of Java (Version 7 Update 25).
Hope you can help me, thanks in advance!
Leave a comment:
-
Originally posted by crazyhottommy View PostI did not quite get it.
I think you mean to say:
10, 12, 8, 10 are values for 4 probes in the same sample ( or in my case TF1)
20,22,18,20 are values for the same 4 probes in my other sample ( the other TF2)
10,12,8,10 are the values for 1 probe across 4 different samples.
20,22,18,20 are the values for a different probe against the same samples.
Per-probe normalisation is a way to be able to compare the relative changes in each probe and removes the effect of the absolute values it has. In this case you'd be interested that sample 2 had the highest value and 3 the lowest even though the absolute values of the probes are quite different.
Leave a comment:
-
Originally posted by simonandrews View PostPer probe normalisation simply adjusts your values by subtracting the median values for that probe across all of your samples from each of the individual probe values. In effect what it shows you is whether the value you have for that probe is higher or lower in this sample than the other samples you're looking at.
Taking some real numbers if you have a probe which has values of 10,12,8 and 10 across 4 samples and you did per-probe normalisation you'd find the median was 10, so the normalised values would be 0,2,-2 and 0. This makes it a lot easier to compare to another data set which might have values of 20,22,18 and 20 for example.
I think you mean to say:
10, 12, 8, 10 are values for 4 probes in the same sample ( or in my case TF1)
20,22,18,20 are values for the same 4 probes in my other sample ( the other TF2)
per-probe normalization is done to compare the value for TF1 and TF2 for the same probe.
I am a little bit confused...
Leave a comment:
-
Originally posted by crazyhottommy View PostWhat does per-probe normalization mean?
Taking some real numbers if you have a probe which has values of 10,12,8 and 10 across 4 samples and you did per-probe normalisation you'd find the median was 10, so the normalised values would be 0,2,-2 and 0. This makes it a lot easier to compare to another data set which might have values of 20,22,18 and 20 for example.
Leave a comment:
-
Originally posted by simonandrews View PostIt's fine to use either value for subsequent analysis. The largest data store option is there so you don't get fractional (or negative if you're on a log scale) read counts, but there's nothing intrinsically wrong with these, it's just that they can look a bit odd on any plot where 0 is supposed to mean something. The magnitude of differences will be exactly the same in either case.
When you're doing a hierarchical cluster you normally do per-probe normalisation anyway, in which case the plot you get will be exactly the same whichever of the earlier normalisation options you selected.
Leave a comment:
-
Originally posted by crazyhottommy View PostWow, prompt answer. Thank you.
another question would be the normalization for the ChIP-seq counts. Seqmonk provides two ways to normalize them.
1. log2 CPM
2. largest data set
my understanding is that normalizing to the largest data set is essentially the same as
log2 CPM but scale up with some factor.
the problem with log2 CPM normalization is that you got some negative values.
Does it still make sense if I use result from the log2 CPM to do the downstream analysis like scatter plot and hierachical cluster? Thanks!
When you're doing a hierarchical cluster you normally do per-probe normalisation anyway, in which case the plot you get will be exactly the same whichever of the earlier normalisation options you selected.
Leave a comment:
-
Originally posted by simonandrews View PostIf you give the two sets of results the same feature name ('peak' for example) then they'll be merged together in the same track and you can use the feature probe generator to make probes from the combined set. You can always separate them later by renaming one of the tracks using the controls in the Annotation Sets folder of the data view.
another question would be the normalization for the ChIP-seq counts. Seqmonk provides two ways to normalize them.
1. log2 CPM
2. largest data set
my understanding is that normalizing to the largest data set is essentially the same as
log2 CPM but scale up with some factor.
the problem with log2 CPM normalization is that you got some negative values.
Does it still make sense if I use result from the log2 CPM to do the downstream analysis like scatter plot and hierachical cluster? Thanks!
Leave a comment:
-
Originally posted by crazyhottommy View PostHi Simon,
Sorry to bother you again. I have two ChIP-seq data, I do a MACS peak call for one, and got some peaks as probes, I then import them as annotations to keep them.
I then did the same thing for my second data set, and also import the probes as the annotations.
How can I pool these two probes (union of these two)? "the combining existing list" only works on active probes, I did not find any function in define probes tab either.
Leave a comment:
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
||
Started by seqadmin, 05-02-2024, 08:06 AM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
05-02-2024, 08:06 AM
|
||
Started by seqadmin, 04-30-2024, 12:17 PM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-30-2024, 12:17 PM
|
Leave a comment: