similar to: Find correlation in Clmm?

Dear list, I'd like to create a visual plot of a clmm() I've fitted using the 'ordinal' package in R. It's possible to do this with a glm() by using the 'effects' package. For example: library(effects) data(BEPS) mod <- lm(political.knowledge ~ age + gender + vote, data=BEPS) eff <- effect("age", mod, default.levels=100) plot(eff,

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

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 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]]

Dear All,   How can I generate multivariate random variable (not multivariate normal )   I am in urgent [[alternative HTML version deleted]]

  I am trying to used "ordgee" from "geepack" for an ordinal dataset. When I write the code it returns  "Warning message:In binomial(link) : use of binomial(link=link) is deprecated" , but the program runs. Even when I run your example for "ohio" and "respdis", it returns the same error. Please guide me [[alternative HTML version deleted]]

Hi all, I am imputing missingness of  90 columns  in a  data frame using mice. But "mice" gives back :  Error in nnet.default(X, Y, w, mask = mask, size = 0, skip = TRUE, softmax = TRUE,  :   too many (1100) weights Any idea to solve this error is welcome, Anera [[alternative HTML version deleted]]

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

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

Hi there, Could anyone please help me to understand what should be done in order not to get this error message: Error: evaluation nested too deeply: infinite recursion / options(expressions=)? Here is my code: determinant<- function(x){det(matrix(c(1.0,0.2,0.5,0.8,0.2,1.0,0.5,0.6,0.5,0.5,0.5,1.0,x,0.8,0.6,x,1.0),ncol=4,byrow=T))} matrix<-

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

Hi, I'm having this Sendmail configuration problem and stuck. I've been doing googling/reading posts and none of the solutions matched my problem. I'm doing some testing, and have been getting this 550 5.1.1 ....User unknown error when I send the email to my company's email address. When I send it to an external email address, such as yahoo, the message delivered successfully.

In this real example (below), all four of the replicates in one treatment combination had zero failures, and this produced a very high standard error in the summary.lm. =20 Just adding one failure to one of the replicates produced a well-behaved standard error. =20 I don't know if this is a bug, but it is certainly hard for users to understand. =20 I would value your comments=20 =20 Thanks =20

Hi folks, Moving the discussion to llvm.dev. None of the changes we talked earlier help. Find attached the C source code that you can use to reproduce the issue. clang --target=aarch64-linux-gnu -c -mcpu=cortex-a57 -Ofast -fno-math-errno test.c -S -o test.s -mllvm -debug-only=licm LICM hoisting to while.body.lr.ph: %21 = load double** %arrayidx8, align 8, !tbaa !5 LICM hoisting to

Can you send me actual LLVM IR or a preprocessed source from using -E? I don't have a machine handy that has headers that target that arch. On Tue Jan 13 2015 at 4:33:29 PM Daniel Berlin <dberlin at dberlin.org> wrote: > Anything other than noalias or mustalias should be getting passed down the > stack, so either that is not happening or CFL aa is giving better answers > and

Yes. I've attached an updated patch that does the following: 1. Fixes the partialalias of globals/arguments 2. Enables partialalias for cases where nothing has been unified to a global/argument 3. Fixes that select was unifying the condition to the other pieces (the condition does not need to be processed :P). This was causing unnecessary aliasing. 4. Adds a regression test to

Hi to all, I initiate in R - Linux and I've some problems to find an editor with R interface as like RWinEdt for WinEdt. Anyone know one? Thanks in advance for your kind cooperation. Best regards Atenciosamente Ana Patricia Martins ------------------------------------------- Serviço Métodos Estatísticos Departamento de Metodologia Estatística Telef: 218 426 100 - Ext: 3210

On 13 January 2015 at 22:11, Daniel Berlin <dberlin at dberlin.org> wrote: > This is caused by CFLAA returning PartialAlias for a query that BasicAA > can prove is NoAlias. > One of them is wrong. Which one? I'm not sure from your description that this is a chaining issue. PartialAlias doesn't chain and isn't supposed to, it's a final answer just like NoAlias and

Oh, sorry, i didn't rebase it when i changed the fix, you would have had to apply the first on top of the second. Here is one against HEAD On Wed, Jan 14, 2015 at 12:32 PM, Ana Pazos <apazos at codeaurora.org> wrote: > Daniel, your patch does not apply cleanly. Are you on the tip? > > The code I see there is no line if (QueryResult == MayAlias|| QueryResult == PartialAlias)

Algo asi? l <- list()l$nombre <- "Ana"l$apellido <- "Pérez"parametro <- "nombre"paste(l[ parametro][[1]], 'OnValidate', sep='')[1] "AnaOnValidate"l[ parametro][[1]][1] "Ana"x <- 1:5AnaOnValidate <- function(x) print(x)if (exists(paste(l[ parametro][[1]], 'OnValidate', sep=''))) do.call(paste(l[