問題タブ [hmisc]

パッケージcut2()の関数を使用して、期間に基づいて因子を作成しようとしています。Hmisc

ここにいくつかのコードがあります：

data.point の各ポイントが時間間隔内に配置されるように、間隔が作成されることを期待していました。しかし、最後にいくつかの NA 値があります。

このオプションminmax=TRUEにより、カットに のすべての値が確実に含まれるようになると期待していましたdata.points。

ここで何が起こっているのか誰でも明確にできますか? cut2 関数を使用して、データ内のすべての値を含む因子を生成するにはどうすればよいですか?

Hmisclatex()関数を使用して植字するために、R データ フレームの行名をグループ化しようとしています。問題は、latex()各行名の前に 2 つのチルダ文字が追加され、これらがドキュメントに表示されることです。

これらの文字を削除するか、表示されないようにするにはどうすればよいですか?

例：

編集：ラテックス関数はknitrチャンク内で発生します。結果の.Rnwファイルは、デフォルトでpdfLatexを使用するknit2pdf関数を介してコンパイルされていると思います。ドキュメント内の他のすべての表/図は、残りの LaTex 構文が表示されることなく、正常にコンパイルされます。

10k 行と 500 列のデータ フレームがあります。列ごとに、行の一意の値ごとにカウントを作成したいと考えています。例えば

生産します：

Hmisc describe 関数はこの種の分析を生成しますが、出力の形式が非常に悪いため役に立ちません。

ありがとう。

ファイル describe.Rnw を次のように作成すると:

そしてそれをコンパイルしてみてください

latex(describe()) の組み合わせで作成されたラテックスに文法上の誤りがあることを示すエラー メッセージが表示されます。

生成されたラテックスが下に貼り付けられているのがわかります-行間が解析を混乱させている可能性がありますか? 私のセットアップ全体は新しく（実際にはテストモードです）、ドキュメントを正常にコンパイルしましたが、knitr は初めてですが、セットアップにも何か問題があるかどうかはわかりません。私の最初の質問は、何か間違ったことをしているのですか、それとももっと一般的な問題ですか??

ここに私のセッション情報があります：

sessionInfo() R バージョン 3.0.1 (2013-05-16) プラットフォーム: x86_64-w64-mingw32/x64 (64 ビット)

編集 - これはknitrによって作成されたdescribe.texドキュメントです。

私の目的は、ある特定の店舗の存在と、それらの店舗を見つけることができる人口密度との間に関連性があるかどうかをピボット テーブルで比較することです。そのために、お店がある地域とない地域の 600 の例を含む CSV ファイルを用意しました。これは 600 行と 2 つの列を持つファイルです: 1/ 1 つの地域の人口密度を表す数字、および 2/ この地域のこの特定の店の数量 (0、1 または 2)。

ピボット テーブルを作成するには、密度をそれぞれ 60 行の 10 グループにグループ化する必要があります (最初のグループでは 60 の大きな密度から、最後のグループでは 60 の小さな密度まで)。そうすれば、密集度が低いか高いかに関わらず、どのくらいの店舗が建てられているかが簡単にわかります。私は理解できますか？:)

本当に難しいことは何もないと思います。しかし、そのために大丈夫な方法（およびパッケージ）がいくつかあります...私は少し迷っています。

私の主な問題:変数をそれぞれ60行の10のグループにグループ化する最も簡単な方法はどれですか? 私は cut()/ cut2()と hist() を試しましたが成功しませんでした。

たとえば（ジャスティンが尋ねたように）。カット（）を使用：

わかりました。確かに「グループ」には、ほぼ同じ行数の 10 個のグループが含まれています。しかし、間隔で示されている特定の値は、私には意味がありません。密度列の最初の行は次のとおりです (徐々にソートされます)。

つまり、最初のグループを見てください。14.9 が最小値であるのに、なぜ 0.492 なのですか? そして、最初の行と値 51.4 の間の行数を手動で数えると、76 行になります。53 行と表示されるのはなぜですか? データフレームが最低から最高に正しくランク付けされていることを正確に示します。

私は確かに何かが恋しいです...しかし、何ですか？

lrmハレルの「rms」パッケージ内の関数を使用して、比例オッズ回帰モデルを適合させました。適合モデルに適用された外部検証データセットの c-index / Somer's D をどのように計算できるのか疑問に思っています。これらの統計は、fitted.model$stats を使用して直接適合モデルに対して計算できますが、外部検証データセットに対してそれらを取得する方法がわかりません。

「rms」パッケージval.prob&のバイナリ応答またはサバイバル応答に対してこれを行うための組み込み関数val.survがあるようですが、序数応答の方法が見つかりません。

R でデータフレームのヒストグラムをプロットしようとしています。Hmisc パッケージを使用して同じことをしています。1 つの問題を除いて、カテゴリ データの非常に優れたプロットが生成されます。確率ではなく頻度を表示します。以下のコードサンプルを見つけてください（パッケージドキュメントから、私のデータはカテゴリカルですが、それは問題ではありません）：

上記のプロットを変更して、正規化された周波数を表示するにはどうすればよいですか?

私はすでに試しhist(dfram, type='density')ましたtype ='density'が、カーネル密度をプロットするための histSpike 関数への引数です。hist(dfram,f=F)andも試しhist(dfram,prob=T)ましたが（基本的には私が理解していることと同じです）、ヒストグラムにはまだ頻度が表示されています。

パッケージを実行latex(rnorm(5),file="")すると、テーブルの両方の内容が表示されます。Hmisc例えば：

ただし、テーブルのプリアンブルとエンディングも表示されます。例えば：

テーブルのプリアンブル/エンディングではなく、フォーマットされたテーブルの内容のみを表示するようにコマンドを強制するにはどうすればよいですか?

ライブラリのaregImpute関数を使用して欠損値を処理するために、R で複数の代入を使用しました。Hmisc分析のために、fit.mult.impute関数を使用して、多重代入プロセスの分散を修正します。

係数と p 値を含むテーブルを作成したいと考えています。fit.mult.impute接尾辞を使用して係数を抽出できますが$coef、p値を抽出することはできません。これどうやってするの？