Rには関数を積分する関数integrate()がある ときに計算不能になる そんなときは、別な方法が逃げ道になるかもしれない(ならないかもしれない)ので、確認しておく install.packages("sfsmisc") library(sfsmisc) help(integrate.xy) x <- 1:4 integrate.xy(x,…

こちらでintegrate()関数を用いて積分をしている 以下に示すように回帰すると、回帰式の係数は結構、異なるかもしれない 回帰曲線自体は同じくらいなのに 推定係数を比較するとずいぶん違うが、積分すると、大差ない 雑な作りなので、バグがあるかもしれない…

この続き 理想値からある程度のばらつきを入れてデータをシミュレーションし、回帰式推定をする 推定係数がどれくらいばらつくかも見てみよう # シミュレーション par(ask=TRUE) Niter<-100 # 係数格納行列２つ Km1<-Km2<-matrix(0,Niter,4) for(ii in 1:Nit…