Displaying 20 results from an estimated 3000 matches similar to: "how to fit with "lme" function"
2006 Jul 30
2
Question about data used to fit the mixed model
Hi everyone,
I would like to ask a question regarding to the data used to fit the mixed
model.
I wonder that, for the response variable data used to fit the mixed model
(either via "spm" or "lme"), we must have several observations per subject
(i.e. Yij, i = 1,..,M, j = 1,.., ni) or it can be just one observation per
subject (i.e. Yi, i = 1,...,M). Since we have to
2012 Aug 10
1
Solving binary integer optimization problem
Hi,
I am new to R for solving optimization problems, I have set of communication
channels with limited capacity with two types of costs, fixed and variable
cost. Each channel has expected gain for a single communication.
I want to determine optimal number of communications for each channel
maximizing ROI)return on investment) with overall budget as constraint.60000
is the budget allocated.
2013 Apr 23
2
Frustration to get help R users group
Dear R users/developers
I requested help to solve the problem of formulating Multivariate Sample
selection model by using Full Information Maximum Likelihood
(FIML)estimation method. I could not get any response. I formulated the
following code of FIML to analyse univariate sample selection problem.
Would you please advise me where is my problem
library (sem)
library(nrmlepln)
Selection
2014 Jul 09
3
error com un archivo
Hola a todos,
Me gustaria pedir vuestra ayuda a encontrar el error que no consigo
encontrar en este archivo. He revisado todo mil veces y probado y no doy
con ello.Adjunto el archivo con Google drive porque es muy grande.
?
monicap_50.csv
<https://docs.google.com/file/d/0B8o2KrPEgG7ATlBMc19lTVk1d3M/edit?usp=drive_web>
?
Este es el script, y lo que no entiendo que pasa es que tengo 592044
2009 May 21
2
Naming a random effect in lmer
Dear guRus:
I am using lmer for a mixed model that includes a random intercept for a
set of effects that have the same distribution, Normal(0, sig2b). This set
of effects is of variable size, so I am using an as.formula statement to
create the formula for lmer. For example, if the set of random effects has
dimension 8, then the lmer call is:
Zs<-
2012 Nov 27
1
Effect of each term in the accuracy of Nonlinear multivariate regression fitting equation
Dear all,
I have a set of data with 4 inputs (independent variables) and one output
(dependent variable). I want to perform a regression analysis in order to
fit these data to a regression model, however due to the non-linearity of
the model I do not have a clue which equation to use. I am thinking of
starting with a very general equation including ^3 terms and interactions
between the variables
2008 Dec 02
3
sampling from data.frame
Hi all,
I have a data frame with "clustered" rows as follows:
Cu1 x1 y1 z1 ...
Cu1 x2 y2 z2 ...
Cu1 x3 y3 z3 ... # end of first cluster Cu1
Cu2 x4 y4 z4 ...
Cu2 x5 y5 z5
Cu2 ... # end of second cluster Cu2
Cu3 ...
...
"cluster"-size is 3 in the example above (rows making up a cluster are
always consecutive). Is there any faster way to sample n clusters
(with
2013 May 03
1
MANOVA summary.manova(m) :" residuals have rank"
Dear All, I am trying to perform MANOVA. I have table with 504 columns(species) and 36 rows) with two grouping (season and location)
Zx <- Z[c(4:504)]
Zxm <- as.matrix(Z)
m<- manova(Zxm~Season*location, data=Z)
when I do summary.aov, I get respond for each species but summary.manova
summary.manova(m) :" residuals have rank" 24<501.
What can it be the reason for this error
2008 Dec 16
1
surface contour plot help
I am trying to do a surface profile plot.
data is
X Y(1) Z(1)
1-jan-02 2002 number
2-jan-02 2002 number
.
.
.
1-jan-03 2003 (Y2) number Z(2)
2-jan-03 2003 (Y2) number Z(2)
.
.
.
until dec 31 2007.
I used the plot3d funtions to build a scatter point plot.
Call rinterface.rrun("library(rgl)")
Call
2009 Nov 06
3
which data structure to choose to keep multile objects?
I have a function called nnmf which takes in one matrix and returns two
matrices. for example,
> X
[,1] [,2] [,3] [,4]
[1,] 1 4 7 10
[2,] 2 5 8 11
[3,] 3 6 9 12
> z=nnmf(X,2)
> z$W
[,1] [,2]
[1,] 0.8645422 0.6643681
[2,] 1.7411863 0.5377504
[3,] 2.6179287 0.4111063
> z$H
[,1] [,2] [,3] [,4]
[1,] 1.14299486
2011 Oct 29
4
[LLVMdev] [llvm-commits] [PATCH] BasicBlock Autovectorization Pass
Ralf, et al.,
Attached is the latest version of my autovectorization patch. llvmdev
has been CC'd (as had been suggested to me); this e-mail contains
additional benchmark results.
First, these are preliminary results because I did not do the things
necessary to make them real (explicitly quiet the machine, bind the
processes to one cpu, etc.). But they should be good enough for
discussion.
2007 May 24
3
Problem with numerical integration and optimization with BFGS
Hi R users,
I have a couple of questions about some problems that I am facing with
regard to numerical integration and optimization of likelihood
functions. Let me provide a little background information: I am trying
to do maximum likelihood estimation of an econometric model that I have
developed recently. I estimate the parameters of the model using the
monthly US unemployment rate series
2011 Oct 29
0
[LLVMdev] [llvm-commits] [PATCH] BasicBlock Autovectorization Pass
On Sat, 2011-10-29 at 12:30 -0500, Hal Finkel wrote:
> Ralf, et al.,
>
> Attached is the latest version of my autovectorization patch. llvmdev
> has been CC'd (as had been suggested to me); this e-mail contains
> additional benchmark results.
>
> First, these are preliminary results because I did not do the things
> necessary to make them real (explicitly quiet the
2013 Apr 17
0
Full Information Maximum Likelihood estimation method for multivariate sample selection problem
Dear R experts/ users
Full Information Maximum Likelihood (FIML) estimation approach is
considered robust over Seemingly Unrelated Regression (SUR) approach
for analysing data of multivariate sample selection problem. The zero
cases in my dependent variables are resulted from three sources:
Irreverent options, not choosing due to negative utility and not used
in the reported time. FIML can
2012 Dec 05
3
data manipulation between vector and matrix
Dear list,
I was curious how to subtract a vector from matrix?
Say, I have
mat <- matrix(1:40, nrow=20, ncol=2)
x <-c(1,2)
I want,
x-mat[1,] and x-mat[2,], and so on... Basically, subtract column elements
of x against column elements in mat.
But x-mat won't do it.
Thanks,
Mike
[[alternative HTML version deleted]]
2004 Apr 11
1
converting lme commands from S-PLUS to R
I'm trying to do some smoothing with lme and am having some difficulty
bringing commands over from S-PLUS to R. I have the following setup
(modified from Ngo and Wand, 2004):
set.seed(1)
x <- runif(200)
y <- sin(3*pi*x) + rnorm(200)*.4
## library(splines)
z <- ns(x, 4)
The following runs without error on S-PLUS
f <- lme(y ~ 1, random = pdIdent(~ -1 + z))
But in R I get
2011 Oct 29
0
[LLVMdev] [llvm-commits] [PATCH] BasicBlock Autovectorization Pass
On Sat, 2011-10-29 at 15:16 -0500, Hal Finkel wrote:
> On Sat, 2011-10-29 at 14:02 -0500, Hal Finkel wrote:
> > On Sat, 2011-10-29 at 12:30 -0500, Hal Finkel wrote:
> > > Ralf, et al.,
> > >
> > > Attached is the latest version of my autovectorization patch. llvmdev
> > > has been CC'd (as had been suggested to me); this e-mail contains
> >
2011 Oct 29
4
[LLVMdev] [llvm-commits] [PATCH] BasicBlock Autovectorization Pass
On Sat, 2011-10-29 at 14:02 -0500, Hal Finkel wrote:
> On Sat, 2011-10-29 at 12:30 -0500, Hal Finkel wrote:
> > Ralf, et al.,
> >
> > Attached is the latest version of my autovectorization patch. llvmdev
> > has been CC'd (as had been suggested to me); this e-mail contains
> > additional benchmark results.
> >
> > First, these are preliminary
2003 Sep 11
1
Customised legend in lattice
Hi List,
Am trying to customize a legend in trellis: Draws 2x5 lines in 5 colors and
2 linetypes. I would like to add two more items to the legend showing the
key for the line types above the colored legend.
Any suggestions welcome - thanks Herry
#############################
#Following example code:
library(gregmisc)
trellis.device(bg="white")
i1=0
i2=-1.89767506
i3=-1.17087085
2002 Mar 01
2
best way to work with cvs?
I use a script to get the lastest samba:
#!/usr/bin/ksh
echo "\npassword: cvs\n"
cvs -d :pserver:cvs@pserver.samba.org:/cvsroot login
cvs -z5 -d :pserver:cvs@pserver.samba.org:/cvsroot co samba
Are there any other commands I could tack on or something I should do
differently to make sure I am getting correct rebuilds?
For example, should I tack on at the end:
cd