私はテスト統計を手で計算しましたが、私の答えが正しいことを確認します。私はRには新しく、これは非常に簡単だと思いますが、ベクトル値を割り当てた後はどこから始めたらいいかわかりません。RでWald/Score/Likelihood Ratio Testを1つのベクトルだけで実行するにはどうすればよいですか？
大きなサンプルの場合、Wald検定統計、スコア検定統計、および尤度比検定統計値は、ほぼchisquare df = 1の分布を持ちます。 P（Y = 1）=πとP（Y = 0）= 1-πとすると、Y = 2項（1、π）観測値Yi、i = 1〜mは、Yとして表示されます。
m = 10の場合、データ（0,1,0,0,1,0,0,0,0,1,0）帰無仮説pi = 0.3を検定する。
免責事項：これは宿題関連ですが、Rパートは宿題の一部ではありません。私は私の答えを確認するためにRでこの問題を実行できるかどうかについて興味があります。
glmの定数項を表すには "1"を入れる必要があります。
x=c(0, 1, 0, 0, 1, 0, 0, 0, 1, 0)
glm(x~1)
これはデータなし、わずか切片でGLMをフィッティングすることと等価です。あなたはlogg-オッズであなたの帰無仮説を表し、オフセット項を設定する必要が
これはpi = 0.5のテストになります –
ありがとう、これは非常に意味があります！ – FettahP
はスケール：
glm(dat ~ offset(rep(log(0.3/0.7),10)), family=binomial)
Call: glm(formula = dat ~ offset(rep(log(0.3/0.7), 10)), family = binomial)
Coefficients:
(Intercept)
1.17e-12
Degrees of Freedom: 9 Total (i.e. Null); 9 Residual
Null Deviance: 12.22
Residual Deviance: 12.22 AIC: 14.22
切片推定値が観測さオッズはtehので値を仮定しているので、本質的に0です。しかし、出来事の数は非常に少なく、結果が何であっても従来のレベルでの統計的有意性は確立されないであろう。 p値は、データが仮説化された値に適合しない確率であるが、この場合はそうである。
'？prop.test'と'？glm' ...を確認しますか？ –
あなた自身でRでこれを行うにはどのような調査をしましたか？あなたはこれをチェックすべきです... http://stackoverflow.com/help/how-to-ask –
私は従属変数と少なくとも1つの独立変数を持っているときに 'glm'を実行することができます。テスト統計。しかし、ベクターが1つしかない場合、どのようにして統計情報を取得できますか？ – FettahP