similar to: How to plot results of clmm()?

Dear List, I am using both the clm() and clmm() functions from the R package 'ordinal'. I am fitting an ordinal dependent variable with 5 categories to 9 continuous predictors, all of which have been normalised (mean subtracted then divided by standard deviation), using a probit link function. From this global model I am generating a confidence set of 200 models using clm() and the

Hi all, I am trying to fit a random effect  model to categorical response variable using package "ordinal" /"clmm". How can I find the correlation between random effects (random intercept and random slope)   Thanks in advance  Ana [[alternative HTML version deleted]]

0down votefavorite <https://stats.stackexchange.com/questions/296569/how-to-obtain-random-slope-random-intercept-plots-for-categorical-response-varia#> I'm trying to generate a random slope random intercept plot after ordinal regression using the clmmfunction from the ordinal package in R. I have aggression levels which are categorical with six levels. Earlier, I made random intercept

Dear Colleagues, We are attempting to create trees using R with our Ruby on Rails application. However, we are running into a problem involving the creation of the graphic. We would like them to be in either jpg or png format so that users can save, but due to a lack of control over our sever we are unable to start X11 server. Is there a way to create these images without using X11 server?

Dear R users, I?d like to pose aquestion about pMCMC and HDP. I have performed a mixed logistic regression by MCMCglmm (a very good package) obtaining the following results: Iterations = 250001:799901 Thinning interval = 100 Sample size = 5500 DIC: 10.17416 G-structure: ~ID_an post.mean l-95% CI u-95% CIeff.samp ID_an 0.7023 0.0001367 3.678 2126 R-structure: ~units post.mean l-95%

Hello, I searched everywhere but I didn't find what I want, that is why I as the question here. Threads discussing this issue on this mailing list are already quite old. Does anybody know of a function in R which allows to estimate ordered probit/logit model with random coefficients. The only mixed effect model I found was clmm of the ordinal package but it only provides random intercepts. I

Dear all, Being a newbie to R, I've trawled through many old posts on this list looking for a solution to my problem, but unfortunately couldn't quite figure it out myself. I'd be very grateful if someone here on this list could perhaps help me out. I have a lattice plot with several panels and would like to add some text next to the y-axis on the right hand side of each row of

I have been working through the book "Applied longitudinal data analysis: modeling change and event occurrence" by Judith D. Singer and John B. Willett. I have been working examples using SAS and also using it as an opportunity for learning to use R for statistical analysis. I ran into some difficulties in chapter 8 which deals with using structural equation modeling. I have tried to

Dear R users, Basically, from the following arbitrary data set: a <- data.frame(id=c(c("A1","A2","A3","A4","A5"),c("A3","A2","A3","A4","A5")),loc=c("B1","B2","B3","B4","B5"),clm=c(rep(("General"),6),rep("Life",4))) > a

Dear R-experts, I need to replace in an expression the character "Cl" by "Cl+beta" But in the following case: form<-expression((Cl-(V *ka) ) +(V *Vm *exp(-(Clm/Vm) *t))) gsub("Cl","(Cl+beta)",as.character(form)) We obtain: [1] "((Cl+beta) - (V * ka)) + (V * Vm * exp(-((Cl+beta)m/Vm) * t))" the character "Clm" has been

Dear list members, I am estimating an ordinal generalized estimating equations (GEE) regression model on repeated measurements data. Is there any way (preferably in R) to test the proportional odds assumption of this model? Thanks in advance. Andreas [[alternative HTML version deleted]]

Hello, I hope this mailing list is correct. I've a cluster pacemaker with a clone resource ocfs2 with ocfs2-tools-1.4.1-25.6.x86_64 ocfs2-tools-o2cb-1.4.1-25.6.x86_64 on Opensuse 11.2 After some network problem on my switch I receive on one of 4 nodes of my cluster the following messages Aug 18 13:12:28 nodo1 openais[8462]: [TOTEM] The token was lost in the OPERATIONAL state. Aug 18 13:12:28

Hi all, I’m trying to integrate Perl’s Sys::Virt into an already existing AnyEvent program. I’m accustomed to writing things like this: use EV; use AnyEvent; use AnyEvent::Handle; my $h = AnyEvent::Handle->new(fh => $fh, …); $h->on_read(sub { … }); EV::run; ## start the event loop I can add some code in the on_read() handler and every time the $fh has something to read, it will fire

I have a big problem with our cluster, if i read my mail with POP3 client (and check the option, keep a copy on server), every time I read my mail, my client retrive ALL messages an the mailbox !!! I test this with: - Apple Mail 1.0 + 2.0 - M$ Outlook 2000 + 2003 - Thunderbird We use: - Dovecot RC1 - Postfix on Debian Sarge (stable) Can you know this bug ? Or have a workaround ? This function

Hi, Sorry if this is a repost. I searched but found no results. I am wondering if it is an easy way to construct the following matrix: r 1 0 0 0 r^2 r 1 0 0 r^3 r^2 r 1 0 r^4 r^3 r^2 r 1 where r could be any number. Thanks. Wen [[alternative HTML version deleted]]

Hello list, I've seen the recent discussions documenting problems with lm. I have encountered the following problem. I use WinXP Pro with service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. > eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP 9.056763e+00 1.815686e+00 WSWOLF 4.703588e+05 1.112832e+05 NPILGRIM 1.017640e+06 2.134335e+05

Hello! I am trying to put two math expressions in the title of a plot. As you can see below, I can place correctly one expression at a time, but not both. Ideally I would like to have them separated by a comma. Any suggestions? > k <- 1 > n.eff <- c(20, 30) > ### this works > plot(0,0, main = substitute(n == k, list(k = k))) > ### this works > plot(0,0, main =

Hello, I am working in a longitudinal study, with a categorical variable as a dependent variable (alcohol consumption: no use, use, abuse and dependence) with repeated measures (baseline, 1 year, 2 years). Besides I have another variable with two groups: control and experimental. I would like to analyze the evolution in each group, I have though in lme, but I am not sure if I can do an lme with

Dear List,   I have some longitudinal data, each patient was followed at times 0, 12, 16, 24 weeks and measure severity of a illness (0-worse, 1-same, 2-better). So, longitudinal response is categorical.  I was wondering whether lmer in R can fit a model for this type of data. If so, how we code? Or any other function in R that can fit this type of longitudinal data? Any suggestion would be

hi, - big and ugly mail. If you don't like them, delete it now :-) - I have collected and classified some information of: http://home.earthlink.net/~rwhron/kernel/bigbox.html And I observed that ext3 performance is worse than previous kernels(2.4.19...). -ac and -aa are here only as reference. Complete information is in the upper URL. dbench: Performance is worse. dbench (Numbers are in