Re: [R] How to detect and exclude outliers in R?

2010-01-19 Thread milton ruser
Hi V.S., Did you search first on r-repositories about this issue prior to ask? May be not. RSiteSearch(outliers) bests milton On Tue, Jan 19, 2010 at 1:08 AM, vikrant vikrant.shi...@tcs.com wrote: Suppose I am reading data from a file and the data contains some outliers. I want to know

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread milton ruser
Dear V.K., 1. http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html 2. it is - in general - good someone that write to the list identify your self. It is a polite way of participate of the list. bests milton On Tue, Jan 19, 2010 at 12:15 AM, vikrant

Re: [R] Data import export zipped files from URLs

2010-01-19 Thread Dieter Menne
Velappan Periasamy wrote: I am not able to import zipped files from the following link. How to get thw same in to R?. mydata - read.csv(http://nseindia.com/content/historical/EQUITIES/2010/JAN/cm15JAN2010bhav.csv.zip;) As Brian Ripley noted in

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread Dieter Menne
vikrant wrote: i am unable to open a file which is saved as .xlsx format in R . The file contains approximately 1,50,000 rows. Use odbcConnectExcel2007 in package RODBC to read these data. So I m not able to save it as csv file. I do not understand why you want to use R to convert

Re: [R] Predict polynomial problem

2010-01-19 Thread Barry Rowlingson
On Tue, Jan 19, 2010 at 1:36 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote: Its the environment thing. I think you want something like this:        models[[i]]=lm( bquote( y ~ poly(x,.(i)) ), data=d) Use        terms( mmn[[3]] ) both with and without this change and        ls( env

[R] read.table with special chars

2010-01-19 Thread Marc Noguera
Hi all, I am just starting with R and I have come across a problem which I guess it is easy to solve. I am reading a table with read.table function. This table contains chars which seem to be problematic when reading them such as ' and #, if I remove these characters fromt the table by hand the

Re: [R] column selection for aggregate()

2010-01-19 Thread Petr PIKAL
Hi If I really wanted aggregate all numerics by all non numerics this is how I would do it my.numerics - which(sapply(zeta, is.numeric)) my.factor - which(sapply(zeta, is.factor)) aggregate(zeta[, my.numerics], zeta[, my.factor], mean) Regards Petr r-help-boun...@r-project.org napsal dne

Re: [R] Help using Cast (Text) Version

2010-01-19 Thread Steve Sidney
Hi Hadley Thanks I have downloaded the intro and the material and will work through it once get a chance Thanks for your interest Regards Steve - Original Message - From: hadley wickham h.wick...@gmail.com To: David Winsemius dwinsem...@comcast.net Cc: Steve Sidney

Re: [R] read.table with special chars

2010-01-19 Thread Jim Lemon
On 01/19/2010 07:41 PM, Marc Noguera wrote: Hi all, I am just starting with R and I have come across a problem which I guess it is easy to solve. I am reading a table with read.table function. This table contains chars which seem to be problematic when reading them such as ' and #, if I remove

Re: [R] read.table with special chars

2010-01-19 Thread Marc Noguera
Thanks Jim, I've tried that but still no luck. Some other suggestions? Thanks again Marc Jim Lemon wrote: On 01/19/2010 07:41 PM, Marc Noguera wrote: Hi all, I am just starting with R and I have come across a problem which I guess it is easy to solve. I am reading a table with read.table

Re: [R] Odp: For loops in R

2010-01-19 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 19.01.2010 02:32:45: Hello Petr. For the random values, I wanted to generate a different random number for each element of my velocity matrix. So will this do it? rmat - matrix(runif(1000), 500,2) rmat2 - matrix(runif(1000), 500,2)

Re: [R] exporting text output to pdf

2010-01-19 Thread Dimitri Shvorob
The modified snippet (with the long paragraph truncated) does not produce anything like the requested document. -- View this message in context: http://n4.nabble.com/exporting-text-output-to-pdf-tp837699p1017332.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] exporting text output to pdf

2010-01-19 Thread Dimitri Shvorob
Thanks a lot! -- View this message in context: http://n4.nabble.com/exporting-text-output-to-pdf-tp837699p1017331.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] column selection for aggregate()

2010-01-19 Thread Ivan Calandra
Not really, I tried without select = - c(MEASUREM, SEL_FACET, SEL_MEAS) and indeed the mean was not computed, but it still appeared in the data, which I didn't want. Thanks a lot for your help Ivan Gabor Grothendieck a écrit : It looks ok except you have both specified the wanted factors and

Re: [R] read.table with special chars

2010-01-19 Thread Jim Lemon
On 01/19/2010 08:26 PM, Marc Noguera wrote: Thanks Jim, I've tried that but still no luck. Some other suggestions? The comment.char change should have worked, as others have reported that it did. The single quote problem came up a few weeks ago and the answer then was to remove the single

Re: [R] Rotating pca scores

2010-01-19 Thread Liviu Andronic
Hello On 1/18/10, francesca.ior...@googlemail.com francesca.ior...@googlemail.com wrote: Does anybody know how I can obtain/calculate rotated PC scores with R? You might want to try principal() in package psych, and see if it does what you need. With this function you can use all the rotations

Re: [R] exporting text output to pdf

2010-01-19 Thread baptiste auguie
Hi, You could play with the splitTextGrob() function from the RGraphics package, string - Lorem ipsum dolor sit amet, consectetur adipiscing elit. Quisque leo ipsum, ultricies scelerisque volutpat non, volutpat et nulla. Curabitur consequat ullamcorper tellus id imperdiet. Duis semper malesuada

Re: [R] call R with un expression (String)?

2010-01-19 Thread baptiste auguie
Hi, ?eval seems like a good candidate HTH, baptiste 2010/1/15 Jiiindo jiin...@yahoo.com: Hello all, I want to call R from java. And I have a expression in Java as a String, example : (variable 1 + variable 2)* variable 3 and i want R calculate this expression. How can I do? ex: Java

Re: [R] Help on using WinBUGS on Mac

2010-01-19 Thread Uwe Ligges
On 19.01.2010 03:23, Wayne (Yanwei) Zhang wrote: Dear all, I had trouble in setting up WinBUGS on my Mac, and I'm seeking for some help here. I followed the instruction by Tom Palmer here: http://www.ruudwetzels.com/MacBUGS/winbugsonmacosx.pdf I installed Darwine 1.1.21, and downloaded

[R] apply command

2010-01-19 Thread marco salvini
Can you please help on the issue? I using the apply command on a matrix below the example: Create a vector x =c(5, 3, 2:4, NA, 7, 3, 9, 2, 1, 5) create a matrix of 2 rows by 6 columns b=matrix(x, 2,6) print(b) [,1] [,2] [,3] [,4] [,5] [,6] [1,]524791 [2,]33

Re: [R] apply command

2010-01-19 Thread Tal Galili
Hello Marco What I would do, is use t to transpose the matrix. Why it is that apply switches the matrix, is beyond my knowledge - and I would love to read more informed replies. Tal Contact Details:--- Contact me:

Re: [R] Help ~

2010-01-19 Thread Dieter Menne
-- View this message in context: http://n4.nabble.com/Help-tp1017274p1017414.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] Macaulay Duration for Group

2010-01-19 Thread Madhavi Bhave
Dear R helpers   I have following csv file which is an input   id   par_value    coupon_rate frequency_coupon   tenure    ytm   1    1000 10  1     5  12   # Here frequency_coupon is coded s.t. 0 means Daily compounding, 1

[R] restricted permutations in permtest()?

2010-01-19 Thread Kay Cecil Cichini
Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan documentation p.24: library(vegan) data(varespec) ## Bray-Curtis

Re: [R] Help ~

2010-01-19 Thread Dieter Menne
Fan Dan wrote: Generate a clustered pattern in [0; 1]2 as follows: (a) Generate n, say 20, independent cluster centers (which can be called parents) that are distributed i.i.d. uniformly in the unit square; (b) then m,say 100, daughters are assigned i.i.d. uniformly to these parents and

[R] PCA scores and loadings

2010-01-19 Thread Paul Dennis
Dear R users group I have performed PCA using the function rda in vegan and then used plot(pcaobject). I have a couple of questions: 1) The default plot shows the individual sites (black) and the variables (red). What I want however is a plot showing the mean of site groups with

Re: [R] How to detect and exclude outliers in R?

2010-01-19 Thread Eik Vettorazzi
fortune(outlier) vikrant schrieb: Suppose I am reading data from a file and the data contains some outliers. I want to know if it is possible in R to automatically detect outliers in a dataset and remove them -- Eik Vettorazzi Institut für Medizinische Biometrie und Epidemiologie

Re: [R] apply command

2010-01-19 Thread baptiste auguie
Hi, I think you could use iapply (search the archives) or the plyr package to save you from transposing the result. HTH, baptiste 2010/1/19 marco salvini marco.salv...@gmail.com: Can you please help on the issue? I using the apply command on a matrix below the example: Create a vector x

Re: [R] how to install spatstat

2010-01-19 Thread Uwe Ligges
On 19.01.2010 03:42, Rolf Turner wrote: On 19/01/2010, at 11:24 AM, Senlin Liang wrote: Hi, I tried to install spatstat using: install.packages(spatstat, dependencies = TRUE, lib = ./R/i486-pc-linux-gnu-library/2.8) but got the following msg: --- Please select a CRAN mirror for use in

Re: [R] ggplot2 histogramm

2010-01-19 Thread Christian Schulz
Hi Dennis, the BaseTheme() are some settings scale the text and title etc. , the plot should work without that. At present i'm experiment with the huge amount of possibilities and try to understand and maybe doing beginner mistakes. Thanks for the hint with scale_fill_gradient and

[R] FW: help with time Series regression please

2010-01-19 Thread snowcat wong
Dear all, I am having difficulty to built a model of quarter sales of spirits data, and deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not appeared right at the end. The data and script is enclosed with this email. I am using the harmonic regression model to

[R] help with time Series regression please

2010-01-19 Thread snowcat wong
Dear all, I am having difficulty to built a model of quarter sales of spirits data, and deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not appeared right at the end. The data and script is enclosed with this email. I am using the harmonic regression model to

Re: [R] apply command

2010-01-19 Thread Linlin Yan
I guess that the matrix dimension changed because matrix in R are filled by columns. Since you try: apply(b, 1, function(y) sort(y, na.last=F)) The second parameter make it scan matrix b row by row but store result by columns, which make the result be a matrix transposed. If you try: apply(b, 2,

[R] Working with text data/text operators

2010-01-19 Thread Mihai.Mirauta
Hello, Could someone tell me, how can I select from a dataframe only those columns whose names contain a certain text? For example, if the column names are Bond1.Creditclass,Bond1.Price,Bond2.Creditclass,Bond2.Price, how do I select only the columns corresponding to Bond1? Thanks a lot,

Re: [R] Working with text data/text operators

2010-01-19 Thread Romain Francois
On 01/19/2010 01:22 PM, mihai.mira...@bafin.de wrote: Hello, Could someone tell me, how can I select from a dataframe only those columns whose names contain a certain text? For example, if the column names are Bond1.Creditclass,Bond1.Price,Bond2.Creditclass,Bond2.Price, how do I select

Re: [R] restricted permutations in permtest()?

2010-01-19 Thread Gavin Simpson
On Tue, 2010-01-19 at 11:39 +0100, Kay Cecil Cichini wrote: Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan

Re: [R] apply command

2010-01-19 Thread Tal Galili
Thank you for the answer Linlin, I am wondering, is there a way to change it so that R will fill matrix's by rows ? Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew)

Re: [R] Data import export zipped files from URLs

2010-01-19 Thread Duncan Temple Lang
Dieter Menne wrote: Velappan Periasamy wrote: I am not able to import zipped files from the following link. How to get thw same in to R?. mydata - read.csv(http://nseindia.com/content/historical/EQUITIES/2010/JAN/cm15JAN2010bhav.csv.zip;) As Brian Ripley noted in

Re: [R] Data import export zipped files from URLs

2010-01-19 Thread Gabor Grothendieck
If you need an example of this look at the yacasInstall function in this file: http://ryacas.googlecode.com/svn/trunk/R/yacasInstall.R from the Ryacas package. It downloads, unzips and installs yacas and associated files for Windows users. On Tue, Jan 19, 2010 at 3:10 AM, Dieter Menne

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread Gabor Grothendieck
See http://wiki.r-project.org/rwiki/doku.php?id=tips:data-io:ms_windows Most of the methods there work with xls files only but a few work with xlsx. On Tue, Jan 19, 2010 at 12:15 AM, vikrant vikrant.shi...@tcs.com wrote: i am unable to open a file which is saved as .xlsx format in R . The file

[R] Macualay Duration code in a Functional Form - Please Help

2010-01-19 Thread Madhavi Bhave
# I have written this code in Notepad++ and copied here. ## ONS - PPA    Duration = function(par_value, coupon_rate, freq_coupon, tenure, ytm) { macaulay_duration  =   NULL modified_duration    =   NULL freq_coupon_new    =   NULL if(freq_coupon = 0) {     freq_coupon_new = 365 }

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread Erich Neuwirth
RExcel mention in the link below can transfer data from Excel 2007 to R. But you have to be on Windows, and you probably have to have Excel 2007. it might work if you have Excel 2003 with the (free) compatibility package installed. On Jan 19, 2010, at 2:09 PM, Gabor Grothendieck wrote: See

Re: [R] restricted permutations in permutest()?

2010-01-19 Thread Kay Cichini
hello gavin, you are right, i didn't get into the documentation to deep and i'm also a beginner, that's why i'm just about to get into the logical part of the syntax. now, the output from perm.disp() says: #No. of permutations: 999 #Permutation type: free #Permutations are unstratified

Re: [R] bootstrapping

2010-01-19 Thread aaron.foley
Thanks all for your help! Aaron Date: Sat, 16 Jan 2010 00:04:59 +0100 From: stephan.kola...@gmx.de To: aaron.fo...@students.tamuk.edu CC: r-help@r-project.org Subject: Re: [R] bootstrapping Hi Aaron, try the argument statistic=mean. Then boot() will give you the mean turn angle

Re: [R] Predict polynomial problem

2010-01-19 Thread Gavin Simpson
On Tue, 2010-01-19 at 08:27 +, Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 1:36 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote: Its the environment thing. I think you want something like this: models[[i]]=lm( bquote( y ~ poly(x,.(i)) ), data=d) Use

[R] Memory usage in read.csv()

2010-01-19 Thread nabble . 30 . miller_2555
I'm sure this has gotten some attention before, but I have two CSV files generated from vmstat and free that are roughly 6-8 Mb (about 80,000 lines) each. When I try to use read.csv(), R allocates all available memory (about 4.9 Gb) when loading the files, which is over 300 times the size of the

Re: [R] Cannot load RInterpreter (SJava)

2010-01-19 Thread Martin Morgan
Jiiindo wrote: Hello, I write a test for call R function from Java by Eclipse. When i run it, raise a error: Loading RInterpreter library Exception in thread main java.lang.UnsatisfiedLinkError: no RInterpreter in java.library.path at

[R] Server hanging despite efforts to correct memory limits

2010-01-19 Thread Nathan Stephens
My group is working with datasets between 100 Mb and 1 GB in size, using multiple log ins. From the documentation, it appears that vsize is limited to 2^30-1, which tends to prove too restrictive for our use. When we drop that restriction (set vsize = NA) we end up hanging the server, which

[R] question on plot in R with mac

2010-01-19 Thread khazaei
-- Hello all My computer is MacBook and I want to draw a plot in R, for example for x - c(1,3,6,9,12) y - c(1.5,2,7,8,15) I use this command plot(x,y). but it dosn't work. Could you please help me? thank you khazaei

Re: [R] Symmetric Matrix classes

2010-01-19 Thread Martin Maechler
Scanning for 'Matrix' in old R-help e-mail, I found GA == Gad Abraham gabra...@csse.unimelb.edu.au on Fri, 27 Nov 2009 13:45:00 +1100 writes: GA Hi, GA I'd like to store large covariance matrices using Matrix classes. GA dsyMatrix seems like the right one, but I want to

Re: [R] Server hanging despite efforts to correct memory limits

2010-01-19 Thread W
I am running into a separate, but related issue. On Linux, one may impose memory limits via the --max-vsize, --max-nsize, and --max-ppsize arguments upon starting R. I do not know if similar arguments are available on Windows. HTH -- View this message in context:

Re: [R] restricted permutations in permutest()?

2010-01-19 Thread Gavin Simpson
On Tue, 2010-01-19 at 06:19 -0800, Kay Cichini wrote: hello gavin, you are right, i didn't get into the documentation to deep and i'm also a beginner, that's why i'm just about to get into the logical part of the syntax. now, the output from perm.disp() says: As I said, you *can't* do

Re: [R] (2nd part) variable name substitution

2010-01-19 Thread Ivan Calandra
Hi again! I feel like I cannot do anything by myself but I would now like to plot for all numeric variables I have (14 of them). I wanted to add a loop then. The code is: -- #defines the function for the plots (as written by Duncan Murdoch) twoplots - function(x, y) { ylab -

Re: [R] (2nd part) variable name substitution

2010-01-19 Thread Carlos Ortega
Hello, You can loop in the subset you need by storing in a variable and looping on that variable with indexes: seq.dat-c(seq(7,10,1), seq(12,17,1)) for( i in 1:length(seq.dat) ) { j-seq.dat[i] with(ssfa, twoplots(TO_POS, ssfa[[j]])) } Regards, Carlos. On Tue, Jan 19, 2010 at 4:53 PM, Ivan

[R] Remove term from formula for predict.lm

2010-01-19 Thread Werner W.
Hi, probably just a quick question: can I somehow change the formula used with predict? E.g., the regression was run on y ~ u + v + w but for the prediction the term v should be removed from the formula contained in the regression object and only y ~ u + w be used. I could use model.matrix

Re: [R] Server hanging despite efforts to correct memory limits

2010-01-19 Thread Thomas Lumley
You should be able to set limits on memory use for a process in the operating system, eg with limits or ulimits under Unix-alike shells. -thomas On Tue, 19 Jan 2010, Nathan Stephens wrote: My group is working with datasets between 100 Mb and 1 GB in size, using multiple log ins. From

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread Henrique Dallazuanna
Try this; mod1 - lm(y ~ u + v + w, data = d) update(mod1, . ~ . -v) On Tue, Jan 19, 2010 at 2:10 PM, Werner W. pensterfuz...@yahoo.de wrote: Hi, probably just a quick question: can I somehow change the formula used with predict? E.g., the regression was run on y ~ u + v + w but for the

Re: [R] (2nd part) variable name substitution

2010-01-19 Thread Ivan Calandra
Thank you for your answer, I got the second part! Ivan Le 1/19/2010 17:03, Carlos Ortega a écrit : Hello, You can loop in the subset you need by storing in a variable and looping on that variable with indexes: seq.dat-c(seq(7,10,1), seq(12,17,1)) for( i in 1:length(seq.dat) ) {

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread Gabor Grothendieck
This recomputes the lm but if that is ok then: mod - lm(y1 ~ x1 + x2 + x3 + x4, anscombe) mod$call$formula - update(as.formula(mod$call$formula), ~ . - x1 - x2) predict(eval(mod$call), list(x3 = 1, x4 = 1)) On Tue, Jan 19, 2010 at 11:10 AM, Werner W. pensterfuz...@yahoo.de wrote: Hi,

Re: [R] question on plot in R with mac

2010-01-19 Thread David Winsemius
On Jan 19, 2010, at 10:21 AM, khaz...@ceremade.dauphine.fr wrote: -- Hello all My computer is MacBook and I want to draw a plot in R, for example for x - c(1,3,6,9,12) y - c(1.5,2,7,8,15) I use this command plot(x,y).

[R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Gavin Simpson
Dear List, A student in the Department where I work would like to produce a graphic similar to this one: http://image.guardian.co.uk/sys-files/Guardian/documents/2009/09/16/Public_spending_160909.pdf Does anyone know if the figure in the pdf can be generated in a specific software application

Re: [R] Model

2010-01-19 Thread Jim Silverton
Hello, I have a barchart. The y-axis represents counts and thex-axis is divided into 10 equal intervals ranging fronm 0 to 0.1, 0.1 to 0.2, ..0.9 to 1.0. Is there a way to model the counts in R? thanks, Jim [[alternative HTML version deleted]]

Re: [R] (2nd part) variable name substitution

2010-01-19 Thread Carlos Ortega
OK. For the names of the variables you can include this code in the loop (variable nv): seq.dat-c(seq(7,10,1), seq(12,17,1)) for( i in 1:length(seq.dat) ) { j-seq.dat[i] nv-names(ssfa)[j] with( ssfa, twoplots(TO_POS, ssfa[[j]], nv) ) } And this modification in the function (nm): #defines

Re: [R] Data import export zipped files from URLs

2010-01-19 Thread Velappan Periasamy
How to unzip this file?. mydata - unzip(http://nseindia.com/content/historical/EQUITIES/2010/JAN/cm15JAN2010bhav.csv.zip;) Warning message: In unzip(http://nseindia.com/content/historical/EQUITIES/2010/JAN/cm15JAN2010bhav.csv.zip;) : error 1 in extracting from zip file

Re: [R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Ista Zahn
You might be able to do that with Rgraphviz or another R package, but if I was doing it I would probably use PGF/tikZ. The homepage is here: http://pgf.sourceforge.net/ -Ista On Tue, Jan 19, 2010 at 4:27 PM, Gavin Simpson gavin.simp...@ucl.ac.uk wrote: Dear List, A student in the Department

Re: [R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Carlos Ortega
MS PowerPoint (version 2007 or beta 2010) although difficult for so dense graphic. Prefearable: MindManager although is $$. Use the Trial. Regards, Carlos. On Tue, Jan 19, 2010 at 5:27 PM, Gavin Simpson gavin.simp...@ucl.ac.ukwrote: Dear List, A student in the Department where I work would

[R] splitting a factor in an analysis of deviance table (negative binomial model)

2010-01-19 Thread Rafael Moral
Dears useRs, I have 2 factors, (for the sake of explanation - A and B), with 4 levels each. I've already fitted a negative binomial generalized linear model to my data, and now I need to split the factors in two distinct analysis of deviance table:  - A within B1, A within B2, A within B3 and A

Re: [R] (2nd part) variable name substitution

2010-01-19 Thread Ivan Calandra
Super, thanks a lot!! I didn't think about using names() Ivan Le 1/19/2010 17:35, Carlos Ortega a écrit : OK. For the names of the variables you can include this code in the loop (variable nv): seq.dat-c(seq(7,10,1), seq(12,17,1)) for( i in 1:length(seq.dat) ) { j-seq.dat[i]

[R] problem with the precision of numbers

2010-01-19 Thread kayj
Hi All, I was wodering if it is possible to increase the precision using R. I ran the script below in R and MAPLE and I got different results when k is large. Any idea how to fix this problem? thanks for your help for (k in 0:2000){ s=0 for(i in 0:k){ s=s+((-1)^i)*3456*(1+i*1/2000)^3000 }

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread Walmes Zeviani
Werner, You could set 0 to that regressors you don't want to consider for prediction. da - expand.grid(A=1:20, B=rnorm(20, 4, 0.2), C=10^seq(1,2,l=20)) da$y - rnorm(da$A, 0, 0.3) m0 - lm(y~A+B+C, data=da) new - da new$C - 0 predict(m0)[1:5] predict(m0, newdata=new)[1:5] At your disposal.

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread werner w
Thanks Gabor and Henrique! Sorry for the imprecise question. I want predict() to use the coefficients estimated by the original regression but to exclude terms from the prediction formula. If I originally estimated y ~ x1 + x2 and got coefficients b0, b1, b2, I would like to remove x2 and

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread David Winsemius
On Jan 19, 2010, at 12:05 PM, werner w wrote: Thanks Gabor and Henrique! Sorry for the imprecise question. I want predict() to use the coefficients estimated by the original regression but to exclude terms from the prediction formula. If I originally estimated y ~ x1 + x2 and got

Re: [R] Data import export zipped files from URLs

2010-01-19 Thread Henrique Dallazuanna
Try this: f - tempfile() download.file(http://nseindia.com/content/historical/EQUITIES/2010/JAN/cm15JAN2010bhav.csv.zip;, f) myData - read.csv(unzip(f)) On Tue, Jan 19, 2010 at 2:56 PM, Velappan Periasamy veepsi...@gmail.com wrote: How to unzip this file?. mydata -

Re: [R] Model

2010-01-19 Thread Walmes Zeviani
Jim, Did you read the posting guide? Did you do a google search, for example, with terms like [R] generalized linear models, [R] count models, [R] poisson regression? I think you should do. Walmes. - ..oooO

Re: [R] Predict polynomial problem

2010-01-19 Thread Charles C. Berry
On Tue, 19 Jan 2010, Charles C. Berry wrote: and the values in those places are different: On Tue, 19 Jan 2010, Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 1:36 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote: Its the environment thing. I think you want something like this: ?

Re: [R] Predict polynomial problem

2010-01-19 Thread Charles C. Berry
and the values in those places are different: On Tue, 19 Jan 2010, Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 1:36 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote: Its the environment thing. I think you want something like this:        models[[i]]=lm( bquote( y ~ poly(x,.(i)) ),

Re: [R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Gabor Grothendieck
See the Rgraphviz package in bioconductor. On Tue, Jan 19, 2010 at 11:27 AM, Gavin Simpson gavin.simp...@ucl.ac.uk wrote: Dear List, A student in the Department where I work would like to produce a graphic similar to this one:

Re: [R] Predict polynomial problem

2010-01-19 Thread Barry Rowlingson
On Tue, Jan 19, 2010 at 5:37 PM, Charles C. Berry cbe...@tajo.ucsd.eduwrote: Note: i - 20 bquote(y ~ poly(x,.(i))) y ~ poly(x, 20) I see it now. bquote(y~poly(x,.(i))) gets it's 'i' there and then, sticks it in the returned expression as the value '20', so any further evaluations get

Re: [R] A model-building strategy in mixed-effects modelling

2010-01-19 Thread Ben Bolker
Stats Wolf stats.wolf at gmail.com writes: Dear all, Consider a completely randomized block design (let's use data(Oats) irrespoctive of the split-plot design it was arranged in). Look: library(nlme) fit - lme(yield ~ nitro, Oats, random = ~1|Block, method=ML) fit2 - lm(yield ~ nitro

[R] Sampling theory

2010-01-19 Thread Christian Hennig
Hi there, are there any R-packages for computations required in sampling theury (such as confidence intervals under random, stratified, cluster sampling; I'd be partoculary interested in confidence intervals for the population variance, which is difficult enough to find even in books)?

Re: [R] problem with the precision of numbers

2010-01-19 Thread Ben Bolker
kayj kjaja27 at yahoo.com writes: Hi All, I was wodering if it is possible to increase the precision using R. I ran the script below in R and MAPLE and I got different results when k is large. Any idea how to fix this problem? thanks for your help for (k in 0:2000){ s=0 for(i in

Re: [R] Remove term from formula for predict.lm

2010-01-19 Thread Gabor Grothendieck
If you want to keep the same coefficients but ignore certain ones just use 0 for the ones you don't want: mod - lm(Sepal.Length ~ Petal.Length + Petal.Width, iris) predict(mod, list(Petal.Length = 3, Petal.Width = 0)) Regarding the problem with the example, the example data was not the best for

Re: [R] Predict polynomial problem

2010-01-19 Thread Thomas Lumley
On Tue, 19 Jan 2010, Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 5:37 PM, Charles C. Berry cbe...@tajo.ucsd.eduwrote: Note: i - 20 bquote(y ~ poly(x,.(i))) y ~ poly(x, 20) I see it now. bquote(y~poly(x,.(i))) gets it's 'i' there and then, sticks it in the returned expression

Re: [R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Rainer M Krug
On Tue, Jan 19, 2010 at 6:56 PM, Ista Zahn istaz...@gmail.com wrote: You might be able to do that with Rgraphviz or another R package, but if I was doing it I would probably use PGF/tikZ. The homepage is here: http://pgf.sourceforge.net/ I second that - gives you really good results.

Re: [R] OT: Software for specific visualisation of data...ideas?

2010-01-19 Thread Steve_Friedman
Whats wrong with Power Point or anyone of its equivalents? Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] Predict polynomial problem

2010-01-19 Thread Peter Ehlers
Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 5:37 PM, Charles C. Berry cbe...@tajo.ucsd.eduwrote: Note: i - 20 bquote(y ~ poly(x,.(i))) y ~ poly(x, 20) I see it now. bquote(y~poly(x,.(i))) gets it's 'i' there and then, sticks it in the returned expression as the value '20', so

Re: [R] Sampling theory

2010-01-19 Thread Thomas Lumley
On Tue, 19 Jan 2010, Christian Hennig wrote: are there any R-packages for computations required in sampling theury (such as confidence intervals under random, stratified, cluster sampling; I'd be partoculary interested in confidence intervals for the population variance, which is difficult

Re: [R] Working with text data/text operators

2010-01-19 Thread Dieter Menne
Mihai.Mirauta wrote: Could someone tell me, how can I select from a dataframe only those columns whose names contain a certain text? For example, if the column names are Bond1.Creditclass,Bond1.Price,Bond2.Creditclass,Bond2.Price, how do I select only the columns corresponding to

Re: [R] problem with the precision of numbers

2010-01-19 Thread Ted Harding
On 19-Jan-10 17:55:43, Ben Bolker wrote: kayj kjaja27 at yahoo.com writes: Hi All, I was wodering if it is possible to increase the precision using R. I ran the script below in R and MAPLE and I got different results when k is large. Any idea how to fix this problem? thanks for your help

Re: [R] problem with the precision of numbers

2010-01-19 Thread Gabor Grothendieck
On Tue, Jan 19, 2010 at 1:41 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: On 19-Jan-10 17:55:43, Ben Bolker wrote: kayj kjaja27 at yahoo.com writes: Hi All, I was wodering if it is possible to increase the precision using R. I ran the script below in R and MAPLE and I got different

Re: [R] Help ~

2010-01-19 Thread Rolf Turner
On 19/01/2010, at 11:40 PM, Dieter Menne wrote: SNIP next time you post homework here, please make sure that you modify the language of the task a bit so that the discrepancy between the task and your helplessness is less evident. This sounds like a fortune to me! How about it

Re: [R] Memory usage in read.csv()

2010-01-19 Thread jim holtman
I read vmstat data in just fine without any problems.  Here is an example of how I do it: VMstat - read.table('vmstat.txt', header=TRUE, as.is=TRUE) vmstat.txt looks like this: date time r b w swap free re mf pi po fr de sr intr syscalls cs user sys id 07/27/05 00:13:06 0 0 0 27755440 13051648

Re: [R] Predict polynomial problem

2010-01-19 Thread Peter Ehlers
Charles C. Berry wrote: On Tue, 19 Jan 2010, Charles C. Berry wrote: and the values in those places are different: On Tue, 19 Jan 2010, Barry Rowlingson wrote: On Tue, Jan 19, 2010 at 1:36 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote: Its the environment thing. I think you

[R] Error compiling R 2.10.1 on AIX

2010-01-19 Thread Mike Waldron
I'm trying to compile R 2.10.1 on AIX 5.3, and am getting the following error: Error in read.dcf(file = descfile) : Line starting 'Package: tools ...' is malformed! Calls: makeLazyLoading ... code2LazyLoadDB - loadNamespace - parseNamespaceFile - read.dcf Execution halted make[3]: *** [all]

Re: [R] restricted permutations in permutest()?

2010-01-19 Thread Kay Cichini
i'll post in r-forge vegan help forum and appreciate your help very much. greetings, kay -- View this message in context: http://n4.nabble.com/restricted-permutations-in-permtest-tp1017422p1017865.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Using the output of strsplit

2010-01-19 Thread James Rome
This suggestion does not work. x seems to have twice the number of entries as spl. x - do.call(rbind, spl) names(x) - c('Date', 'quarter') x$Date - as.Date(x$Date) Error in x$Date : $ operator is invalid for atomic vectors x$quarter - as.numeric(x$quarter) Error in x$quarter : $ operator is

Re: [R] Using the output of strsplit

2010-01-19 Thread jim holtman
'x' is a matrix and not a dataframe. You should be doing colnames(x) - c(Date, quarter) x[,Date] - as.Date(x[,Date]) It would help if you took a look at the structure you were using to understand how to access. 'names' applied to a vector would give you the output for 13000 more entries.

Re: [R] problem with the precision of numbers

2010-01-19 Thread Ted Harding
On 19-Jan-10 18:48:47, Gabor Grothendieck wrote: On Tue, Jan 19, 2010 at 1:41 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: On 19-Jan-10 17:55:43, Ben Bolker wrote: kayj kjaja27 at yahoo.com writes: Hi All, I was wodering if it is possible to increase the precision using R. I ran the

Re: [R] Memory usage in read.csv()

2010-01-19 Thread Gabor Grothendieck
You could also try read.csv.sql in sqldf. See examples on sqldf home page: http://code.google.com/p/sqldf/#Example_13._read.csv.sql_and_read.csv2.sql On Tue, Jan 19, 2010 at 9:25 AM, nabble.30.miller_2...@spamgourmet.com wrote: I'm sure this has gotten some attention before, but I have two

[R] Number of download.

2010-01-19 Thread Christophe Genolini
Hi the list Is there a way to know how many times an R package (on CRAN) has been download ? Christophe __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

  1   2   >