Displaying 20 results from an estimated 100 matches similar to: "How to fit bivaraite longitudinal mixed model ?"
2005 Sep 12
1
Glmm for multiple outcomes
Dear All,
I wonder if there is an efficient way to fit the generalized linear mixed model for multivariate outcomes.
More specifically, Suppose that for a given subject i and at a given time j we observe a multivariate outcome Yij = (Y_ij1, Y_ij2, ..., Y_ijK).
where Y_ijk is a binomial(n_ijk, p_ijk).
One way to jointly model the data is to use the following specification:
g(p_ijk) =
2005 Jul 14
1
A statistical modeling problem
Hi,
This is not an R related question and I apologize for that, but given the
brain power of the R community it is hard for me to resist posting this
here.
I have a problem where each participant is shown a series of visual cues
(displayed on a computer screen in a random order) and asked to respond by
pressing a button (from a finite number of buttons) that corresponds to the
correct
2004 Apr 18
2
lm with data=(means,sds,ns)
Hi Folks,
I am dealing with data which have been presented as
at each x_i, mean m_i of the y-values at x_i,
sd s_i of the y-values at x_i
number n_i of the y-values at x_i
and I want to linearly regress y on x.
There does not seem to be an option to 'lm' which can
deal with such data directly, though the regression
problem could be algebraically
2005 Jun 04
1
can R do Fixed-effects (within) regression (panel data)?
i want to ask 2 questions.
1) can R do Random-effects GLS regression which i can get from Stata?
the following result is frome Stata.can I get the alike result from R?
xtreg lwage educ black hisp exper expersq married union, re
Random-effects GLS regression Number of obs = 4360
Group variable (i) : nr Number of groups = 545
R-sq:
2006 Nov 21
3
Fitting mixed-effects models with lme with fixed error term variances
Dear R users,
I am writing to you because I have a few question on how to fix
the error term variances in lme in the hope that you could help me. To
my knowledge, the closest possibility is to fix the var-cov structure,
but not the whole var-cov matrix. I found an old thread (a few years
ago) about this, and it seems that the only alternative is to write the
likelihood down and use optim or a
2011 Nov 22
1
Generate Simulation
Hallo everybody,
I'm new in r and I"ll appreciate some help!
I have a matrix of nrow=30 and ncoll=54,and I would like to generate 50
simulations with tha same size of the matrix!!!That is to say that I want
to generate 50 matrices -for my 50 simulations - with the same dimensions!
I took my 1st matrix according to the formula that I want to implement:
D<-mean_m + U_i*mat_DELTA
2001 Oct 09
1
PROC MIXED user trying to use (n)lme...
Dear R-users
Coming from a proc mixed (SAS) background I am trying to get into
the use of (n)lme.
In this connection, I have some (presumably stupid) questions
which I am sure someone out there can answer:
1) With proc mixed it is easy to get a hold on the estimated
variance parameters as they can be put out into a SAS data set.
How do I do the same with lme-objects? For example, I can see the
2010 Jun 09
1
equivalent of stata command in R
Dear all,
I need to use R for one estimation, and i have readily available stata command, but i need also the R version of the same command.
the estimation in stata is as following:
1. Compute mean values of relevant variables
. sum inno lnE lnM
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
2008 May 07
0
Help with Mixed effect modeling in R
Hi everyone,
I want to fit the following mixed effect model
Y_ij = b_0i + b_1i * (t_ij*grp_ij == 1) + b_2i * (t_ij*grp_ij == 2) +
v_0i + v_1i*t_ij + e_ij
with a different covariance matrix of random effects for each group.
(Y is the response
t is time
grp is the group indicator
b 's are fixed effects
v 's are random effects)
I know that this is possible in SAS
2010 Jun 09
1
equivalent of stata command in R
From: saint-filth@hotmail.com
To: saint-filth@hotmail.com
Subject: RE:
Date: Wed, 9 Jun 2010 09:53:20 +0000
OK! sorry thats my fault,
here the translations of the stata commands
1st step is to get the mean values of the variables, well that doesnt need explanation i guess,
2nd step is to estimate the model on panel data estimation method
which is:
2010 Feb 05
3
metafor package: effect sizes are not fully independent
In a classical meta analysis model y_i = X_i * beta_i + e_i, data
{y_i} are assumed to be independent effect sizes. However, I'm
encountering the following two scenarios:
(1) Each source has multiple effect sizes, thus {y_i} are not fully
independent with each other.
(2) Each source has multiple effect sizes, each of the effect size
from a source can be categorized as one of a factor levels
2006 Oct 17
1
About compositional data analysis
The compositional data xi=(x_i1,x_i2,...,x_in), for each fixed i , xij>0,
and sum(xij)=1;
I want to compare the mean( u_i) of several groups
i.e.
H0: u_1=u_2=...=u_N
or
H0: u_11=u_21=...=u_N1
Are there any ANOVA tpye tools to do this work in R?
Thanks,
WEN S Q
[[alternative HTML version deleted]]
2011 Aug 26
2
How to generate a random variate that is correlated with a given right-censored random variate?
Hi,
I have a right-censored (positive) random variable (e.g. failure times subject to right censoring) that is observed for N subjects: Y_i, I = 1, 2, ..., N. Note that Y_i = min(T_i, C_i), where T_i is the true failure time and C_i is the censored time. Let us assume that C_i is independent of T_i. Now, I would like to generate another random variable U_i, I = 1, 2, ..., N, which is
2011 Oct 31
1
Question on estimating standard errors with noisy signals using the quantreg package
Dear all,
My question might be more of a statistics question than a question on R,
although it's on how to apply the 'quantreg' package. Please accept my
apologies if you believe I am strongly misusing this list.
To be very brief, the problem is that I have data on only a random draw, not
all of doctors' patients. I am interested in the, say, median number of
patients of
2004 Jul 10
1
Exact Maximum Likelihood Package
Dear R users,
I am a mathematics postdoc at UC Berkeley. I have written a package
in a Computational Algebra System named Singular
http://www.singular.uni-kl.de
to compute the Maximum Likelihood of a given probability distribution over
several discrete random variables. This package gives exact answers to the
problem. But more importantly, it gives All MLE solutions.
My understanding is that
2009 May 07
2
lasso based selection for mixed model
Dear useRs (called Frank Harrell, most likely),
after having preached for years to my medical colleagues to be cautious
with stepwise selection procedures, they chanted back asking for an
alternative when using mixed models.
There is a half dozen laXXX packages around for all types of linear models,
but as far I see there is none for mixed models such as lme. Even
boot.stepAIC (which I
2002 Aug 12
0
help with pseudo-random numbers
Dear People,
I have a vexing problem related to pseudo-random number generation, and
would appreciate any help and advice. This problem is not directly related
to R, and the only reason I am posting it to this list is that my
implementation is using R. Let me describe my problem by giving an
example, that is close to what I am trying to do.
Suppose we are given a stream of pseudo-random numbers,
2012 Nov 29
0
constrOptim
Dear R users,
I am using the function "constrOptim" to minimize the -1*log-likelihood
where \beta_i>=0 i=1,...,p and \beta_0 is unconstrained.
I construct u_i as
0 0 0 ... 0
0 1 0 ... 0
0 0 1 ... 0
. . . ... 0
. . . ... 0
.
2008 Aug 30
1
Exiting ssh when MaxSessions=0
Hi,
I've been experimenting with MaxSessions=0 in the sshd_config and have
encountered one unfortunate problem. Once the client authenticates to
the server, it ceases to respond to keyboard input.
At first glance, it looks like the client is in a hung state and does
not time out. If port forwarding was requested on the command-line and
the server accepts the request, that continues to work.
2004 Apr 02
0
Hessian in constrOptim
Dear R-users,
In the function constrOptim there is an option to get an approximation
to the hessian of the surrogate function R at MLE by declaring
hessian=TRUE in the calls to the function optim. I would like to ask
if it is advisable to get an approximate hessian for the funcrion f as
follows:
f''(theta)=R''(theta|theta_k)-B''(theta)
where