[R] cox regression warning/error messages

2008-12-28 Thread Bob Green
Hello, I am hoping for some advice regarding warning/error messages I received when running a Cox regression # message 1 - obtained while creating a plot of residuals plot (NV.zph, main = groupNUSM - UNFIT, var= 'groupNUSM') Warning messages: 1: In approx(xx, xtime, seq(min(xx), max(xx),

Re: [R] how to generalize the arguments for lm() (r-h...@lists.r-project.org)

2008-12-28 Thread Gabor Grothendieck
Using the built in CO2 data frame this regresses uptake on all the other columns: lm(uptake ~., CO2) On Sun, Dec 28, 2008 at 2:28 AM, Math Girl mathmathg...@yahoo.com wrote: Hello, How could I generalize the following statement for an arbitrary number of columns instead of 7?

[R] how to calculate DWT maximum decomposition level

2008-12-28 Thread mauede
Given a time series of length N, I am trying to figure out its maximum DWT decomposition level. With reference to dwt function of R package wavelets, running the provided example I get the following and wonder how the maximum decomposition level (which is not an integer number) is calculated.

[R] how to calculate the DWT maximum decomposition level

2008-12-28 Thread mauede
Given a time series of length N, I am trying to figure out its maximum DWT decomposition level. With reference to dwt function of R package wavelets, running the provided example I get the following and wonder how the maximum decomposition level (which is not an integer number) is calculated.

[R] 2 player Nash Equilbrium game

2008-12-28 Thread tao . wang
Dear all, I am looking to compute a 2-player Nash Equibrium. Has anybody looked into similar things before? I would like to implement it in R instead of doing it in Matlab. Thanks. Tao [[alternative HTML version deleted]] __

[R] model.matrix and missing values

2008-12-28 Thread Jarrod Hadfield
Hi, Does anyone know an easy way of retaining rows in a model.matrix where missing values are present in the predictors. Ideally I'd be able to retain these rows as zeros. Thanks, Jarrod -- The University of Edinburgh is a charitable body, registered in Scotland, with registration

Re: [R] model.matrix and missing values

2008-12-28 Thread John Fox
Dear Jarrod, If you use na.action=na.exclude in fitting the model, then you can produce rows of NAs in the model matrix with na.resid(); e.g., for a model m, mm - naresid(m$na.action, model.matrix(m)) Then if you really want to convert the rows of NAs in mm to 0s, you could use

[R] how to get degree according to the name of the node?

2008-12-28 Thread Weijia You
Hi all, I have two networks for the same group of the users. I want to compare individual's degree in different networks. But how could I get the degrere of the nodes according to its name? When I use degree(g1), I could only get a vector of the degree of each node. But when I turn to g2, I

Re: [R] Using a constant scale across X-Y plots

2008-12-28 Thread Dieter Menne
Lisa lschwei at mac.com writes: I am working off an example from Deepayan Sarkar's Lattice:Multivariate Data Visualiization with R. I am trying to create Figure 5.6, essentially, but I would like to be able depict different metro areas. These of course have different lat/longs, so I

[R] Line graphs with NA

2008-12-28 Thread Chris Poliquin
Hi, I have sets of three points provided by subjects that I want to graph as lines over a specific range, but the subjects were split into two groups and provided different points. The subjects provided a y-value for the given x-value, for example: Subject: 1 2 3 4

Re: [R] Line graphs with NA

2008-12-28 Thread Gabor Grothendieck
Try this. na.approx fills in missing values. See ?na.approx, ?approx and ?plot.zoo and the three zoo vignettes. Lines - 'Subject: 1 2 3 4 1NA3 2 NA 24NANA 4 366.5 6 5.5 47NA NA7 5

[R] month-day-year to Date

2008-12-28 Thread Terry Therneau
I have a data from with 3 numeric variables, and wish to create a Date object. The old date library had a function mdy.date to do this. I've chased all of the See Also sections I can find for the Date class, and didn't find anything comparable. Yes, I can use as.Date(paste(data$year+1900,

Re: [R] month-day-year to Date

2008-12-28 Thread Gabor Grothendieck
I would use paste but as.Date(ISOdate(...)) is another possibility although it uses paste internally itself. On Sun, Dec 28, 2008 at 2:11 PM, Terry Therneau thern...@mayo.edu wrote: I have a data from with 3 numeric variables, and wish to create a Date object. The old date library had a

[R] Spatial Statistics e-book.

2008-12-28 Thread Marcus Vinicius
Dear R user´s, Is there anyone that may send me *Spatial Statistics* e-book from *Brian Ripley* ? Thanks a lot. Best regards. -- Marcus Vinicius P. de Souza Juiz de Fora - Minas Gerais Brasil [[alternative HTML version deleted]] __

[R] naming rows and columns in a matrix

2008-12-28 Thread Math Girl
Hello, I have a 3 x 7 matrix. I want to assign a name to each row and each column. I tried using dimnames, but for some reason I cannot retrieve the names of the rows or the names of the columns. Is there another way to name the rows and columns and refer to them?

Re: [R] naming rows and columns in a matrix

2008-12-28 Thread jim holtman
?rownames ?colnames x - matrix(nrow=7, ncol=3) rownames(x) - paste('name', 1:7) rownames(x) [1] name 1 name 2 name 3 name 4 name 5 name 6 name 7 On Sun, Dec 28, 2008 at 12:43 PM, Math Girl mathmathg...@yahoo.com wrote: Hello, I have a 3 x 7 matrix. I want to assign a name to each row

[R] Logistic regression with rcs() and inequality constraints?

2008-12-28 Thread Stephan Kolassa
Dear guRus, I am doing a logistic regression using restricted cubic splines via rcs(). However, the fitted probabilities should be nondecreasing with increasing predictor. Example: predictor - seq(1,20) y - c(rep(0,9),rep(1,10),0) model - glm(y~rcs(predictor,n.knots=3),family=binomial)

[R] Random coefficients model with a covariate: coxme function

2008-12-28 Thread Asier Rodríguez
Dear R users: I'm new to R and am trying to fit a mixed model Cox regression model with coxme function. I have one two-level factor (treat) and one covariate (covar) and 32 different groups (centers). I'd like to fit a random coefficients model, with treat and covar as fixed factors and a

[R] help on ylab of xyplot.zoo

2008-12-28 Thread Rosa Trancoso
Hello, I have a zoo object that I would like to plot with lattice, because I need the legend outside the plots. However, I don't want to use the strips, because these figures are for publication. I noticed that if I choose strip=FALSE but add the ylab argument, the xyplot.zoo functions plots all

Re: [R] help on ylab of xyplot.zoo

2008-12-28 Thread Gabor Grothendieck
On Sun, Dec 28, 2008 at 3:06 PM, Rosa Trancoso ar...@ist.utl.pt wrote: Hello, I have a zoo object that I would like to plot with lattice, because I need the legend outside the plots. However, I don't want to use the strips, because these figures are for publication. I noticed that if I

Re: [R] help on ylab of xyplot.zoo

2008-12-28 Thread Rosa Trancoso
Dear Gabor, Thank you for the quick answer. I will use plot.zoo and strech the x-axis with xlim to fit the legend on the right side of the plot. Best regards, Rosa On Sun, Dec 28, 2008 at 10:52 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Sun, Dec 28, 2008 at 3:06 PM, Rosa

Re: [R] Conditional operation on multiple columns from two data frames

2008-12-28 Thread tsippel
The suggestion below was made. df1$Date - as.Date(df1$Date) df2$Date - as.Date(df2$Date) ifelse(df1$ID==df2$ID df1$Date-df2$Date0.5,df1$y-df2$y, NA) However, because my dataframe rows do not align, I need the conditionals to be tested on every combination of cells. I'm starting to think I

Re: [R] Conditional operation on multiple columns from two data frames

2008-12-28 Thread Peter Dalgaard
tsippel wrote: The suggestion below was made. df1$Date - as.Date(df1$Date) df2$Date - as.Date(df2$Date) ifelse(df1$ID==df2$ID df1$Date-df2$Date0.5,df1$y-df2$y, NA) However, because my dataframe rows do not align, I need the conditionals to be tested on every combination of cells. I'm

Re: [R] Conditional operation on multiple columns from two data frames

2008-12-28 Thread jim holtman
Merge will give you an output like below and you should be able to find the value you need: merge(df1, df2, by='ID', all=TRUE) ID Date.x y.x x.x Date.y y.y x.y 1 1NANANA 2007-05-31 23:00:00 -20.82907 217.7022 2 1