Skip to content

cbind and mean by week

2 messages · Felipe Carrillo, jim holtman

Felipe Carrillo
# 
Hi:
 I have been able to finally crunch my data by
importing it by week(thank you all for your help),but
here we go again..
 Now I'am trying to do it for the the whole year.
Since the
 dataset is huge I'm only making a 3 weeks dataframe.
- I want to get the mean of pd by week
- I want to count the number of days by week and bind
it to the existing dataframe(x)
- I want to remove the NA's by week and bind another
column counting the rows after the NA's are removed.
 Don't think that I'am crazy but I need those fields
 for my annual report. I was trying to use some of
 the functions below but nothing seems to work..One
 of these days I'll be an R pro(we'll see about that)
 Thanks
 
week<-c(28,28,28,28,28,28,28,29,29,29,29,29,29,29,30,30,30,30,30,30,30)
td<-c(0.015,0.012,NA,0.015,NA,0.014,0.014,0.013,0.013,0.013,0.013,NA,0.010,0.013,0.015,0.011,0.014,0.014,0.014,0.014,0.016)
pd<-c(0,0,NA,80,NA,45,223,228,224,488,525,NA,80,417,413,373,914,213,651,521,979)
 x<-data.frame(week,td,pd)
 x
  pmean<- tapply(x[,3], x$week, mean)
    pdmean
    pdmean<- data.frame(pmean)
    pdmean
 xweek <- aggregate(x, list(x$week), mean)
xweek

 Felipe D. Carrillo
  Fishery Biologist
  Department of the Interior
  US Fish & Wildlife Service
  California, USA



      ____________________________________________________________________________________
Be a better friend, newshound, and
jim holtman
# 
Is this what you want?  You have to add "na.rm=TRUE" to the calculations:
week    td  pd
1    28 0.015   0
2    28 0.012   0
3    28    NA  NA
4    28 0.015  80
5    28    NA  NA
6    28 0.014  45
7    28 0.014 223
8    29 0.013 228
9    29 0.013 224
10   29 0.013 488
11   29 0.013 525
12   29    NA  NA
13   29 0.010  80
14   29 0.013 417
15   30 0.015 413
16   30 0.011 373
17   30 0.014 914
18   30 0.014 213
19   30 0.014 651
20   30 0.014 521
21   30 0.016 979
28       29       30
 69.6000 327.0000 580.5714
pmean
28  69.6000
29 327.0000
30 580.5714
Group.1 week     td       pd
1      28   28 0.0140  69.6000
2      29   29 0.0125 327.0000
3      30   30 0.0140 580.5714
On Tue, Mar 25, 2008 at 10:18 PM, Felipe Carrillo
<mazatlanmexico at yahoo.com> wrote: