similar to: write.table error

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,

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 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

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 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" &

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 would appreciate help translating the following lme model to an lmer function. lme(lognrms ~ Group*Rotation*muscle*side*support*arms, random=~1|Subject/Stratum2/rep, data=Data) Many thanks Ross Darnell r.darnell at uq.edu.au

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

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 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 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

I am introducing the scan() function to my class. Consider the following file (Scanexamp.Rnw ) \documentclass[12pt]{article} \begin{document} <<>>= height = scan() 64 62 66 65 62 69 72 72 70 part = scan(what = character(0)) "Soprano" "Soprano" "Soprano" "Alto" "Alto" "Tenor" "Tenor" "Bass"

> From: Ross Darnell > > Liaw, Andy wrote: > >>From: Liaw, Andy > >> > >> > >>>From: Ross Darnell > >>> > >>>A good point but what is the value of storing a large set of > >>>predicted > >>>values when the values of the explanatory variables are lost > >>>(predicted >

-----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

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

I would like to "keep" the call formulae from the stepAIC sequence. Does anyone know of an example showing how to use the keep argument? Help would be much appreciated. Ross Darnell Email: <r.darnell at shrs.uq.edu.au>

Sweavers As the title suggests I would appreciate any help to include S code in ordinary paragraph mode. I can use the textsl font but is isn't the same. Thanks Ross Darnell

I have not been able to solve a problem trying to put a legend onto an image. After generating the image I try to add a legend with legend(locator(1),c("0","0.2","0.4","0.6","0.8","1"),horiz=T, fill=c("white",rev(heat.colors(5)))) If I place the legend on top of the image the text and fill boxes are shown, but

R^2 for a model is usually defined as 1-RSS/TSS where TSS is the SS about the mean and RSS is the residual SS from the model. Consider the model in R z <- runif(20) y <- z+rnorm(20) my.model <- lm(y~offset(z)) summary(my.model)$r.squared Here the RSS is equivalent to the TSS and gives 0 when it should (IMHO and a few others perhaps) be 1 - RSS/TSS(corrected for the mean only) RSS =

I cannot find the glmmNQ function in the MASS package (or anywhere else I have tried) mentioned on page 296 of MASS4. I would appreciate directions. Thanks -- Ross Darnell