similar to: survival:: plotting survfit with two predictors

Hi, I'm working on fitting a glm model to my data using Gamma error structure and reciprocal link. I've been using dropterm (MASS) in the model simplification process, but the F values from analysis of deviance tables reported by dropterm and anova functions are different - sometimes significantly so. However, the reported residual deviances, degrees of freedom, etc. are not different.

I am trying to run the following code in R on a Linux cluster. I would like to use the full processing power (specifying cores/nodes/memory). The code essentially runs predictions based on a GAM regression and saves the results as a CSV file for multiple sets of data (here I only show two). Is it possible to run this code using HPC packages such as Rmpi/snow/doParallel? Thank you!

I've noticed a change in the way grep() behaves between the 1.9.0 release and a recent R-patched. On 1.9.0 I get the following output: > x <- dget(file = url("http://www.biostat.jhsph.edu/~rpeng/names.R")) > length(grep("^l\\w+tmean", x, perl = TRUE, value = TRUE)) [1] 84 And on R-patched (2004-06-11) I get > x <- dget(file =

I have a set of time-series climate data with missing entries. I need to add rows for these missing entries to this data set. The only way I know to do this is unsing a foreloop, but this won't work on a list. I've tried to convert the list to a data frame, but that won't happen, either. I want to fill rows in this table: > newtest[10:15,] yrmos yearmo snow.sum snow.mean

Dear everyone, I'm encountering an unknown error message when using gamm function: > fitoutput <- gamm(cvd~as.factor(dow)+pm10+s(time,bs="cr",k=15,fx=TRUE)+s(tmean,bs="cr",k=7,fx=TRUE) + ,correlation=corAR1(form=~1|city),family=poisson,random=list(city=~pm10),data=mimp) Maximum number of PQL iterations: 20 iteration 1 iteration 2 iteration 3 iteration 4

Dear List, I'm trying to transfer a script that uses S-plus functions for the kernel and integration calculations to R. This it what it looks like in S-plus: qrm(cpt,don,ect,h,intctk,nbcol,nbl1,nbl2,nbl3,nlc,nlgn,overlap,Results,tdon,tdon1,tdon11,tdon2,tdon21) don<-guiDisplayFileDialog() import.data (DataFrame="tdon",FileName=don,FileType="EXCEL")

I have one data frame with a column of dates and I want to fill another data frame with one column of dates, one of years, one of months, one of a unique combination of year and month, and one of days, but R seems to have some problems with this. My initial data frame looks like this (ignore the NAs in the other fields): > mans[1:10,] date loc snow.new prcp tmin snow.dep tmax 1

All, I am having trouble with a "read.table()" function that is inside of another function. But if I call the function by itself, it works fine. Moreover, if I run the script on a Mac OS X (with the default Mac OS X version of R installed, rev 2.8), it works fine. But it does not work if I run it on windows vista (also default Windows version of R, rev. 2.8). Again, both

Dear list, For educational purposes I have been working with the script below. I have a observation: line 31 #CScal[i] = (amo^2) # IT IS WRONG, I KNOW, BUT IT MAKE R TO CRASHES! I'm thinking this is a possible bug in the R! And I have a couple of doubts: 1) line 51 #curve(dchisq(x, n-1), add = T, col = 'red' I think that it is correct, but the function no

Message: 63 Date: Wed, 26 Jan 2005 04:28:51 +0000 (UTC) From: Gabor Grothendieck <ggrothendieck at myway.com> Subject: Re: [R] chron: parsing dates into a data frame using a forloop To: r-help at stat.math.ethz.ch Message-ID: <loom.20050126T052153-333 at post.gmane.org> Content-Type: text/plain; charset=us-ascii Benjamin M. Osborne <Benjamin.Osborne <at> uvm.edu>

Is there a package in R that can do a variability plot? A variability plot is a kind of categorized dot plot. (If there is a lot of data in each category, box plots are used rather than dot plots.) Usually, the categories are factor level combinations. All the dot plots appear in the same window; below the x-axis a hierarchy of factors shows which dot plot corresponds to which factor-level

--oyUTqETQ0mS9luUI Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In case you want this - and feel free to delete this - here are the data for the bug I just reported concerning missing data in stripchart. You need to run hrp2.R. Please excuse the extraneous junk. This was written mostly by a grad student who was learning R in the process (Andrea Gurmankin). She