Re: [R] Error: C stack usage

2014-06-16 Thread Prof Brian Ripley
On 17/06/2014 04:27, Mohan Radhakrishnan wrote: 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 ? No. It is an error in the C code in a package you are using which has corrupted R internals and

[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 t

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 e

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

2014-06-16 Thread Charles Berry
peter dalgaard gmail.com> writes: > > > On 16 Jun 2014, at 05:22 , array chip 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

[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, comb

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 Cc: "r-help@r-project.org" Sent: Monday, June 16, 2014 2:22 AM Subject: Re: [R] prediction based on conditional logistic regression clogit > Hi, I am using clogit

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 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 having this error show up

[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 co

[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] 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 tim

[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 was

[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 a

[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 po

[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 unde

[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 ar(x,

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

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 equal

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 wrote: > Hi, > > > Looking for and function which produces the minimum r (pearson > cor

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 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 of the 10-o'cloc

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 : > 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 Cerone wr

[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 https://stat.ethz.ch/mailman/listinfo

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

2014-06-16 Thread peter dalgaard
On 16 Jun 2014, at 05:22 , array chip 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 <- data.frame(set=rep(1:50,each=

[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 <- xts

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) keys.list=list(Agree=c(-1,2:5),Conscientious=c(6:8,-9,-10),Extraversion=c(-11,-12,