similar to: Frailty models and omnibus test

Dear all, I'm quite new in R and especially with linear mixed effects models and I'm not completely sure to read the lmer table in the right way. for example: head(march.f) fam subjID Cond Code reg total first second log.total log.second cat 3 f 30 an fDan1 3 1.2304688 0.6679688 0.56250000 0.20739519 0.44628710 f

Relatively new to R, I'm trying to do a relatively simple task. I have data set that has several variables arranged by SubjID and visit, with multiple observations for that combination. I do linear regression on those multiple observations, then generated a set of interpolated values from the regression at fixed intervals along "x". I now want to average each of those across all the

Hi, I have a multilevel situation where subjects are nested within clinics, and each subject has multiple measurements. For simplicity, suppose there 4 clinics, 3 subjects per clinic, and each subject has 3 repeated measures. Outcome is continuous. I am trying to implement this model with lmer function in lme4 library.

Dear all, I'm quite new in R and I have a problem with the function which. When I use it to select a subset of a dataframe it works well but somewhere R takes trace of the past dataframe and this creates problems with following operations. For example: sentences <- read.xls("frasi.tot.march.3.xls", header=TRUE) head(sentences) fam subjID Cond Code reg total first

Dear All, I met problems when doing contrast and now really need some help in the model below: Fit=gam(y~treat+SEQUENCE+PERIOD+SEX+s(x),data=dat, random=list(SUBJID=~1),correlation=corAR1(form=~1|SUBJID)) And error message keeps coming out when I want to compare the differences between treatments: Diff=contrast(Fit, list(treat=treatment[-placebo.pos]),list(treat="Placebo"),

gem install camping-omnibus --source code.whytheluckystiff.net Considering the number of optional dependencies now in Camping, I''ve put together an eigengem which will load all of the recommended libraries that most devoted campers will want: ActiveRecord, Mongrel, SQLite3, RedCloth and the acts_as_versioned gem. (Those last two are used in the examples.) Stuff like mosquito and

I am looking for an R package with which one can calculate an effect size for a set of contrasts given an omnibus chi-square test statistic (more than 1 degree of freedom). Is there such a package? Presumably, it would implement the procedure (or something like it) described by Roznow and Rosenthal 1996 Psychological Methods 1: 331-340. many thanks! Steven Orzack Fresh Pond Research Institute

If one estimates a model using multinom, is it possible to perform the omnibus LR test ( the analogue to omnibus F in linear models ) using the output from multinom ? The residual deviance is there but I was hoping I could somehow pull out the deviance based on just using an intercept ? Sample code is below from the CAR book but I wasn't sure how to do it based on that example. Thanks

I am trying to optimize some code to take advantage of R loop-avoiding capabilities when working on vectors/arrays that contain time intervals. The calculation involves adding (for each time interval) the time portion (of events defined by their start and end times) that elapsed during time intervals. Any advice on how to improve this code. I searched the email archive and looked at the MASS

I''ve noticed that the copy of Mongrel installed by the camping-omnibus gem doesn''t work with Ruby 1.8.6. Or to be more specific, cgi_multipart_eof_fix (which Mongrel is dependent upon) doesn''t work: > $ sudo gem install mongrel --source http://code.whytheluckystiff.net > > ERROR: Error installing mongrel: > cgi_multipart_eof_fix requires Ruby version

R 2.9.1 Windows XP I am fitting a random effects ANOVA with two factors Group which has two levels and Time which has three levels: fita<-lme(Post~Time+factor(Group)+factor(Group)*Time, random=~1|SS,data=blah$alldata) I want to get the omnibus significance tests for each factor and the interaction. I believe I can get the omnibus test for the interaction by running the model:

Hello, I've been running some HLMs using the lme function quite happily; it does what I want and I'm pretty sure I understand it. The issue is that I'm currently trying to estimate a model with a 14-level "nusiance" factor as an independent variable...which makes the output quite ugly. All I'm really interested in is the question of whether these factor as a whole

Dear Prof Azzalini, You have an interesting example of recursive use of uniroot(). [re-cited at the end] However, note that R currently has the same problem as S-plus: Uniroot() doesn't work reliably, recursively. When you found it to be better, you were just lucky. The relevant file in R's source, src/main/optimize.c says /* WARNING : As things stand, these routines should

Daniel, When Group is entered as a factor, and the factor has two levels, the ANOVA table gives a p value for each level of the factor. What I am looking for is the omnibus p value for the factor, i.e. the test that the factor (with all its levels) improves the prediction of the outcome. You are correct that normally one could rely on the fact that the model

Hi guys, I am trying to compile a 64bit version of R 2.7.0 on AIX 5.3. And I am running into some troubles. I also try the latest version 2.7.1 and get the very similar errors. The machine info: ----------------------- $uname -srvp AIX 3 5 powerpc ---------------------- My configuration is this: ---------------------------------------------------------- ./configure \

I am working on a nlme model that has multiple fixed effects (linear and nonlinear) with a nonlinear (asymptotic) random effect. asymporig<-function(x,th1,th2)th1*(1-exp(-exp(th2)*x)) asymporigb<-function(x,th1b,th2b)th1b*(1-exp(-exp(th2b)*x)) mod.vol.nlme<-nlme(fa20~(ah*habdiv+ads*ds+ads2*ds2+at*trout)+asymporig(da.p,th1,th2)+ asymporigb(vol,th1b,th2b),

Hello, This morning I downloaded, unzipped, and compiled the latest version of R for unix, in my HOME/APPLICATIONS directory. My current unix machine is a sparc-sun-solaris2.9, running SunOS 5.9. Here are the last few lines of the output following ./configure: R is now configured for sparc-sun-solaris2.9 Source directory: . Installation directory: /usr/local C compiler:

Hi R users, I want to build R-2.8.1 on AIX5.3, but I got following error message: Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared library '/rnd/homes/jixu/tmp/R-2.8.1/library/stats/libs/stats.so': rtld: 0712-001 Symbol d1mach was referenced from module /rnd/homes/jixu/tmp/R-2.8.1/library/stats/libs/stats.so(), but a runtime definition of

On 04-Jan-16 23:24, Michael Felt wrote: > The bulk is on my forums - the final post for today is: > > Results to date: > > A. It looks like I am going to need a newer compiler for C - xlc/xlC > V11 apparently does not understand this code: > > "/data/prj/cran/R-3.2.3/src/main/memory.c", line 2149.31: 1506-046 (S) > Syntax error. > > I will have to check

I cannot seem to get a L'abbe plot to work on R. I do not understand what the X coordinates, or alternatively an object of class metabin, is supposed to mean. What is a class of metabin? Institute of Behavioral Genetics University of Colorado, Boulder Whitney.Melroy at Colorado.EDU