r - issues with groupedData() -

-

i'm following tutorial on mixed effects models. tutorial uses egsingle dataset mlmrev package. part of tutorial, author uses groupeddata() as:

egsingle <- groupeddata(math ~ year | schoolid/childid, data = egsingle)

could me understand "schoolid/childid" refers to?

please note schoolid , childid both factors!

also, later in tutorial, author takes sample of size 50 , uses lmlist() fit ols regression each subject using:

egsingle <- groupeddata(math ~ year | schoolid/childid, data = egsingle) samp <- sample(levels(egsingle$childid), 50) level2.subgroup <- subset(egsingle, childid %in% samp)  # fitting separate ols regression line each student level2 <- lmlist(math ~ year | childid, data = level2.subgroup) plot(augpred(level2))

when run lmlist command above, these errors:

error in eval(expr, envir, enclos) : object 'childid' not found in addition: warning messages: 1: in lmlist(math ~ year | childid, data = level2.subgroup) : lmlist not (yet) work correctly on groupeddata objects 2: in ops.factor(schoolid, childid) : ‘/’ not meaningful factors

could me figure out why these errors?

  1. as roman luštrik comments, schoolid/chilidid means "school id" , "child id nested within school id" both grouping variables. nesting format formally constructs interaction between higher , lower levels; heuristically, lets computer know "child 1 in school 1" , "child 1 in school 2" different individuals.

  2. you're having problem conflicts between versions of lmlist in nlme , lme4 packages. if run exactly these lines clean r session:

## load data without loading package & dependencies data(egsingle, package="mlmrev") library("nlme")  egsingle <- groupeddata(math ~ year | schoolid/childid, data = egsingle) samp <- sample(levels(egsingle$childid), 50) level2.subgroup <- subset(egsingle, childid %in% samp)  # fitting separate ols regression line each student level2 <- lmlist(math ~ year | childid, data = level2.subgroup) plot(augpred(level2))

it should work fine. should work if library("mlmrev") before load nlme (so nlme before lme4 in search path), or if explicitly specify nlme::lmlist don't accidentally pick lme4::lmlist.


Comments

Post a Comment

Popular posts from this blog

qt - Using float or double for own QML classes -

-
i've created components , helper methods using c++ , qt used in qml gui. of methods calculate gui elements (size, transformations, scene graph items creation). target arm 32 bit using opengl / eglfs render output. opengl using float , qml using real (defined double ). what should use in c++ do not lose precision? if use qreal ( double ) shouldn't lose precision. if git grep float in qtbase.git , can find out bit more: dist/changes-5.2.0-**************************************************************************** dist/changes-5.2.0-* important behavior changes * dist/changes-5.2.0-**************************************************************************** dist/changes-5.2.0- dist/changes-5.2.0- - qt compiled qreal typedef'ed double on dist/changes-5.2.0: platforms. qreal float on arm chipsets before. guarantees more dist/changes-5.2.0- consistent behavior between platforms qt supports, binary dist/chan...
Read more

Create Outlook appointment via C# .Net -

-
i need create outlook appointment using microsoft.office.interop.outlook , while can work locally on own workstation , works if run web application via server's browser, not work when connect server externally. hence, think might permission issue. i changed application pool identity " localsystem " don't access denied error. unfortunately, doesn't work. the app behaves if appointment created, appointment doesn't appear in outlook. here include file on top of aspx.cs page using outlook = microsoft.office.interop.outlook; here code using pop appointment. outlook.application apptapp = new outlook.application(); outlook.appointmentitem appt = apptapp.createitem(outlook.olitemtype.olappointmentitem) outlook.appointmentitem; appt.subject = txtfirstname.text + " " + txtlastname.text; appt.body = txtcomment.text; appt.alldayevent = false; appt.start = datetime.parse(txtreminderdate.text + " 8:00 am"); appt.end = datetime.p...
Read more

ios - Swift Array Resetting Itself -

-
i have 2 arrays: pointhistory , datehistory . have nstimer , , when timer runs out, action called. in action appending int pointhistory , , current date string datehistory . here timer action: func timeraction() { let steppervalue = int(stepper.value) employee.numberofpoints = steppervalue networking.saveemployee(employee, completion: { (error) -> void in if error != nil { let alert = uialertcontroller(title: "error", message: "there error updating \(self.employee.name)'s points", preferredstyle: .alert) let alertaction = uialertaction(title: "ok", style: .cancel, handler: nil) alert.addaction(alertaction) self.presentviewcontroller(alert, animated: true, completion: nil) } }) var dateformatter = nsdateformatter() dateformatter.dateformat = "mm-dd-yyyy" dateformatter.timezone = nstimezone.localtimezone() var date = nsdate() var ...
Read more