Re: [R] creating dummy variables

R is for dummies (like me, but I don't use dummy variables) or for the non-Dummies like all experts who help us all the time@@.. so dummy variables are not needed! :) QED... On Sat, Apr 20, 2013 at 6:16 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote: On 21/04/13 10:56, Eva Prieto Castro wrote:

[R] cluster gene list

Hi, I have created a heatmap using heatmap.2 having 7 clusters. I would like to extract the list of genes that are in these 7 clusters. Is there any function that can be used to extract genes for each cluster? Cheers, Sudhir -- __

Re: [R] Welcome to the R-help mailing list

I have a question and please help. In SPSS, we normally perform data screening to eliminate bad items before doing factor analysis, for instance, those negatively affect the reliability of the scale. do we have to do the same with R, if yes, how? thank you very much! Subject: Welcome to the

Re: [R] Welcome to the R-help mailing list

On Sun, Apr 21, 2013 at 4:04 AM, 李究 li...@hotmail.com wrote: You must know your password to change your options (including changing the password, itself) or to unsubscribe without confirmation. It is: [[*snip*]] Not an answer to your question, but I'd like to point out that you just sent

[R] How to set frequncy?

Hi Boss, i have data of one stock price from 2 January 2006 to 31 December 2012.But each year had different markert day, I want to calculate daily returns and display the series,so how to set about the frequency? If I set frequency=260,is that OK? [[alternative HTML version deleted]]

Re: [R] Excluding observations

Hi Farnoosh, set.seed(25) Data.All.Var- data.frame(ID= sample(LETTERS[1:10],25,replace=TRUE), value=rnorm(25),stringsAsFactors=FALSE) set.seed(28) Y.Out- data.frame(ID=sample(LETTERS[1:5],25,replace=TRUE),value=rnorm(25),stringsAsFactors=FALSE) unique(Y.Out$ID) #[1] A C E D B  

[R] Fwd: Bug in VGAM z value and coefficient ?

Dear, I dont see my topic in the daily summary : Thanks -- Forwarded message -- From: Olivier Merle oliviermerl...@gmail.com Date: 2013/4/17 Subject: Bug in VGAM z value and coefficient ? To: r-help@r-project.org Dear, When i multiply the y of a regression by 10, I would

[R] Pls help to prevent my post from being indexed on google

Hi is it possible to add this line to my earlier post meta name=robots content=noindex to prevent it from being indexed by google? the post is at https://stat.ethz.ch/pipermail/r-help//2013-April/350857.html From: boon_lo...@hotmail.com To: r-help@r-project.org Subject: Help for bootstrappingþ

[R] cedta decided 'igraph' wasn't data.table aware

Hi, what does this mean? --8---cut here---start-8--- graph - graph.data.frame(merged[!v,], vertices=ve, directed=FALSE) cedta decided 'igraph' wasn't data.table aware cedta decided 'igraph' wasn't data.table aware cedta decided 'igraph' wasn't data.table aware

Re: [R] cedta decided 'igraph' wasn't data.table aware

cedta (Calling Environment is Data.Table Aware) is a private function in package:data.table. It looks like it only prints those comments if option(datatable.verbose = TRUE) has been set. Ask the data.table maintainer for more details maintainer(data.table) [1] M Dowle

Re: [R] Pls help to prevent my post from being indexed on google

Hi, No that is not really possible for several reasons. One is that these posts are not archived on just one web server. There are multiple archives managed by different people/places. Another is that they are actually scrubbed to plain text, so my guess is that bots would miss that tag anyway

[R] double exponential regression R

Hello all! I have a problem with a double exponential equation. This are my data's structure(list(proc = c(1870.52067384719, 766.789388745793, 358.701545859122, 237.113777545511, 43.2726259059654, 148.985133316262, 92.6242882655781, 88.4521557193262, 56.6404686159112, 27.0374477259404,

[R] how to import several files every day

Hi All I want to import several .dat files every day of the week from the same folder. Let say, on day 1, I have about 100 files in the folder. By using this code, everything works perfectly (maybe there is a more efficient way to do it): filenames -list.files(path=pathtofile, full.names=TRUE)

Re: [R] double exponential regression R

Dear Catalin, I think the package dcemriS4 has this capability. Hope it helps, Jose From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of catalin roibu [catalinro...@gmail.com] Sent: 21 April 2013 20:34 To: r-help@r-project.org

Re: [R] how to import several files every day

On 04/22/2013 05:36 AM, Martin Lavoie wrote: Hi All I want to import several .dat files every day of the week from the same folder. Let say, on day 1, I have about 100 files in the folder. By using this code, everything works perfectly (maybe there is a more efficient way to do it):

[R] lsoda question from deSolve package

Dear List, Wonder if you have some thoughts on the following question using lsoda in desolve: I have the following data and function: require(deSolve) times - c(0:24) tin  - 0.5 D - 400 V    - 26.3 k -0.056 k12  - 0.197118 k21  - 0.022665 yini - c(dy1 = 0,dy2 = 0)  events - data.frame(var

Re: [R] how to import several files every day

Hi Jim thanks for the help. I think I understand all the steps you suggested. I tried an example with your code but there is something not working I think. See below. all.filenames-list.files(**path=pathtofile, full.names=TRUE) *#this step works* old.filenames-read.table(**old.filenames.tab)

Re: [R] double exponential regression R

Hi CR, Assuming your data are stored in d, perhaps something like this: m - nls(log(proc) ~ k + a*(b^cls), start = list(a = 2, b = .25, k = 0), data = d) summary(m) # gives ## Formula: log(proc) ~ k + a * (b^cls) ## Parameters: ## Estimate Std. Error t value Pr(|t|) ## a 5.246170.50905

Re: [R] Excluding observations

Thanks a lot. That worked:-) Cc: R help r-help@r-project.org Sent: Sunday, April 21, 2013 7:32 AM Subject: Re: Excluding observations Hi Farnoosh, set.seed(25) Data.All.Var- data.frame(ID= sample(LETTERS[1:10],25,replace=TRUE),

Re: [R] Pls help to prevent my post from being indexed on google

See also: fortune(celebrity) cheers, Rolf Turner On 22/04/13 07:32, Joshua Wiley wrote: Hi, No that is not really possible for several reasons. One is that these posts are not archived on just one web server. There are multiple archives managed by different

Re: [R] matrix of size 30^5

Benjamin Caldwell btcaldwell at berkeley.edu writes: Dear R helpers Reproducible example: #warning - this causes a hard freeze on the machines I've tried it on matrix.holder- matrix(rnorm(150), nrow=30, ncol=5) Out= expand.grid(matrix.holder[,1],matrix.holder[,2],matrix.holder[,3],

Re: [R] Reshape or Plyr?

Hi Bruce, From your second email, I think you solved the problem.  But, still there is confusion due to the wordings. res1-mutate(ddply(dat1,.(SPEC_CODE),function(x) colSums(x[,-c(1:3,6)])),RA=10*(AI/Survey_Time))  res1    SPEC_CODE Survey_Time AI RA 1   Buzz    72.8  8 

Re: [R] Reshape or Plyr?

HI, dat1- structure(list(Location_name = c(079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna, 079-f2p1-Acetuna,

[R] Using copulas with user-defined marginal functions

I am trying to make a loglikelihood function using copulas. I am trying to use mvdc to find the density function. When I run this I got the error that the pdf and cdf of my function tobit doesn't exist. Can somebody guide me where my mistake is? dtobit - function(beta,sigma, x, y) {ifelse(y0,

Re: [R] Reshape or Plyr?

Hi Bruce It sounds like the aggregate function is what you need. Here's an example using the ChickWeight dataset in R: chick.agg-aggregate(ChickWeight[,c('weight', 'Time')], by = ChickWeight[c('Chick','Diet')], sum) weight.time-chick.agg$weight/chick.agg$Time

Re: [R] Reshape or Plyr?

Hi Bruce, May be this helps. library(plyr)  res-mutate(ddply(dat1,.(SPEC_CODE,Location_name),function(x) colSums(x[,-c(1:3,6)])),RA=10*(AI/Survey_Time))  res #   SPEC_CODE    Location_name Survey_Time AI RA #1   Buzz 079-f2p1-Acetuna    72.8  8  1.0989011 #2 Eumspp

Re: [R] NAMESPACE and imports

Dear Rxperts, I am trying to use pdflatex from Hmisc, I get the error message. pdflatex is not available... The version of R is version platform i386-w64-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 15.2 year

Re: [R] NAMESPACE and imports

Please don't hijack a thread... start a new one with a new subject. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] Using copulas with user-defined marginal functions

Readers of this list should have been advised that this question has also been posted to StackOverflow: (And the unnamed package from which dMvdc is assumed to be coming from is 'copula'.) -- David On Apr 21, 2013, at 4:12 PM, Elisa Taveras Pena wrote: I am trying to make a loglikelihood

Re: [R] lsoda question from deSolve package

Hi Andreas, adding a varying amount of something over time is even easier than implementing an immediate event. It can be implemented as a forcing, see ?forcings and example below. Hope it helps, Thomas require(deSolve) pkmod - function(t, y, p) { inp - forc(t) if (t tin) R -