ダメ出し：シミュレーションの条件設定 - 裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学 http://d.hatena.ne.jp/foo22222/20120115 Nt <- 1000 df1 <- 1 set.seed(101) ts1 <- rchisq(Nt, df1) set.seed(101) # set.seed()で生成される乱数を同じものにする。 ts2 <- rchisq(Nt, df1) t1 <- sample(ts1, 500, replace = FALSE) t2 <- sample(ts2, 500, replace = FALSE) p1 <- pchisq(t1, df1, lower.tail = FALSE) p2 <- pchisq(t2, df1, lower.tail = FALSE) plot(p1, p2) としている