経時データ解析の検出力を計算する必要があったので調べました 同一個体内での繰り返し間で変化する因子の効果を調べる場合には（時間に対する傾きや、投与前後での繰り返し測定など）、繰り返し間の相関(rho)が少ないほうが有利です。 alpha = 0.05 power =…

一般線形最小二乗法を使うために行列を固有値分解して、-0.5乗などの行列を求める必要があり、このあたりをどうするかを調べました。 Rには、eigenというFunctionがあります。 > A = matrix(c(1,2,3,2,3,4,3,4,5), ncol = 3) > A [,1] [,2] [,3] [1,] 1 2 3 …