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"

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

Dear R community: it is possible to replace NA?s in a data frame with zeroes? what should I do? Thanks in advance Juan Pablo _________________________________________________________________ MSN Photos es la manera m?s sencilla de compartir e imprimir sus fotos: http://photos.latam.msn.com/Support/WorldWide.aspx -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

> 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