Re: [R] Model selection: On the use of the coefficient determination(R2) versus the frequenstist (AIC) and Bayesian (AIC) approaches

2013-04-13 Thread Bert Gunter
This is off topic here-- it has nothing to do with R, per se. Post on a statistics list such as stats.stackexchange.com instead. -- Bert On Sat, Apr 13, 2013 at 5:41 PM, Kaptue Tchuente, Armel wrote: > Dear all, > > I'm modeling growth curve of some ecosystems with respect to their > rainfall

Re: [R] Comparison of Date format

2013-04-13 Thread farnoosh sheikhi
Thank you so much. You are the BEST:-) Now, I'm gonna learn all these codes and understand.   Cc: R help Sent: Saturday, April 13, 2013 5:18 PM Subject: Re: Comparison of Date format Hi, DataA<- structure(list(ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,

Re: [R] CCA report "species environment correlation"

2013-04-13 Thread SRuhl
Thank you thank you thank you! In hindsight this seems so easy I cannot believe I couldn´t find it anywhere, Thanks :) -- View this message in context: http://r.789695.n4.nabble.com/CCA-report-species-environment-correlation-tp4664089p4664144.html Sent from the R help mailing list archive at Na

[R] Model selection: On the use of the coefficient determination(R2) versus the frequenstist (AIC) and Bayesian (AIC) approaches

2013-04-13 Thread Kaptue Tchuente, Armel
Dear all, I'm modeling growth curve of some ecosystems with respect to their rainfall-productivity relationship using a simple linear regression (ANPP(t)=a+b*Rain(t)) and a modified version of the Brody Model ANPP(t)=a*(1-exp(-b*rain(t))) I would like to know why the "best model" is function o

[R] help on smoothing volatility surface..

2013-04-13 Thread C
This script below pulls yahoo data via a function in quantmod, then massages the data around to forumalate a 3D graph with RGL library, attached is a ggplot to show the data i'm trying to create a surface with in separate line geoms . the issue is that the 3D graph looks very ugly and cut up becaus

Re: [R] Comparison of Date format

2013-04-13 Thread arun
Hi, DataA<- structure(list(ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L), Status = c("A", "B", "A", "B", "A", "B", "A", "A", "B", "B", "A", "A", "A",

Re: [R] how can I convert a result (text) in table format in R?

2013-04-13 Thread David Winsemius
On Apr 13, 2013, at 12:36 PM, David Winsemius wrote: > > On Apr 13, 2013, at 3:34 AM, Kristi Glover wrote: > >> Hi R user, >> Could you please give me some hints on how I can convert text in table >> format in R? >> I was doing model assessment using dismo package >> for example: >> bg <- ra

Re: [R] how to add a row vector in a dataframe

2013-04-13 Thread arun
Hi, Using S=1000 and simdata <- replicate(S, generate(3000)) #If you want both "m1" and "m0" #here the missing values are 0 res1<-sapply(seq_len(ncol(simdata.psm1)),function(i) {x1<-merge(simdata.psm0[,i],simdata.psm1[,i],all=TRUE); x1[is.na(x1)]<-0; x1}) res1[,997:1000] #  [,1] [,

Re: [R] Reshaping Data for bi-partite Network Analysis [SOLVED]

2013-04-13 Thread sylvain willart
Wow ! so many thanks Arun and Rui works like a charm problem solved 2013/4/13 arun > Hi, > Try this; > library(reshape2) > res<-dcast(Input,people~place,value.var="time") > res[is.na(res)]<-0 > res > # people beach home school sport > #1Joe 53 0 1 > #2 Marc 04

Re: [R] Reshaping Data for bi-partite Network Analysis

2013-04-13 Thread Rui Barradas
Hello, With the following the order of both rows and columns will be different than the order of your example output, but the table is basically the same. xtabs(time ~ people + place, data = Input) Hope this helps, Rui Barradas Em 13-04-2013 22:03, sylvain willart escreveu: Hello I have

Re: [R] Reshaping Data for bi-partite Network Analysis

2013-04-13 Thread arun
Hi, Try this; library(reshape2) res<-dcast(Input,people~place,value.var="time") res[is.na(res)]<-0  res #  people beach home school sport #1    Joe 5    3  0 1 #2   Marc 0    4  2 0 #3   Mary 0    0  4 0 #or  xtabs(time~.,Input) #  place #people beach home s

[R] Reshaping Data for bi-partite Network Analysis

2013-04-13 Thread sylvain willart
Hello I have a dataset of people spending time in places. But most people don't hang out in all the places. it looks like: > Input<-data.frame(people=c("Marc","Marc","Joe","Joe","Joe","Mary"), + place=c("school","home","home","sport","beach","school"), + time=c(2,4,3,1,

Re: [R] Windows RT or Windows 8

2013-04-13 Thread Leandro Marino
Dear Levke, R works perfect with windows 8. I didn't update to the 3.0 version but I guess that will work too. Windows Rt differ from Windows 8 in a lot of things really don't know if it will work. Kindly, Leandro *LEANDRO MARINO* | Showmetech Estatístico | Fotó

Re: [R] how can I convert a result (text) in table format in R?

2013-04-13 Thread David Winsemius
On Apr 13, 2013, at 3:34 AM, Kristi Glover wrote: > Hi R user, > Could you please give me some hints on how I can convert text in table format > in R? > I was doing model assessment using dismo package > for example: > bg <- randomPoints(current, 500) I doubt very much that 'randomPoints' is

Re: [R] CCA report "species environment correlation"

2013-04-13 Thread Gavin Simpson
On Fri, 2013-04-12 at 11:00 -0700, SRuhl wrote: > Hi everyone, > I did a CCA with R in the vegan package and got all the outputs I need/want > to report except for one, that is the "species environment correlation" > values. I only know of them because of one of the sources I read up on about > CCA

[R] Windows RT or Windows 8

2013-04-13 Thread Levke Godbersen
Dear R Experts, will R run on Windows RT or Windows 8?  I understand these platforms to be significantly different from former Windows versions. Thank you for your time and answer.   Kind regards Levke  [[alternative HTML version deleted]]

[R] how can I convert a result (text) in table format in R?

2013-04-13 Thread Kristi Glover
Hi R user, Could you please give me some hints on how I can convert text in table format in R? I was doing model assessment using dismo package for example: bg <- randomPoints(current, 500) pvtest <- data.frame(extract(current, occtest)) avtest <- data.frame(extract(current, bg)) e2 = evaluate(m

Re: [R] Solving an integral in R gives the error “The integral is probably divergent”

2013-04-13 Thread peter dalgaard
On Apr 13, 2013, at 06:47 , Thomas Lumley wrote: > On Sat, Apr 13, 2013 at 6:31 AM, peter dalgaard wrote: > But is it supposed to be t^{-3/2} or t^{-0.5}?? The formula has the former > and the code the latter, and the integral is clearly divergent with the > former. > > > I didn't look at th

Re: [R] model frame and formula mismatch in model.matrix()

2013-04-13 Thread Eva Prieto Castro
Hi Julien, >From my point of view this error may be produced by the fact you are working >with a lot of variables, so the number of characters if you take them into >account as character is bigger than the system can support. Have you tried to test an example with 50 variables, for example?. S