Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

Spencer: Hmmm Well, I don't get what's going on here -- as.character.default is internal -- but could you method-ize a simple synonym: asChar- function(e,...)UseMethod(asChar) asChar.call - function(e,...)deparse(e,...) asChar.default - function(e,...)as.character(e,...) asChar(xDy) [1]

Re: [R] group by and merge two dataframes

yes thanks, that's correct! here a slight variation inspired by your solution: a cartesian product restricted to non duplicated records to get the logical vector i to be used in the next natural join i-!duplicated(merge(df1$id,df1$item, by=NULL)) merge(df1[i,],df2) thanks Il 08/05/2014

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

On 09/05/2014, 2:41 AM, Bert Gunter wrote: Spencer: Hmmm Well, I don't get what's going on here -- as.character.default is internal -- but could you method-ize a simple synonym: See ?InternalMethods: For efficiency, internal dispatch only occurs on objects, that is those for which

Re: [R] problem in r-code

Wrong list. This is an R list not Bugs. You may want to consult Bugs materials: http://www2.mrc-bsu.cam.ac.uk/bugs/weblinks/webresource.shtml On 8 May 2014 11:36, thanoon younis thanoon.youni...@gmail.com wrote: dear all members is there anyone explain to me the code below and how can i

Re: [R] How can I make this nested loop faster?

Your code is not re-producable. Can you provide a working example using a standard dataset from R? But, you could first try to use compiler package, see ?enableJIT. Another option would be to use doMC/foreach packages if you can run your assignment in the nested loop in parallel, see %dopar%. On

[R] How to examine the parameter values during iteration (nls)

Hi everyone, I am using nls() to run a simple sigmoid curve type of regression and would like to see how the parameter values change through iterations. How can I see those values? Or ideally, can I even extract those values? For example, nls(Response~E0+Emax*Conc/(EC50+Conc), data=data,

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

Ahhh. Thanks Duncan. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is not knowledge. And knowledge is certainly not wisdom. H. Gilbert Welch On Fri, May 9, 2014 at 2:41 AM, Duncan Murdoch murdoch.dun...@gmail.com wrote: On

Re: [R] How to examine the parameter values during iteration (nls)

Read ?nls and note the trace argument. Does nls(...,trace = TRUE) not do what you want? By using ?capture.output you could then capture the trace I would think, but I haven't tried. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information.

Re: [R] DESeq vs DESeq2 different DEGs results

You're on the wrong list. This is more appropriate on the bioconductor mailing list. On Mon, May 5, 2014 at 9:42 AM, Catalina Aguilar Hurtado cata...@gmail.comwrote: Hi, I want to compare DESeq vs DESeq2 and I am getting different number of DEGs which I will expect to be normal. However,

Re: [R] Tukey no factors

HI, May be the below links help: http://r.789695.n4.nabble.com/TukeyHSD-troubles-td1570205.html http://r.789695.n4.nabble.com/TukeyHSD-error-td3051059.html A.K. Hi I'm trying to run a Tukey test on mortality data, where I want to test whether mortality is influenced by the amount of copper

[R] Resquest

Dear Sir, I need Your help,I am working on the comparison of l moments and trimmed l moments .,Now i face a problem in R programming of trimmed l moments. for this purpose i need your help .for example let suppose i have the following data vector abcc

Re: [R] Error during working with wgcna and R

I am also working on co-expression analysis. It seems like there is no way to use TOMsimilarityFromExpr for large datasets. The option 'maxBlockSize' exists for module detection but not for TOMsimilarity? The only solution seems to reduce the dataset On Sunday, July 8, 2012 1:02:54 PM UTC+2,

[R] Pearson III distribution

Please I need your help.I'm intrested to know if there is any R-package for fit a Pearson III extreme value distribution to data.Thanks a lot Noelia [[alternative HTML version deleted]] __

Re: [R] Error during working with wgcna and R

WGCNA maintainer here. When working with a large data set, you have a few options. 1. Without being snarky, the best option is to get (or get access to) a computer with large-enough RAM. Many universities, departments, and other research institutes have computer clusters with nodes with at least

Re: [R] Pearson III distribution

On May 9, 2014, at 8:17 AM, NOELIA LEGAL wrote: Please I need your help.I'm intrested to know if there is any R-package for fit a Pearson III extreme value distribution to data.Thanks a lot You are requested to do some searching yourself before posting. If you have failed to find what you

Re: [R] How to write a loop in which i has to be in POSIX format?

i is an element in zn so replace zn[i] with just i for (i in zn){ treat$su[treat$Zugnacht==as.POSIXct(i, UTC)] - min(treat$Vollzeit[treat$Zugnacht==as.POSIXct(zni, UTC)]) treat$sa[treat$Zugnacht==as.POSIXct(zni, UTC)] - max(treat$Vollzeit[treat$Zugnacht==as.POSIXct(zni, UTC)]) }

[R] transition from depends to imports

Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages uses 'wilcoxsign_test' from the 'coin' package. And this is the only function from the coin package that it uses (directly). (Everything works fine, as long as my package

Re: [R] Help with passing in name of XTS file into write.table as a variable.

It is a simple change. If you don't save the results of getSymbol() into a variable, then it's not available for write.table() to write. Try mydata - getSymbols(ticker,from='1990-01-01') write.table( mydata, {etc}) The help for write.table says that x (the first argument to write.table),

Re: [R] transition from depends to imports

On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages uses 'wilcoxsign_test' from the 'coin' package. And this is the only function from the coin package that it uses (directly).

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

Hi, Duncan: Thanks very much. I used to think that everything in R was a object. Now I know that is.object(quote(x)) is FALSE. (A decade ago, S-Plus asked me if I wanted to save changes to history. I thought, Wow! Do I get to change history? Hadley's Advanced R book

[R] adding rows

Dear useRs, I have a matrix, say el of 30 rows and 10 columns, as el-matrix(sample(1:300),ncol=10) I want to sum up various sets of three rows of each column in the following manner sum(el[c(1,2,3),]) ##adding row number 1, 2 and 3 of each column sum(el[c(2,3,4),])##adding row number 2, 3 and 4

Re: [R] transition from depends to imports

Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch murdoch.dun...@gmail.com writes: On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages uses

Re: [R] transition from depends to imports

On 09/05/2014, 5:39 PM, Andreas Leha wrote: Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch murdoch.dun...@gmail.com writes: On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from depends to imports in one of my

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

Beware of the is.* functions: * is.object() does not test the usual definition of objects * is.vector() does not test the usual definition of vectors * is.numeric() does not work the same way as is.character() or is.integer() * is.Date() doesn't exist * is.nan() doesn't return TRUE for some NaNs

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

Dear Hadley: Thanks for that. Digits are not numbers. Numbers are not data. Data is not information. Information is not intelligence. Intelligence is not knowledge. Knowledge is not wisdom. And your is. warnings are more useful than my trivia here. Spencer On 5/9/2014

Re: [R] transition from depends to imports

Duncan Murdoch murdoch.dun...@gmail.com writes: On 09/05/2014, 5:39 PM, Andreas Leha wrote: Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch murdoch.dun...@gmail.com writes: On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a

Re: [R] adding rows

Hello, Try the following. sapply(1:(30 - 2), function(i) sum(el[i:(i+2), ])) but with number of rows instead of 30. Hope this helps, Rui Barradas Em 09-05-2014 22:35, eliza botto escreveu: Dear useRs, I have a matrix, say el of 30 rows and 10 columns, as el-matrix(sample(1:300),ncol=10) I

Re: [R] adding rows

Dear Rui and Murphy, Thanks for your help. Eliza Date: Fri, 9 May 2014 22:55:27 +0100 From: ruipbarra...@sapo.pt To: eliza_bo...@hotmail.com; r-help@r-project.org Subject: Re: [R] adding rows Hello, Try the following. sapply(1:(30 - 2), function(i) sum(el[i:(i+2), ])) but with

Re: [R] Lattice Histogram with Normal Curve - Y axis as percentages

Just an afterthought if any one really needs to do it again. The crux of the matter is the different limits in the 2 plot x and y scales. It may be easier to accomplish this with a prepanel function to get the limits eg panel.loess Duncan -Original Message- From: