Hi,
Try this:
res1<-do.call(cbind,t(sapply(res,c))[,1]) #get the acf values only.? Assuming
that is what you wanted.
?write.csv(res1,"AMSacf.csv",row.names=FALSE)
A.K.
________________________________
From: Zilefac Elvis <zilefacelvis at yahoo.com>
To: arun <smartpink111 at yahoo.com>
Sent: Monday, June 24, 2013 1:00 PM
Subject: Re: [R] Apply acf to data frame containing 'NA'
Hi AK,
How can I write the 'res' output to file? Attached is my data.
dat1<-AMS1[,-1]
pdf("AMS1_acf.pdf")
res<-lapply(dat1,function(x)acf(x,na.action=na.pass))
dev.off()
I tried something like this: ?data<-as.data.frame(t(sapply(res,c))) but it
did not work.
I would like to export the file and store in column mood.
Thanks,
Atem.
________________________________
From: arun <smartpink111 at yahoo.com>
To: Zilefac Elvis <zilefacelvis at yahoo.com>
Sent: Sunday, June 23, 2013 2:33 PM
Subject: Re: [R] Apply acf to data frame containing 'NA'
HI Atem,
No problem.
Arun
________________________________
From: Zilefac Elvis <zilefacelvis at yahoo.com>
To: arun <smartpink111 at yahoo.com>
Sent: Sunday, June 23, 2013 4:30 PM
Subject: Re: [R] Apply acf to data frame containing 'NA'
Hi AK,
Thanks so much for the fruitful correction.
AE.
________________________________
From: arun <smartpink111 at yahoo.com>
To: Zilefac Elvis <zilefacelvis at yahoo.com>
Cc: "barron.hr at gmail.com" <barron.hr at gmail.com>
Sent: Sunday, June 23, 2013 1:11 PM
Subject: Re: [R] Apply acf to data frame containing 'NA'
HI Atem,
________________________________
From: Hannah Barron <hannahbarron at trentu.ca>
To: Zilefac Elvis <zilefacelvis at yahoo.com>
Cc: R help <r-help at r-project.org>
Sent: Sunday, June 23, 2013 2:04 PM
Subject: Re: [R] Apply acf to data frame containing
'NA'
HI,
This makes more sense.? Sorry, I didn't check the ?acf() before I send the
reply.
lapply(dat1,function(x)acf(x,na.action=na.pass))
A.K.
Hi,
The acf function requires that you specify how to deal with NAs in your
data. By default, it uses na.fail ("Error in *na.fail.default*(as.ts(x)) :
missing values in object).
Try using na.action=na.pass as an argument.
While this will pass over your NAs, make sure that you understand how that
will affect your results. See the acf R documentation for further
explanation.
Good luck,
H
On Sun, Jun 23, 2013 at 1:51 PM, Zilefac Elvis <zilefacelvis at
yahoo.com>wrote:
> 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.
>
> ______________________________________________
> 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.
>
>
??? [[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.