Search Archives

Thanks to the Davids! I will brush up on my reading then. > The best description of this that I have found is the article by > Heagerty and Zeger in Statistical Science (2000, I think). For those interested, this article is...

Dear R users, I have not found anything on this in the archives. Does anyone know whehther the parameter use= is not functioning in cor or enlighten me what it is supposed to do? My R version is "R version...

Hi, I have a question regarding the creation of new variables on the basis of existing ones in R. I have two variables containing information on parents' educational degree (e.g. 1 'high school degree', 2 'college degree', etc.). I...

Dear all, Thanks for the recent comments on my segmentation fault problem when using loess. I realized that I should have mentioned two things. Firstly, I am already working with R --vsize 100 and still get a crash at n...

See help file of glmmPQL (the first thing you should be doing): glmmPQL(fixed, random, family, data, correlation, weights, control, niter = 10, verbose = TRUE, ...) ^^^^^^^^^^ Arguments ... niter maximum number of iterations. ... Thus, you can increase the value of niter to say...

I don't think, this has been answered: > I'm trying to run a 3-way within-subject anova in lme with 3 > fixed factors (Trust, Sex, and Freq), but get stuck with handling > the random effects. As I want...

Hi everyone, I have a question with reference to "glm.cluster" from the package "miceadds" and hope that someone can help me. I am trying to calculate cluster-robust standard errors for a glm-model with multiply imputed datasets. Everything...

> Now I want to evaluate GroupCov as a covariate to Treat. I > can then start with either m1 or m2 as base, but what is most > correct when GroupCov has only one value for each Group? > > m3 <- lm(Yield ~ Treat...

> This graph -> > http://www.math.hope.edu/~tanis/dallas/images/disth36.gif > is an example I found at > http://www.math.hope.edu/~tanis/dallas/disth1.html > created by Maple. > > Does anybody know how to create something similar in R...

> > If I understand correctly, you want to include the interactions > > between the random and fixed terms? > > Yes that is exactly I wanted to model. > > > This is done like: > > > > model.lme <- lme(Beta ~ Trust*Sex*Freq, > > random = ~Trust*Sex*Freq|Subj...

Joris, Basically, the I() function, as it is used in a regression model, allows the user to perform arithmetic operations on a variable that would otherwise be interpreted by the formula. It is not trapped as a special function as...

> I need to do multiple comparisons following nlme analysis (Compare > the effects of different treatments on a response measured > repeatedly over time; > fixed = response ~ treat*time). If you have an interaction it does not really make sense to conduct a...

Dear all, is there a standard way to calculate a nested, unbalanced ANOVA with R? I am glad for pointers, tips, ideas. The problem that I am trying to evaluate is a comparison of (average) group sizes. The explanatory variables...

Hi Tine, > I have been trying to look for documentation on making a 3D > histogram in R (so a histogram on which there are several distributions > are plotted), but I cannot find anything on it. Is it possible to > construct...

There is something I'm really missing here. The function show is a standardGeneric function, so the correct way to write it as method like this: setMethod("show", "Person", function(object) { for an object of class "Person" for example. Dave...

Hi everyone, I am trying to perform a multiple imputation with mice on a dataset of about 13000 observations and 178 variables. I can start an "empty" imputation ("imp_start <- mice(data, maxit=0)"), but after a few minutes R...

> The experimental design was suppose to consist of 4 > treatments replicated 3 time, Source 1 and applied at 10 cm > and source 2 applied at 20 cm. During the construction of the > treatmetns the depths vary considerably so i can...

Dear all, dear Hadley, thank you Hadley for proposing the (in hindsight) quite obvious > devtools::install_github("dmbates/RePsychLing", build_vignettes = TRUE) which brought me on the right track. After fiddling a bit with installing correctly all necessary R-packages...

Did you verify that the correct percentages were above/below the regression lines? I did a quick check and for example did not consistently get 50% of the observed response values greater than the tau=.5 line. I did when...

Dear all, I am still trying to find out what goes wrong when I use ace (I get an 'NAs in foreign fuction call (arg3)' but when I go through it with debug there seem no data to be missing...