Search Archives

Browse[1]> j <- c(1,2,NA) Browse[1]> j[j==1][-1] [1] NA Browse[1]> j[j==1][-2] [1] 1 Browse[1]> j[j==1][-2] <- 2 Error during wrapup: NAs are not allowed in subscripted assignments...

That is correct: clmm in package ordinal cannot fit random slope models (or other types of vector-valued random-effects). Perhaps MCMCglmm will be able to assist? Cheers, Rune On 4 October 2012 23:29, Eiko Fried <torvon at gmail...

Duncan Temple Lang wrote (Mon Nov 7 22:35:22 CET 2005): >R is not yet thread safe. >We are working on it, and I hope to make some progress before >the end of the year. (This one even!) > > D...

Simon Urbanek skrev: > Rune, > > On Aug 15, 2008, at 8:50 AM, bugs at philosof.dk wrote: > >> Full_Name: Rune Philosof >> Version: 2.6.2 and 2.7.1 >> OS: Ubuntu 8.04 >> Submission from: (NULL) (130.226.87.177...

On 16 February 2012 02:35, David Duffy <David.Duffy at qimr.edu.au> wrote: > On Wed, 15 Feb 2012, Felipe Nunes wrote: > >> glmer(dummy ~ x + time + (time | subject), data, family=binomial(link >> = "logit"), REML=T, verbose=T) Will data...

When I execute mclapply it creates the needed processes, but these processes never begin computing anything, they just wait indefinitely. I recently upgraded to version 2.9.1, which might have caused the problem. I use Ubuntu 8.10. > sessionInfo...

I am trying to optimize some calculations using the package multicore. I have a function that returns three variables of interest in a list. Normally I would do a for loop over (i in 1:n) but the calculations are...

On 7 June 2018 at 13:00, Michael Dewey <lists at dewey.myzen.co.uk> wrote: > > One thing which occurs to me though is whether the maintaner of lmerTest > would accept it into the package? If it has a different...

If you are willing to leave the GEE world and use a mixed effects version of the proportional odds model, then you can test the proportional odds assumption with the clmm function in package ordinal. You simply include those variables...

Yes; see clm and clmm2 (mixed effects) in the ordinal package for fitting proportional odds models. See section 3 of http://cran.r-project.org/web/packages/ordinal/vignettes/clm_tutorial.pdf to see how to test the proportional odds...

Hi, I am R beginner. I've tried googling and reading, but this might be too simple to be found in the documentation. I have a dissimilarity index (symmetric matrix) from which I have extracted the unique values using the...

On 26 November 2014 at 17:55, Charlotte Whitham <charlotte.whitham at gmail.com> wrote: > Dear Rune, > > Thank you for your prompt reply and it looks like the ordinal package could be the answer I was looking for! > > If you...

Try library(lmerTest) fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy) summary(fm1) Linear mixed model fit by REML Formula: Reaction ~ Days + (Days | Subject) Data: sleepstudy AIC BIC logLik deviance REMLdev 1756 1775 -871.8 1752 1744 Random effects: Groups Name...

On 5 May 2011 01:08, Chris Gast <cmgast at gmail.com> wrote: > I'm having trouble replying directly to the list, so hopefully this gets > through as intended. > > Multinomial mixed models can be programmed in ADMB (www.admb-project...

Dear Markus Example code: library(multicore) mclapply(1:10, function(i) i^2) Markus Schmidberger skrev: > Dear Rune, > > I am using Ubuntu 9.04, R 2.9.1 and multicore 0.1-4. > This combination is working great. Try to...

On 20 March 2018 at 18:34, Douglas Bates <bates at stat.wisc.edu> wrote: > Kind of looks like SPSS went for bug-for-bug compatibility with SAS on this > one. In SAS PROC MIXED, "REPEATED" and "RANDOM" are two...

Peter Dalgaard skrev: > Rune Schjellerup Philosof wrote: > >> Hi >> >> I want to use a function (update) that in its body uses >> eval(call, parent.frame()) >> >> I would like to use this function in a function that does not contain >> the variables...

On 5 January 2015 at 21:08, Ben Bolker <bbolker at gmail.com> wrote: > Roger Coppock <rcoppock <at> cox.net> writes: > >> >> When will "R" implement the "se.fit" option to the >> predict.nls() function? Is there some schedule? >> > > I think...

This is to announce the new R-package ?ordinal? that implements cumulative link (mixed) models for ordinal (ordered categorical) data (http://www.cran.r-project.org/package=ordinal/). The main features are: - scale (multiplicative) as well as location (additive) effects...

Aurore, I don't know if car::Anova is able/should be able to produce anova tables for clmm objects; I usually use drop1() (and sometimes add1) to test terms in CLMMS: > library(ordinal) > fm1 <- clmm(rating ~ temp + contact + (1...