Re: [R] Kolmogorov-Smirnov test

2011-09-24 Thread rommel
Dear Dr. Snow, nbsp; Thank you for your reply. nbsp; 1. Are you doing the 2 sample KS test? Comparing if 2 samples come from the same distribution? -Yes, I am doing 2-sample KS test nbsp; 2. With 3,000 points you will still likely have power to find meaningless differences, what exactly

[R] ggplot2: expression() in legend labels?

2011-09-24 Thread Casper Ti. Vector
Is there any way to use expression() in legend labels with ggplot2? It seems that things like scale_shape_manual(value = c( x = expression(italic(x)), y = expression(italic(y)) )) don't work. Thanks very much :) -- Using GPG/PGP? Please get my current public key (ID: 0xAEF6A134,

Re: [R] R-help Digest, Vol 103, Issue 23

2011-09-24 Thread Richard Valliant
I will be out of the office Sep 32-28, 2011, with limited email access. For immediate help, please call the JPSM main number, 301-314-7911. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] R-help Digest, Vol 103, Issue 23

2011-09-24 Thread mihalicza . peter
Szeptember 12-től 26-ig irodán kívül vagyok, és az emailjeimet nem érem el. Sürgős esetben kérem forduljon Kárpáti Edithez (karpati.e...@gyemszi.hu). Üdvözlettel, Mihalicza Péter I will be out of the office from 12 till 26 September with no access to my emails. In urgent cases please contact

[R] Assessing prediction performance of SVM using e1071 package

2011-09-24 Thread Jokel Meyer
Dear R-Users! I am using the svm function (e1071 package) to classify two groups using a set of 180 indicator variables. Now I am confused about the cross-validation procedure. (A) On one hand I use the setting cross=10 in the svm function to run 10 cross-validation iterations and to get an

[R] 45 degress line in zoo plot

2011-09-24 Thread Manta
Dear all, I am not able to draw a 45 degrees line in a zoo plot. The time index is a sequence of 15 min bins from 7 in the morning to 18:45 in the evening. Thanks, Marco -- View this message in context: http://r.789695.n4.nabble.com/45-degress-line-in-zoo-plot-tp3839182p3839182.html Sent from

[R] help

2011-09-24 Thread Mathew Brown
Mathew Brown Institute of Bioclimatology University of Göttingen Büsgenweg 2 37077 Göttingen, Germany t: +49 551 39 9359 mathew.br...@forst.uni-goettingen.de On 9/24/2011 6:00 PM, r-help-requ...@r-project.org wrote: Send R-help mailing list submissions to r-help@r-project.org To

Re: [R] help

2011-09-24 Thread peter dalgaard
Please read instructions CAREFULLY. This is the R-help mailing list, not R-help-request. People on this list are not interested in an empty mail, much less one that is followed by an entire digest message! On Sep 24, 2011, at 13:56 , Mathew Brown wrote: Mathew Brown Institute of

Re: [R] ncdf - install error

2011-09-24 Thread David William Pierce
2011/9/24 Muhammad Rahiz muhammad.ra...@ouce.ox.ac.uk: Dear all, I'm having issues with the installation of the ncdf package. It returns a non-zero exit status. Can anyone suggest what I should do next? FYI, I do not have problems installing other packages. [...] /usr/bin/ld:

[R] Help with Integration Output

2011-09-24 Thread ChrisQ
Hi, I need to do a calculation in R which involves adding the estimate of a numerical integration with another number - however, if I add the number to the output of the intergrate function (which when run on its own provides me with a number with error), I get the message 'non-numeric argument

Re: [R] comparing mixed binomial model against the same model without random effect

2011-09-24 Thread lincoln
Any answer to this? I really need to compare a mixed model with binomial error against the same model without the random effect. I would use anova() but I don't know how to specify both models in order to make them comparable. Thanks for any answer Simone -- View this message in context:

Re: [R] Help with Integration Output

2011-09-24 Thread B77S
It is not simply the answer, it is a list str(integrate(dnorm, -1.96, 1.96)) List of 5 $ value : num 0.95 $ abs.error : num 1.05e-11 $ subdivisions: int 1 $ message : chr OK $ call: language integrate(f = dnorm, lower = -1.96, upper = 1.96) - attr(*, class)= chr

[R] Supporting R/Membership

2011-09-24 Thread Joshua Wiley
Dear R Users, I know issues like this have come up in the past. I am wondering how many people would be interested in an online mechanism for donating to the R Foundation or for becoming supporting members. Right now, the link for it is somewhat buried on the R home page, and once you get

Re: [R] Supporting R/Membership

2011-09-24 Thread David Cross
For my part, knowing all of my (mostly bad) habits, I would only donate if I could do it online. And I would be happy to donate! Cheers David Cross d.cr...@tcu.edu www.davidcross.us On Sep 24, 2011, at 10:43 AM, Joshua Wiley wrote: Dear R Users, I know issues like this have come up in

[R] two ecdf in the same figure

2011-09-24 Thread Adel ESSAFI
Hello, is ot possible to draw two ecdf of vectors (say s1 and s2) on the same figire with R. plot function draws a new plot and there is no function like points or lines to draw a second ecdf on the figure. Regards -- *PhD candidate in Computer Science Address 3 avenue lamine, cité ezzahra,

Re: [R] Cross Spectrum : Conversion of 2-D spectrum into a single complex array

2011-09-24 Thread Joseph Park
Perhaps this question is inappropriate or posted to the wrong list? Any guidance would be appreciated. Thanks. --- On Fri, 9/23/11, Joseph Park jpark...@att.net wrote: From: Joseph Park jpark...@att.net Subject: Cross Spectrum : Conversion of 2-D spectrum into a single complex array To:

Re: [R] two ecdf in the same figure

2011-09-24 Thread R. Michael Weylandt
What version of R are you using? x = ecdf(rnorm(50)); y = ecdf(rnorm(50)) plot(x); lines(y, col=2) works for me on R 2.13.1. Specifically, there is a lines.stepfun method (albeit non-visible) which handles this. Michael On Sat, Sep 24, 2011 at 12:23 PM, Adel ESSAFI adel.s...@imag.fr wrote:

[R] Can't reliably use RefClass methods in Snowfall

2011-09-24 Thread henry
Greetings, I am trying to use Reference Class methods in Snowfall, using R 2.12.1 on Ubuntu Natty. Using then directly seems to work (stanza 2 below), but using them indirectly does not (stanza 3 below). I get an attempt to apply non-function error. In addition to exporting the instance

Re: [R] Cross Spectrum : Conversion of 2-D spectrum into a single complex array

2011-09-24 Thread Ben Bolker
Joseph Park jpark.us at att.net writes: Perhaps this question is inappropriate or posted to the wrong list? Any guidance would be appreciated. Thanks. [snip] This seems like a perfectly reasonable question. Sometimes people just don't bite. Ben Bolker

Re: [R] Tobit Fixed Effects

2011-09-24 Thread Felipe Nunes
Hi Arne, my problem persists. I am still using censReg [version - 0.5-7] to run a random effects model in my data (50,000 cases), but I always get the message. tob7 - censReg(transfers.cap ~ pt.pt + psdb.pt + pt.opp + pt.coa + psdb.coa + pib.cap + transfers.cap.lag + pib.cap + ifdm +

Re: [R] Tobit Fixed Effects

2011-09-24 Thread Joshua Wiley
Can you fit a simpler model successfully? I would play around with all your RHS variables to try to see if there is any particular one that causes the estimation problems. Two ways I would do this would be to drop them one at a time, or to start from an intercept only model (i.e., transfers.cap

Re: [R] two ecdf in the same figure

2011-09-24 Thread Joshua Wiley
Or just ask let the second plot know to add it to the existing plot: set.seed(2) plot(ecdf(rnorm(50))) plot(ecdf(rnorm(50,,2)), add = TRUE, pch = 2) legend(-2.9, .95, legend = c(First ECDF, Second ECDF), pch = c(16, 2)) note the limits are set by the first plot, so watch out for that. You may

[R] maptools::sunriset() daylight savings to stardard time change

2011-09-24 Thread Folkes, Michael
Hello all, After beating my head against the wall for a few hours, I give up. I have two problems. 1. If I use seq() to generate a POSIXct series that crosses a daylight savings to stardard time change I get two values for the first day of change. The time change is Nov 6,2011. The value for

Re: [R] comparing mixed binomial model against the same model without random effect

2011-09-24 Thread Billy
Hello Simone, Given that your response variable is binary and, consequently, you should use generalized models, just occurs to me a comparison between a Generalized Linear Model (the model without the random effect) and a Generalized Linear Mixed Model (the model with the random effect). You

Re: [R] function censReg in panel data setting

2011-09-24 Thread felipnunes
Hi guys, did you figure out a way to solve your problem, Igor? I'm still trying to understand what's happening to my estimations. My data set is also big (50,000) and I get the same error: Error in maxNRCompute(fn = logLikAttr, fnOrig = fn, gradOrig = grad, hessOrig = hess, : NA in the

Re: [R] Tobit Fixed Effects

2011-09-24 Thread Felipe Nunes
That was a great idea, Josh! I found a problem even when I tried a simple model. Do you have any idea why? tob6 - censReg(transfers.cap ~ 1, left=0, right=Inf, method=BHHH, nGHQ=20, iterlim=1, data = pdata2) Error in uniqval[as.character(effect), , drop = F] : incorrect number of

Re: [R] Tobit Fixed Effects

2011-09-24 Thread Joshua Wiley
On Sat, Sep 24, 2011 at 4:55 PM, Felipe Nunes felipnu...@gmail.com wrote: That was a great idea, Josh! I found a problem even when I tried a simple model. Do you have any idea why? I forgot you were using a pdata.framecensReg hands off estimation in part to plm, which wants between and

Re: [R] function censReg in panel data setting

2011-09-24 Thread Arne Henningsen
On 24 September 2011 23:58, felipnunes felipnu...@gmail.com wrote: Hi guys, did you figure out a way to solve your problem, Igor? I'm still trying to understand what's happening to my estimations. My data set is also big (50,000) and I get the same error: Error in maxNRCompute(fn =

Re: [R] Tobit Fixed Effects

2011-09-24 Thread Arne Henningsen
On 25 September 2011 00:16, Felipe Nunes felipnu...@gmail.com wrote: Hi Arne, my problem persists. I am still using censReg [version - 0.5-7] to run a random effects model in my data (50,000 cases), but I always get the message. tob7 - censReg(transfers.cap ~ pt.pt + psdb.pt + pt.opp + pt.coa

Re: [R] maptools::sunriset() daylight savings to stardard time change

2011-09-24 Thread Jeff Newmiller
Sys.setenv(TZ=Etc/GMT+8) will use standard time only for time calculations. --- Jeff Newmiller The . . Go Live... DCN:jdnew...@dcn.davis.ca.us Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research