similar to: a question about nlme

Hi, I am using Splus to run a multiphase mixed-effects model. The quations of the models are as below: gf[ij]=b0[i]+b1[i]*age[ij]+b2[i]*max(0,(age[ij]-tau[i]))^2+e[ij] b0[i]=b00+e[i0] b1[i]=b10+e[i1] b2[i]=b20+e[i2] tau[i]=tau+e[i3] i: 1,2,...,100 subjects j: 1,2,...,6 occasions The main scripts of Splus is: simu1<-groupedData(gf~age|id) simu.nlme<-nlme(gf~(b0 + b1 * age + b2 *

Hi everyone, I have posted a similar question to this list, but I don't get a reply. I really want to solve this problem, so I post it again... I am trying to use R to fit some mixed-effects models for a nested data. The data is a simulated data with 111 subjects. Each subject has 6 waves' data. Below are the first two subjects' data : > simu1[1:12,] Grouped Data: gf ~ age |

Dear all, I need to read in 4 sets of tab delimited data in a loop. The 4 data sets are called "simu1.dat", "simu2.dat" and so on. I know what I need on the righthand side of the read.table expression but I can't the left hand side of it to work (see the line in bold below). Can you kindly help? Many thanks. simu1 <- matrix(0,30,3) simu2 <- matrix(0,30,3) simu3

I'm hoping that this is a relatively easy question for someone familiar with the lme4 package. I'm accustomed to using HLM software and writing a simple 2 level [null] equation like this: L1 - Yij = b0 + e L2 - b0 = B00 + u0 The following command in R provides results that are identical to the HLM program. results <- lmer( Y ~ 1 |id , PanelData4) I can't

Hello, everyone! I picked up R several months ago and have adopted it as my choice for statistical programming. Coming from a Java background, I can honestly say that R is not only free, it is better tha S-plus: the lexical scope in R makes it very simple to simulate Java's object model. For this, I encourage everyone to read the artcle: Robert Gentleman and Ross Ihaka (2000) "Lexical

Dear R users, We are happy to announce that the R package 'forestinventory: Design-Based Global and Small-Area Estimations for Multiphase Forest Inventories' is now on CRAN (https://cran.r-project.org/web/packages/forestinventory/). The aim of our package is to provide global- and smallarea estimators for twophase and threephase forest inventories under simple and cluster sampling. The

Dear R users, We are happy to announce that the R package 'forestinventory: Design-Based Global and Small-Area Estimations for Multiphase Forest Inventories' is now on CRAN (https://cran.r-project.org/web/packages/forestinventory/). The aim of our package is to provide global- and smallarea estimators for twophase and threephase forest inventories under simple and cluster sampling. The

Dear all, I have a question regarding setIs and method dispatch in S4 classes: Consider the following example: ##################################################### ## A02 "contains" A01 by setIs: setClass("A01", representation(a="numeric",b="numeric")) setClass("A02", representation(a="numeric",d="numeric"))

It runs more than twice as slowly using 8 core than using a single core in inverting large matrix. Tested on 8 core Windows XP 64 machine. > n = 1000 > n.simu = 100 > func1 = function() + { + x = rnorm(n*n) + dim(x)=c(n,n) + y = solve(x) + } > > setMKLthreads(1) > system.time(for(i in 1:n.simu) func1()) user system elapsed 69.48 2.42 71.91

On Wed, Jul 31, 2013 at 7:40 PM, Daniel O'Connor <doconnor at gsoft.com.au> wrote: > > On 01/08/2013, at 24:37, Gene Cumm <gene.cumm at gmail.com> wrote: >> The Perl version should also be in the source tree >> (utils/isohybrid.in) and built in the archives (util/isohybrid.pl) >> which I'd recommend using everything from the 6.01 binary/source >>

To be fair I don't really understand what the error is telling me so I don't know how to correct it. I'm trying to reshape my data so that I do repeated measures ANOVA on it. The warning message: "Error in d[, timevar] <- times[1L] : subscript out of bounds" The program: library(QRMlib) library(Hmisc) ss<-5 mu<- c(0,2,3,12) mu2<- c(0,2,10,20)

Dear R listers, I know I'm breaking the rules by asking a "homework" related question-- I hope you'll forgive me. I am a social psychology graduate student, and the only one in my department who uses R. I successfully completed my multiple regression and structural equation modeling courses using R (John Fox's car and sem packages were a big help, as was his book).

Dear R developers, I was looking for a function analogous to base::diag() for getting and setting bands of a matrix. The closest I could find was Matrix::band(), but this was not exactly what I wanted for two reasons. Firstly, Matrix::band() returns a matrix rather than just the specified band. Secondly, Matrix::band() cannot be used for setting the values for a matrix band. Setting or

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. The details of the model are: Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij}) y_{ij} = mu_{ij} + w_{ij} * *with* *logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. I'm so sorry. In the last email, I forgot to say that W is also a unknown parameter in the mixed beta regression model. In any case, here I send you the correct formulation. ** Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~

Hi All, I a C function code, which has been successed in integrating with Splus on Unix. Now I want to move it on windows. So I want to integrate this C function with R on windows XP. I download tools.zip and saved it in c:\bin directory, download MinGW and saved it in c:\MinGw directory, download Perl and saved it in c:\Perl directory. My path is set as:

Dear r-help members, Could any of you help me with this model, please? This model gives error when some value touch whatever border and I do not know how to correct it. The 80% of the seeds produced by a plant will fell into the parent cell, the 15% in the first ring according to the king movement (in chess), and a 5% in the second ring defined by the queen2 matrix. Someone said me the functions

Dear all, As our previous email did not get any response, we try again with a reformulated question! We are trying to do something which needs an efficient loop over a huge array, possibly functions such as apply and related (tapply, lapply...?), but can't really understand syntax and examples in practice...i.e. cant' make it work. to be more specific: we are trying to apply a mask

I am out of my league with this question. The following code starts the java imaging program ImageJ from within R, and displays an image (assuming ImageJ is installed on your computer). library(RImageJ) img <- IJ$openImage( file.choose() ) #pick an available .tif file img$show() # make the image object visible # An image is now displayed # find out about the objects involved >

Hello. =) I''m having a bit of trouble with my nfs xen installation and have spent much time trying to research/figure out this problem with no luck, so I''m hoping someone might be able to help me here. =) I''ve successfully installed xen 3.02 on Linux RH4 on my Sun single proc opteron v20z server. I''ve installed what I believe to be a base Linux OS on a