Re: [R] Summarizing data based on Date

2015-06-11 Thread Shivi82
Thank you John for spending time on this query and helping out. It really helped me and finally i am able to achieve the desired results. Thanks a ton to all others as well to spending time and furbishing solution. Regards, Shivi -- View this message in context:

Re: [R] Summarizing data based on Date

2015-06-11 Thread PIKAL Petr
Hi I (wrongly) understood that Shivi82 wanted to summarise on month values. Therefore format(test$CR_DT,%m) shall give you month number and list is required by aggregate. All the problem was in test$CHG_WT which seems to be a factor (for whatever reason) Cheers Petr -Original

Re: [R] Issue with mcapply

2015-06-11 Thread AROONALOK PYNE
Thanks a lot Rainer On Thu, Jun 11, 2015 at 4:24 PM, Rainer M Krug rai...@krugs.de wrote: Again: please keep this on r-help! AROONALOK PYNE aroonalok.p...@gmail.com writes: Hi, My system has 4GB memory. On running htop command on linux terminal, for c(999), a() rug on

Re: [R] Issue with mcapply

2015-06-11 Thread Rainer M Krug
Again: please keep this on r-help! AROONALOK PYNE aroonalok.p...@gmail.com writes: Hi, My system has 4GB memory. On running htop command on linux terminal, for c(999), a() rug on single cpu core showed stats for 100% cpu utilization. On reaching b(), the entire system hanged (no

Re: [R] Issue with mcapply

2015-06-11 Thread Rainer M Krug
AROONALOK PYNE aroonalok.p...@gmail.com writes: Please check this code : library(parallel) workerFunc - function(n) { return(n^2) } a - function(){ CurrentTime - Sys.time() res - lapply(values, workerFunc) TimeTaken - Sys.time() - CurrentTime print(TimeTaken) } b - function(){

Re: [R] Issue with mcapply

2015-06-11 Thread Rainer M Krug
Please keep this on the r-help mailing list. AROONALOK PYNE aroonalok.p...@gmail.com writes: R version 3.0.2 (2013-09-25) Large value : 999 I rerun the code as c(1000) for which your machine works fine but my code still hangs on reaching b(). I am running it from Linux Terminal.

[R] Fwd: Issue with mcapply

2015-06-11 Thread AROONALOK PYNE
-- Forwarded message -- From: AROONALOK PYNE aroonalok.p...@gmail.com Date: Thu, Jun 11, 2015 at 2:25 PM Subject: Re: Issue with mcapply To: Rainer M Krug rai...@krugs.de R version 3.0.2 (2013-09-25) Large value : 999 I rerun the code as c(1000) for which your machine

[R] Fwd: Issue with mcapply

2015-06-11 Thread AROONALOK PYNE
-- Forwarded message -- From: AROONALOK PYNE aroonalok.p...@gmail.com Date: Thu, Jun 11, 2015 at 3:49 PM Subject: Re: Issue with mcapply To: Rainer M Krug rai...@krugs.de Hi, My system has 4GB memory. On running htop command on linux terminal, for c(999), a() rug on

Re: [R] Different behavior of model.matrix between R 3.2 and R3.1.1

2015-06-11 Thread Therneau, Terry M., Ph.D.
Frank, I'm not sure what is going on. The following test function works for me in both 3.1.1 and 3.2, i.e, the second model matrix has fewer columns. As I indicated to you earlier, the coxph code removes the strata() columns after creating X because I found it easier to correctly create

Re: [R] Quantile regression model with nonparametric effect and interaction

2015-06-11 Thread Roger Koenker
The main effect trend seems rather dangerous, why not just estimate the f’s in a loop? url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] Problem with Volcanoplot

2015-06-11 Thread John Kane
Do you already have an R package that will do a Spearman Correlation Volcano plot ? What do the data look like? Have a look at http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example and http://adv-r.had.co.nz/Reproducibility.html for some suggestions on how to

Re: [R] Boxplot function error-help required

2015-06-11 Thread MacQueen, Don
In addition to the other answers, I would suggest that the next time you get the could not find function message, try like this: help.search('Boxplot') Among the output from that you should see graphics::boxplot Box Plots which should lead you to boxplot instead of

Re: [R] Boxplot function error-help required

2015-06-11 Thread John Kane
Well it might have worked for your supervisor but I don't see how. As was mentioned it is boxplot , not Boxplot and the rest of the syntax looks dodgy to say the least. Try boxplot(Acc_S$Subj ~ Acc_S$Acc) I don't see how label = will work , ?boxplot says it should be names = and as the

Re: [R] powerTransform Convergence erro

2015-06-11 Thread John Fox
Dear Brittany, There is an essentially perfect linear dependency among the variables in your data (note the last eigenvalue, which is 0 within rounding error): eigen(cor(problem.data.boxcox[,-1]), only.values=TRUE) $values [1] 3.644257e+00 1.821582e+00 1.712152e+00 1.205091e+00

[R] Getting XML Data From Sharepoint

2015-06-11 Thread David Wagle
I'm trying, and so far failing, to extract data from a very large SharePoint list. I have a URL for getting XML from sharepoint: http:// site/_vti_bin/owssvr.dll?Cmd=DisplayList={GUID}Query=*XMLDATA=TRUE When I do something like the following: require(XML) require(RCurl) URL - http://

Re: [R] Boxplot function error-help required

2015-06-11 Thread David Winsemius
On Jun 11, 2015, at 2:19 PM, John Kane wrote: Well it might have worked for your supervisor but I don't see how. As was mentioned it is boxplot , not Boxplot and the rest of the syntax looks dodgy to say the least. There is a Boxplot function in package 'car' and it has a labels

Re: [R] for loop incorrect row count

2015-06-11 Thread David Winsemius
On Jun 11, 2015, at 12:25 PM, Kevin Kowitski wrote: Hey, I am having an issue with a for loop that is intended to read index values by row and column so that it can pull out the valuable information. My issue is that I am using a data.frame(which(df==1, arr.ind=TRUE)) That would be

Re: [R] Boxplot function error-help required

2015-06-11 Thread David Winsemius
On Jun 11, 2015, at 2:44 PM, MacQueen, Don wrote: In addition to the other answers, I would suggest that the next time you get the could not find function message, try like this: help.search('Boxplot') Spencer Graves uses RSiteSearch() as the underlying function for sos::findFn --

Re: [R] powerTransform Convergence erro

2015-06-11 Thread Brittany Demmitt
Hi John, That does help, thanks! Brittany On Jun 11, 2015, at 4:02 PM, John Fox j...@mcmaster.ca wrote: Dear Brittany, There is an essentially perfect linear dependency among the variables in your data (note the last eigenvalue, which is 0 within rounding error):

[R] for loop incorrect row count

2015-06-11 Thread Kevin Kowitski
Hey,    I am having an issue with a for loop that is intended to read index values by row and column so that it can pull out the valuable information.  My issue is that I am using a data.frame(which(df==1, arr.ind=TRUE)) to find the index of the values in my data frame that are equal to 1.  

[R] GAMM - error when using spline smoother, only allowed to use tensor product

2015-06-11 Thread Hannah Linder
Hello, I am familiar with the basics of statistical regression models, including GAMs, but I am stumped on a particular implementation issue. I am constructing a GAMM to fit to data that is autocorrelated. On of my covariates (covariate 4) is a sin function the represents time of day. It is

Re: [R] accessing the ellements in an array and factors.

2015-06-11 Thread David Winsemius
Cross-posting to SO and Rhelp is deprecated. On Jun 11, 2015, at 8:40 PM, Bogdan Tanasa wrote: Dear all, please could you please with a simple question : I do have an array of 32 elements, where each element is indexed by a name : eg : list_triplet_wells -c(A1:A2:A3, A4:A5:A6 ,

Re: [R] Finding index of specific values in a data.frame

2015-06-11 Thread Kevin Kowitski
I did not realize that there is a coursers assignment similar to this. I am running this for data analysis at work, not for coursers. However I will look through the link you provided and see if it is applicable. Thanks, Kevin Sent from my iPhone On Jun 10, 2015, at 3:21 PM, David Winsemius

Re: [R] graphs, need urgent help (deadline :( )

2015-06-11 Thread Rosa Oliveira
Dear Don, thank you very much. I really wasn’t being able to figure the problem. You were a big (huge) help. Seeing the graphs, I think I’ll try to put the 3 settings (sample size) in different graphs. I’ll try to use trellis graphs :) using sample size as the “factor” Thank you very much

Re: [R] Finding index of specific values in a data.frame

2015-06-11 Thread Kevin Kowitski
Oh I see, I'm sorry I just plopped it in GitHub for ease of help, I didn't notice I put it under coursera work. This task is not related to coursera, I will separate it out. -Kevin Sent from my iPhone On Jun 10, 2015, at 3:21 PM, David Winsemius dwinsem...@comcast.net wrote: On Jun 10,

[R] Quantile regression model with nonparametric effect and interaction

2015-06-11 Thread Waltl, Sofie (sofie.wa...@uni-graz.at)
Dear all, I would like to estimate a quantile regression model including a bivariate nonparametric term which should be interacted with a dummy variable, i.e., log p ~ year + f(a,b):year. I tried to use Roger Koenker's quantreg package and the functions rqss and qss but it turns out that

Re: [R] graphs, need urgent help (deadline :( )

2015-06-11 Thread Don McKenzie
Thanks John! My eyes aren't good enough to see that. I actually checked (I thought). This was the default window on Mac console, for others who might care. Sent from my iPad On Jun 10, 2015, at 6:17 PM, John Kane jrkrid...@inbox.com wrote: You have curly quotes rather than plain ones here

[R] accessing the ellements in an array and factors.

2015-06-11 Thread Bogdan Tanasa
Dear all, please could you please with a simple question : I do have an array of 32 elements, where each element is indexed by a name : eg : list_triplet_wells -c(A1:A2:A3, A4:A5:A6 , A7:A8:A9, A10:A11:A12 ) xxx -array(0, dim=4) dimnames(xxx) = list(list_triplet_wells) From another script, I

Re: [R] Boxplot function error-help required

2015-06-11 Thread John Kane
Thanks Don, I suspected there was a Boxplot() out there by was too lazy to look. I still don't see how the original code would work if label only had onevalue but I must admit what Boxplot() is actually doing is still confusing me. John Kane Kingston ON Canada -Original Message-

Re: [R] for loop incorrect row count

2015-06-11 Thread Bert Gunter
Oh, Swami, gazing into the crystal ball one can see ... ;-} Cheers, Bert Bert Gunter Data is not information. Information is not knowledge. And knowledge is certainly not wisdom. -- Clifford Stoll On Thu, Jun 11, 2015 at 4:48 PM, David Winsemius dwinsem...@comcast.net wrote: On Jun 11,

Re: [R] for loop incorrect row count

2015-06-11 Thread John Kane
Hi Kevin, I don't even pretend to try to read it. If David is having a problem I am going for tea. I think we need some sample data (see ?dput for the best way to supply it) and an succinct description in English of what the problem is and what you need to get out of the data. I suspect

Re: [R-es] Duda glmer

2015-06-11 Thread Guillermo.Vinue
Muchas gracias por tu respuesta y el código Olivier, me es de mucha ayuda. Un saludo, Guillermo Guillermo, me temo que en tu simulación, el enfoque multinivel carezca de sentido. Ten en cuenta que en este tipo de modelo la agrupación de los datos, es decir el segundo nivel en la jerarquía,

[R] Issue with mcapply

2015-06-11 Thread AROONALOK PYNE
Please check this code : library(parallel) workerFunc - function(n) { return(n^2) } a - function(){ CurrentTime - Sys.time() res - lapply(values, workerFunc) TimeTaken - Sys.time() - CurrentTime print(TimeTaken) } b - function(){ CurrentTime - Sys.time() numWorkers - detectCores()

[R] Boxplot function error-help required

2015-06-11 Thread Kris Singh
Dear Sirm/Madam, Just wondering if someone could help me. I've tried running a code on R and the code includes the following: Boxplot(~Acc_S$Acc, label=Acc_S$Subj) But I receive the following error message: *Error: could not find function Boxplot* I have tried installing all the packages

Re: [R] graphs, need urgent help (deadline :( )

2015-06-11 Thread Jim Lemon
Rosa Oliveira wrote: Dear Jim, when I run your code (even the one you send me, not in my data), I get: Don't know how to automatically pick scale for object of type function. Defaulting to continuous Error in data.frame(x = c(0.1, 0.2, 0.1, 0.2, 0.1, 0.2, 0.1, 0.2, 0.1, : arguments

Re: [R] Quantile regression model with nonparametric effect and interaction

2015-06-11 Thread Waltl, Sofie (sofie.wa...@uni-graz.at)
The idea is to move from regional dummies interacted with time dummies (model 1) to a smooth spline (defined on longitudes and latitudes) interacted with time dummies (model 2), i.e., Model 1: Log p ~ X\beta + REGION*YEAR Model 2: Log p ~ X\beta + f(long, lat)*YEAR Estimating the f's in a loop

[R] poLCA - Latent Class Analysis - How long should analysis take?

2015-06-11 Thread Steven King
This is the script i am working from. library(poLCA) f - cbind(bq70,bq72_1,bq72_2,bq72_3,bq72_4,bq72_5,bq72_6,bq72_7,bq73a_1,bq73a_2,bq73a_3,bq73a_4)~ zq88+zq89+dm_zq101_2+dm_zq101_3+dm_zq101_4+dm_zq101_5+dm_zq101_6+dm_zq101_7+dm_zq101_8+dm_zq101_9 for(i in 2:14){max_II--100 min_bic-10

[R] Problem with Volcanoplot

2015-06-11 Thread Pierlot Gabin
Hi all,  I have a data frame composed by 25 numerical variables. I want to do a Spearman Correlation Volcano plot (i. e. x = correlation coefficient and y = -log10(p value)) I'm a begginer in R, so how can I do this ? PS : Sorry for my English, this is not my mother tongue. Thank you !

[R] Problem in ksvm prediction

2015-06-11 Thread bharath
Dear all, I am using ksvm in R for SVM classification. I trained the SVM classifier with prob.model=TRUE. I have the two following outputs when I predicted. 1) label-predict(model, TestData) I have classification accuracy about 93%, where as in another mode 2) Probvalues-predict(model,

Re: [R] Boxplot function error-help required

2015-06-11 Thread Richard M. Heiberger
R is case sensitive. try boxplot not Boxplot On Thu, Jun 11, 2015 at 7:20 AM, Kris Singh kris.si...@research.uwa.edu.au wrote: Dear Sirm/Madam, Just wondering if someone could help me. I've tried running a code on R and the code includes the following: Boxplot(~Acc_S$Acc,

Re: [R] Problem with Volcanoplot

2015-06-11 Thread Bert Gunter
Start by going through an R tutorial or two? You need to do some minimal homework BEFORE posting here. Cheers, Bert Bert Gunter Data is not information. Information is not knowledge. And knowledge is certainly not wisdom. -- Clifford Stoll On Thu, Jun 11, 2015 at 8:05 AM, Pierlot Gabin