我需要在fd对象中使用MBD(修改后的波段深度),但它变得越来越复杂。我一直试图从我的fd对象中提取数据,但是虽然在应用某些函数时得到了结果,但我找不到它们。问题是我有两个fd数据库,我需要迭代粘贴它们,所以fd数据没用。有谁知道如何将fd对象转换为fData对象。
非常感谢你。
答案 0 :(得分:0)
右。 fd和fData是功能数据的对象,我在fd对象中有两个功能数据集,比如X和Y.我需要计算这些数据的修改带深度(MBD),包括对另一个的每个观察。我的意思是,对于X中的每条曲线,我必须创建一个功能数据集,将Y中的所有观察结果与X中的观察结合起来。我需要多次这样做,问题是MBD函数仅适用于fData对象。我无法将fd转换为fData。
SE_hat_U<-fd(basisobj=BSpl)
for(i in 1:61){
SE_hat_U$coefs<-cbind(SE_hat_U$coefs, sm[[i]]$fd$coefs)
}
SE_hat_U$fdnames$reps<-(1:61)
X<-SE_hat_U[as.numeric(NT$Id)]
Y<-SE_hat_U[!datos$Id%in%(NT$Id)]
这是我的代码,但我不知道如何将SE_hat_U转换为fData对象而不是fd对象。