search for: childid

...ve a data set with over 7000 students with about 4 observations over time per student. I want to examine the within-group fits of a random sample of this group as it takes forever to compute and draw all 7000 regressions. Here is the code I have used so far. >group<-groupedData(math~year|childid, data=scores) >group.list<-lmList(group) >plot(augPred(group.list)) How might I select a random sample of say 100 students so that I can visually examine their trajectories? Thank you ------ Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington St...

...; firefox windows open, but I don't know why I have three. > Looking at this now .. here is how it appears to work: there is a wrapper .. /usr/lib64/firefox/firefox That is one process .. and it starts the rest ========================================== then there is a process called: -childID 1 That seems to be the main broswer process .. all the things besides the actual TABS. =========================================== Then each tab you open spawns a seperate process .. so the first TAB (that is your home page) spawns a process with '-childID 2' ===========================...

Dear List: I was able to run the following in nlme successfully, but the same model and code (same dataset) failed to run in lme4 and gave me the error message below. Any thoughts? lme(math~year, data=egsingle, random=~year|schoolid/childid) Error in lme(formula = math ~ year, data = egsingle, random = structure(list( : Unable to invert singular factor of downdated X'X Harold [[alternative HTML version deleted]]

Johnny Hughes wrote: > On 07/13/2018 02:36 AM, Nicolas Kovacs wrote: > >> Le 12/07/2018 ? 16:36, mark a ?crit : >> >>> I just updated the "critical" firefox update, and it is *seriously* >>> buggy. >> <snip> > Well .. it is working like the upstream one, that's why its released. > If we knew about any bugs other than ones that

...Rev package has repeated measures on students and those students are nested in schools. We can fit a model with random intercepts and slopes for students and students within schools. The first model is fit using lme and the second is fit using lmer. fm3 <- lme(math ~ year, random=~year|schoolid/childid, egsingle) fm4 <- lmer(math ~ year +(year|schoolid:childid) +(year|schoolid), egsingle, control=list(gradient = FALSE, niterEM = 0)) Both result in parameter estimates that are exactly the same. The newest release of lme4 allows for an even easier transition between lme and lmer. Model syntax i...

On 07/13/2018 02:36 AM, Nicolas Kovacs wrote: > Le 12/07/2018 ? 16:36, mark a ?crit : >> I just updated the "critical" firefox update, and it is *seriously* >> buggy. > > I just upgraded Firefox on a sandbox installation, and here's what I > noticed so far. > > 1. On startup, the whole window is black for 1 - 2 seconds, and then > Firefox displays

...evel 3, but not at the other levels. In HLM this is rather easy to do. However, how can one specify a covariate in R for only one of the levels? I have a sample data set with the structure as described above. I fit the unconditional model in R as model1<-lme(math~year, random=~year|schoolid/childid, data=datafile) Now, if I want to enter "female" as a covariate at level 2 only, how might I modify the code to accomplish this? ------ Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria, Virginia 22314 703.647.1628...

Dear list: I am trying to merge two files together from output I get based on the coef() command. Here is what I am running into. I have two simple linear mixed models > mod1.lme<-lme(math~year, data=sample, random=~year|childid/schoolid) > mod2.lme<-lme(math~year, data=sample, random=~year|childid) I then call the coefficients and store them in the following objects using > mod1.coef<-coef(mod1.lme, level=2) > mod2.coef<-coef(mod2.lme) The problem is that the IDs are different for the two dataframes...

Paul: I may have found the issue (which is similar to your conclusion). I checked using egsingle in the mlmRev package as these individuals are strictly nested in this case: library(mlmRev) library(nlme) fm1 <- lme(math ~ year, random=~1|schoolid/childid, egsingle) fm2 <- lmer(math ~ year +(1|schoolid:childid) + (1|schoolid), egsingle) Checking the summary of both models, the output is exactly the same w.r.t. the fixed effects, variance components, standard errors etc. The prior two lines of code fit the same models. However, the following does...

...or from GEE output. -----------sample output using list (geemodel)------------------------ Link: Identity Variance to Mean Relation: Gaussian Correlation Structure: Exchangeable Call: gee(formula = days.sick1 ~ bmi + age + gender + surveyround2 + surveyround3, id = childid, data = dat, family = gaussian, corstr = "exchangeable") Number of observations : 1110 Maximum cluster size : 3 Coefficients: (Intercept) bmi age gender -1.12830236 -0.05216691 -0.01230798 0.14646632 Estimated Scale Parameter: 2.620455 Number of It...

* Please cc me if you reply as I am a digest subscriber * Hi, I am wondering how I can run a multilevel survival model in R? Below is some of my data. > head(bi0.test) childid famid lifedxm sex age delta 1 22.02 22 CONTROL MALES 21.36893 0 2 13.02 13 MAJOR MALES 21.18001 0 3 64.02 64 CONTROL MALES 20.09377 0 4 5.02 5 CONTROL FEMALES 19.94524 0 5 183.02 183 BIPOLAR FEMALES 19.87406 0 6 37.02 37 BIPOLAR M...

...- "TEST0000853" library("RODBC") query <- function(qu){ con <- odbcConnect ("Development ",uid="username",pwd="password",believeNRows=FALSE,rows_at_time=50) sqlQuery(con,qu) } qu <- paste("m.parentid as SET from setmap m where m.childid in (select t.experiment from experimentsets t where t.setname='",inSET,"')",sep="") ExperimentSet <- as.character(query(qu)[,1]) The error only occurs when I call my "query" function. Does anyone know if there is a way to run sql queries in th...

...fitted values at different levels of a multilevel model. My data set is a series of student test scores over time with a total of 7,280 observations, 1,720 students nested witin 60 schools. The data set is not balanced. The model was fit using eg.model.1<-lme(math~year, random=~year|schoolid/childid, data=single). When I call the random effects at all levels using EB.1<-data.frame(ranef(eg.model1, level=1)) and EB.2<-data.frame(ranef(eg.model1, level=2)), I get the shrinkage estimators that I expect. That is, I get 2 random effects for each child (1 intercept and 1 slope) and 2 for...

Hi list I have a question related to an has_many :through association. Basicly I have a form that updates the entire has_many :through collection AND the parent object. from the form I have: params[:parentobject] (all the attributes related to the parent object) and params[:childid] an array of all the destination ID of the has many through association How can I cleanly (and maybe atomicly since I''m using postgres) destroy and create the needed records in the join table within the update statement of the parent? Thanks Paolo

...I did not have a problem creating the graphic when I used the Orthodont data frame or 2 other data sets when there are no missing values. My data has missing values. Do I need to remove the missing values before the lmList function will work? for a small example: > a schoolid spring score childid 1 1 0 550 345 2 1 1 568 345 3 1 0 560 456 4 1 1 NA 456 5 2 0 540 32 6 2 1 562 32 7 2 0 579 34 8 2 1 599 34 > a1=groupedData(score~spring|schoolid,d...

...ng each outcome and with time variables for each outcome. The model was fit by removing the overall intercept, but creating fixed effects and random effects for each using the following: mult2.lme<-lme(fixed=score~-1+read+math+time.m+time.r, data=mult.samp, random=~-1+read+math+time.r+time.m|childid) This worked and seemed to produce reasonable estimates. I then ran a model using only a single outcome (reading) and found that the estimates are very similar, so I am relatively confident in the results of the model. Now, I have a couple of questions regarding the DM LME: 1) If I wanted to...

...9-1 Dear All, ? I was trying to use the lmList function to get the lmList graphic but? have a problem creating the graphic.? My data has missing values and?an error occurred? as stated below. ################################################### library(nlme) > a ??schoolid spring score childid 1 1 0 550 345 2 1 1 568 345 3 1 0 560 456 4 1 1 NA 456 5 2 0 540 32 6 2 1 562 32 7 2 0 579 34 8 2 1 599 34 (lmlis1 <- lmList(score ~ childid | spring, data=a, na.action=T)) Error in UseMethod("getGroups") : no applicable method for "getGroups" ###################...

This is a first draft of a subvolume quota implementation. It is possible to limit subvolumes and any group of subvolumes and also to track the amount of space that will get freed when deleting snapshots. The current version is functionally incomplete, with the main missing feature being the initial scan and rescan of an existing filesystem. I put some effort into writing an introduction into