學生回報的問題
library("fda")
data=t(cbind(matrix(growth$hgtm,31,39),matrix(growth$hgtf,31,54)))
ff=supc1(data, rp=0.2, t="dynamic")
Center 的距離很近
dist(ff$centers)
## 1 2 3 4 5
## 2 1.467968e-04
## 3 4.535352e+01 4.535345e+01
## 4 2.667317e+01 2.667321e+01 3.182091e+01
## 5 5.324208e+01 5.324212e+01 9.736920e+01 7.777192e+01
## 6 8.455447e+01 8.455441e+01 3.976392e+01 6.571311e+01 1.370650e+02
## 7 9.341955e+01 9.341955e+01 5.960144e+01 6.810036e+01 1.456070e+02
## 6
## 2
## 3
## 4
## 5
## 6
## 7 4.040460e+01
其實和tolerance有關
ff2 <- supc1(data, rp = 0.2, t = "dynamic", tolerance = 1e-5)
dist(ff2$centers)
## 1 2 3 4 5
## 2 45.35349
## 3 26.67319 31.82091
## 4 53.24210 97.36920 77.77192
## 5 84.55444 39.76392 65.71311 137.06499
## 6 93.41955 59.60144 68.10036 145.60696 40.40460