similar to: Meta-analysis question

Hi all: I have a question about lm on t-test. data(sleep) I wanna perform t-test to test the difference between the 2 groups: I can use: t.test(extra~group) The t.test result shows that:t = -1.8608; mean1=0.75,mean2=2.33 But I still wanna use: summary(lm(extra~group)) Intercept=0.75,which is mean1,just the same as t.test. group2=1.58 means the difference of the 2 groups,so

Dear all! After hours of trying around, I gave up: I have a 2-dimensional array, and I know how to split it into its rows and how to get the mean for every row using 'sapply'. But what I want is to calculate the mean over the first n rows, and then the second n rows, etc., so that I get a vector like: v == mean1(row 1:5), mean2(row6:10),... (trivial, you might say. I find it rather

I want to obtain a barplot with error bars and a mean plot with error bars with other scale on the left in the same plot. I need help to get the same x axis (centered) when overlay two plots (barplot2 and plotCI (or lineplot.CI) with errorbars), with diferent y axis. I use par (new=T) but the X axis names are not centered with the figures. An example of my data:

I have a rather basic question. How can I get R to generate a ANOVA table and a f-value for a hypothesis test such as: Data: group1 values: 5.2 4.5 6.0 6.1 6.7 5.8 group2 values: 6.5 8.0 6.1 7.5 5.9 5.6 ... H0: mean1 = mean2 = mean3 = mean4 HA: at least two means different where I want to evaluate using a f test statistic? F = MSTr/MSE I'd like a table similar to one that

Dear All I want to generate variable with Bivariate Normal Distribution by use mean1 = a, variance1 = b, mean2 = c, variance2 = d, rho = e. How I can do this. Many Thanks. IRD [[alternative HTML version deleted]]

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100404/6f9d26c0/attachment.pl>

Hi, I am trying to write a function; I want to subtract the mean of each class in level 2 from the mean of each class in level 1 and square the answer, eg..... level.1 level.2 observation 1 1 0.5 1 1 0.2 1 2 0.6 1 2 0.4 2 3

Hi, I am a Ph.D. student from Québec, Canada. I’m a beginner with R and Bioconductor. Until now the only experience I have is in analyzing microarray data using affy and limma packages. Now I am trying to analyze Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth moderated t test on those arrays. Since no cdf official package is available for those arrays, after reading many

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20080422/9bb74193/attachment.pl

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100501/c1302c2e/attachment.pl>

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20111202/b4d64bba/attachment.pl>

Hello All, I posted a similar question before, but the direction was driven to whether my case is suitable for a wilcoxon test. After research about the appropriateness, I am pretty sure that a wilcoxon test is the right tool for my case. But how to compute the power of the test is still an unanswered question bothering me. The basic stats of my two paired samples are: mean1 = 0.0032, sd1 =

Dear all, I found a discrepancy while performing a power calculation for a two sample t-test in R and S-PLUS, respectively. For given values of sample number (5 each), sd (0.2) , significance level (0.01), and a desired power (80%) I looked for the difference in means. These values differ: 0.5488882 in R and 0.4322771 in S-PLUS (see dump below). Did I overlook any detail or confuse some

Please help me in writing the R code for this problem. I've been solving this for 4 days. It was hard for me to solve it. It's a simulation problem in R. The problem is My true model is a normal mixture which is given as 0.5 N(-0.8,1) + 0.5 N(0.8,1). This model has two components. I will get a random sample of size 100 from this model. I will do this 300 times. That means, I will have

Dear all, I could not draw a forest plot for meta-analysis under ramdon models using the rmeta package. The rmeta has a default function for MH (fixed-effect) model. Has the rmeta package been updated for such a function? Or someone revised it and kept a private code? I would appreciate it if you could provide some information on this question. Thanks, Andrew This email is intended

CRAN (and crantastic) updates this week New packages ------------ Updated packages ---------------- New reviews ----------- This email provided as a service for the R community by http://crantastic.org. Like it? Hate it? Please let us know: cranatic at gmail.com.

Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed? > library(rmeta) > data(cochrane) > cochrane name ev.trt n.trt ev.ctrl n.ctrl 1 Auckland 36 532 60 538 2 Block 1 69 5

Dear Colleagues, does anybody know how to aggregate a data.frame using different functions for different columns? Sincerely ___________________ Markus Preisetanz Consultant Client Vela GmbH Albert-Roßhaupter-Str. 32 81369 München fon: +49 (0) 89 742 17-113 fax: +49 (0) 89 742 17-150 mailto:markus.preisetanz@clientvela.com

Dear R-list, I'd like to do a meta-analysis plot similar to install.packages('rmeta') require(rmeta) data(catheter) a <- meta.MH(n.trt, n.ctrl, col.trt, col.ctrl, data=catheter, names=Name, subset=c(13,6,5,3,7,12,4,11,1,8,10,2)) summary(a) plot(a) (see attached file) by using my own OR (Odds Ratio) and 95% Confidence Interval data set, which looks like

Hello, I am facing a problem with optimization in R from 2-3 weeks. I have some Gaussian mixtures parameters and I want to find the maximum in that *Parameters are in the form * mean1 mean2 mean3 sigma1 sigma2 sigma3 c1 c2 c3 506.8644 672.8448 829.902 61.02859 9.149168 74.84682 0.1241933 0.6329082 0.2428986 I have used optima and optimx to find the