similar to: sem package and growth curves

I am trying to work my way through the book "Singer, JD and Willett, JB, Applied Longitudinal Data Analysis. Oxford University Press, 2003" using R. I have the SAS code and S-Plus code from the UCLA site (doesn't include chapter 8 or later problems). In chapter 8, there is a structural equation/path model which can be specified for the sem package as follows S <- cov(al2)

Dear all, I'm writing code to do automatically several calculations and gather the results in a matrix. I have the following vectors and I just want to make a matrix with them: VD1 VD2 VD3 ... VD12 *In this case, this series of vectors stops there, in 12, but it could be a higher number or smaller. (for example: VD1, VD2 and VD3; only) (or for example: VD1, VD2, ... till VD78) All

Hello, I'm running mixed models in GAMM4 with 2 (non-nested) random intercepts and I want to include a spline term for one of my exposure variables. However, when I include a spline term, I always get reported degrees of freedom of less than 1, even when I know that my spline is using more than 1 degree of freedom. For example, here is the code for my model: >

Dear list, I'd like to create a visual plot of a clmm() I've fitted using the 'ordinal' package in R. It's possible to do this with a glm() by using the 'effects' package. For example: library(effects) data(BEPS) mod <- lm(political.knowledge ~ age + gender + vote, data=BEPS) eff <- effect("age", mod, default.levels=100) plot(eff,

All - I'm having an issue with trying to plot a model derived from nls() onto a simple plot.? I have included a sample data set and the code that I've been using. ?? year month day?????? date location mileage? cost gallon????? cpg ? mpg????????? x 2009???? 1?? 4?? 1/4/2009????? BZN? 124585 19.39? 14.37 1.349339 10.71677 2009-01-04 2009???? 1? 15? 1/15/2009????? BZN? 124888? 23.2? 16.12

I have considered opening a bug report on this, but wanted to get some feedback and make sure I am not missing something in the way of a simple work around. What is the scenario in which this impacts your implementation? Ours is the desire to use the same realtime SIP database for many asterisk servers, and route the call based on a "home server" value in the realtime database. The

Howdy, I have a strange issue, first some back story, I failed to upgrade from Debian 6.0 to 7.0 and had to do a fresh install of 7.0 I had a few PV''s and a Windows 7 VM on this system. they still work for the most part for the router I have the following issue that I am not able to reproduce on the other PV''s root@router:/home/kd5lct# iptables --list libkmod: ERROR

Dear All, I am working with linear mixed-effects models using the lme4 package in R. I created a model with the lmer function including some main effects, a two-way interaction and a random effect. Now I am searching how I could incorporate an interaction between the random effect and one of the fixed effects. I tried to express the interaction in:

I'm using 1.6.2 on a win 2k box. (I know I'm due for an upgrade.) I've used brute force to label a series of vectors, and I'm wondering if there's a better way to do this. I wanted to create a biplot of the 2nd & 3rd PCs, and did this: test <- edit(loadings(cv.prc.spr)) test.1 <- test[,1] test.2 <- test[,2] test.3 <- test[,3] x0 <- 0 y0 <- 0 i <-

Hello all, In lme4 if you want to model two non-nested random effects you code it like this: mod1 <- lmer(y~x + (1|randomvar1) + (1|randomvar2)) How would you go about to model something similar in nlme? In my database I have two variables for which I have repeated measures, lets call them "individual" and "year". But none of the "individuals" were measured in

(Warning: There is some potential that this message is redundant.) I think that i have spotted an issue with R CMD check that does not trace back to a rare case, so it must have been reported already or i am doing something stupid. However, to be sure and because I did not find any comments on this while searching the R mailing lists, here is what i did: R CMD check . from within a package

It appears that file.copy() does not protect against the case when one copies one directory to a subdirectory of itself resulting in it trying to create and endless recursive copy of itself. REPRODUCIBLE EXAMPLE: src <- "foo" dir.create(src); file.copy(src, src, recursive=FALSE) # ok file.copy(src, src, recursive=TRUE) # PROBLEM

Take an example where there is two sip users defined in sip.conf as follows; [peer1] Host=192.168.1.1 ... [peer2] Host=dynamic Secret=password ... [Peer3] Config not relevant ... The intention is to accept calls from peer1 without authentication (ip address authentication only), but require authentication from peer2 If by chance a SIP invite comes "From"

Hi, I too have found this matching to be frustrating. I would like it to behave as you describe. Doug -- Doug Meredith 506-854-7997 ext. 801 ________________________________ From: asterisk-users-bounces@lists.digium.com [mailto:asterisk-users-bounces@lists.digium.com] On Behalf Of Damon Estep Sent: Thursday, January 04, 2007 1:50 AM To: Asterisk Users Mailing List - Non-Commercial

I am fitting individual growth models using nlme (multilevel models with repeated measurements nested within the individual), and I am trying to calculate the Pseudo R-squared for the models (an overall summary of the total outcome variability explained). Singer and Willett (2003) recommend calculating Pseudo R-squared in multilevel modeling by squaring the sample correlation between observed and

Hi I was asked if lme can use FIML (Full Information Maximum Likelihood) instead of REML or ML but I don't know the answer. Does anybody know if this is implemented in R? Thanks Francisco

Greetings, Rowland Penny! > Yes, don't post the output of 'testparm -v', just post the output of > 'testparm'. Also don't post a 'diff', post two separate smb.conf files, > better still post them into a post here. > I refuse to try and wade through all those parameters, when most of > them are DEFAULT settings. I wasn't expecting anything

I am trying to write the following code in R. The code works in S+ and i am trying to do the program in R. x=discrim(admit~gpa+gmat,prior=c("uniform"),data=data.mm) i wrote the following in R: x=lda(admit~gpa+gmat,data=data.mm) i could not figure out how to write prior=c("uniform") in R. I would get an error every time. I think that it has something to do with

I am trying to model data in which subjects are followed through time to determine if they fall, or do not fall. Some of the subjects fall once, some fall several times. Follow-up time varies from subject to subject. I know how to model time to the first fall (e.g. Cox Proportional Hazards, Kaplan-Meir analyses, etc.) but I am not sure how I can model the data if I include the data for those

Dear R-help list members I gathered longitudinal data on fish behaviour which I try to analyse using a multi level model for change. Mostly, I am following Singer & Willett (2003), who provide also the S/R code for their examples in the book (e.g. http://www.ats.ucla.edu/stat/Splus/examples/alda/ch4.htm). Of course I am interested in change over time, but I am also very much interested in