My stupdity I made a late edit for clarity and forgot to run it to be sure I
had changed everything.
It should read:
tabx <- ddply(meltx, .(Seamount, variable), summarize, mean = mean(value),
sd = sd(value))
My appologies.
John Kane
Kingston ON Canada
> -----Original Message-----
> From: smartpink...@yahoo.com
> Sent: Sat, 8 Sep 2012 11:14:11 -0700 (PDT)
> To: jrkrid...@inbox.com
> Subject: Re: [R] Producing a table with mean values
>
> Hi John,
>
> I am getting error messages with your solution. tabx <- ddply(nn,
> .(Seamount, variable), summarize, mean = mean(value), sd = sd(value))
> #Error in empty(.data) : object 'nn' not found
>
> A.K.
>
>
>
> ----- Original Message -----
> From: John Kane <jrkrid...@inbox.com>
> To: Tinus Sonnekus <tsonne...@gmail.com>; r-help@r-project.org
> Cc:
> Sent: Saturday, September 8, 2012 1:19 PM
> Subject: Re: [R] Producing a table with mean values
>
> x <- " Seamount Pico Nano Micro Total_Ch
> Off_Mount1 0.0691 0.24200 0.00100 0.31210
> Off_Mount1 0.0938 0.00521 0.02060 0.11961
> Off_Mount1 0.1130 0.20000 0.06620 0.37920
> Off_Mount1 0.0864 0.15900 0.22300 0.46840
> Off_Mount1 0.0262 0.04570 0.00261 0.07451
> Off_Mount2 0.0314 0.17400 0.12800 0.33340
> Off_Mount2 0.0314 0.17400 0.12800 0.23340
> Off_Mount2 0.0414 0.17400 0.02800 0.23340"
>
> xx <- read.table(textConnection(x), header=TRUE, as.is=TRUE)
>
> library(reshape)
> meltx <- melt(xx)
>
> tabx <- ddply(nn, .(Seamount, variable), summarize, mean = mean(value),
> sd = sd(value))
>
> tabx
>
> John Kane
> Kingston ON Canada
>
>
>> -----Original Message-----
>> From: tsonne...@gmail.com
>> Sent: Fri, 7 Sep 2012 22:49:55 +0200
>> To: r-help@r-project.org
>> Subject: [R] Producing a table with mean values
>>
>> Hi All,
>>
>> I have a data set wit three size classes (pico, nano and micro) and 12
>> different sites (Seamounts). I want to produce a table with the mean and
>> standard deviation values for each site.
>>
> > Seamount Pico Nano Micro Total_Ch
>> 1 Off_Mount 1 0.0691 0.24200 0.00100 0.31210
>> 2 Off_Mount 1 0.0938 0.00521 0.02060 0.11961
>> 3 Off_Mount 1 0.1130 0.20000 0.06620 0.37920
>> 4 Off_Mount 1 0.0864 0.15900 0.22300 0.46840
>> 5 Off_Mount 1 0.0262 0.04570 0.00261 0.07451
>> 6 Off_Mount 2 0.0314 0.17400 0.12800 0.33340
>>
>> I tried the following script but get an error message
>>
>> *Error in results[i, "u.Pico", "u.Nano", "u.Micro"] <- sapply(z, mean) :
>> *
>> * incorrect number of subscripts *
>>
>> The code I used:
>>
>> *SChla <- read.csv("SM_Chla_data.csv")*
>> *sm <- as.character(unique(SChla$Seamount))*
>> *
>> *
>> *results <-
>> matrix(NA,nrow=length(sm),ncol=6,dimnames=list(sm,c("u.Pico","u.Nano","u.Micro","sd.Pico","sd.Nano","sd.Micro")))
>> *
>> *
>> *
>> *for (i in sm){*
>> *z <- subset(SChla, Seamount==i, select=c(Pico, Nano, Micro))*
>> *results[i,"u.Pico","u.Nano","u.Micro"] <- sapply(z, mean)*
>> *results[i,"sd.Pico","sd.Nano","sd.Micro"] <- sapply(z, sd)*
>> *}*
>> *
>> *
>> *print(results)*
>>
>> Please can some one advise me how to fix the error or maybe have an
>> alternative solution I will appreciate it.
>>
>> Thank you.
>> Tinus
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help@r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ____________________________________________________________
> GET FREE SMILEYS FOR YOUR IM & EMAIL - Learn more at
> http://www.inbox.com/smileys
> Works with AIM®, MSN® Messenger, Yahoo!® Messenger, ICQ®, Google Talk™
> and most webmails
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
____________________________________________________________
FREE 3D EARTH SCREENSAVER - Watch the Earth right on your desktop!
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
