I have data that has been subjected to one of two treatments (solution added or not) and is in a random block design. I have 8 blocks (numbered 1-8) but in the aov, df=1. R didn't recognise there were 8, so I lettered them (B1, B2 etc). Now there is no F value or p value in the two-way aov, although now the df=7.

Originally:

Output

Now:

Output

RI を使用して Python で p 値を計算したいのですが、このパッケージ rpy2 を使用しています。スクリプトを実行すると、Python が予期せず終了し、次のエラー メッセージが表示されます。

「エラー: 'rho' は NULL でない環境でなければなりません: C レベルの eval で検出されました 起動中 - 警告メッセージ:

中止トラップ: 6"

データは以下のとおりです。

以下は私のコードです：

値のペアでPythonでp値を計算するための提案、またはオプションはありますか?

binom オブジェクトが計算されます。

データ: c(94L, 107L) 成功回数 = 94、試行回数 = 201、p 値 = 0.3974 対立
仮説: 真の成功確率は 0.5 と等しくない 95% 信頼区間: 0.3971286 0.5391627
サンプル推定:
成功確率
0.4676617

ただし、p 値の抽出中に、次のエラーが発生します。

ファイル "/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/rpy2/robjects/vectors.py"、233 行目、getitem res = super(Vector, self) 内。getitem (i) TypeError: 'str' オブジェクトはインデックスとして解釈できません

ロジスティック回帰の従属変数 (私のデータ セットでは: dat$admit) を、使用可能なすべての変数、ペア、およびトリオ (3 つの独立変数) で実行したいと思います。各回帰は、異なる独立変数と従属変数で行われます。私が取得したい結果は、各回帰要約のリストです: coeff,p-value ,AUC,CI 95%. 以下に提出されたデータセットを使用すると、7 つの回帰があるはずです。

サンプル データ セットを次に示します (dat$admit はロジスティック回帰従属変数です)。

@marek のコメントによると、出力は次のようになります (女性のみの場合、および female & apcalc から): # Intercept Estimate P-Value (Intercept) P-Value (Estimate) AUC 0.5

女性+apcalc 0.000000e+00 0.000000e+00 1 1 0.5

@David Arenburg が書いた、統計を生成する良いコードがありますが、ペアとトリオのモデル作成はありません。そのため、モデル作成を追加する方法を知りたいです。ここに David Arenburg のコードがありますか?

このリストを作成する方法はありますか? ありがとう、ロン

R では、返される配列の最初の値が p 値であるphyper () 関数を使用して、超幾何分布の p 値を計算できます。

Go / Golang にこの計算を完全に Go 内で実行できるパッケージがあるかどうか疑問に思っていましたか?

私は多くの回帰を実行していますが、特定の変数の係数と p 値への影響のみに関心があります。したがって、私のスクリプトでは、glm サマリーから p 値を抽出できるようにしたいと考えています (係数自体を取得するのは簡単です)。私が知っている p 値を表示する唯一の方法は、summary(myReg) を使用することです。他の方法はありますか？

例えば：

行列として扱ってみましfit$coefficientsたが、p 値を単純に抽出することはまだできません。

これを行うことは可能ですか？

ありがとう！

自由度 n の t 分布の P 値を決定する方法を教えてください。

このテーマに関する調査は、このスタック交換の回答を示しています: https://stackoverflow.com/a/17604216

np.abs(tt) が T 値であると仮定しますが、自由度でどのように機能するのですか?それは n-1 ですか?

前もって感謝します

corと を使用して、2 つのラスター ブリック間のピクセル単位の相関と有意性 (p 値) を取得しようとしていますcor.test。私のデータはここにあります：

• レンガ 1
• レンガ 2

どちらもかなり小さく、全体で 2MB 未満です。

StackOverflow と r-sig-geo に関する以前の議論から、次の 2 つのコード (両方とも Robert Hijmans から) を見つけました。

どちらのコードも、cor相関グリッドを生成する関数で期待どおりに機能します。ただし、 p 値を抽出するために を にcor置き換えてみました。cor.test

次のエラーが発生しました（RStudioのトラックバックを使用）：

以前の r-sig-geo ディスカッションで、ユーザーがこのエラーについて質問しましたが、回答がありませんでした。corそのため、もう一度質問したところ、行列を入力できて入力できないと指摘されたがcor.test、データを数値ベクトルに変換した後でも、問い合わせに対する 1 つの応答を受け取りました。

次のエラーに直面しています：

誰かがこれを手伝ってくれるかどうか疑問に思っていましたか？

私の sessionInfo() は以下の通りです:

ありがとう！

から p 値を抽出する方法についての投稿がたくさんあることは知っていますaov。ただし、数千のサンプルを含むリストがあります。2 つの異なる処理を比較するために各サンプルの aov を実行しましたが、1 つずつコピーできないため、すべての p 値を含むリストを取得する方法を探しています..

これは可能ですか？によって作成されたp値に対してこれを行うのに問題はありませんでしたttest：

データはインポートした .csv であり、[1:3] は列 [4:6] と比較される列を示します

それは実際には問題ではありませんでしたが、aov に対して同様のことを行うことはできないようです:

すべての p 値 (Pr(>F) と呼ばれる) のリストを取得できません。