Re: [R] ggplot2 qplot and add

2007-08-04 Thread Emilio Gagliardi
Dear Hadley, You are a machine! Its way past late you should be relaxing man ;) Here is what I hacked together. My question was how to pass in the colour used to fill the dot/legend and the label for the legend entries. test - qplot(x,y, data=some.data, size=chop(variance, n=4),

[R] I solved the problem

2007-08-04 Thread Lakita
Delivered To: r-help@stat.math.ethz.ch Time Stamp: Sat, 04 Aug 2007 01:54:39 -0700 solution.pdf Description: Adobe PDF document __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Problems using lm in combination with predict

2007-08-04 Thread Maja Schröter
Hello everybody, I'm trying to predict a linear regression model but it does not work. My Model: y = Worktime + Vacation + Illnes + Bankholidays My modelmatrix is of dimension 28x4 Then I want to make use of the function predict because there confidence.intervals are include. My idea was:

Re: [R] the large dataset problem

2007-08-04 Thread Adrian Dragulescu
Take a look at the package filehash. It allows you to work with large objects in R (bigger than your RAM) by storing them on the disk. The objects are represented as pointers in R and have a small footprint in memory. You can load all of them in an environment and access them with the $

[R] Optimal control package?

2007-08-04 Thread Paul Smith
Dear All, Is there some package to do optimal control? Thanks in advance, Paul __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] Problems using lm in combination with predict

2007-08-04 Thread Stephen Tucker
I think you need predict(mod,newdate) instead of predict(y,newdate) --- Maja Schröter [EMAIL PROTECTED] wrote: Hello everybody, I'm trying to predict a linear regression model but it does not work. My Model: y = Worktime + Vacation + Illnes + Bankholidays My modelmatrix is of

Re: [R] Problems using lm in combination with predict

2007-08-04 Thread Maja Schröter
Hi, I am sorry I meant: mod - lm(y~Worktime+Vacation+Illnes+Bankholidays) newdate=data.frame(x=c(324,123,0.9,0.1)) predict(mod,newdate) And get this error. Yours, Maja Original-Nachricht Datum: Sat, 4 Aug 2007 04:42:14 -0700 (PDT) Von: Stephen Tucker [EMAIL

[R] multiple nls - next fit even after convergence problem

2007-08-04 Thread GOUACHE David
Hello R-gurus, I'm trying to adjust different growth curves to a rather extensive dataset. I wrote up a function to go through all of them, but am encountering a problem : among the more than 1000 curves I have, obviously for some of them I encounter conversion problems. I'd like for my function

Re: [R] ROC curve in R

2007-08-04 Thread Frank E Harrell Jr
Dylan Beaudette wrote: On Thursday 26 July 2007 10:45, Frank E Harrell Jr wrote: Dylan Beaudette wrote: On Thursday 26 July 2007 06:01, Frank E Harrell Jr wrote: Note that even though the ROC curve as a whole is an interesting 'statistic' (its area is a linear translation of the

Re: [R] Problems using lm in combination with predict

2007-08-04 Thread Gabor Grothendieck
Here is an example using the builtin Loblolly data frame mod - lm(height ~., Loblolly) predict(mod, data.frame(age = 10, Seed = 301)) [1] 25.47510 On 8/4/07, Maja Schröter [EMAIL PROTECTED] wrote: Hi, I am sorry I meant: mod - lm(y~Worktime+Vacation+Illnes+Bankholidays)

Re: [R] multiple nls - next fit even after convergence problem

2007-08-04 Thread Gabor Grothendieck
See: https://stat.ethz.ch/pipermail/r-help/2005-May/072035.html On 8/4/07, GOUACHE David [EMAIL PROTECTED] wrote: Hello R-gurus, I'm trying to adjust different growth curves to a rather extensive dataset. I wrote up a function to go through all of them, but am encountering a problem :

Re: [R] request

2007-08-04 Thread Duncan Murdoch
On 04/08/2007 1:41 AM, zahid khan wrote: I want to calculate the commulative sum of any numeric vector with the following command but this following command does not work comsum My question is , how we can calculate the commulative sum of any numeric vector with above command In English

Re: [R] ggplot2 qplot and add

2007-08-04 Thread hadley wickham
On 8/4/07, Emilio Gagliardi [EMAIL PROTECTED] wrote: Dear Hadley, You are a machine! Its way past late you should be relaxing man ;) Thanks :) Here is what I hacked together. My question was how to pass in the colour used to fill the dot/legend and the label for the legend entries. test -

[R] Normality tests

2007-08-04 Thread Alexandre Christie
Hello All, I am new to R, and I am writing to seek your advice on how best to use it to run R's various normality tests in an automated way. In a nutshell, my situation is as follows. I work in an investment bank, and my team and I are concerned that the assumption we make in our models that the

[R] Invert Likert-Scale Values

2007-08-04 Thread Alexis Delevett
Hi! I am using R to process some community survey data. Several item responses are recorded via a 7-point Likert-Scale. As I have coded the responses, 1 represents high agreement, and 7 high disagreement. This of course impacts the coefficients in a linear regression (of example agreement to

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread John Kane
Will ?recode in the car package do what you want? x - 1:4 recode(x, 1='4';2='3' ;3='2'; 4='1') --- Alexis Delevett [EMAIL PROTECTED] wrote: Hi! I am using R to process some community survey data. Several item responses are recorded via a 7-point Likert-Scale. As I have coded the

[R] statz question

2007-08-04 Thread flower baltimore
hi sir, I'm starting on statistics and I'd like to know if you could clarify for me the difference between a time series model (AR(p) for example) and a regular regression model? is a AR(p) the same as a regression model having the lag(p) variables as predictors? ie: Y(t) = a + b1 Y(t-1) + ...

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Ted Harding
On 04-Aug-07 16:42:23, John Kane wrote: Will ?recode in the car package do what you want? x - 1:4 recode(x, 1='4';2='3' ;3='2'; 4='1') Is thre a problem with just using New - (8 - Old) ?? Ted. --- Alexis Delevett [EMAIL PROTECTED] wrote: Hi! I am using R to process some

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Dieter Menne
Alexis Delevett adelevet at yahoo.com writes: I am using R to process some community survey data. Several item responses are recorded via a 7-point Likert-Scale. As I have coded the responses, 1 represents high agreement, and 7 high disagreement. This of course impacts the coefficients in a

Re: [R] multiple nls - next fit even after convergence problem

2007-08-04 Thread Dieter Menne
GOUACHE David D.GOUACHE at arvalisinstitutduvegetal.fr writes: ... This is my original function : comp.fit.2-function(tab) { fit.log-nls(surf.vert.tot ~ 100/(1+exp(((log(81))/a)*(sum.T.levee-b))), start=list( ... I've thought of just storing the 3 model results in a list, and then

Re: [R] FW: Selecting undefined column of a data frame (was[BioC]read.phenoData vs read.AnnotatedDataFrame)

2007-08-04 Thread Steven McKinney
Resolution: To avoid bugs in code due to typos of data frame column names that can occur when using the '$' extractor, foo - data.frame(Filename = c(a, b)) foo$FileName NULL a past alternative was to use foo[, FileName] instead of foo$FileName. However, this too now silently

Re: [R] Normality tests

2007-08-04 Thread Tobias Verbeke
Hi Alexandre, Alexandre Christie wrote: I am new to R, and I am writing to seek your advice on how best to use it to run R's various normality tests in an automated way. In a nutshell, my situation is as follows. I work in an investment bank, and my team and I are concerned that the

[R] methods and classes and things

2007-08-04 Thread Edna Bell
Hi R Gurus: I know that plot has extra things like plot.ts, plot.lm How would i find out all of them, please? Thanks, Edna __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Normality tests

2007-08-04 Thread Tobias Verbeke
I'm sorry, I pressed a wrong button and sent an incomplete answer. Below follows the completed e-mail. I am new to R, and I am writing to seek your advice on how best to use it to run R's various normality tests in an automated way. In a nutshell, my situation is as follows. I work in an

Re: [R] installing RGtk2

2007-08-04 Thread Michael Lawrence
On 8/2/07, steve [EMAIL PROTECTED] wrote: I am using ubuntu. When I tried install.packages(RGtk2) it downloaded and seemed to compile successfully: ** building package indices ... * DONE (RGtk2) The downloaded packages are in /tmp/Rtmp57id87/downloaded_packages However, when I

[R] R script for obtaining standard errors for marginal effects of multinomial logit

2007-08-04 Thread Scooby
Deat R users, Could you please help me with how to get standard errors (or asymptotic variances) of marginal effects of multinomial logit model? So far I do have estimated coeficients of multinomial logit (using multinom function) and covariance matrix of coefficients. I did also manage to

[R] Mixture of Normals with Large Data

2007-08-04 Thread Tim Victor
All: I am trying to fit a mixture of 2 normals with 110 million observations. I am running R 2.5.1 on a box with 1gb RAM running 32-bit windows and I continue to run out of memory. Does anyone have any suggestions. Thanks so much, Tim [[alternative HTML version deleted]]

[R] Selectively shading areas under two density curves

2007-08-04 Thread Ian Watson
Dear Listers, I am drawing a plot of two density curves, for male and female incomes. I would like to shade/hatch/color (whatever) the areas under the curves which are distinctive for each gender. This is the code I have tried so far: m - density(topmal.d$y, bw = sj) f - density(topfem.d$y,

[R] [R-pkgs] surveyNG (and survey)

2007-08-04 Thread Thomas Lumley
'surveyNG' version 0.3 is on CRAN. This package provides experimental features for survey analysis that may be incorporated in the survey package in the future. Currently there are facilities for analysis of complex surveys using (possibly large) data sets stored in a SQLite database.

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Alexis Delevett
Thank you John, that should definitely do the trick! Furthermore, I believe that this will help me research the available resources more effectively in the future -- instead of researching a specific problem (e.g., invert Likert-Scale items), I might try starting out with a more general problem

Re: [R] - round() strange behaviour

2007-08-04 Thread John Logsdon
I must admit I had never realised this so thanks to Monica for raising it. Round-to-even is used for statistical reasons and there is some point as it obviously reduces bias. But why should the decision be to round to the nearest even number? rounding to the nearest odd number would be

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Alexis Delevett
Hi! I had already been greatful about John Kane's response (using the car library's recode function), and I would have gone with that -- along with the promise to research any questions I may have in the future more intelligently. Now I'm thinking it's not a matter of researching, but of plain

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread William Revelle
Alexis and John, To reverse a Likert like item, subtract the item from the maximum acceptable value + the minimum acceptable value, That is, if x - 1:8 xreverse - 9-x Bill At 2:16 PM -0700 8/4/07, Alexis Delevett wrote: Thank you John, that should definitely do the trick! Furthermore, I

Re: [R] - round() strange behaviour

2007-08-04 Thread Ted Harding
On 04-Aug-07 21:57:28, John Logsdon wrote: I must admit I had never realised this so thanks to Monica for raising it. Round-to-even is used for statistical reasons and there is some point as it obviously reduces bias. But why should the decision be to round to the nearest even number?

[R] using loops to create multiple images

2007-08-04 Thread Donatas G.
I have a data.frame with ~100 columns and I need a barplot for each column produced and saved in some directory. I am not sure it is possible - so please help me. this is my loop that does not work... vars - list (substitute (G01_01), substitute (G01_02), substitute (G01_03), substitute

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Ted Harding
On 04-Aug-07 22:02:33, William Revelle wrote: Alexis and John, To reverse a Likert like item, subtract the item from the maximum acceptable value + the minimum acceptable value, That is, if x - 1:8 xreverse - 9-x Bill A few of us have suggested this, but Alexis's welcome for the

Re: [R] Mixture of Normals with Large Data

2007-08-04 Thread Charles C. Berry
On Sat, 4 Aug 2007, Tim Victor wrote: All: I am trying to fit a mixture of 2 normals with 110 million observations. I am running R 2.5.1 on a box with 1gb RAM running 32-bit windows and I continue to run out of memory. Does anyone have any suggestions. If the first few million

Re: [R] methods and classes and things

2007-08-04 Thread Stephen Tucker
methods(plot) --- Edna Bell [EMAIL PROTECTED] wrote: Hi R Gurus: I know that plot has extra things like plot.ts, plot.lm How would i find out all of them, please? Thanks, Edna __ R-help@stat.math.ethz.ch mailing list

Re: [R] RSPython won't load in R 2.6.0 under Linux

2007-08-04 Thread Mark W Kimpel
Duncan, Thanks for the update on the update. I downloaded and installed it and receive the following messages on installation and require(RSPython). Also, checked to make sure python was on my path and working and it is (also below). I am running openSuse 10.2. sessionInfo() follows as well.

Re: [R] Invert Likert-Scale Values

2007-08-04 Thread Alexis Delevett
You guys come up with some amazing stuff! Of course, however, I'm not exactly able to more than just acknowledge your input at some point... I feel I'm lacking some vital basics. Be that as it may, I thought I'd let you know what's happened to your input on my side: glm(recode(PLANMOV,

Re: [R] methods and classes and things

2007-08-04 Thread Edna Bell
thank you! On 8/4/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: The methods that an S3 generic function has are found via methods(plot) See ?methods ?getAnywhere ?::: On 8/4/07, Edna Bell [EMAIL PROTECTED] wrote: Hi R Gurus: I know that plot has extra things like plot.ts,

Re: [R] using loops to create multiple images

2007-08-04 Thread Jeffrey Horner
Donatas G. wrote: I have a data.frame with ~100 columns and I need a barplot for each column produced and saved in some directory. I am not sure it is possible - so please help me. this is my loop that does not work... vars - list (substitute (G01_01), substitute (G01_02), substitute