[R] Help needed!

[Peter Ehlers]

On 2011-04-20 14:33, Shuangyan Xiong wrote:
> Hi everyone,
>
> I have a question. Now I am reading the resource code of the package
> "ssfcov". The resource code is as following. I cannot find the resource
> code of the function "myss2d" anywhere in the package. Can anyone give
> me a hint how to find it in the package. Thanks a lot!!bv
>
>
>   >  ssfcov
> function (time, x, subject, nbasis = 5, centered = FALSE, noDiag = TRUE)
> {
>       if (!centered) {
>           fit<- smooth.spline(time, x)
>           x<- x - fitted(fit)
>       }
>       gg<- NULL
>       for (zz in unique(subject)) {
>           if (sum(subject == zz)>  1) {
>               tt<- time[subject == zz]
>               xx<- x[subject == zz]
>               g<- expand.grid(t1 = tt, t2 = tt)
>               scov<- xx %*% t(xx)
>               if (noDiag)
>                   scov<- scov + diag(rep(Inf, length(xx)))
>               g$z<- matrix(scov, ncol = 1)
>               gg<- rbind(gg, g[g$z<  Inf, ])
>           }
>       }
>       nobs<- nrow(gg)
>       tt<- min(time) + (max(time) - min(time)) * (1:nbasis)/(nbasis +
>           1)
>       g<- expand.grid(t1 = tt, t2 = tt)
>       g$z<- 0
>       gg<- rbind(gg, g)
>       fit<- myss2d(z ~ t1 * t2, data = gg, id.basis = ((nobs +
>           1):(nobs + nbasis * nbasis)))
>       class(fit)<- "rfcovObj"
>       return(fit)
> }
> <environment: namespace:ssfcov>

Maybe I'm just slow tonight, but where did you get
that package? I can't see on CRAN or R-Forge.

You might also consider a more informative subject
line. "Help needed" is true for all questions (not
answers) on R-help. Just think, why is this list
called R-***help***?

Peter Ehlers

>
>
> Best,
> Shuangyan
>   >
>
> ______________________________________________
> 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.