[R] rms::cr.setup and Hmisc::fit.mult.impute

I have fitted a proportional odds model, but would like to compare it to a continuation ratio model. However, I am unable to fit the CR model _including_ imputated data. I guess my troubles start with settuping the data for the CR model. Any hint is appreciated! Christian library(Hmisc)

Re: [R] lme4 and incomplete block design

Many thanks, Bill and Emmanuel! Christian Emmanuel Charpentier schrieb: Le dimanche 08 novembre 2009 à 00:05 +0100, Christian Lerch a écrit : Dear list members, I try to simulate an incomplete block design in which every participants receives 3 out of 4 possible treatment. The outcome

Re: [R] correlated binary data and overall probability

To answer my own question: The package mvtBinaryEP was helpful! Christian Lerch schrieb: Dear All, I try to simulate correlated binary data for a clinical research project. Unfortunately, I do not come to grips with bindata(). Consider corr-data.frame(ID=as.factor(rep(c(1:10), each=5

[R] lme4 and incomplete block design

Dear list members, I try to simulate an incomplete block design in which every participants receives 3 out of 4 possible treatment. The outcome in binary. Assigning a binary outcome to the BIB or PBIB dataset of the package SASmixed gives the appropriate output. With the code below, fixed

[R] data.frame extracting data row-wise

Dear All, I am struggling with extracting data from a data frame: x=data.frame(a=1:11,b=100:110) What I want is a list/vector in this sence: 1 100 2 101 3 102... For single rows, this works fine: as.matrix(x)[1,] For, say 2 rows, this works fine: z=c(as.matrix(x)[1,],as.matrix(x)[2,]) But

Re: [R] data.frame extracting data row-wise

Nice workaround :-). Thank you. Best, Christian Original-Nachricht Datum: Fri, 30 Oct 2009 10:54:49 +0100 Von: Karl Ove Hufthammer k...@huftis.org An: r-h...@stat.math.ethz.ch Betreff: Re: [R] data.frame extracting data row-wise On Fri, 30 Oct 2009 10:42:12 +0100 Christian

[R] correlated binary data and overall probability

Dear All, I try to simulate correlated binary data for a clinical research project. Unfortunately, I do not come to grips with bindata(). Consider corr-data.frame(ID=as.factor(rep(c(1:10), each=5)), task=as.factor(rep(c(1:5),10))) [this format might be more appropriate:]