similar to: Translating lme model call to lme4

I am relatively new to R, so this may be a much simpler question than it seems to me. I am trying to create a function that includes two summations, and can't figure out how. I am attaching the equation as a pdf file. This function will then be optimized over a chosen range of values, but right now I just need help with the function. Thank you. Zack Darnell Duke University -- Zack

I have fitted a logistic regression model > failed.lr2$call lrm(formula = failed ~ Age + task2 + Age:task2, data = time.long, na.action = na.omit) using the Design package functions and would like to generate a nomogram from this model. the datadist information is generated and stored in > ddist time.long$Age time.long$task2 Low:effect 45

I would like to create a page of two graphs (2 rows by 1 col) and then draw vertical lines (segments?) on both graphs from the minimum values to the corresponding maximum value. So I have tried # > y <- rnorm(3000) > par(mfrow=c(2,1)) > plot(y,type="l") > plot(cumsum(y),type="l") > segments(1000,min(cumsum(y)),1000,max(cumsum(y))) > par(mfg=c(1,1)) >

I would like to define the order of the levels of a factor. The relevel function would work but since I have 20 levels I would prefer to declare the order explicitly. Using a smaller example levels(oldfactor) "b1" "b2" "r1" "r2" nufactor <- order(oldfactor,order=c("b1","r1","b2","r2")) # my fabricated function

Maybe a useful addition to the predict functions would be to return the values of the predictor variables. It just (unless there are problems) requires an extra line. I have inserted an example below. "predict.glm" <- function (object, newdata = NULL, type = c("link", "response", "terms"), se.fit = FALSE,

Using R (version 0.63.3) for MS windows, I try the following command > shell(paste("cd ",getenv("RHOME"),sep="")) which replies with an error message Too many parameters - FILES\RW0633 which appears to suggest that the space if the path name is causing difficulties to the cd command. getenv returns > getenv("RHOME") RHOME

John, I spoke with Level(3) last week regarding SIP termination. They quoted $0.01/minute, with an 11 Million Minute / Month minimum. Ugh! -dg -------------------------------------------------------------- Darnell Gadberry President binaryMedia darnell AT binmedia DOT com ------------ Date: Fri, 19 Dec 2003 21:12:22 -0500 To: asterisk-users@lists.digium.com From: John Todd

I am using the Windows menu functions below which will work on the first pass, but if I repeat the same script I cannot get the WinMenuAddItem to work. This is a problem if I change the menu structure and reread the source code I am forced to quit and restart Rgui. "try.menu" <- function(){ OS <- .Platform$OS.type GUI <- .Platform$GUI if (!(OS == "windows" &

-----Original Message----- >From: Tommy Christensen [mailto:tommy.christensen@tpack.net] >Sent: Wednesday, October 08, 2003 10:09 AM >To: Christian Darnell >Cc: 'Linux 802.1Q VLAN'; Bart De Schuymer; netdev@oss.sgi.com; bridge >Subject: Re: [Bridge] RE: [VLAN] Re: [PATCH/RFC] Let {ip, arp}tables >"see" bridged VLAN tagged{I,AR}P packets > > > >This

Full_Name: Ross Darnell Version: 2.4.1 OS: Linux Submission from: (NULL) (130.102.133.33) rdarnell at pc:~/Data$ ls -al .RData -rwxrwxrwx 2 root root 13551 2006-12-06 08:58 .RData rdarnell at pc:~/Data$ R R version 2.4.1 (2006-12-18) Copyright (C) 2006 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

The FAQ says in [5.2 How can add-on packages be installed?] to install a package to a private tree, use [$ R INSTALL -l lib pkgdir_1 ... pkgdir_n] where lib gives the path to the library tree to install to. which, for me, returns Package '-l' does not exist. Has there been a change? Thank you Ross -- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~. |Ross

Is it possible to produce a histogram directly using the hist() function with the common borders removed? It can be done by plotting the histogram object using type 's'teps. my.hist <- hist(x,plot=FALSE) plot(my.hist$breaks,c(0,my.hist$counts),type='s') I would appreciate help Ross Darnell -- University of Queensland, Brisbane QLD 4067 AUSTRALIA Email: <r.darnell at

I apologise if this type of question has been answered before, however how do I count the number of TRUE's in the upper half (non-symmetric) matrix of logicals. Thank you Ross Darnell University of Newcastle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I am trying to adjust the sigma symbol down in the y direction with the following command but R recycles the sigma and puts another at -0.1 in the x axis. Any suggestions would be appreciated. mtext(expression(sigma),1,outer=TRUE,adj=c(0.5,-0.1)) Regards Ross Darnell -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I am trying to place a legend onto an image but the legend to clipped/nonexitent outside the boundary of the image. I am wanting to place the legend where the main title normally would appear. I would appreciate someone telling me why I cannot do this. (Which part of the help file I haven't read most likely) Regards Ross Darnell

The Rows function which is called from plot.compareFits in the nlme package is not found. > plot(compareFits(coef(bp.model3),coef(bp.model3M))) Error in plot.compareFits(compareFits(coef(bp.model3), coef(bp.model3M))) : couldn't find function "Rows" > Can I find it elswhere? Have I missed a required package? Thanks Ross Darnell > library(help=nlme) nlme

I am getting the following when trying to write a dataframe to a file. I would appreciate any advice to where the problem lies. Thanks Ross Darnell > write.table(all,file="agree.csv",row.names=F) Error in "rownames<-"(*tmp*, value = row.names(x)) : length of dimnames[1] not equal to array extent > is.data.frame(all) [1] TRUE > dim(all) [1] 854 18 >

I think I have a data.frame called "statlab". > is.data.frame(statlab) [1] TRUE and I can perform summary, dim, ... but when I try > attach(statlab) Error in attach(statlab) : attach only works for lists and data frames I am using > version _ platform i686-unknown-linux arch i686 os linux system i686,

I am wanting to generate panels showing scatterplots with the linear fitted line for two groups within each panel superimposed. I have two conditioning factors, "variable" and "Group" and I want separate panels for each level of "variable" with different symbols and "lmline"s for each level of "Group". However all observations for a group are

The following is evidence of what is surely an undesirable feature. The issue is the handling, in calls to model.frame(), of an explanatory variable that has been derived as an unclassed factor. (Ross Darnell drew this to my attention.) ## Data are slightly modified from p.191 of MASS > worms <- data.frame(sex=gl(2,6), Dose=factor(rep(2^(0:5),2)), +