我怎么才能把方程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)