search for: calpha

Hej, Calculation of confidence intervals for means based on a model fitted with lmer using the package multcomp - yields results for calpha=adjusted_calpha - NA's for calpha=univariate_calpha Example: library(lme4) library(multcomp) ### Generate data set.seed(8) d<-expand.grid(treat=1:2,block=1:3) e<-rnorm(3) names(e)<-1:3 d$y<-rnorm(nrow(d)) + e[d$block] d<-transform(d,treat=factor(treat),block=factor...

...arded message ---------- > Date: Sat, 3 Sep 2011 23:56:35 +0200 > From: Ulrich Halekoh <Ulrich.Halekoh at agrsci.dk> > To: "r-help at r-project.org" <r-help at r-project.org> > Subject: [R] glht (multcomp): NA's for confidence intervals using > univariate_calpha > > Hej, > > Calculation of confidence intervals for means > based on a model fitted with lmer > > using the package multcomp > > - yields results for calpha=adjusted_calpha > - NA's for calpha=univariate_calpha > > > Example: > library(lm...

...vec[1]-1)*varvec[1]+(nvec[2]-1)*varvec[2])/error.df diff12estimate=ybar[1]-ybar[2] stderr=sqrt(pooled.var*((1/nvec[1])+(1/nvec[2]))) tratio=diff12estimate/stderr twosidedP=2*(1-pt(abs(tratio),error.df)) tcrit=qt(1-alpha/2,error.df) lower=diff12estimate-tcrit*stderr upper=diff12estimate+tcrit*stderr calpha=1-alpha out=data.frame(diff12estimate,stderr,tratio,twosidedP,lower,upper,alpha) names(out)=c("Estimator","SE","T","P-value","Lower CI","Upper CI","Confidence") out }

Dear R-users, I'm currently using the multcomp package to produce plots of means with 95% confidence intervals i.e. mult<-glht(lm(response~treatment, data=statdata), linfct=mcp(treatment="Means")) plot(confint(mult,calpha = sig)) Unfortunately the y-axis on the plot appears to be fixed and hence if the labels on the y-axis (treatment levels) are too long, then they are not displayed in full on the plot. Of course I could always make the labels shorter but I was wondering if there was a way to make the position of t...

...vec[1]-1)*varvec[1]+(nvec[2]-1)*varvec[2])/error.df diff12estimate=ybar[1]-ybar[2] stderr=sqrt(pooled.var*((1/nvec[1])+(1/nvec[2]))) tratio=diff12estimate/stderr twosidedP=2*(1-pt(abs(tratio),error.df)) tcrit=qt(1-alpha/2,error.df) lower=diff12estimate-tcrit*stderr upper=diff12estimate+tcrit*stderr calpha=1-alpha out=data.frame(diff12estimate,stderr,tratio,twosidedP,lower,upper,alpha) names(out)=c("Estimator","SE","T","P-value","Lower CI","Upper CI","Confidence") out } I need to alter the above function twosamp so that its outpu...

...-exclude .svn --exclude "*~" --delete --delete-excluded CONTRIBUTORS COPYRIGHT Makefile arrangements src lib .gitignore AEIArrangements AEIPhysics bbhfactory bin cactus.config carpet carpet-stable carpet-stable-2 carpet-stable-3 kranc parfiles eschnett@zeke.loni.org:/work/default/eschnett/Calpha After checking about 8000 files, I receive the error sending incremental file list rsync: connection unexpectedly closed (794 bytes received so far) [sender] rsync error: error allocating core memory buffers (code 22) at io.c(635) [sender=3.0.2] With two additional --verbose options, the last sc...