2016-10-03 2 views
0

私は散布図にプロットするために2列のデータを持っていますが、その点はあたかも軸がシフトされているかのように正しい場所にないようです並べ替えられています。 データはより大きなセットに由来し、計画は "線"を使用してプロットに輪郭線を追加することでした。私は長年Rを使用しておらず、QuickRウェブサイトの例の折れ線グラフに従っています。データはすでにビニングされており、分布カーブを作成する必要があります。計画は、動作している場合には、Rコードをknime内で使用することです。データポイントが私のRプロットの間違った場所にあります

これは、これは1.5、5時時点の私のプロット distribution scatter plot 1 sample from R

ノートで自分のデータとコード

1     | 0.2 |     ? 
2     | 0.3 |0.9040357526438414 
3     | 0.4 |5.411174917770767 
4     | 0.5 |6.680564306410919 
5     | 0.6 |6.209330779324634 
6     | 0.7 |6.064339545114229 
7     | 0.8 |6.689841006370736 
8     | 0.9 |8.755978393214562 
9     | 1.0 |13.11661304772278 
10     | 1.1 |17.093411510054928 
11     | 1.2 |15.479560872141883 
12     | 1.3 |8.420296969382726 
13     | 1.4 |3.5346457261075566 
14     | 1.5 |1.6256413200367157 
15     | 1.6 |2.2848049576096026 
16     | 1.7 |1.2394427974974978 
17     | 1.8 |0.28169014084507044 
18     | 1.9 |     ? 
19     | 2.0 |     ? 
20     | 2.1 |     ? 
21     | 2.2 |     ? 
22     | 2.3 |     ? 
23     | 2.4 |     ? 
24     | 2.5 |     ? 
25     | 2.6 |     ? 
26     | 2.7 |     ? 
27     | 2.8 |     ? 
28     | 2.9 |     ? 

nwells=c(2, 6, 10, 14, 18) #to take the columns I want from the larger dataset 
plot(3, 20, type="n", xlim=c(0, 3), ylim=c(0, 26), xlab="Intensity", 
   ylab="Proportion", xaxs="i", yaxs="i") 
colors <- rainbow(length(nwells)) 
linetype <- c(1:length(nwells)) 
plotchar <- seq(nwells) 

# add lines 

   well <- data1[2] 
    bin <-data1[1] 
    data<-data.frame(bin, well) 
   lines(data, type="p", lwd=1.5, 
     lty=1, col=colors[n], pch=plotchar[n], grid()) 

である - それは、このデータセットの中から来たのでしょうか?

誰でも私が間違っていることは何ですか?

答えて

0

私が原因あなたのデータ列2の要因に数値変数から正しい変換にいくつかの問題があると思い、次はそれを修正する必要があります(ポイント1.5、5はもうありませんではない):

data[,2] <- as.numeric(as.character(data[,2])) 
well <- data1[2] 
bin <-data1[1] 
data<-data.frame(bin, well) 
lines(data, type="p", lwd=2.5, 
     lty=1, col=colors, pch=plotchar, grid()) 

enter image description here

+0

おかげで、私は、数値/因子問題を割れたと思いました。あなたのソリューションはエラーを出します: - >データ[、1] < - as.numeric(as.character(data [、1])) 警告メッセージ: NACは強制で導入されましたが、これは、「?」で表される値が欠落しているためです。 'a < - read.csv(" sorted table.txt "、header = T、sep =" \t "、dec ="。 "、na.strings =("? "))'が問題を解決しました。ありがとう – tech133

関連する問題