Re: [R] Latex no where to be seen

2010-08-20 Thread Barry Rowlingson
On Fri, Aug 20, 2010 at 5:18 AM, Donald Paul Winston satchwins...@yahoo.com wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But ?latex() gives me an error. The package manager does not have it. The package installer can't find it. Where is

[R] How to make 'step' faster?

2010-08-20 Thread sambit rath
Dear all, I am fairly new to R. I would like to perform a step-wise logit regression aiming to select a model on the basis of AIC. I am using some large datasets (up to a million rows and 97 variables). It is taking the 'step' function just too long to complete a single routine. Now, I have

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread r.ookie
From my experience, *.pdf is the best of the three file formats. On Aug 19, 2010, at 8:32 PM, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph that I generate from R code and save it into word document. Which format is better? pdf, jpeg or tiff? Thank you.

Re: [R] Latex no where to be seen

2010-08-20 Thread Joshua Wiley
On Thu, Aug 19, 2010 at 9:18 PM, Donald Paul Winston satchwins...@yahoo.com wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But ?latex() gives me an error. The package Just as a side note, using ?foo() will always return an error. If there

Re: [R] Rserve (Anyone?)

2010-08-20 Thread Romain Francois
Questions about Rserve typically go to this mailing list: http://mailman.rz.uni-augsburg.de/mailman/listinfo/stats-rosuda-devel Le 20/08/10 06:25, Donald Paul Winston a écrit : REXP has an asBytes() method. Will this capture the output of an R plot function if a proper graphics device is used?

Re: [R] Latex no where to be seen

2010-08-20 Thread Tim Gruene
Dear (person who did not put her name under her message unless it's --), maybe whoever told you referred to the system latex command. You could also issue '??tex' at the R-prompt. When I then search for 'TeX' I quickly find 'tools::texi2dvi Compile LaTeX Files' Is that what you are

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Tim Gruene
Dear , PDF format has the advantage that it is resolution independent (unless it contains an embedded bitmap, but I guess for a graph from R it won't), so that's one great advantage of pdf. However, as far as I know, MS product are a little ignorant of PostScript and PDF, and if that's the case

Re: [R] Latex no where to be seen

2010-08-20 Thread Joshua Wiley
On Fri, Aug 20, 2010 at 12:25 AM, Philippe Grosjean phgrosj...@sciviews.org wrote: On 20/08/10 09:11, Joshua Wiley wrote: On Thu, Aug 19, 2010 at 9:18 PM, Donald Paul Winston satchwins...@yahoo.com  wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Cedric Laczny
Hi Rosalina, I created simple barplots with R and saved them as jpeg oder png. I wanted to use them in a LaTex-beamer presentation. The results were simply awful. They looked like if I had copied them out of some other document, thus losing most of the resolution. However, saving the plot as

[R] How to use FitAR with independent variables (like xreg for arima)?

2010-08-20 Thread Steffen Wagner
Hello, I would like to fit a linear model with autoregressive terms and independent variables. Since the largest AR-Term I am interested in is lag.max=336, the arima-function (package stats) does not work (maximum for arima: lag.max=100). The FitAR function (package FitAR) does allow model

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Duncan Murdoch
On 20/08/2010 1:58 AM, Jeff Newmiller wrote: I use Rgui and copy/paste special windows metafile. You can also use win.metafile() as the graphics device and write directly to a file. I haven't used Word in years, but this used to be the recommended way to include graphics there. Duncan

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Joshua Wiley
On Fri, Aug 20, 2010 at 1:07 AM, Cedric Laczny cedric.lac...@gmx.de wrote: Hi Rosalina, I created simple barplots with R and saved them as jpeg oder png. I wanted to use them in a LaTex-beamer presentation. The results were simply awful. They looked like if I had copied them out of some other

Re: [R] How to make 'step' faster?

2010-08-20 Thread Gavin Simpson
On Fri, 2010-08-20 at 12:25 +0530, sambit rath wrote: Dear all, I am fairly new to R. I would like to perform a step-wise logit regression aiming to select a model on the basis of AIC. I am using some large datasets (up to a million rows and 97 variables). It is taking the 'step' function

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Gavin Simpson
On Thu, 2010-08-19 at 20:32 -0700, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph that I generate from R code and save it into word document. Whichformat is better? pdf, jpeg or tiff? Not that I have used word for a while myself, but when I did EPS files were my

Re: [R] How to make 'step' faster?

2010-08-20 Thread sambit rath
Thank you Gavin! I am aware of the lasso regularization routine. But, in this case, my brief was to perform a stepwise AIC procedure. I guess, subsetting the data and cross validating it over the rest of the data is the only answer. sambit On 20 August 2010 14:43, Gavin Simpson

[R] Seeking advice

2010-08-20 Thread Stephen Liu
Hi folks, Has any folks read before; Basic Statistics Rand R.Wilcox http://www.ebooksx.com/Basic-Statistics-Understanding-Conventional-Methods-and-Modern-Insights_303718.html It is good for a layman in statistics? TIA B.R. Stephen L __

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Ted Harding
On 20-Aug-10 09:24:17, Gavin Simpson wrote: On Thu, 2010-08-19 at 20:32 -0700, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph that I generate from R code and save it into word document. Whichformat is better? pdf, jpeg or tiff? Not that I have used word for a

[R] Deviance Residuals

2010-08-20 Thread Iasonas Lamprianou
Dear all, I am running a logistic regression and this is the output: glm(formula = educationUniv ~ brncntr, family = binomial) Deviance Residuals: Min 1Q Median 3Q Max # αυτά είναι τα υπόλοιπα -0.8825 -0.7684 -0.7684 1.5044 1.6516 Coefficients:

Re: [R] how to set chart output size in rgl (surface3d)?

2010-08-20 Thread Karl Brand
Hi Mandy and esteemed R users, I'm searching for the functions you used to produce the shaded bounding box (and axis for that matter) of the plot image* attached to your orginal post. Would you or anyone else be patient enough to post an example of such code or point me to the relevant

Re: [R] Latex no where to be seen

2010-08-20 Thread Philippe Grosjean
On 20/08/10 09:11, Joshua Wiley wrote: On Thu, Aug 19, 2010 at 9:18 PM, Donald Paul Winston satchwins...@yahoo.com wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But ?latex() gives me an error. The package Just as a side note, using

[R] Determining the length of unique items in a vector

2010-08-20 Thread Ron Michael
Dear all, let suppose I have following vector:   dat1 - c(rep(asd, 5), rep(xyz, 12), rep(erd, 17)) dat1 - dat1[sample(1:length(dat1), length(dat1), replace=F)] dat1  [1] erd xyz erd asd asd erd xyz asd erd erd asd xyz erd asd xyz xyz erd xyz erd [20] erd erd xyz xyz erd erd erd erd xyz xyz xyz

Re: [R] Pass By Value Questions

2010-08-20 Thread Matthew Dowle
To: r-help Cc: Jeff, Matt, Duncan, Hadley [ using Nabble to cc ] Jeff, Matt, How about the 'refdata' class in package ref. Also, Hadley's immutable data.frame in plyr 1.1. Both allow you to refer to subsets of a data.frame or matrix by reference I believe, if I understand correctly.

Re: [R] Determining the length of unique items in a vector

2010-08-20 Thread Barry Rowlingson
On Fri, Aug 20, 2010 at 11:01 AM, Ron Michael ron_michae...@yahoo.com wrote: Dear all, let suppose I have following vector: dat1 - c(rep(asd, 5), rep(xyz, 12), rep(erd, 17)) dat1 - dat1[sample(1:length(dat1), length(dat1), replace=F)] dat1  [1] erd xyz erd asd asd erd xyz asd erd erd asd xyz

Re: [R] Determining the length of unique items in a vector

2010-08-20 Thread Alain Guillet
Hi, You can try sapply(levels(as.factor(dat1)),nchar) Alain On 20-Aug-10 12:01, Ron Michael wrote: Dear all, let suppose I have following vector: dat1- c(rep(asd, 5), rep(xyz, 12), rep(erd, 17)) dat1- dat1[sample(1:length(dat1), length(dat1), replace=F)] dat1 [1] erd xyz erd asd asd erd

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Gavin Simpson
On Fri, 2010-08-20 at 10:54 +0100, ted.hard...@manchester.ac.uk wrote: On 20-Aug-10 09:24:17, Gavin Simpson wrote: On Thu, 2010-08-19 at 20:32 -0700, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph that I generate from R code and save it into word document.

Re: [R] Determining the length of unique items in a vector

2010-08-20 Thread Ivan Calandra
Hi, I thought you were looking for table(), but the other answers gave you something really different; I might have wrongly understood your question. HTH, Ivan Le 8/20/2010 12:32, Alain Guillet a écrit : Hi, You can try sapply(levels(as.factor(dat1)),nchar) Alain On 20-Aug-10 12:01,

[R] paired samples, matching rows, merge()

2010-08-20 Thread Cecilia Carmo
Hi everyone! I'm matching two samples to create one sample that have pairs of observations equal for the k1 variable. Merge() doesn't work because I dont't want to recycle the values. x - data.frame(k1=c(1,1,2,3,3,5), k2=c(20,21,22,23,24,25)) x y - data.frame(k1=c(1,1,2,2,3,4,5,5),

Re: [R] Determining the length of unique items in a vector

2010-08-20 Thread Alain Guillet
Hi Ivan, Now I read the other answers I also think I misunderstood the question... The good thing is that one of use certainly gave the right answer to the question ;-) Alain On 20-Aug-10 12:37, Ivan Calandra wrote: Hi, I thought you were looking for table(), but the other answers

[R] error heatmap and stack overflow

2010-08-20 Thread michy
Hello, Im trying to create a heatmap with a dataset (38 x 15037) but get the error below: Error: protect(): protection stack overflow Execution halted or Error: C stack usage is too close to the limit Execution halted I tried to increase the stack size by changing: extern uintptr_t

[R] plotting moving range control chart with qcc. . .

2010-08-20 Thread jim smith
This regards an old post that posed the question: Tom Hodgess wrote: The problem is the (apparent?) inability to produce moving range process behavior (a.k.a. control) charts with individuals data in the package qcc (v. 2.0). I have also struggled with the same limitation in package IQCC (v.

Re: [R] Latex no where to be seen

2010-08-20 Thread Duncan Murdoch
Philippe Grosjean wrote: On 20/08/10 09:11, Joshua Wiley wrote: On Thu, Aug 19, 2010 at 9:18 PM, Donald Paul Winston satchwins...@yahoo.com wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But ?latex() gives me an error. The package

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Ted Harding
See in-line below. On 20-Aug-10 10:37:50, Gavin Simpson wrote: On Fri, 2010-08-20 at 10:54 +0100, ted.hard...@manchester.ac.uk wrote: On 20-Aug-10 09:24:17, Gavin Simpson wrote: On Thu, 2010-08-19 at 20:32 -0700, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Scionforbai
Which format is better? pdf, jpeg or tiff? png and jpeg are bitmap formats, and you can trigger the resolution playing with options width, height (in pixels) and res (nominal DPI) till you find what is suitable for your purpose. Png is definitely the way to go for graphs. Pdf is the overall best

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Barry Rowlingson
On Fri, Aug 20, 2010 at 12:32 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: Now that you remind me, I do recall sending EPS (not EPSI) files to someone and they could see them in Word (recent version). So perhaps Office has moved on (slightly)! Like you, I use Linux and try to avoid

Re: [R] Pass By Value Questions

2010-08-20 Thread Jens Oehlschlägel
Jeff, R has 'environments' as a general mechanism to pass around objects by reference. However, that does not help with most functions like 'apply' which take arguments other than environments. I'm familiar with FF and BigMemory, but are there any packages/tricks which allow for passing

Re: [R] Deviance Residuals

2010-08-20 Thread John Fox
Dear Iasonas, -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Iasonas Lamprianou Sent: August-20-10 5:55 AM To: r-help@r-project.org Subject: [R] Deviance Residuals Dear all, I am running a logistic regression and this

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Wu Gong
Hi Cecília, Assuming that you want to pair two samples as more as possible, let's relabel the k1 indices. The x$k1 c(1,1,2,3,3,5) could be relabeled as c(1.01,1.02,2.01,3.01,3.02,5.01), so as y$k1. x$k1 - x$k1+sequence(rle(x$k1)$lengths)/100 y$k1 - y$k1+sequence(rle(y$k1)$lengths)/100

Re: [R] Deviance Residuals

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 5:54 AM, Iasonas Lamprianou wrote: Dear all, I am running a logistic regression and this is the output: glm(formula = educationUniv ~ brncntr, family = binomial) Deviance Residuals: Min 1Q Median 3Q Max # αυτά είναι τα υπόλοιπα -0.8825 -0.7684

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 6:44 AM, Cecilia Carmo wrote: Hi everyone! I'm matching two samples to create one sample that have pairs of observations equal for the k1 variable. Merge() doesn't work because I dont't want to recycle the values. When there is more than one possible match in either y

Re: [R] Deviance Residuals

2010-08-20 Thread Iasonas Lamprianou
that's enough, thank you now I need to do some more background reading... cheers Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-713178 Fax:

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 6:44 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: Hi everyone! I'm matching two samples to create one sample that have pairs of observations equal for the k1 variable. Merge() doesn't work because I dont't want to recycle the values. x - data.frame(k1=c(1,1,2,3,3,5),

Re: [R] Deviance Residuals

2010-08-20 Thread Iasonas Lamprianou
Thanks, I'll try to put my hands on the reference By the way, would it be easier if I just checked out the code by which the glm function computes the residuals? Or maybe this is not a very good idea. And if it is, how can I check out the source, I never really found out! jason Dr. Iasonas

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Kay Cichini
hello cecilia, i tried: yn-y[y$k1%in%x$k1,] xn-x[x$k1%in%y$k1,] x1st-unique(match(yn$k1,xn$k1)) y1st-unique(match(xn$k1,yn$k1)) new-cbind(k1=intersect(y$k1,x$k1),k2.x=xn[x1st,2],k2.y=yn[y1st,2]) giving: k1 k2.x k2.y [1,] 1 20 10 [2,] 2 22 12 [3,] 3 23 14 [4,] 5 25 16

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Kay Cichini
..gabor gave the solution during i was typing - so please disregard this. yours, kay - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context:

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 5:54 AM, (Ted Harding) wrote: On 20-Aug-10 09:24:17, Gavin Simpson wrote: On Thu, 2010-08-19 at 20:32 -0700, Roslina Zakaria wrote: Hi, I need some opinion. I would like to use graph that I generate from R code and save it into word document. Whichformat is better? pdf,

[R] Shifting of Principal amount while calculating Present Value

2010-08-20 Thread Sarah Sanchez
Dear R Helpers I have following data - cash_flow = c(7, 7, 107)  # 107 = Principle 100 + interest of 7% t = c(1,2,3) and zero rate table as rating year1   year2   year3 AAA    3.60    4.17  4.73 AA  3.65    4.22

Re: [R] Latex no where to be seen

2010-08-20 Thread Philippe Grosjean
On 20/08/10 13:22, Duncan Murdoch wrote: Philippe Grosjean wrote: On 20/08/10 09:11, Joshua Wiley wrote: On Thu, Aug 19, 2010 at 9:18 PM, Donald Paul Winston satchwins...@yahoo.com wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Cecilia Carmo
This is what I need, but my dataframe has many rows and it returns to me the following message: Error: cannot allocate vector of size 120 Kb In addition: There were 18 warnings (use warnings() to see them) Could you help me, Thanks Cecília Em Fri, 20 Aug 2010 08:30:16 -0400 Gabor

[R] Assign statistically relevant groups following multimodal distribution of data

2010-08-20 Thread Sbarrato, Thomas
Dear R mailing list members,   I am trying to find a way to assign samples into groups depending on their mathematical distribution in a computational way (that is, looped many times). My input is by row, a series of 34 samples per row with n rows. The goal would be to compute the distribution

Re: [R] dimnames

2010-08-20 Thread Jimmy Söderly
Thanks Joshua, David and Peter. 2010/8/18 Peter Dalgaard pda...@gmail.com On 08/18/2010 06:20 PM, Joshua Wiley wrote: matrix(1:4, nrow = 4, ncol = 1, dimnames = list(NULL, true value)) or even cbind(true value=1:4) -- Peter Dalgaard Center for Statistics, Copenhagen Business School

Re: [R] Deviance Residuals

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 8:25 AM, Iasonas Lamprianou wrote: Thanks, I'll try to put my hands on the reference By the way, would it be easier if I just checked out the code by which the glm function computes the residuals? Or maybe this is not a very good idea. And if it is, how can I check out

Re: [R] logistic regression tree

2010-08-20 Thread Kay Cichini
hello gavin achim, thanks for responding. by logistic regression tree i meant a regression tree for a binary response variable. but as you say i could also use a classification tree - in my case with only two outcomes. i'm not aware if there are substantial differences to expect for the two

[R] how to interpret KS test

2010-08-20 Thread Izidine Pinto
Dear R users I am using KS test to compare two different distribution for the same variable (temperature) for two different time periods. H0: the two distributions are equal H1: the two distributions are different ks.test (temp12, temp22) Two-sample Kolmogorov-Smirnov test data: temp12

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Cecilia Carmo
The rule is not important to me. I'm selecting a sample that must have one important feature: the same number of obs from x and from y with the same k1. Thanks Cecília Em Fri, 20 Aug 2010 08:28:24 -0400 David Winsemius dwinsem...@comcast.net escreveu: On Aug 20, 2010, at 6:44 AM, Cecilia

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 9:07 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: This is what I need, but my dataframe has many rows and it returns to me the following message: Error: cannot allocate vector of size 120 Kb In addition: There were 18 warnings (use warnings() to see them) Assuming

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Cecilia Carmo
It wasn't the merge command. It doesn't create the variable x3. Cecília Em Fri, 20 Aug 2010 09:22:29 -0400 Gabor Grothendieck ggrothendi...@gmail.com escreveu: On Fri, Aug 20, 2010 at 9:07 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: This is what I need, but my dataframe has many rows and

Re: [R] Latex no where to be seen

2010-08-20 Thread Erik Iverson
On 08/19/2010 11:18 PM, Donald Paul Winston wrote: I'm experimenting using R as a report writer. I'm told LaTex is the destination for my quest. But ?latex() gives me an error. The package manager does not have it. The package installer can't find it. Where is it? It amazes me that there's not

[R] handling recursion relation

2010-08-20 Thread Shant Ch
Hi, I wanted to compute the value of the function ifn at certain values of n. But I am receiving the following error when I was using the following code(given at the end). Error: evaluation nested too deeply: infinite recursion / options(expressions=)? I feel that since the function Grx is

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 9:27 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: It wasn't the merge command. It doesn't create the variable x3. Cecília How about: x$k3 - with(x, unlist(tapply(k1, k1, seq_along))) y$k3 - with(y, unlist(tapply(k1, k1, seq_along)))

[R] Reshape package and variable assignment

2010-08-20 Thread Hosack, Michael
R Experts, I would like to create a series of variables without having to assign a separate line of code for each new variable. My abbreviated dataframe (DF) contains two groups of linked variables (SP1:SP3) and (CAU1:CAU3). Within SP1:SP3 are character codes (full dataframe contains 26

Re: [R] how to interpret KS test

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 9:19 AM, Izidine Pinto wrote: Dear R users I am using KS test to compare two different distribution for the same variable (temperature) for two different time periods. H0: the two distributions are equal H1: the two distributions are different ks.test (temp12, temp22)

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 9:33 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Fri, Aug 20, 2010 at 9:27 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: It wasn't the merge command. It doesn't create the variable x3. Cecília How about: x$k3 - with(x, unlist(tapply(k1, k1,

Re: [R] how to interpret KS test

2010-08-20 Thread Cedric Laczny
Hi Izidine, On Friday, 20. August 2010 15:19:24 Izidine Pinto wrote: Dear R users I am using KS test to compare two different distribution for the same variable (temperature) for two different time periods. H0: the two distributions are equal H1: the two distributions are different

Re: [R] logistic regression tree

2010-08-20 Thread Achim Zeileis
On Fri, 20 Aug 2010, Kay Cichini wrote: hello gavin achim, thanks for responding. by logistic regression tree i meant a regression tree for a binary response variable. but as you say i could also use a classification tree - in my case with only two outcomes. i'm not aware if there are

Re: [R] Latex no where to be seen

2010-08-20 Thread David Hajage
It amazes me that there's not a built in report function that can produce the same kinds of reports that every report writer and data analysis software in the whole word can do. (see SAS, Crystal Reports, SPSS, Oracle Reports, Actuate, Hyperion, Cognos, ..etc) To do some reports, most of R users

Re: [R] handling recursion relation

2010-08-20 Thread jim holtman
Look at your first function: Grx-function(x,r) { G0=Grx(x,0) G0=1; fu1-function(t){exp(-t^2/(2*r*(r+1)))*Grx(x,r-1)} return(integrate(fu1,0,x)$value) } It immediately calls itself for infinite recursion. How did you expect to stop it? I would have expected some type of test before the call. I

Re: [R] Latex no where to be seen

2010-08-20 Thread jim holtman
I will ask a variation on my favorite question that is on my signature line: Tell me what you want to do, not how you want to do it. On Fri, Aug 20, 2010 at 9:48 AM, David Hajage dhajag...@gmail.com wrote: It amazes me that there's not a built in report function that can produce the same kinds

Re: [R] how to interpret KS test

2010-08-20 Thread Cedric Laczny
On Friday, 20. August 2010 15:19:24 Izidine Pinto wrote: Dear R users I am using KS test to compare two different distribution for the same variable (temperature) for two different time periods. H0: the two distributions are equal H1: the two distributions are different ks.test (temp12,

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Cecilia Carmo
The last one worked! The other one, not. Thank you very much! Another question about merge(): sometimes I'm merging two dataframes and the merged dataframe has much more rows than the two merged? I think this should not happen,does it? Cecília Em Fri, 20 Aug 2010 09:37:26 -0400 Gabor

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 9:59 AM, Cecilia Carmo cecilia.ca...@ua.pt wrote: The last one worked! The other one, not. Thank you very much! Another question about merge(): sometimes I'm merging two dataframes and the merged dataframe has much more rows than the two merged? I think this should

Re: [R] probabilities from predict.svm

2010-08-20 Thread Steve Lianoglou
Awesome! Good news, James. Thanks for letting us know. Glad you were able to sort this out. -steve On Thu, Aug 19, 2010 at 5:00 PM, Watling,James I watli...@ufl.edu wrote: Hi Steve-- I spent some more time tuning the model with alternative gamma and cost values, but still kept coming back

Re: [R] How to make 'step' faster?

2010-08-20 Thread Martin Maechler
GS == Gavin Simpson gavin.simp...@ucl.ac.uk on Fri, 20 Aug 2010 10:13:38 +0100 writes: GS On Fri, 2010-08-20 at 12:25 +0530, sambit rath wrote: Dear all, I am fairly new to R. I would like to perform a step-wise logit regression aiming to select a model on the basis

Re: [R] Deviance Residuals

2010-08-20 Thread Renzo Tascheri
Dear Dr. Iasonas, perhaps the following reference could also be useful: Dunn, P.K., Smyth, G.K., 1996. Randomized quantile residuals. J. Comput. Graph. Stat. 5, 236–244. An introduction to this residual definition and a link to the paper could be find here:

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Stuart Luppescu
On Fri, 2010-08-20 at 01:30 -0700, Joshua Wiley wrote: I usually save them from R as a PDF or postscript file, rasterize them in GIMP (free answer to Photoshop) at the desired resolution, and finally choose the desired format/compression (jpeg, png, bitmap, tiff, etc.) to save it as from

[R] Problem Installing R-Commander

2010-08-20 Thread Stephen P Molnar
I am trying to install R Commander in my R installation on a laptop running the 64 bit version of Windows 7 Professional. I was running R as the Administrator. The downloads proceed normally via CRAN, but then I get the following error messages: package 'slam' successfully unpacked and

Re: [R] Latex no where to be seen

2010-08-20 Thread Michael Sumner
Fortune alert: It's actually slightly more complicated: ?ls() really is parsed as ? of ls(), not (?ls)(). On Fri, Aug 20, 2010 at 11:53 PM, jim holtman jholt...@gmail.com wrote: I will ask a variation on my favorite question that is on my signature line: Tell me what you want to do, not how

[R] Date Inconsistencies? Buglets?

2010-08-20 Thread ivo welch
The treatment of dates seems to be a little inconsistent in R 2.11.1 (2010-05-31): [1] The choice of origins? as.integer(as.Date(1970-01-01)) works and assumes as origin 1970-01-01. However, as.Date(1) does not work. It requires an origin (as.Date(1, origin=1970-01-01)). If we set a

Re: [R] Aggregate Help

2010-08-20 Thread Hall, Ken (CDC/OSELS/NCPHI)
The first one worked great. Thanks for the quick response. The second one, with the with statement, gave me the following error: Error in aggregate.date.frame(as,dat.frame(x), ...): arguments must have same length Ken PS. What is the proper etiquette? Should I respond to you individually and

Re: [R] Aggregate Help

2010-08-20 Thread Hall, Ken (CDC/OSELS/NCPHI)
Both of these worked great. Thanks for the quick response. Question: which of the three codes (these two and Phil's) would be preferred from a performance in R perspective? They all ran extremely fast. Ken -Original Message- From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]

Re: [R] Date Inconsistencies? Buglets?

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 10:35 AM, ivo welch ivo.we...@gmail.com wrote: The treatment of dates seems to be a little inconsistent in R 2.11.1 (2010-05-31): [1] The choice of origins?   as.integer(as.Date(1970-01-01)) works and assumes as origin 1970-01-01.  However,   as.Date(1) The zoo

[R] Problem with POSIXct in ave

2010-08-20 Thread Sasha Hafner
Hi, I am having trouble using the ave function with a POSIXct object. For example: x-Sys.time()+0:9*3600 dat-data.frame(id=rep(c('a',' b','c'),each=10),dt=rep(x,3),i=rep(1:10,3)) dat # This is what I want to do: dat$time.elapsed-unsplit(lapply(split(dat$dt,dat$id),function(x) x-x[1]),f=dat$id)

Re: [R] Aggregate Help

2010-08-20 Thread William Dunlap
Ken, Your error message Error in aggregate.date.frame(as,dat.frame(x), ...): arguments must have same length would arise if you called the 'Date' variable 'date'. R is a case-sensitive language: 'date' is a function object (with length 1) and 'Date' is the variable in your data.frame

Re: [R] logistic regression tree

2010-08-20 Thread Frank Harrell
It would be good to tell us of the frequency of observations in each category of Y, and the number of continuous X's. Recursive partitioning will require perhaps 50,000 observations in the less frequent Y category for its structure and predicted values to validate, depending on X and the

[R] Surface generated as product of two curves

2010-08-20 Thread mau...@alice.it
I have the following data table: [,1] [,2] [,3] [,4] [,5] [1,] 2575 1927 1754 581 354 [2,] 1156 810 730 541 237 [3,] 417 297 199 125 110 [4,] 281 132 132 47 58 [5,] 152 84 87 54 19 Each column represents a time interval expressed in nano-seconds Each row

Re: [R] 64 bit RSQLite

2010-08-20 Thread Seth Falcon
Hi, On 8/14/10 11:43 PM, Stephen Liu wrote: Could you please explain a little bid further whether after connecting the website I can select either 32 or 64 bit version? R packages come in pre-compiled binary flavors for some OS/architectures and always as source packages. If you are on

Re: [R] Problem with POSIXct in ave

2010-08-20 Thread David Winsemius
On Aug 20, 2010, at 10:56 AM, Sasha Hafner wrote: Hi, I am having trouble using the ave function with a POSIXct object. For example: x-Sys.time()+0:9*3600 dat-data.frame(id=rep(c('a',' b','c'),each=10),dt=rep(x,3),i=rep(1:10,3)) dat # This is what I want to do:

Re: [R] aggregate with cummax ??

2010-08-20 Thread skan
Hi Just in case somebody could be interested I've found that unlist(lapply(split(mydata,as.Date), cummax)) could make the job but I don't know if it's the proper way if I want to use zoo objects. Any suggestion?? -- View this message in context:

Re: [R] Aggregate Help

2010-08-20 Thread Hall, Ken (CDC/OSELS/NCPHI)
David pointed out that it would have helped to include the results of str(pilot) in my original post. I ran the str(pilot) and, sure enough, it revealed the problem (my error), that the column name should be date1. Also, the editor in Outlook capitalized Date, so that was another problem. Now,

[R] R and Java

2010-08-20 Thread Allan Freitas
Hi folks, Maybe that is not the best place to ask, but I 'd like some of you could help me... I'm using jri to run R commands under Java, by its examples I can eval expressions, returning single values or vectors to my Java variables. But, I would like to know how to assign the values of a Java

Re: [R] Date Inconsistencies? Buglets?

2010-08-20 Thread Gavin Simpson
On Fri, 2010-08-20 at 10:35 -0400, ivo welch wrote: The treatment of dates seems to be a little inconsistent in R 2.11.1 (2010-05-31): snip / [2] How do dates on axes work? plot( c(as.Date(1:20, origin=1970-01-01)), 1:20 ) axis( side=3, c(as.Date(1:20, origin=1970-01-01))) The

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread Joshua Wiley
On Fri, Aug 20, 2010 at 7:30 AM, Stuart Luppescu s...@ccsr.uchicago.edu wrote: On Fri, 2010-08-20 at 01:30 -0700, Joshua Wiley wrote: I usually save them from R as a PDF or postscript file, rasterize them in GIMP (free answer to Photoshop) at the desired resolution, and finally choose the

Re: [R] R and Java

2010-08-20 Thread Romain Francois
Le 20/08/10 17:56, Allan Freitas a écrit : Hi folks, Maybe that is not the best place to ask, No it is not, this is: http://mailman.rz.uni-augsburg.de/mailman/listinfo/stats-rosuda-devel but I 'd like some of you could help me... I'm using jri to run R commands under Java, by its

Re: [R] aggregate with cummax ??

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 11:38 AM, skan juanp...@gmail.com wrote: Hi Just in case somebody could be interested I've found that unlist(lapply(split(mydata,as.Date), cummax)) could make the job but I don't know if it's the proper way if I want to use zoo objects. Any suggestion?? See

Re: [R] Surface generated as product of two curves

2010-08-20 Thread Bert Gunter
Essentially, you just want to fit a smooth surface over your grid of x,y values. Lots of ways to do this, but check out ?loess for a simple base R approach. The mgcv package does tensor product splines. Also search on thin plate splines, 2d spatial and the like for other approaches. There are

[R] U value from wilcox.test

2010-08-20 Thread Chloe
Dear all, I want to compare the efficiency of 2 methods in extracting proteins from algal samples. I collected 6 independant algal samples and I extracted 3 by the method 1 and 3 others by the method 2. So I have 2 groups of 3 samples, that are not paired. I would like to know if the results

Re: [R] which one give clear picture-pdf, jpg or tiff?

2010-08-20 Thread kees duineveld
So do I. You can even convert these to MS graph format and edit some. It is sharp, because this is also vector format. On Fri, 20 Aug 2010 07:58:01 +0200, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: I use Rgui and copy/paste special windows metafile. Roslina Zakaria zrosl...@yahoo.com

[R] differecing a zoo series

2010-08-20 Thread steven mosher
A quick question x - as.yearmon(2000 + seq(0, 23)/12) x [1] Jan 2000 Feb 2000 Mar 2000 Apr 2000 May 2000 Jun 2000 Jul 2000 Aug 2000 Sep 2000 Oct 2000 Nov 2000 Dec 2000 Jan 2001 [14] Feb 2001 Mar 2001 Apr 2001 May 2001 Jun 2001 Jul 2001 Aug 2001 Sep 2001 Oct 2001 Nov 2001 Dec 2001

Re: [R] differecing a zoo series

2010-08-20 Thread Gabor Grothendieck
On Fri, Aug 20, 2010 at 2:19 PM, steven mosher mosherste...@gmail.com wrote: A quick question  x - as.yearmon(2000 + seq(0, 23)/12)  x  [1] Jan 2000 Feb 2000 Mar 2000 Apr 2000 May 2000 Jun 2000 Jul 2000 Aug 2000 Sep 2000 Oct 2000 Nov 2000 Dec 2000 Jan 2001 [14] Feb 2001 Mar 2001 Apr 2001

Re: [R] U value from wilcox.test

2010-08-20 Thread Cedric Laczny
Hi Chloe, first of all, I want to note, that you should be careful using the WMW-test. Even though it is often reported to be some sort of a swiss-army-knife for comparing two distributions, recent research on this test has revelaed that it is crucial what hypotheses you consider. Also the

Re: [R] U value from wilcox.test

2010-08-20 Thread Charles C. Berry
Issues like that in this thread can often be resolved by reading the help page for the relevant function. From: ?wilcox.test Note The literature is not unanimous about the definitions of the Wilcoxon rank sum and Mann-Whitney tests. The two most common definitions correspond to the

  1   2   >