科学・IＴ・登山の話題
投稿日：2014年1月31日 更新日：
平均が異なるランダムデータを作成する
> x <- rnorm(12, mean=rep(1:3,each=4), sd=0.2) > y <- rnorm(12, mean=rep(c(1,2,1), each=4), sd=0.2)
現在ある分布と回帰分析の結果を利用したランダムデータの作成
newGalton <- data.frame(parent=rep(NA,1e6),child=rep(NA,1e6))
// NAを使ってデータフレームの箱を作る
newGalton$parent <- rnorm(1e6,mean=mean(galton$parent),sd=sd(galton$parent))
// galtonの親データを1e6個作成する
// データはオリジナルデータの平均と標準偏差の正規分布から取得される。
newGalton$child <- lm1$coeff[1] + lm1$coeff[2]*newGalton$parent + rnorm(1e6,sd=sd(lm1$residuals))
執筆者：admin
関連記事
マニュアル conjointパッケージ caFactorialDesign Rでconjointパッケージを利用した例 Rでコンジョイント分析 Rでコンジョイント分析 |極めて個人的なメモ コンジョイ …
Rでヒストグラムはhist()関数を使って描画する。この時与えるデータによりブレーク引数が異なることがわかった。 まず下記がデフォルトのhist()関数のヘルプである。 hist(x, breaks …
ベクトルからパラメータにもっとも近い値を探す方法 > x <- 1:10 > num <- which( abs(x-5.2) == min(abs(x-5.2)) ) & …
factorと数値型の変換でとりあえず区切り文字付き数字の変換方法についてめどは立ったが、そもそも区切り文字付き数字がなぜfactorになるのかがわからない。文字列にするにはパラメータを指定する必要が …
データセットならstrを使う > str(iris) ‘data.frame’: 150 obs. of 5 variables: $ Sepal.Length: num 5.1 4. …
newGalton <- data.frame(parent=rep(NA,1e6),child=rep(NA,1e6)) のように「箱を作る」必要はありません。
newGalton <- rnorm(1e6,mean=mean(galton$parent),sd=sd(galton$parent))
newGalton$child <- predict(lm1, newGalton)
で十分です。
[…] Rでのランダムデータ作成 – 上級編にコメントをいただいたので調べてみた。結果を対比するためにプロットも入れた。 […]