Ana Marija
2019-Oct-23 01:17 UTC
[R] how to compare two distributions and calculate p value?
Hello, I would like to calculate a p value from two distributions, one looks like this:> head(b)gene_id number_of_eqtles_per_gene 1: ENSG00000237683.5 5 2: ENSG00000225972.1 267 3: ENSG00000225630.1 97 4: ENSG00000237973.1 257 5: ENSG00000240409.1 19 6: ENSG00000248527.1 41> dim(b)[1] 31385 2 and the other one (call it "a") looks very similar to "b" only has 103 entries. How would I do say t test considering these two have different number of entries? t.test(a,b...) Thanks Ana
Bert Gunter
2019-Oct-23 01:34 UTC
[R] how to compare two distributions and calculate p value?
"How would I do say t test considering these two have different number of entries?" Read and follow ?t.test . Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Oct 22, 2019 at 6:18 PM Ana Marija <sokovic.anamarija at gmail.com> wrote:> Hello, > > I would like to calculate a p value from two distributions, one looks like > this: > > head(b) > gene_id number_of_eqtles_per_gene > 1: ENSG00000237683.5 5 > 2: ENSG00000225972.1 267 > 3: ENSG00000225630.1 97 > 4: ENSG00000237973.1 257 > 5: ENSG00000240409.1 19 > 6: ENSG00000248527.1 41 > > dim(b) > [1] 31385 2 > and the other one (call it "a") looks very similar to "b" only has 103 > entries. > > How would I do say t test considering these two have different number > of entries? > t.test(a,b...) > > Thanks > Ana > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
Ana Marija
2019-Oct-23 02:24 UTC
[R] how to compare two distributions and calculate p value?
ok, I would have paired = FALSE option. What other options you would recommend? On Tue, Oct 22, 2019 at 8:34 PM Bert Gunter <bgunter.4567 at gmail.com> wrote:> > "How would I do say t test considering these two have different number > of entries?" > > Read and follow ?t.test . > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along and sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > > On Tue, Oct 22, 2019 at 6:18 PM Ana Marija <sokovic.anamarija at gmail.com> wrote: >> >> Hello, >> >> I would like to calculate a p value from two distributions, one looks like this: >> > head(b) >> gene_id number_of_eqtles_per_gene >> 1: ENSG00000237683.5 5 >> 2: ENSG00000225972.1 267 >> 3: ENSG00000225630.1 97 >> 4: ENSG00000237973.1 257 >> 5: ENSG00000240409.1 19 >> 6: ENSG00000248527.1 41 >> > dim(b) >> [1] 31385 2 >> and the other one (call it "a") looks very similar to "b" only has 103 entries. >> >> How would I do say t test considering these two have different number >> of entries? >> t.test(a,b...) >> >> Thanks >> Ana >> >> ______________________________________________ >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code.