Re: [R] plot in package psych with function error.bars.by

2014-06-16 Thread Tham Tran
Hi William, I've just updated your latest package psych_1.4.6.11.zip from server personality-project/r/src/contrib/. One time the updating process was finished, i tried to run based samples code: require(psych)

[R] Aggregating 15 minute xts sequence to hourly

2014-06-16 Thread Costas Vorlow
Dear all, Why aggregation of 15 minute xts data happens on the 45th (3rd quarter) and not the exact hour close (i.e., 00) time? For example, temp below is an xts sequence with 15-minute frequency: quarters - ISOdatetime(2012,05,02,9,0,0) + seq(0:39)*15*60; set.seed(42); observation -

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread peter dalgaard
On 16 Jun 2014, at 05:22 , array chip arrayprof...@yahoo.com wrote: Hi, I am using clogit() from survival package to do conditional logistic regression. I also need to make prediction on an independent dataset to calculate predicted probability. Here is an example: dat -

[R] correlation given p value and sample size

2014-06-16 Thread Witold E Wolski
Hi, Looking for and function which produces the minimum r (pearson correlation) so that H0 (r=0) can be rejected, given sample size and p-value? Witold -- Witold Eryk Wolski __ R-help@r-project.org mailing list

Re: [R] Defining default method for S3, S4 and R5 classes

2014-06-16 Thread Luca Cerone
thanks Suzen 2014-06-15 2:34 GMT+02:00 Suzen, Mehmet msu...@gmail.com: There is a nice tutorial on this: http://adv-r.had.co.nz/OO-essentials.html For an in depth guide, have a look at the book from John Chambers, Software for data analysis programming with R. On 13 June 2014 12:20, Luca

Re: [R] Aggregating 15 minute xts sequence to hourly

2014-06-16 Thread Joshua Ulrich
On Mon, Jun 16, 2014 at 3:41 AM, Costas Vorlow costas.vor...@gmail.com wrote: Dear all, Why aggregation of 15 minute xts data happens on the 45th (3rd quarter) and not the exact hour close (i.e., 00) time? The 00 time is the beginning of the hour, not the end. E.g., 10:00:00 is the beginning

Re: [R] correlation given p value and sample size

2014-06-16 Thread peter dalgaard
There's a simple relation t = r / sqrt(1 - r^2) * sqrt(n - 2) r = t / sqrt(n - 2 + t^2) where t has a t distribution on n-2 df. Insert t = +-qt(p/2, n-2). -pd On 16 Jun 2014, at 11:23 , Witold E Wolski wewol...@gmail.com wrote: Hi, Looking for and function which produces the minimum r

Re: [R] Help with SEM package - model significance

2014-06-16 Thread John Fox
Dear Bernardo, The df for the LR chisquare over-identification test come not from the number of observations, but from the difference between the number of observable variances and covariances, on the one hand, and free parameters to estimate, on the other. In your case, these numbers are

Re: [R] Aggregating 15 minute xts sequence to hourly

2014-06-16 Thread Costas Vorlow
Dear Joshua, Thanks for your reply. As I see, the solution you suggest aligns the time stamps as required but leaves the aggregation results as is. Hence, the last quarter data of every hour are not aggregated still... Am I right or am I understanding something wrongly? I tried to move ahead

[R] Determination lag order - problem with daily data and AR / ARIMA

2014-06-16 Thread serena1234
Hello, I am trying to determine a lag order for my data with the help of AIC and/ or BIC in order to conduct further tests. It is about prices measured at a daily frequency (weekends and holidays excluded). My first approach was to approximate the process with an AR model using the function

[R] xml package - free document / memory leak

2014-06-16 Thread Loos, Martin
Dear R helpers I have a problem of releasing memory after having read .xml-files using the xml package (version 3.98-1.1, R 3.0.2, Windows 7, 64bit). The problem has appeared previously and several solutions/bug fixes have been proposed. I went through many of these and have also read (and

[R] Power graph for two Proportion

2014-06-16 Thread Suganthie Jeyaganth
Dear R mailing listers, I am try to find the different power calculation for p1=0.2 and p2=0.4 , with significant level=0.05  (one sided test) I would like to  have a graph   y -axis as a power and x-axis as a sample size . I run this command for different value of power. and get the n and

[R] Determination lag order - problem with daily data and AR / ARIMA

2014-06-16 Thread serena1234
Hello, I am trying to determine a lag order for my data with the help of AIC and/ or BIC in order to conduct further tests. It is about prices measured at a daily frequency (weekends and holidays excluded). 1) My first approach was to approximate the process with an AR model using the function

[R] Question on JAVA and R

2014-06-16 Thread Luis Eduardo Castillo Méndez
As I discuss today we are trying to do with my partner Steven Penaloza and a small program of supervised classification of images with unconventional using R methods, connecting libraries and own code of R to an environment of JAVA, in our case we are using NETBEANS to interface, the connection

[R] Any refit function available for 'car' package?

2014-06-16 Thread Gang Chen
Suppose that I need to run a multivariate linear model Y = X B + E many times with the same model matrix X but each time with different response matrix Y. Is there a function available in 'car' package similar to refit() in lme4 package so that the model matrix X would not be reassembled each

[R] R, Rserve logging

2014-06-16 Thread Hema Seshadri
I am running R on daemon mode Rserve (to connect to java). It spits out a constant flow output. Is there a way to turn it off or set it to rotate etc. In other words I want to find out how R, Rserve handles logging. [[alternative HTML version deleted]]

[R] glm.fit: fitted probabilities numerically 0 or 1 occurred for a continuous variable?

2014-06-16 Thread Nwinters
I have gotten the this error before: glm.fit: fitted probabilities numerically 0 or 1 occurred and the problem was usually solved by combining one or more categories were there were no observations. I am now having this error show up for a variable that is continuous (not categorical). What

Re: [R] glm.fit: fitted probabilities numerically 0 or 1 occurred for a continuous variable?

2014-06-16 Thread Marc Schwartz
On Jun 16, 2014, at 2:34 PM, Nwinters nicholas.wint...@mail.mcgill.ca wrote: I have gotten the this error before: glm.fit: fitted probabilities numerically 0 or 1 occurred and the problem was usually solved by combining one or more categories were there were no observations. I am now

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread array chip
Thank you Peter. Any other suggestions are absolutely welcome!! John From: peter dalgaard pda...@gmail.com Cc: r-help@r-project.org r-help@r-project.org Sent: Monday, June 16, 2014 2:22 AM Subject: Re: [R] prediction based on conditional logistic regression

[R] Error in merge [negative length vectors are not allowed]

2014-06-16 Thread Kate Ignatius
Hi All, I'm trying to merge two files together using: combinedfiles - merge(comb1,comb2,by=c(Place,Stall,Menu)) comb1 is about 2 million + rows (158MB) and comb2 is about 600K+ rows (52MB). When I try to merge using the above syntax I get the error: Error in merge.data.frame(comb1, comb2, by

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread Charles Berry
peter dalgaard pdalgd at gmail.com writes: On 16 Jun 2014, at 05:22 , array chip arrayprofile at yahoo.com wrote: Hi, I am using clogit() from survival package to do conditional logistic regression. I also need to make prediction on an independent dataset to calculate predicted

Re: [R] Error in merge [negative length vectors are not allowed]

2014-06-16 Thread Prof Brian Ripley
On 16/06/2014 23:21, Kate Ignatius wrote: Hi All, I'm trying to merge two files together using: combinedfiles - merge(comb1,comb2,by=c(Place,Stall,Menu)) comb1 is about 2 million + rows (158MB) and comb2 is about 600K+ rows (52MB). When I try to merge using the above syntax I get the error:

[R] Error: C stack usage

2014-06-16 Thread Mohan Radhakrishnan
Hi, I am not relating this to shiny but this error was thrown by my shiny server code. Does it have anything to do with low memory settings ? I am spawning a JVM using rJava. That JVM is using the attach API to connect to another JVM Thanks, Mohan Error: C stack usage 140730070087404 is