Re: [R] EDCR build from source with namespace problem

Thank you very much Roy. I will the wait Luis Tito Le 19/04/2013 15:48, Roy Mendelssohn - NOAA Federal a écrit : Several changes in R and in Java broke the code in EDCR. We have in beta a new version that fixes these problems. I will be certain to let you know when it is ready. -Roy M.

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 <

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 fu

Re: [R] NAMESPACE and imports

Please don't hijack a thread... start a new one with a new subject. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go...

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] 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 079-f2p1-

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 chick.agg.df<-data.frame(chick.

[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(y>0, dn

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-A

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  1.098

Re: [R] matrix of size 30^5

Benjamin Caldwell 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], ma

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 people/place

Re: [R] Excluding observations

Thanks a lot. That worked:-) Cc: R help 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), value=rnorm(25),stringsAsFactors=FALSE) set.seed

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] 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"

[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.fr

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): filenames

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 Sub

[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=TRU

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

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

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 " Bill Dunlap

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

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

Hi is it possible to add this line to my earlier post 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þ Date: Thu, 4 Apr 2013 15:14:05 +080

[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 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 expect that the coefficient

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"  unique

[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] Welcome to the "R-help" mailing list

On Sun, Apr 21, 2013 at 4:04 AM, 李究 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 your passwo

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