曲線回帰
こちらでintegrate()関数を用いて積分をしている 以下に示すように回帰すると、回帰式の係数は結構、異なるかもしれない 回帰曲線自体は同じくらいなのに 推定係数を比較するとずいぶん違うが、積分すると、大差ない 雑な作りなので、バグがあるかもしれない…
この続き 理想値からある程度のばらつきを入れてデータをシミュレーションし、回帰式推定をする 推定係数がどれくらいばらつくかも見てみよう # シミュレーション par(ask=TRUE) Niter<-100 # 係数格納行列2つ Km1<-Km2<-matrix(0,Niter,4) for(ii in 1:Nit…