Re: [R] Is it possible to coerce R to continue proceeding t he next command in a loop after an error message ?

2005-07-18 Thread Dieter Menne
Yimeng Lu yl2058 at columbia.edu writes: In a loop, if a function, such as nls, gives an error, is it possible to coerce R to continue proceeding the next command with the same loop? You may also have a look at nlsList in package(nlme). In handles errors returns a nice list of result.

Re: [R] Nested ANOVA with a random nested factor (how to use the lme function?)

2005-07-18 Thread Simon Blomberg
At 01:59 PM 18/07/2005, Addison, Prue wrote: Hi, I am having trouble using the lme function to perform a nested ANOVA with a random nested factor. My design is as follows: Location (n=6) (Random) Site nested within each Location (n=12) (2 Sites nested within each Location) (Random) Dependent

Re: [R] Contingency-Coefficient, Variance

2005-07-18 Thread Christoph Buser
Hi You can easily implement it by: Mag. Ferri Leberl writes: Dear Everybody! Excuse me for this newbie-questions, but I have not found answers to these Question: - Is there a command calculating the variance if the whole group is known (thus dividing by n instead of n-1)?

[R] dataframes of unequal size

2005-07-18 Thread Renuka Sane
I have two dataframes C and C1. Each has three columns viz. state, psu and weight. The dataframes are of unequal size i.e. C1 could be 2/25/50 rows and C has 42000 rows. C1 is the master table i.e. C1$state, C1$psu and C1$weight are never the same. ThisA. P., Urban, 0 is not so for C. For

Re: [R] rmpi in windows

2005-07-18 Thread Uzuner, Tolga
Hi Uwe, Thanks for your kind response. best, Tolga Please follow the attached hyperlink to an important disclaimer http://www.csfb.com/legal_terms/disclaimer_europe.shtml -Original Message- From: Uwe Ligges [mailto:[EMAIL PROTECTED] Sent: 16 July 2005 14:29 To: Uzuner, Tolga Cc:

[R] Quantile Regression, S-Function Rreg

2005-07-18 Thread Stefan Hoderlein
I have the following problem: I would like to do a nonparametric quatile regression. Thus far I have used the quantreg package and done a local quadratic, but it does not seem to work well. Alternatively, I have tried with an older S version I have the function rreg, and used

Re: [R] Proportion test in three-chices experiment

2005-07-18 Thread NOEL Yvonnick
Rafael, when testing binomial hypotheses with both repeated measures and inter-group factors, you should make explicit your model on the intra-subject part of the data. You can't do Chi-square comparisons on count data that mix independent and dependent measures. But you can define a

[R] about 3d surface plot

2005-07-18 Thread 吴 昊
Hi I have a data format like this: x coordinates y coordinates z value 3.77E+002 7.13E+002 0,0 1.27E+003 5.52E+002 2,756785261 1.06E+003 4.76E+002 2,583918174 3.86E+002 7.15E+002 0,158626133 3.60E+002 1.77E+002

[R] column-wise deletion in data-frames

2005-07-18 Thread jhainm
Hi, I have a huge dataframe and like to delete all those variables from it that that have NAs. The deletion of vars should be done column-wise, and not row-wise as na.omit would do it, because I have some vars that have NAs for all rows thus using na.omit I would end up with no obs. Is there a

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Frank E Harrell Jr
Stephen wrote: Hi All, I am currently conducting some survival analyses. I would like to extract coefficients at each level of the IVs. I read on a previous posting that dummy regression using coxph was not possible. I'm not sure what that means. Therefore I though, hey why

Re: [R] column-wise deletion in data-frames

2005-07-18 Thread Chuck Cleland
[EMAIL PROTECTED] wrote: Hi, I have a huge dataframe and like to delete all those variables from it that that have NAs. The deletion of vars should be done column-wise, and not row-wise as na.omit would do it, because I have some vars that have NAs for all rows thus using na.omit I would

Re: [R] column-wise deletion in data-frames

2005-07-18 Thread Peter Dalgaard
Chuck Cleland [EMAIL PROTECTED] writes: data - as.data.frame(cbind(X1,X2,X3,X4,X5)) So only X1, X3 and X5 are vars without any NAs and there are some vars (X2 and X4 stacked in between that have NAs). Now, how can I extract those former vars in a new dataset or remove all those

Re: [R] about 3d surface plot

2005-07-18 Thread Huntsinger, Reid
Have a look at the rgl package on CRAN. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Monday, July 18, 2005 8:24 AM To: r-help@stat.math.ethz.ch Subject: [R] about 3d surface plot Hi I have a data format like

[R] levels() deletes other attributes

2005-07-18 Thread Heinz Tuechler
Dear All, it seems to me that levels() deletes other attributes. See the following example: ## example with levels f1 - factor(c('level c','level b','level a','level c'), ordered=TRUE) attr(f1, 'testattribute') - 'teststring' attributes(f1) levels(f1) - c('L-A', 'L-B', 'L-C') attributes(f1) If

[R] how to change bar colours in plot.stl

2005-07-18 Thread Michael Townsley
Dear helpeRs, Is it possible to change the shading colour of the range bars in the plot generated by plot.stl? By default they are grey, but I would prefer them white (I am preparing some graphics for a powerpoint presentation so I'm inverting all colours). As far as I can see plot.stl

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Stephen
Hi 1. To clarify: There is a posting saying that dummy regression using the coxph function is not possible... That posting may be outdated 2. Q. You say 'Make sure that eventbefore is a pre time zero measurement' please explain: Do you mean that if someone who is not left censored and has no

Re: [R] column-wise deletion in data-frames

2005-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2005, Peter Dalgaard wrote: Chuck Cleland [EMAIL PROTECTED] writes: data - as.data.frame(cbind(X1,X2,X3,X4,X5)) So only X1, X3 and X5 are vars without any NAs and there are some vars (X2 and X4 stacked in between that have NAs). Now, how can I extract those former vars

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Prof Brian Ripley
This is almost unreadable. On Mon, 18 Jul 2005, Stephen wrote: Hi 1. To clarify: There is a posting saying that dummy regression using the coxph function is not possible... That posting may be outdated Clarification needs 1) an explanation of what you mean by `dummy regression' and 2) a

Re: [R] levels() deletes other attributes

2005-07-18 Thread Heinz Tuechler
At 09:29 18.07.2005 -0500, Frank E Harrell Jr wrote: Heinz Tuechler wrote: Dear All, it seems to me that levels() deletes other attributes. See the following example: ## example with levels f1 - factor(c('level c','level b','level a','level c'), ordered=TRUE) attr(f1, 'testattribute') -

Re: [R] how to change bar colours in plot.stl

2005-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2005, Michael Townsley wrote: Dear helpeRs, Is it possible to change the shading colour of the range bars in the plot generated by plot.stl? By default they are grey, but I would prefer them white (I am preparing some graphics for a powerpoint presentation so I'm inverting

Re: [R] column-wise deletion in data-frames

2005-07-18 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: On Mon, 18 Jul 2005, Peter Dalgaard wrote: Chuck Cleland [EMAIL PROTECTED] writes: data - as.data.frame(cbind(X1,X2,X3,X4,X5)) So only X1, X3 and X5 are vars without any NAs and there are some vars (X2 and X4 stacked in between that

Re: [R] levels() deletes other attributes

2005-07-18 Thread Heinz Tuechler
At 16:53 18.07.2005 +0200, Heinz Tuechler wrote: At 09:29 18.07.2005 -0500, Frank E Harrell Jr wrote: Heinz Tuechler wrote: Dear All, it seems to me that levels() deletes other attributes. See the following example: ## example with levels f1 - factor(c('level c','level b','level a','level

[R] colnames

2005-07-18 Thread Gilbert Wu
Hi, I have a matrix with column names starting with a character in [0-9]. After some matrix operations (e.g. copy to another matrix), R seems to add a character 'X' in front of the column name. Is this a normal default behaviour of R? Why has it got this behaviour? Can it be changed? What

Re: [R] Quantile Regression, S-Function Rreg

2005-07-18 Thread Stefan Hoderlein
Dear Brian, thanks for your mail. For other reasons I need a local polynomial. The nonparametric regression code is very scetchy, but I have used it as base anyway. Best Stefan __ R-help@stat.math.ethz.ch mailing list

[R] Picking a subset of explanatory variables corresponding to a subset or a transformed subsert of the response

2005-07-18 Thread Luwis Tapiwa Diya
I have a dataframe in which I have used a cerain procedure to select the response a subset of the response and obtained a vector y.Now I need to fit a glm of the response y given x(the covariates) but now I get an error that the variables are not of equal length.How can I come up with a way to

[R] read large amount of data

2005-07-18 Thread Weiwei Shi
Hi, I have a dataset with 2194651x135, in which all the numbers are 0,1,2, and is bar-delimited. I used the following approach which can handle 100,000 lines: t-scan('fv', sep='|', nlines=10) t1-matrix(t, nrow=135, ncol=10) t2-t(t1) t3-as.data.frame(t2) I changed my plan into using

Re: [R] read large amount of data

2005-07-18 Thread Thomas Lumley
On Mon, 18 Jul 2005, Weiwei Shi wrote: Hi, I have a dataset with 2194651x135, in which all the numbers are 0,1,2, and is bar-delimited. I used the following approach which can handle 100,000 lines: t-scan('fv', sep='|', nlines=10) t1-matrix(t, nrow=135, ncol=10) t2-t(t1)

Re: [R] colnames

2005-07-18 Thread Adaikalavan Ramasamy
This normally happens to me when I read in a table where the rownames will be appended by an X. Read help(make.names) for more information. Remember that R is primarily a statistical software and thus likes colnames classes to be characters. mat1 - matrix( 1:12, nc=3, dimnames=list(NULL,

Re: [R] column-wise deletion in data-frames

2005-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2005, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: On Mon, 18 Jul 2005, Peter Dalgaard wrote: Chuck Cleland [EMAIL PROTECTED] writes: data - as.data.frame(cbind(X1,X2,X3,X4,X5)) So only X1, X3 and X5 are vars without any NAs and there are some vars

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Frank E Harrell Jr
Stephen wrote: Hi 1. To clarify: There is a posting saying that dummy regression using the coxph function is not possible... That posting may be outdated That does not make sense. 2. Q. You say 'Make sure that eventbefore is a pre time zero measurement' please explain: Do you mean

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Frank E Harrell Jr
Stephen wrote: Hi 1. Right perhaps this should clarify. I would like to extract coefficeints for different levels of the IVs (covariate). So for instance, age of onset I would want Hazards etc for every 5 years and so on... The approach I took was to categorize the variables (e.g., age of

Re: [R] read large amount of data

2005-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2005, Thomas Lumley wrote: On Mon, 18 Jul 2005, Weiwei Shi wrote: Hi, I have a dataset with 2194651x135, in which all the numbers are 0,1,2, and is bar-delimited. I used the following approach which can handle 100,000 lines: t-scan('fv', sep='|', nlines=10)

[R] definition of index.array and boot.return in the code for boot

2005-07-18 Thread Obrien, Josh
Dear R friends, I am reading the code for the function boot in package:boot in an attempt to learn how and where it implements the random resampling used by the non-parametric bootstraps. The code contains two (apparent) functions - 'index.array' and 'boot.return' - for which I can find no

Re: [R] Time Series Count Models

2005-07-18 Thread Brett Gordon
Thanks for the suggestion. Is such a model appropriate for count data? The library you reference seems to just be form standard regressions (ie those with continuous dependent variables). Thanks, Brett On 7/16/05, Spencer Graves [EMAIL PROTECTED] wrote: Have you considered lme in

[R] listing datasets from all my packages

2005-07-18 Thread Elizabeth Purdom
Hi, I am using R 2.1.0 on Windows XP and when I type data() to list the datasets in R, there is a helpful hint to type 'data(package = .packages(all.available = TRUE))' to see the datasets in all of the packages -- not just the active ones. However, when I do this, I get the following message:

[R] New functions supporting GIF file format in R

2005-07-18 Thread Tuszynski, Jaroslaw W.
Hi, A minor announcement. I just added two functions for reading and writing GIF files to my caTools package. Input and output is in the form of standard R matrices or arrays, and standard R color-maps (palettes). The functions can read and write both regular GIF images, as well as,

[R] Package Vegan: species accumlumation functions

2005-07-18 Thread Steven K Friedman
Hi everyone, I am working with a data frame consisting of 1009 sampling locations, 138 species incidence and abundance data, and eight forest community types. My goal is to develop species acumulation curves and extrapolated estimates for each community type. I am using the following

[R] how to get dissimilarity matrix

2005-07-18 Thread Baoqiang Cao
Hello All, I'm learning R. Just wonder, any package or function that I can use to get the dissimilarity matrix? Thanks. Best regards, Baoqiang Cao __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] definition of index.array and boot.return in the code for boot

2005-07-18 Thread Spencer Graves
Excellent question. Try 'getAnywhere(index.array)'. It's hidden in namespace:boot. Ditto for boot.return. spencer graves Obrien, Josh wrote: Dear R friends, I am reading the code for the function boot in package:boot in an attempt to learn how and where it

Re: [R] Time Series Count Models

2005-07-18 Thread Paul Johnson
Dear Brett: There are books for this topic that are more narrowly tailored to your question. Lindsey's Models for Repeated Measurements and Diggle, et al's Analysis of Longitudinal Data. Lindsey offers an R package on his web site. If you dig around, you will find many modeling papers on

Re: [R] Time Series Count Models

2005-07-18 Thread Spencer Graves
We are leveraging too far on speculation, at least from what I can see. PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html;. In particular, try the simplest example you can find that illustrates your question, and explain your concerns to us in terms of a

Re: [R] Time Series Count Models

2005-07-18 Thread Brett Gordon
Paul, Thank you so much for your thoughtful reply. I agree - there are many possible descriptions for my data, and I realize that I don't want to get bogged down with figuring out the 'best' model if something simple will work well. For me, I think the difficulty is going to be handling the

[R] svmlight running error

2005-07-18 Thread Luke
Dear R Users, When I used svmlight, I got below error: my command is: foo - svmlight(y~., data= myData) the results: Error in file(con, r) : unable to open connection In addition: Warning messages: 1: svm_learn not found 2: cannot open file '_model_1.txt' myData[1:2,] y X1 X2 X3 X4 X5 X6

[R] Obtaining argument name within a function

2005-07-18 Thread Francisco J. Zagmutt
Dear all How can I obtain the name of the argument passed in a function? Here is a simplistic example of what I would like to obtain: myfunction= function(name) { print(paste(The parameter name was,unknownFunction(name)) } myfunction(myobject) [1] The parameter name was myobject

Re: [R] Obtaining argument name within a function

2005-07-18 Thread Sebastian Luque
Francisco J. Zagmutt [EMAIL PROTECTED] wrote: Dear all How can I obtain the name of the argument passed in a function? Here is a simplistic example of what I would like to obtain: myfunction= function(name) { print(paste(The parameter name was,unknownFunction(name)) }

Re: [R] Obtaining argument name within a function

2005-07-18 Thread John Sorkin
Francisco, I had exactly the same question a few days ago. Try the following: z-function(x) { xName - deparse(substitute(x)) cat(xName) cat(The parameter name was ,xName,\n) } Let me know if this works. n.b. the \n in the cat function is a signal to print the next piece of data on the next

[R] help: how to change the column name of data.frame

2005-07-18 Thread wu sz
Hello, I have a data frame with 15 variables, and want to exchange the data of 4th column and 6th column. First I append a column in the data frame, copy the 4th column data there, then copy the 6th column data to 4th column, and copy the appended column data to 6th column, but the names of the

Re: [R] Survival dummy variables and some questions

2005-07-18 Thread Stephen
Many thanks I follow you what you say You can request predicted values at any sequence of ages - I guess there are plenty of postings on how to do that Regards, Stephen - Original Message - From: Frank E Harrell Jr To: Stephen Cc: Prof Brian Ripley ; Sent: Monday, July 18, 2005