confidence-interval

7熱

1答えて

bbmleのパッケージmle2のコマンドを使用しようとしています。私はBolkerの "bbmleパッケージでの最尤推定と解析"のp2を見ています。どういうわけか私は正しい開始値を入力することができません。 l.lik.probit <-function(par, ivs, dv){ Y <- as.matrix(dv) X <- as.matrix(ivs) K <-ncol(X) b

4熱

1答えて

シミュレーションデータからのプロットに信頼区間を追加するR

私は尤度関数とシミュレーションに基づいてプロビットシミュレーションを作成しました。これらはすべて以下のコードで複製できます。これは尤度関数である： my.probit <- function(y,x) { # use OLS to get start values par <- lm(y~x)$coefficients ytilde <- 2*y-1 # Run optim res <

3熱

1答えて

ブートストラップの信頼区間R

私は新しいRユーザーです。ブートパッケージの使用に問題があります。 x <- rnorm(100, 1, .5) 任意のヒント：私がやりたいすべてのような、数のベクトルの平均の周りに信頼区間を生成するためにブートストラップを使用しているのですか？

1熱

1答えて

信頼区間の値を計算し、それに応じて曲線をプロットする

私はプロットしている[x、y]データのペアを持っています。私は、最小二乗回帰モデルを使ってカーブのための最良のフィットラインを計算し、それに応じてy値を予測しています。これで、95％信頼区間の値を生成し、それをプロットする必要があります。誰かが信頼区間のy値を計算するための公式を私に助けることができたら、それは素晴らしいことです。私はここに投稿する前にグーグルで試した。しかし、サイトのほとんどは

2熱

1答えて

r/sciplot：lineplot.CIのウィスカーが重複する

sciplotライブラリでlineplot.CIを使用して相互作用プロットを作成すると、エラーバーが複数のグループにまたがることがあります。例えば、 data = c(1,5,3,7,3,7,5,9) grp1 = c(1,1,1,1,2,2,2,2) grp2 = c(1,1,2,2,1,1,2,2) lineplot.CI(grp1, data, grp2) グループがグループ化変数

4熱

1答えて

既定で分数回帰のブートストラップによって信頼区間を得る方法

私の統計クラスでは、私はStataを使用しています。私はRユーザーなので、Rで同じことをしたいのですが、正しい結果が得られました。信頼区間と同じくらい単純なものを得るのはやや厄介なやり方です。ここに私の粗溶液です： library(quantreg) na = round(runif(100, min=127, max=144)) f <- rq(na~1, tau=.5, data=ds)

2熱

1答えて

正規分布の標本から複数の信頼区間を生成するR

私は統計学の学生であり、割り当ての一部としてランダムに生成された正規分布の標本に対して複数の信頼区間を生成しようとするR初心者（年の過小評価）です。 Iは、（二重マトリックス中）N（50、6^2）分布から25個のサイズのサンプルをN = 20を生成する機能を data <- replicate(25, rnorm(20, 50, 6)) を用います。私の質問は、この配布のサンプルごとに95％