[R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
Thanks, It seems for loop spends less time ;) with dim(DataFrame) [1] 338 70 For loop has user system elapsed 0.012 0.000 0.012 and apply has user system elapsed 0.020 0.000 0.021 Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ On Wed, Dec

Re: [R] rle with data.table - is it possible?

2014-12-31 Thread Jeff Newmiller
I do not understand the value of using the rle function in your description, but the code below appears to produce the table you want. Note that better support for the data.table package might be found at stackexchange as the documentation specifies. x - read.table( text= Dad Mum Child Group

Re: [R] Problem in R2WinBUGS

2014-12-31 Thread Berend Hasselman
You should reply to the R-help list and not only to me . yo in your answer to me does not correspond with the original mail. It’s not possible to make sense of your code. Berend On 31-12-2014, at 10:14, thanoon younis thanoon.youni...@gmail.com wrote: thank you for your response N-2000;

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Martin Morgan
On 12/31/2014 12:22 AM, Karim Mezhoud wrote: Thanks, It seems for loop spends less time ;) with dim(DataFrame) [1] 338 70 For loop has user system elapsed 0.012 0.000 0.012 and apply has user system elapsed 0.020 0.000 0.021 The timings are so short that the answer

Re: [R] Moving average

2014-12-31 Thread Achim Zeileis
On Tue, 30 Dec 2014, jim holtman wrote: Try this: x - read.csv(text = Date,Open,High,Low,Close,Volume,Adj Close + 2014-12-26,162.27,163.09,162.01,162.34,1912200,162.34 + 2014-12-24,162.88,162.99,161.61,161.82,1868100,161.82 + 2014-12-23,162.23,162.90,161.61,162.24,4043300,162.24 +

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
Many Many Many thanks! it is a demonstrative lesson. I need time to test all examples :) Thank you for your time and support. Happy and Healthy New Year Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ On Wed, Dec 31, 2014 at 2:38 PM, Martin Morgan

[R] Saving an edited R function - RStudio and R

2014-12-31 Thread John Sorkin
Windows 7 Colleagues, I used the fix() function to edit an existing function when using RStudio. After editing the function, I am given the option to SAVE the modified function. I would like to know (1) where the modified function is stored (the save button does not have an option to specify

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
Concretely I request cbioportal through cgsdr package. Depending of Cases and Genetic profiles I receive in general data.frame with heterogeneous structure. The bad one if the returned data.frame is composed by numeric and character columns. in this case numeric columns are considered as factor.

Re: [R] Saving an edited R function - RStudio and R

2014-12-31 Thread Sarah Goslee
Hi, On Wed, Dec 31, 2014 at 11:24 AM, John Sorkin jsor...@grecc.umaryland.edu wrote: Windows 7 Colleagues, I used the fix() function to edit an existing function when using RStudio. After editing the function, I am given the option to SAVE the modified function. I would like to know (1)

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
Yes the last one this the best. But I need to test if returned data.frame is with factor or character: cidx - sapply(df, is.factor) or cidx - sapply(df, is.character) Thanks Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ On Wed, Dec 31, 2014 at 5:24 PM, Karim

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
for both cidx - !(sapply(df, is.numeric)) df[cidx] - lapply(df[cidx], as.numeric) Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ On Wed, Dec 31, 2014 at 5:51 PM, Karim Mezhoud kmezh...@gmail.com wrote: Yes the last one this the best. But I need to test if

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread William Dunlap
But this heterogeneity comes even with only supposed numeric data.frame (gene expression). here an example ibrary(cgdsr) GeneList - c(DDR2, HPGDS, MS4A2,SSUH2,MLH1 ,MSH2, ATM ,ATR, MDC1 ,PARP1) cgds-CGDS(http://www.cbioportal.org/public-portal/;) str(getProfileData(cgds,GeneList,

Re: [R] Fwd: which is faster for or apply

2014-12-31 Thread Karim Mezhoud
Thanks, please find what I got: str(getProfileData(cgds,GeneList, stad_tcga_methylation_hm27,stad_tcga_methylation_hm27)) 'data.frame':48 obs. of 10 variables: $ ATM : num NA NA NA NA NA NA NA NA NA NA ... $ ATR : num NA NA NA NA NA NA NA NA NA NA ... $ DDR2 : num 0.714 0.857 0.549

[R] R-Hierarchical binomial proportion modelling

2014-12-31 Thread zaid bhatti
Hi, i am running a model without covariates, along with shape file and neighbourhood matrix, i am facing difficulty to edit this code for 34 provinces, can any help me to edit this below given code. FinalModels_WINScriptToRunModel_AllMetrics v2 final_20141205 inits1-list(precu =

[R] Memory usage problem while using nlm function

2014-12-31 Thread Kushal
Hi, I am trying to do nonlinear minimization using nlm() function, but for large amount of data it is going out of memory. Code which i am using: f-function(p,n11,E){ sum(-log((p[5] * dnbinom(n11, size=p[1], prob=p[2]/(p[2]+E)) + (1-p[5]) * dnbinom(n11, size=p[3],

[R] Help with finding tutors for Linux, R, Perl, Python MATLAB and/or Cytoscape for yeast microarray analysis, next generation sequencing and constructing gene interaction networks

2014-12-31 Thread thomas hahn
​ Help with finding tutors for Linux, R, Perl, Python MATLAB and/or Cytoscape for yeast microarray analysis, next generation sequencing and constructing gene interaction networks Hi I am a visually impaired bioinformatics graduate student using microarray data for my master’s thesis aimed

[R] rounding down with as.integer

2014-12-31 Thread Mike Miller
This is probably a FAQ, and I don't really have a question about it, but I just ran across this in something I was working on: as.integer(1000*1.003) [1] 1002 I didn't expect it, but maybe I should have. I guess it's about the machine precision added to the fact that as.integer always

[R] Can not save files

2014-12-31 Thread John Sorkin
When I try to save a file, I get Permission denied. Can someone let me know why this happens? I am running under Windows 7, Rstudio save(TheResults,file=c:\\data) Error in gzfile(file, wb) : cannot open the connection In addition: Warning message: In gzfile(file, wb) : cannot open

Re: [R] Can not save files

2014-12-31 Thread Duncan Murdoch
On 31/12/2014 6:29 PM, John Sorkin wrote: When I try to save a file, I get Permission denied. Can someone let me know why this happens? I am running under Windows 7, Rstudio save(TheResults,file=c:\\data) Error in gzfile(file, wb) : cannot open the connection In addition: Warning

Re: [R] rounding down with as.integer

2014-12-31 Thread David Winsemius
On Dec 31, 2014, at 3:24 PM, Mike Miller wrote: This is probably a FAQ, and I don't really have a question about it, but I just ran across this in something I was working on: as.integer(1000*1.003) [1] 1002 I didn't expect it, but maybe I should have. I guess it's about the machine

Re: [R] Can not save files

2014-12-31 Thread David Winsemius
On Dec 31, 2014, at 3:29 PM, John Sorkin wrote: When I try to save a file, I get Permission denied. Because your OS does not allow you to save things there. Since you are specifying an absolute path, R is not using your working directory. You may want to read the help pages for:

Re: [R] Can not save files

2014-12-31 Thread Jeff Newmiller
Because the root directory is not a place for saving your files in. Put your files in your Documents directory. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.us

Re: [R] rounding down with as.integer

2014-12-31 Thread Zaid Bhatti
How can I unsubscribe to not receive loop e mails? Sent from my Huawei Mobile David Winsemius dwinsem...@comcast.net wrote: On Dec 31, 2014, at 3:24 PM, Mike Miller wrote: This is probably a FAQ, and I don't really have a question about it, but I just ran across this in something I was

Re: [R] rounding down with as.integer

2014-12-31 Thread David Winsemius
Read the message at the bottom of every message from rhelp. -- David. On Dec 31, 2014, at 8:09 PM, Zaid Bhatti wrote: How can I unsubscribe to not receive loop e mails? Sent from my Huawei Mobile David Winsemius dwinsem...@comcast.net wrote: On Dec 31, 2014, at 3:24 PM, Mike Miller

Re: [R] Can not save files

2014-12-31 Thread John Sorkin
David, I checked getwd() and tried to save to the directory. It did not work: getwd() [1] C:/Users/John/Documents save(TheResults,file=c:\\users\\john\\documents) Error in gzfile(file, wb) : cannot open the connection In addition: Warning message: In gzfile(file, wb) : cannot open compressed

[R] normalmixEM does not produce the same results when run twice using the same data.

2014-12-31 Thread John Sorkin
Windows 7 I am running normalmixEM to fit two normal curves to my data. For some reason, that I don't understand, the results of running the function twice, on the same data, results in two different (although similar) sets of values. I would expect the function run twice on the same data

[R] please un-subscribe my e-mail id to not recieve loop e-mails

2014-12-31 Thread Zaid Bhatti
Best, Zaid Bhatti Sr. Data Analyst Department of Paediatrics and Child Health The Aga Khan University Stadium Road, Karachi-74800, Pakistan. Tel: +92 213 2005466 Cell: +92 (0) 345 2573155 --- This e-mail may contain information

Re: [R] please un-subscribe my e-mail id to not recieve loop e-mails

2014-12-31 Thread Berend Hasselman
As you were told before: read the message at the very end of each message to R-help: R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help Berend On 01-01-2015, at 06:20, Zaid Bhatti zaid.bha...@aku.edu wrote: Best,

Re: [R] Can not save files

2014-12-31 Thread Jeff Newmiller
You are attempting to save a file on top of a directory, which is fortunately protected from your mistake. Try adding a filename to the path you are saving to. save(TheResults,file=c:\\users\\john\\documents\\test.rda)

Re: [R] rle with data.table - is it possible?

2014-12-31 Thread Kate Ignatius
Is it possible to add the following code or similar in data.table: childseg-0 x:=sumchild -0 span-rle(x)$lengths[rle(x)$values==TRUE childseg[x]-rep(seq_along(span), times = spanLOH) childseg[childseg == 0]-'' I was hoping to do this code by SNPEFF_GENE_NAME for mum, dad and child. The problem

Re: [R] rle with data.table - is it possible?

2014-12-31 Thread Kate Ignatius
correct code: childseg-0 x:=sumchild -0 span-rle(x)$lengths[rle(x)$values==TRUE childseg[x]-rep(seq_along(span), times = span) childseg[childseg == 0]-'' On Thu, Jan 1, 2015 at 1:56 AM, Kate Ignatius kate.ignat...@gmail.com wrote: Is it possible to add the following code or similar in

[R-es] survival1

2014-12-31 Thread jbetancourt
Estimados Esto es segunda parte del mensaje anterior, aqui logro que salgan en el output todos los gráficos pero no salen los resultados numericos Saludos Jose cmd_sur.R Description: Binary data surF.R Description: Binary data ___ R-help-es