Re: [R] Score Test Function

2011-06-11 Thread Bill.Venables
The score test looks at the effect of adding extra columns to the model matrix. The function glm.scoretest takes the fitted model object as the first argument and the extra column, or columns, as the second argument. Your x2 argument has length only 3. Is this really what you want? I would h

[R] Score Test Function

2011-06-11 Thread Tyler Rinker
Greeting R Community, I'm trying to learn Logistic Regression on my own and am using An Introduction to Logistic Regression Analysis and Reporting (Peng, C., Lee, K., & Ingersoll, G. ,2002). This article uses a Score Test Stat as a measure of overall fit for a logistic regression model. The

Re: [R] How to do Clustered Standard Errors for Regression in R?

2011-06-11 Thread kevingoulding
Hi Bo, See the following blog post: http://thetarzan.wordpress.com/2011/06/11/clustered-standard-errors-in-r/ http://thetarzan.wordpress.com/2011/06/11/clustered-standard-errors-in-r/ where I calculate clustered standard errors in R. The result is the same as the "cluster" command in Stata. I

Re: [R] How to compute the P-value for a mixture of chi^2 distributions

2011-06-11 Thread Tiago Pereira
thank you very much for your advice, Thomas! Extremely helpful! Cheers! Tiago -- View this message in context: http://r.789695.n4.nabble.com/How-to-compute-the-P-value-for-a-mixture-of-chi-2-distributions-tp3591276p3591365.html Sent from the R help mailing list archive at Nabble.com. _

[R] smoothScatter function (color density question) and adding a legend

2011-06-11 Thread Clayton K Collings
Dear R experts, I am resending my questions below one more time just in case someone out there could help but missed my email. I don't think my questions are too hard. I am most concerned about the transformation function. See below. Thanks, Clayton Hello, I have a few questions, regarding th

Re: [R] How to compute the P-value for a mixture of chi^2 distributions

2011-06-11 Thread Thomas Lumley
On Sun, Jun 12, 2011 at 12:44 PM, Tiago Pereira wrote: > The test I am working on has an asymptotic 0.5*chi2(1)+0.5*chi2(2) > distribution, where numbers inside parenthesis stand for the degrees of > freedom. Is is possible to compute quickly in R the cumulative distribution > of that distributio

[R] How to compute the P-value for a mixture of chi^2 distributions

2011-06-11 Thread Tiago Pereira
Hello, The test I am working on has an asymptotic 0.5*chi2(1)+0.5*chi2(2) distribution, where numbers inside parenthesis stand for the degrees of freedom. Is is possible to compute quickly in R the cumulative distribution of that distribution? Thanks in advance. Tiago -- View this message in co

Re: [R] Factor Analysis/Inputting Correlation Matrix

2011-06-11 Thread Joshua Wiley
Hi Matt, Did you try reading the documentation for factanal()? You can pull it up by typing: help("factanal") These give basically identical results using the raw data, the covariance matrix, and the correlation matrix. factanal(x = mtcars, factors = 3) factanal(factors = 3, covmat = cov(mtcars

[R] Factor Analysis/Inputting Correlation Matrix

2011-06-11 Thread Matt Stati
Can someone please direct me to how to run a factor analysis in R by first inputting a correlation matrix? Does the function "factanal" allow one to read a correlation matrix instead of data vectors? Thanks, Matt. [[alternative HTML version deleted]]

Re: [R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Mike Marchywka
> Date: Sat, 11 Jun 2011 19:57:47 +0200 > Subject: Re: [R] Amazon AWS, RGenoud, Parallel Computing > From: lui.r.proj...@googlemail.com > To: marchy...@hotmail.com > CC: r-help@r-project.org > > Hello Mike, > [[elided Hotmail spam]] > Best to my k

[R] Is there an implementation of loess with more than 3 parametric predictors or a trick to a similar effect? [re-posting as plain text to pass char-set filter]

2011-06-11 Thread Dr. D. P. Kreil (Boku)
Dear R experts, I have a problem that is a related to the question raised in this earlier post     https://stat.ethz.ch/pipermail/r-help/2007-January/124064.html My situation is different in that I have only 2 predictors (coordinates x,y) for local regression but a number of global ("parametric")

Re: [R] Help with Median test and Coxon-Mann-Whittney test

2011-06-11 Thread Ethan Brown
Hi Rob, This list is primarily intended for questions about how to do things in R, so you're more likely to get a helpful response elsewhere. You might want to try some place like the Cross-Validated web site ( http://stats.stackexchange.com/) for general statistics and data analysis questions. B

Re: [R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Lui ##
Hello Mike, thank you very much for your response! Best to my knowledge the sort algorithm implemented in R is already "backed by C++" code and not natively written in R. Writing the code in C++ is not really an option either (i think rGenoud is also written in C++). I am not sure whether there re

[R] Is there an implementation loess with more than 4 parametric predictors or a trick to similar effect?

2011-06-11 Thread Dr. D. P. Kreil (Boku)
I have revised my earlier question to and would be grateful for any comments! -- Subject: Is there an implementation of loess with more than 3 parametric predictors or a trick to a similar effect? Dear R experts, I have a problem that is a related to the question raised in this earlier

Re: [R] running R commands asynchronously

2011-06-11 Thread Martin Morgan
On 06/10/2011 02:29 PM, Richard M. Heiberger wrote: I am interested in running R commands asynchronously. My first choice is in the same R session that I am currently in. Here, the goal would be to run something like RunAsynchSameSession(myfunction(), "outputname.rda") Once RunAsynchSame

Re: [R] Inspecting C code in an R package

2011-06-11 Thread Barry Rowlingson
On Sat, Jun 11, 2011 at 5:31 PM, Layman123 wrote: > Thank you! Of course, I will read the posting guidelines. A subscriber helped > me via e-mail telling me to use the grep-command, that is type in: grep > "somename" *.c. For Windows users it's: findstr "somename" *.c. The problem here is that t

Re: [R] plot of tree

2011-06-11 Thread John
On Saturday, June 11, 2011 04:26:44 AM Shuiwang Ji wrote: > friends, > > I need some tool to visualize the results of hierarchical clustering. > Specifically, I want to plot it as a radial plot as a phylogenetic tree. In > addition, I want to specify the color to each leaf node. I search all > phy

Re: [R] Inspecting C code in an R package

2011-06-11 Thread Layman123
Thank you! Of course, I will read the posting guidelines. A subscriber helped me via e-mail telling me to use the grep-command, that is type in: grep "somename" *.c. For Windows users it's: findstr "somename" *.c. -- View this message in context: http://r.789695.n4.nabble.com/Inspecting-C-code-in

Re: [R] running R commands asynchronously

2011-06-11 Thread Greg Snow
Good points, I was going to include something along the lines of don't change any data in the console that is involved in the asynchronous process, but figured Richard was smart enough to not do something like that. But your points add to that and make it clear for others who read this that don

Re: [R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Mike Marchywka
> Date: Sat, 11 Jun 2011 13:03:10 +0200 > From: lui.r.proj...@googlemail.com > To: r-help@r-project.org > Subject: [R] Amazon AWS, RGenoud, Parallel Computing > > Dear R group, > > [...] > I am a little bit puzzled now about what I could do... It seems

Re: [R] Inspecting C code in an R package

2011-06-11 Thread Jeff Newmiller
Presumably in the source code of the package (something.tar.gz). Follow the posting guidelines if you want further assistance. --- Jeff Newmiller The . . Go Live... DCN: Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead:

Re: [R] Inspecting C code in an R package

2011-06-11 Thread Layman123
Thank you Sarah for your quick answer! I've just downloaded the source package, but now I don't know in which file the C-Code is stored ".C("somename",)" is calling - there is no file with the name "somename". How could one figure that out? -- View this message in context: http://r.789695.n4.

Re: [R] Inspecting C code in an R package

2011-06-11 Thread Sarah Goslee
The easiest thing to do is download the source package from your local CRAN mirror. That will contain all the R and other code. Sarah On Sat, Jun 11, 2011 at 10:50 AM, Layman123 wrote: > Hello everyone, > > Trying to comprehend code of an R package, I encountered the problem that > the interesti

[R] Inspecting C code in an R package

2011-06-11 Thread Layman123
Hello everyone, Trying to comprehend code of an R package, I encountered the problem that the interesting part of the function I'm inspecting is written in C-Code and called by ".C("somename", )". Now I can't inspect the C-Code the function is calling since I can't find it in the folder of the

[R] plot of hierarchical clustering results

2011-06-11 Thread Shuiwang Ji
friends, I need some tool to visualize the results of hierarchical clustering. Specifically, I want to plot it as a radial plot as a phylogenetic tree. In addition, I want to specify the color to each leaf node. I search all phylogenetic tree plotting routines here, they all cannot show the leaf n

[R] plot of tree

2011-06-11 Thread Shuiwang Ji
friends, I need some tool to visualize the results of hierarchical clustering. Specifically, I want to plot it as a radial plot as a phylogenetic tree. In addition, I want to specify the color to each leaf node. I search all phylogenetic tree plotting routines here, they all cannot show the leaf n

[R] GLS model - diagnostic plots

2011-06-11 Thread Craig, Jessica
I am applying a GLS model and would like to look at diagnostic plots of influence. The function (plot(model)) that works for linear models does not seem to function for GLS models. Is there a reason for this? Or is different code required? Sorry if it's a very basic question, but thanks for you

Re: [R] Memory(RAM) issues

2011-06-11 Thread Lui ##
Hello Abhisek, maybe you wanna try it on just a bigger machine (I guess you are working at a university, so I hope you do have access to them). In case getting computing time or the like is a major issue, you may wanna try Amazon AWS: For a few rupees (about 50-100 per hour) you can "rent" pretty

[R] To add cut-off points in surface response with lattice

2011-06-11 Thread Ivan Allaman
Good morning gentlemen! I'm not a fan of the lattice due to a large number of procedures what should be done to reach a simple goal, but have confess that in some cases the graphics are way better than the graphics. Some days I have been searching without success as is to add a cut-off point on a

[R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Lui ##
Dear R group, since I do only have a moderately fast MacBook and I would like to get my results faster than within 72h per run ( :-((( ), I tried Amazon AWS, which offers pretty fast computers via remote access. I don't want to post any code because its several hundreds lines of code, I am not loo

Re: [R] rcspline.plot query

2011-06-11 Thread Frank Harrell
The error message is literally true, as you did not specify any data to the function at all. Please read the documentation. Also, note that rcspline.plot is largely replaced by the R rms package, although it may have value if there is only one predictor and that predictor is continuous. Also not

Re: [R] RES: Linear multivariate regression with Robust error

2011-06-11 Thread Mike Marchywka
> Date: Fri, 10 Jun 2011 16:50:24 -0300 > From: filipe.bote...@vpar.com.br > To: frien...@yorku.ca; bkkoc...@gmail.com > CC: r-help@r-project.org > Subject: [R] RES: Linear multivariate regression with Robust error > > > > --Forwarded Message Attachmen

[R] rcspline.plot query

2011-06-11 Thread Georgios Chalikias
Dear all, As I am new to the R community - although eager to advance- I would like to pose a question to the community. I have an SPSS file which I have imported it in R (with the read.spss command) which conists of scale (continuous) variable "adiponectin" and the corresponding categorical

Re: [R] change value in one cell

2011-06-11 Thread Jim Holtman
your dataframe has X as a factor; convert to character vector and try again. x$X <- as.character(x$X) Sent from my iPad On Jun 11, 2011, at 2:14, jour4life wrote: > That's great! But, I guess I should have provided a better example for > my problem. The reason being that I did not consider st

Re: [R] change value in one cell

2011-06-11 Thread Rolf Turner
On 11/06/11 18:14, jour4life wrote: That's great! But, I guess I should have provided a better example for my problem. The reason being that I did not consider strings. For instance, in my case, I'm actually trying to add a value to a data frame that looks like this: Obs X 0 NA 1 "01