問題タブ [plotrix]

次の極座標プロットの例を検討してください。

角度 30 と 330、および 150 と 210 (中心から外側) に線を追加したいと思います。line 関数を試してみましたが、うまくいきませんでした。

次の形式のような時系列パネル データセットがあります。ここで、V1 は文字 (ここでは A、B、C...) で、V2 は数値 (5,8,6...) です。V1 には 65 の異なるアイテムがあります (合計約 50000 の奇妙な観測値)。V1=A または B などの条件で V2 をプロットする方法を見つけようとしています。私が簡単にできることはplot(V1[1:1065])、plot([V1:1066:2085])などですが、これを行うための堅牢な方法を探していました。次のようなことを試しました：

しかし、R の 'if' はベクトルを受け入れないため、エラー メッセージが表示されます。"the condition has length > 1 and only the first element will be used"

誰でも何か提案してもらえますか？

テイラー図に複数のモデルをプロットしようとしていますが、コードに少し苦労しています。図を作成することはできましたが（画像を参照）、軸が大きすぎるため軸を縮小し、軸を標準化して1、2、3、4のラベルを付け、相関に目盛りを追加する方法を理解できません-目盛り付き0.1 ごとにメジャー ティック、最大 0.9 までは 0.05 ごとにマイナー ティックを設定したいと考えていました。その後、この時点で 0.95 にメジャー ティックを設定し、0.01 ごとにマイナー ティックを設定しようとしました (これが理にかなっている場合)。上記のヘルプ/アドバイスが役立ちます。「plotrix」パッケージ内で「taylor.diagram」を使用し（ガイドを読みましたが、Rの経験は比較的ありません）、これまでに（やや基本的な）コードを添付しましたが、プロットはかなり乱雑に見えます。ありがとう

私は最近、JESS3 による肥満と飢餓の関係の視覚化を発見しました: http://jess3.com/obesity-and-hunger-infographic/

情報/空間比率が気に入ったので、R で同様のグラフを作成しようとしました (つまり、中央に散布点があり、中央の端にある極座標プロットで、積み上げ棒グラフは 0 から新しい軸で始まります)しかし、plotrix または ggplot2 で異なる polarplot バージョンのスタックを実装する方法が見つかりませんでした。

非常に基本的なバージョンで 2 つの個別のグラフがどのように見えるかを説明するために、mtcars パッケージを使用して簡単な例を作成しました。

可能であれば、新しい軸が開始され、バーが散布図の外に追加されるようにするにはどうすればよいですか?

事前にどうもありがとう、すべてのヒントは大歓迎です!

RStudio で素敵なプロットを生成するスクリプトの作成に成功しました。

abscissa = nrow(LoopVariable) library(plotrix) plotCI(1:abscissa, LoopVariable[ ,13], uiw=LoopVariable[ ,"F_corrected_normed_error"], err="y", main=(LoopVariable[1,"Samples..Sample_ID"]), xlab="", ylab="correct Y axis", pch=19) abline(h=LoopVariable[1,"mean_fraction-modern"], lwd=2, col="red") mtext("no additional error added")

ここまでは順調ですね。次に、このコードを書きました：

png(filename=paste0("I:/C14Data/QualityControl/Johannes/code_output/", RNumberList[i], "_", date, " (R#_yyyy-mm-dd_hh-mm-ss).png"), type="cairo", units="cm", width=abscissa, height=10, pointsize=12, res=96) plotCI(1:abscissa, LoopVariable[ ,13], uiw=LoopVariable[ ,"F_corrected_normed_error"], err="y", main=(LoopVariable[1,"Samples..Sample_ID"]), xlab="", ylab="fraction modern", pch=19) dev.off() これにより、希望どおりにpngファイルBUTがエクスポートされます。次の2行をエクスポート関数に追加して、前のコードとまったく同じに見えるようにするにはどうすればよいですか?: abline(h=LoopVariable[1,"mean_fraction-modern"], lwd=2, col="red") mtext("no additional error added")

「反対の積み上げ棒グラフ」を作成しようとしていますpyramid.plotが、plotrixパッケージからはうまくいくようです。（一部のユーザーにとっては頼りになるソリューションになると思いますが、このグラフィックスにggplot2固執することを望んでいます。）base

残念ながら、制限を非整数値に設定しようとすると、x 軸で奇妙なことをしているようです。制限を自動的に定義させると、それらは整数になり、私の場合は空白が多すぎます。しかし、それらを として定義するとxlim=c(1.5,1.5)、以下の奇妙な結果が生成されます。

ドキュメントから正しく理解している場合、たとえば軸を抑制して後で追加したり、目盛りを定義したりするなど、追加のグラフィカルパラメーターを渡す方法はありません。より柔軟にする方法はありますか?

以下のプロットを作成するために使用される最小限の作業例を次に示します。

他の人が興味を持っている場合に備えて、人口ピラミッドではなく、負の値の積み上げ棒グラフを試しています。上記のコードには、両側に異なる数の棒のセットを持つことを可能にするために使用する「トリック」が含まれています。つまり、マトリックスに空の列を追加することです。最小限でした！

私は次のデータを持っています、

SampleID Pos Dep Pvalues
sample_1 849 62 0.02755358
sample_1 859 63 0.07406833
sample_1 864 63 0.00351564
sample_1 883 60 0.02780868
sample_1 893 58 0.00451450
sample_1 895 58 0.03600795
sample_2 54 66 0.11864407
sample_2 55 67 0.01515152
sample_2 71 91 0.02712367
sample_2 78 97 0.00077325

各バーの上に頻度値を持つ P 値のヒストグラムを生成しました。以下、コードです

pval_at_site <- read.table("samples.pval")
s <- hist(pval_at_site$Pvalues, xlab="Pval",cex=0.8)
text(s$mids,s$counts,s$count,srt=90,位置 = 3、オフセット = 1、cex = 0.6)

ここで、私がやりたいことは、P 値の頻度とともに、各バーの上にサンプル数を追加することです。

たとえば、最初の間隔に 1000 個のデータポイントがあり、これらの値が 20 個の一意のサンプルから取得された場合、プロットの最初のバーの上に「1000,20」と表示する必要があります。

これについてどうすればよいか教えてください。私が自分自身を明確にしたことを願っています。

ありがとう。

私はかなり長い間、plotrix を使用して極座標プロットでエラー バーを作成しようと試みてきました。散布関数は問題ではありませんでしたが、エラー バーを追加できませんでした。たとえば、ラジアン モードで 2 つ目の極座標プロットを使用して線を追加すると、それらをプロットの中心から離して散布から開始することができないためです。ポイント...私の最善の推測は、散布図ポイントのラジアン+エラーバーを使用して2番目の極座標プロットを上に追加し、その上に散布図ポイント-エラーバーを含む白い線を追加することでした。下限までのエラーバー...これは明らかにコードだけでなくグラフでも非常に醜いので、それに対するより良い解決策があるかどうか疑問に思っていました。これまでの私のコード (無意味な例として mtcars を使用) は次のようになります。

カスタマイズされた x 軸を使用して積み上げ面グラフを作成したいと考えています。具体的には、x 軸の目盛りに cex.axis=0.8 のラベルを付け (すべての値に合わせて小さくします)、赤色を使用して x 値の一部を強調表示したいと思います (たとえば、週末平日とは対照的に日付）。

plotrix パッケージの stackpoly は、プロットの作成に優れています。しかし、cex.axis を制御したり、x 軸のプロットを抑制したりすることはできません。

これを実現するための私の2つの試みは次のとおりです。