Search found 17 matches

by MannyGomez
Mon Oct 05, 2015 11:59 am
Forum: R2MLwiN user forum
Topic: New syntax for runMLwiN: mixed hierarchical model
Replies: 3
Views: 3220

Re: New syntax for runMLwiN: mixed hierarchical model

Thanks Chris. First of all, just wanted to understand a bit better the new syntax. In particular, the level-2 random component. Are the numbers in square brackets 1[1] + 1[2]|USUBJID referring to the continuous outcomes (1 and 2)? The R2MLwiN guide suggests that the square brackets are only used whe...
by MannyGomez
Thu Oct 01, 2015 4:37 pm
Forum: R2MLwiN user forum
Topic: New syntax for runMLwiN: mixed hierarchical model
Replies: 3
Views: 3220

New syntax for runMLwiN: mixed hierarchical model

Hi Chris, I'm using runMLwiN in R to estimate a mixed hierarchical model (2 continuous and 2 binary). The data is multilevel with patients (USUBJID) clustered within groups (STUDYID). In the past I have been using successfully this code: runMLwiN(Formula ="c(walk1, QOL1, probit(death12mo,denomb...
by MannyGomez
Fri Aug 01, 2014 3:12 pm
Forum: R2MLwiN user forum
Topic: R2MLwiN and starting values for MCMC estimation
Replies: 7
Views: 7916

Re: R2MLwiN and starting values for MCMC estimation

Hi, After a helpful discussion with Chris, here are some useful notes than may be of interest to other users: - when prior variance is not definite positive, one may have to fix some of the covariance elements of the covariance matrices stored in C1096. - This can be done via R2MLwiN by including th...
by MannyGomez
Tue Jul 29, 2014 3:26 pm
Forum: R2MLwiN user forum
Topic: R2MLwiN and starting values for MCMC estimation
Replies: 7
Views: 7916

Re: R2MLwiN and starting values for MCMC estimation

Thanks Chris.
The GUI window now opens but the MLwiN is still crashing when trying to estimate the model by MCMC.
I will email you an example dataset.
Manuel
by MannyGomez
Tue Jul 29, 2014 10:41 am
Forum: R2MLwiN user forum
Topic: R2MLwiN and starting values for MCMC estimation
Replies: 7
Views: 7916

Re: R2MLwiN and starting values for MCMC estimation

Hi Chris, Thanks for the suggestion. When I included 'debugmode=TRUE', it did not open a MLwiN window? It should, right? Anyway, I tried to estimate the model manually using MLwiN. I first estimated the model by IGLS to generate the starting values for the MCMC estimation, then went to C1096 and man...
by MannyGomez
Mon Jul 28, 2014 2:30 pm
Forum: R2MLwiN user forum
Topic: R2MLwiN and starting values for MCMC estimation
Replies: 7
Views: 7916

R2MLwiN and starting values for MCMC estimation

Hi, I'm estimating a multilevel mixed (2 continuous, 2 binary) model using R2MLwiN. The model runs fine when estimated by RIGLS. However, when I try estimating the same model using MCMC, I get: MCMC Error 0315: Prior variance matrix is not positive definite The most common suggestion here in the for...
by MannyGomez
Fri May 09, 2014 3:13 pm
Forum: R2MLwiN user forum
Topic: Using R2MLwiN to write BUGS code
Replies: 3
Views: 4156

Re: Using R2MLwiN to write BUGS code

Bill,
Many thanks for your detailed response.
I will compare this approach with the latent Normal variable (Albert and Chib, 1993) and I will let you know.
Best,
Manny
by MannyGomez
Tue May 06, 2014 9:27 pm
Forum: R2MLwiN user forum
Topic: Using R2MLwiN to write BUGS code
Replies: 3
Views: 4156

Using R2MLwiN to write BUGS code

Hi, I'm trying to use R2MLwiN to write some BUGS code of a bivariate mixed model, but I'm getting the error below: mod<- runMLwiN(Formula ="c(y1, probit(y2,cons)) ~ (0s|cons+treat) + (1s|cons.y1) + (2s|cons)", D=c("Mixed", "Normal", "Binomial"), levID = c(&quo...
by MannyGomez
Sat Apr 05, 2014 9:56 pm
Forum: R2MLwiN user forum
Topic: extracting covariance between coefficients
Replies: 2
Views: 4066

Re: extracting covariance between coefficients

Thanks Chris!
Apologies for having missed the FP.cov component included in the object. :oops:
by MannyGomez
Wed Apr 02, 2014 8:17 pm
Forum: R2MLwiN user forum
Topic: extracting covariance between coefficients
Replies: 2
Views: 4066

extracting covariance between coefficients

Dear Richard, After estimating an hierarchical model using R2MLwiN, is there a way of directly extracting the covariance between two estimated (fixed) coefficients? I am looking into the object from runMLwiN, say 'estimates' ----> estimates["estIGLS"] and check the 2nd column (_FP_v) which...