Displaying 3 results from an estimated 3 matches for "reglm".
Did you mean:
realm
2006 Jun 28
1
lme - Random Effects Struture
...is my SAS code
proc mixed data=Reg;
class ID;
model y=Time Time*x1 Time*x2 Time*x3 /S;
random intercept Time /S type=UN subject=ID G GCORR V;
repeated /subject = ID R RCORR;
run; **
(Type =UN for random effects)
The eqivalent lme statement I am using is :
reglme <- lme(y ~ Time+Time*x1+Time*x2+Time*x3, data=Reg, random = ~ Time |
ID)
When I compare the results, the values differ by considerable margin; I
suppose this is due to the Random effects covariance structure. R output
tells me that the structure is
"Structure: General positive-defin...
2002 Apr 01
2
writing a package for generalized linear mixed modesl
Happy new month, everyone!
I am planning to write a NIH grant proposal to study ways to speed
Monte Carlo based maximum likelihood algorithm for hierarchical models
with a focus on generalized linear mixed models (GLM with random
effects). I thought it would be nice and also increase the chance of
funding if I could produce an R package in the process. I understand
that Prof. Pinheiro ang Bates
2002 Apr 12
1
summary: Generalized linear mixed model software
...ne issue of
Bioemtrika. REML estimation has the same asymptotic property as
standard MLE but less bias for finite samples. Of course, one can argue
that Bayesian formulation automatically generates REML for variance
components.
Now the software.
PQL and variations: GLME in beta by José Pinheiro, reglm by Gordon
Smyth, glmmPQL in package MASS by Venables & Ripley
Bayesian method: GLMMgibbs and BUGS (not a R or S-plus package).
Standard MLE: glmm (only a random intercept) in one of Jim Lindsey's
packages?
Ok, this still leaves us without a good way of doing standard MLE. I
have prog...