科学・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
関連記事
因子分析とは 因子分析では複雑な現象を単純な要因で理解するための手法。集められた個別のデータは結果である目的変数とその目的変数の原因である説明変数に分かれる。目的変数を観察できた結果とし、目的変数に影 …
練習問題-2 勉強時間と定期試験の点数についてデータを作成する。 study_time <- c(1, 3, 10, 12, 6, 3, 8, 4, 1, 5) point <- c(20 …
Structure of Data Analysis 質問を定義する 理想的なデータセットの定義 アクセスできるデータの検索 データの取得 クレンジング Exploratory 統計的予測・モデリング …
rvestでスクレイピングをしようとしたら、%>%という記述が大量頻出。これ何かということで調べてみた。 https://stackoverflow.com/questions/24536154 …
グループ毎に度数分布を作成する by(housing$VAL, housing$SERIALNO, length) グループで平均値を計算するための式 DT[,mean(pwgtp15),by=SEX …
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でのランダムデータ作成 – 上級編にコメントをいただいたので調べてみた。結果を対比するためにプロットも入れた。 […]