R的帮助页面的例子运行成功了,但是当我把参数中“y”代入矩阵时,就提示下标出界,就是用帮助页面中的那组数据也不行。
R帮助页面中使用数据是hgtf。提供的代码是:
age <- c( seq(1, 2, 0.25), seq(3, 8, 1), seq(8.5, 18, 0.5))rng <- c(1,18)
norder <- 6
nage <- length(age)
nbasis <- nage + norder - 2
wbasis <- create.bspline.basis(rng, nbasis, norder, age)
cvec0 <- matrix(0,nbasis,1)
Wfd0 <- fd(cvec0, wbasis)
Lfdobj <- 3
lambda <- 10^(-0.5)
growfdPar <- fdPar(Wfd0, Lfdobj, lambda)
wgt <- rep(1,nage)
result <- smooth.monotone(age, hgt1, growfdPar, wgt, conv=0.1)
最后一步里hgt1换成hgtf[,1],hgtf[,5]都可以运行并给出结果。
但是直接用hgtf,或者hgtf[,1:2]都会提示“错误于coef0[, icurve] : 下标出界”。
刚开始学FDA,希望有前辈指导一下,谢谢


雷达卡


京公网安备 11010802022788号







