[R] Apply acf to data frame containing 'NA'

[arun]

HI,
dat1<- read.csv("AMS3.csv",sep="\t",row.names=1)
pdf("AMS3_acf.pdf")
 lapply(dat1,function(x)acf(na.omit(x)))
dev.off()

A.K.






________________________________
From: Zilefac Elvis <zilefacelvis at yahoo.com>
To: arun <smartpink111 at yahoo.com> 
Cc: R help <r-help at r-project.org> 
Sent: Sunday, June 23, 2013 1:51 PM
Subject: Re: [R] Apply acf  to data frame containing 'NA'



Hi,
Thanks so much for the quick reply.
It worked on your data but failed on mine.
Attached is my data. 
I received this error:
"Error in acf(na.omit(x)) : 'lag.max' must be at least 0
In addition: Warning message:
In split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) :
  data length is not a multiple of split variable"

AE.


________________________________
From: arun <smartpink111 at yahoo.com>
To: Zilefac Elvis <zilefacelvis at yahoo.com> 
Cc: R help <r-help at r-project.org> 
Sent: Sunday, June 23, 2013 11:41 AM
Subject: Re: [R] Apply acf  to data frame containing 'NA'


Hi,
May be this helps:
set.seed(24)
A <- matrix(sample(c(NA,rnorm(1500)),1500,replace=FALSE),nrow=500)
par(mfrow=c(ncol(A),1))
 lapply(split(A,col(A)),function(x) acf(na.omit(x)))
A.K.




----- Original Message -----
From: Zilefac Elvis <zilefacelvis at yahoo.com>
To: "r-help at r-project.org" <r-help at r-project.org>
Cc: 
Sent: Sunday, June 23, 2013 1:29 PM
Subject: [R] Apply acf  to data frame containing 'NA'

Hi,
I have a data frame with each column representing a separate site.
Following this code, I can apply acf to all columns in A. 'A' contains randomly generated data.

A <- matrix(rnorm(1500),nrow=500)
par(mfrow=c(ncol(A),1))

lapply(split(A,col(A)), acf)

#OR
lapply(split(A,col(A)), function(snow) acf(snow, lag.max=5))


Problem: Some of my data contains 'NA' in some of the columns. I get this error
"Error in na.fail.default(as.ts(x)) : missing values in object"
How can I resolve this problem?

Thanks.
AE.
    [[alternative HTML version deleted]]

______________________________________________
R-help at 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.