我怎么才能把方程calculation计算得到的pxt数据导出来,放到 45×2的矩阵当中,我只能得到最后sum(s/(1.035596^n))的结果
market_force=c(1:45)
market_force_replicate<-as.data.frame(matrix(market_force,nrow=45,ncol=2))
a<-as.data.frame(matrix(NA,nrow=45,ncol=1))
calculation <- function (data){
pxt<-matrix(NA,nrow=45,ncol=1)
s<-as.data.frame(matrix(NA,nrow=45,ncol=1))
pxt<-as.data.frame(exp(-data^(1-(2*rbeta(45,1000,1000)-1))))
s[1,]<-pxt[1,]
for (i in 2:45)
{s[i,]<-s[i-1,]*pxt[i,]
}
sum(s/(1.035596^n))
}
result_100010000<-sapply(market_force_replicate,calculation)