Hello,

r1$df and r2$df don't exist.

Regards,
Pascal



2013/8/13 Ingo Wardinski <ingo@gfz-potsdam.de>

> G'day
> I try do compute some F-statistics of a singular spectrum analysis of a
> timeseries sv
> I run:
> > require(Rssa)
> > s <- ssa(sv)
> > summary(sv)
>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
>  -4.238   2.761   6.594   6.324  10.410  15.180
> > r1  <- reconstruct(s,groups = list(1:5))
> > r2  <- reconstruct(s,groups = list(1:6))
> > SSE_M1 <- sum(residuals(r1)^2)
> > SSE_M2 <- sum(residuals(r2)^2)
> > df.num <- r1$df - r2$df
> > df.den <- r2$df
> >  F <- ((SSE_M2 - SSE_M1) / df.num) / (SSE_M1 / df.den)
> and eventually
> > p.value <- 1 - pf(F, df.num, df.den)
> Error in pf(F, df.num, df.den) :
>   Non-numeric argument to mathematical function
>
> > summary(df.num)
>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
>
> > summary(df.den)
> Length  Class   Mode
>      0   NULL   NULL
> > summary(F)
>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
>
> I need to compute the p.value, but something is going wrong, and I can't
> see what.
> Any help would be very much appreciated
> ingo
>
> ______________________________**________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
> PLEASE do read the posting guide http://www.R-project.org/**
> posting-guide.html <http://www.R-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]

