similar to: hierarchical linear models, mixed models and lme

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

Hi, I''m trying to understand the lme output and procedure. I''m using the Crawley''s book. I''m try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :125.0 Min. :1 Min. :1.0 Min. :1 1st Qu.:135.8

Hello R friends, I'm wondering why I get funny qqnorm() results. It seems that they should all be reflected in the normal qqline(). For instance: if I qqnorm() bimodal or uniform data I get a sigmoidal in which the qqnorm() points lie above the qqline() at -ve theoretical quantiles, and the qqnorm() points lie below the qqline() at +ve theoretical quantiles. Yet I expect such platykurtic

Hi, I'm trying to figure out how anova works in R by translating the examples in Sokal And Rohlf's (1995 3rd edition) Biometry. I've hit a snag with planned comparisons, their box 9.4 and section 9.6. It's a basic anova design: treatment <- factor(rep(c("control", "glucose", "fructose", "gluc+fruct",

This question is mainly aimed at Kurt Hornik as author of the ctest package, but I'm cc'ing it to r-help as I suspect there will be other valuable opinions out there. I have been attempting 2 sample Kolmogorov-Smirnov tests using the ks.test function from the ctest package (ctest v.0.9-15, R v.0.63.3 win32). I am comparing fish length-frequency distributions. My main reference for the

Hola, tengo una serie de datos datExpr, al usar cor() : cor(datExpr ,method = "pearson", use ="pairwise.complete.obs") me da el siguiente error allocMatrix: too many elements specified Trate con "complete.obs", "na.or.complete", y el resto de las opciones para "use", pero siempre me da algun error.  ¿Alguna idea de como puedo hacer que cor() lea

Hi All, I would like to ask a question on fixed and random effecti in lme. I am fiddlying around Mick Crawley dataset "rats" : http://www.bio.ic.ac.uk/research/mjcraw/statcomp/data/ The advantage is that most work is already done in Crawley's book (page 361 onwards) so I can check what I am doing. I am tryg to reproduce the nested analysis on page 368:

Dear r-helpers, Prior to analysis of variance, I ran the Boxcox function (MASS library) to find the best power transformation of my data. However, reading the Boxcox help file, I cannot figure out if this function (through its associated log-likelihood function) corrects for * normality only * or if it also induces * homogeneity of variances *. I found in Biometry (Sokal and Rohlf, p. 419)

Ronaldo, Further to my previous posting on your Glycogen nested aov model. Having read Douglas Bates' response and Reflected on his lmer analysis output of your aov nested model example as given.The Glycogen treatment has to be a Fixed Effect.If a 'treatment' isn't a Fixed Effect what is ? If Douglas Bates' lmer model is modified to treat Glycogen Treatment as a purely

I frequently want to test for differences between animal size frequency distributions. The obvious test (I think) to use is the Kolmogorov-Smirnov two sample test (provided in R as the function ks.test in package ctest). The KS test is for continuous variables and this obviously includes length, weight etc. However, limitations in measuring (e.g length to the nearest cm/mm, weight to the nearest

Hello list, I'm trying to figure out how exactly the specification of nested random effects works in the lmer function of lme4. To give a concrete example, consider the rat-liver dataset from the R book (rats.txt from: http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/ ). Crawley suggests to analyze this data in the following way: library(lme4) attach(rats) Treatment <-

Dear list, Could someone show me where can I find the detailed formula on how to calculate the two way anova with unbalanced design? Say, if I have 2*2 design with 10,20,30,40 samples in each of the 2*2 cells. Most of the places I've googled only show how to calculate using software such as R, but not clear the detailed formula for calculating this. Thanks, Jack [[alternative HTML version

Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the same output in all 3 cases below (forgetting about the convergence problem). Would there be a reason to prefer labeling the levels one way or

Hello all. This may be a trivially simple question to answer, but I'm a little bit stumped with respect to the calculation of the F statistics in nested anovas in R. If I understand correctly, the F statistic for the among-subgroups but within groups hypothesis is calculated as MS_subgroups/MS_error, while the F statistic for the factor is calculated as MS_factor/MS_subgroups (I'm

Dear professors and collegues I need to perform a analysis of dates from a nested experimental design. From "Bioestatical Analysis" of Zar "Bimetry of Sokal" & Rohlf "Design and Analysis of Experiments" of Montgomery I have: Sum (mean(x)_i - mean(x)_T)2 / (a-1) -> var(epsilon) + n sigma2_B + n b (sum alfa_i)2 / (a-1) Sum (mean(x)_ij - mean(x)_i)2 /

A common point made in discussion of contrasts, type I, II, III SS etc is that for sensible comparisons one should use contrasts that are 'orthogonal in the row-basis of the model matrix' (to quote from http://finzi.psych.upenn.edu/R/Rhelp02/archive/111550.html) Question: How would one check, in R, that this is so for a particular fitted linear model object? Steve Ellison

Dear R-listers, Does anybody know what is the correct source of "rats" dataset in survival package? The help gives the following information: Rat data from survival5 Description: 48 rats were injected with a carcinogen, and then randomized to either drug or placebo. The number of tumors ranges from 0 to 13; all rats were censored at 6 months after randomization.

Hi, I am trying to run the command "forward.sel.par," however I receive the error message: "Error: could not find function 'simpleRDA2'." I have the vegan library loaded. The documentation on "varpart" has not helped me to understand why I cannot call this function. Maybe I am missing something obvious because I am still an 'R' novice. Below is a

Hi, I am new to R and AIC scores but what I get from coxme seems wrong. The AIC score increases as p-values decrease. Since lower AIC scores mean better models and lower p-values mean stronger effects or differences then shouldn't they change in the same direction? I found this happens with the data set rats as well as my own data. Below is the output for two models constructed with the rats

I am using the following: library(RODBC) chan = odbcConnectExcel("rats-lda") rats.lda = sqlFetch(chan, "data") close(chan) And getting the following error message: > library(RODBC) Error in library(RODBC) : there is no package called ?RODBC? > chan = odbcConnectExcel("rats-lda") Error: could not find function "odbcConnectExcel" > rats.lda =