search for: zcor

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070905/6d1002c1/attachment.pl

...both sub-group and meta-regression. In subgroup-analyses, I have stratified the database to create a separate csv file just for European women from the original database and conducted the following: women_west<-read.csv("women_west.csv") print(women_west) dat<-escalc(measure="ZCOR",ri=Pearson,ni=N,data=women_west,append=TRUE) res<-rma(yi,vi,data=dat) is.factor(dat$year) forest(res,transf=transf.ztor) In meta-regression, I used the original database, but used categorical moderators for sex (=women), and ethnicity (=european) to find the effect specifically in Europea...

...ing for residual spatial correlation . Several questions: 1. Are GEE’s the right choice? 2. How can I specify a user-defined spatial (exponential) correlation structure in geeglm or geese? I did not really get this from the help-sites. I guess I need to set corstr = "userdefined" and zcor = Zcor, with Zcor being in some way the specified correlation matrix – but how? I earlier specified spatial exponential correlation structures in glmmPQL, but it seems that this does not work in the same way for geeglm or geese. Do I need to use the function genZcor? 3. Makes it sense to use conti...

...e of work, Wolfgang! Thanks heaps.) I specified the "rma.uni' function so that it looks like this: MAergebnis<-rma.uni(ri=PosOutc, ni=N, mods=NULL, data=dm, intercept=TRUE, slab=c(dm$Article, dm$StudyNo), subset=NULL, measure="ZCOR", add=1/2, to="only0", vtype="LS", method="DL", weighted=TRUE, level=95, digits=4, btt=NULL, tau2=NULL, knha=FALSE, control=list()) With PosOutc being the column containing the raw correlations and N being the column with the...

...d fine: gee1<-geese(Activity~Water, id=Individual, data=dataF, family=gaussian) but when I use other variables e.g gee2<-geese(Activity~Hgt+Size+Water, id=Individual, data=dataF, family=gaussian) I get the error message "Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor, corp, : nrow(zsca) and length(y) not match" which I am not able to understand. The same problem has been reported in the list in 2006, but I have not found any response to it. Any suggestion? Giacomo

...ilar error. Can anyone help explain the error? > cor<-c(-0.3018, 0.667, -3.8002, -0.607, -0.4885, -3.8002, -0.0701, 0.1348, -0.9505, -0.5709, -0.6127, -1.2419, -0.1511, -0.1054)> n<-c(3,4,3,3,3,3,16,36,30,9,3,3,30,4)> library(meta)> metacor(cor, n, data=NULL, subset=NULL, sm="ZCOR", level=0.95, level.comb=level, comb.fixed=TRUE, comb.random=TRUE, hakn=NULL, method.tau="DL", tau.preset=NULL, TE.tau=NULL, method.bias="linreg", title="title", complab="comparison", outclab="outcome") Error in data.frame(subset = NULL, comb.f...

Hi I have 2 data.frames each of the same number of rows (approximately 30000 or more entries). They also have the same number of columns, lets say 2. One column has the date, the other column has a double precision number. Let the column names be V1, V2. Now I want to calculate the correlation of the 2 sets of data, for the last 100 days for every day available in the data.frames. My code looks

..."exchangeable",data=dat,alphainit=0.05) However, I am making a mistake with: library(geepack) fit.geese <- geese(score ~ chem + time, id=id, family=gaussian,corstr="exch",data=dat) I obtain the following error: Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor, corp, : nrow(zsca) and length(y) not match Could someone tell me what I have done incorrectly. Thanks for your time, Juliet. Data Below: id treat time1 time2 time3 time4 chem1 chem2 chem3 chem4 time score chem 1 1 20 18 15 15 1000 1100 1200 1300 0 20 1000 1 1 20 18 15 15 1000 1100 1200 1300 2...