2013-03-10 144 views
1

我无法弄清楚为什么当我绘制它们时,我的价值观不合理。防爆。 x轴0显示值为0.2,x轴0.2显示值为0.6,x轴0.6显示1.2和x轴的值1.2显示值为0.阴谋破坏

R版本2.15.2(2012-10-26) - “捣蛋”在Mac OS X上运行10.6.8

方法秒杀ELISA试剂盒鱼谱浓CONC2

1 Mekebri CNTRL AbraxisRBT abraxis rbt 450 0.0900000  0.09 
2 Mekebri CNTRL AbraxisRBT abraxis rbt 450 0.0700000  0.09 
3 Mekebri CNTRL AbraxisRBT abraxis rbt 450 0.0700000  0.08 
4 Mekebri CNTRL EnviroRBT enviro rbt 450 0.0900000  0.09 
5 Mekebri CNTRL EnviroRBT enviro rbt 450 0.0700000  0.09 
6 Mekebri CNTRL EnviroRBT enviro rbt 450 0.0700000  0.08 
7 Mekebri 0.2 AbraxisRBT abraxis rbt 450 0.1100000  0.12 
8 Mekebri 0.2 AbraxisRBT abraxis rbt 450 0.3000000  0.32 
9 Mekebri 0.2 AbraxisRBT abraxis rbt 450 0.1000000  0.10  
10 Mekebri 0.2 EnviroRBT enviro rbt 450 0.1100000  0.12 
11 Mekebri 0.2 EnviroRBT enviro rbt 450 0.3000000  0.32 
12 Mekebri 0.2 EnviroRBT enviro rbt 450 0.1000000  0.10  
13 Mekebri 0.6 AbraxisRBT abraxis rbt 450 0.1600000  0.16 
14 Mekebri 0.6 AbraxisRBT abraxis rbt 450 0.1800000  0.18 
15 Mekebri 0.6 AbraxisRBT abraxis rbt 450 0.1700000  0.17 
16 Mekebri 0.6 EnviroRBT enviro rbt 450 0.1600000  0.16 
17 Mekebri 0.6 EnviroRBT enviro rbt 450 0.1800000  0.18 
18 Mekebri 0.6 EnviroRBT enviro rbt 450 0.1700000  0.17 
19 Mekebri 1.2 AbraxisRBT abraxis rbt 450 0.9680557 0.963486175 
20 Mekebri 1.2 AbraxisRBT abraxis rbt 450 0.6040148 0.622156567 
21 Mekebri 1.2 AbraxisRBT abraxis rbt 450 0.5665602 0.5849501 
22 Mekebri 1.2 EnviroRBT enviro rbt 450 0.9680557 0.963486175 
23 Mekebri 1.2 EnviroRBT enviro rbt 450 0.6040148 0.622156567 
24 Mekebri 1.2 EnviroRBT enviro rbt 450 0.5665602 0.5849501 

plot(c(0,0.2,0.6,1.2), 
with(mc, tapply(conc2, list(kit,spike), mean,na.rm=T))[1,], 
type="b",lwd=2, 
ylim=c(0,1), 
xlab=expression(paste("Spike, ",mu,"g/L")), 
ylab=expression(paste(mu,"g/L")), 
col="blue") 
points(c(0,0.2,0.6,1.2), 
with(mc, tapply(conc2, list(kit,spike), mean,na.rm=T))[2,], 
type="b",lwd=2, 
ylim=c(0,1), 
col="red") 
legend(0.2,0.8,lty=1,lwd=2, 
col=c("blue","red"),legend=c("Abraxis","EnviroLogix")) 

dput(head(mc)) 
structure(list(method = structure(c(3L, 3L, 3L, 3L, 3L, 3L), .Label = c("Geis", 
"Mag", "Mekebri"), class = "factor"), spike = structure(c(4L, 
4L, 4L, 4L, 4L, 4L), .Label = c("0.2", "0.6", "1.2", "CNTRL"), class = "factor"), 
elisa = structure(c(2L, 2L, 2L, 2L, 2L, 1L), .Label = c("AbraxisBlank", 
"AbraxisRBT", "EnviroBlank", "EnviroRBT"), class = "factor"), 
kit = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Abraxis", 
"Enviro"), class = "factor"), fish = structure(c(2L, 2L, 
2L, 2L, 2L, 1L), .Label = c("blank", "rbt"), class = "factor"), 
spectral = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("450", 
"600", "scan"), class = "factor"), conc = structure(c(1L, 
1L, 1L, 8L, 1L, 2L), .Label = c("> DL", ">DL", "0", "0.01", 
"0.02", "0.03", "0.04", "0.05", "0.06", "0.07", "0.08", "0.09", 
"0.1", "0.11", "0.12", "0.13", "0.14", "0.15", "0.16", "0.17", 
"0.18", "0.19", "0.2", "0.21", "0.22", "0.23", "0.24", "0.25", 
"0.26", "0.27", "0.28", "0.29", "0.294871066", "0.3", "0.308253804", 
"0.31", "0.32", "0.33", "0.34", "0.344304772", "0.35", "0.350277282", 
"0.353189188", "0.359024076", "0.36", "0.360435916", "0.37", 
"0.370993533", "0.37858631", "0.378888547", "0.38", "0.384568909", 
"0.39", "0.4", "0.401289641", "0.41", "0.42", "0.43", "0.44", 
"0.45", "0.46", "0.47", "0.473548535", "0.48", "0.489942496", 
"0.49", "0.5", "0.51", "0.52", "0.53", "0.54", "0.55", "0.56", 
"0.566560247", "0.57", "0.58", "0.5849501", "0.6", "0.604014755", 
"0.61", "0.62", "0.622156567", "0.64", "0.65", "0.66", "0.67", 
"0.69", "0.7", "0.71", "0.73", "0.75", "0.77", "0.78", "0.79", 
"0.8", "0.81", "0.82", "0.83", "0.84", "0.85", "0.87", "0.88", 
"0.9", "0.91", "0.92", "0.963486175", "0.968055663", "0.97", 
"0.98", "1", "1.01", "1.02", "1.03", "1.04", "1.1", "1.11", 
"1.2"), class = "factor"), conc2 = c(0, 0, 0, 0.05, 0, 0)), .Names = c("method", 
"spike", "elisa", "kit", "fish", "spectral", "conc", "conc2"), row.names = c(NA, 
6L), class = "data.frame") 
+2

请重复此操作。看到这里:http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example – 2013-03-10 23:13:59

+3

当这种情况发生在我身上时,通常是因为其中一列出乎意料地是类的因素'。你可以把你的问题粘贴到'dput(head(mc))'的结果中,这样我们就可以检查了吗? – Ben 2013-03-10 23:16:25

+0

希望我现在增加了足够的信息... – Ellen 2013-03-10 23:55:37

回答

3

它看起来像你想的CNTRL水平的0​​变量,以对应于零,但这不是因素水平的顺序。尝试

mc <- transform(mc,spike=factor(spike,levels=c("CNTRL","0.2","0.6","1.2"))) 
+0

谢谢!成功:) – Ellen 2013-03-11 01:07:49

+0

很高兴它的工作。您可以(但不一定要)点击答案旁边的复选标记以将其标记为接受的答案... – 2013-03-11 03:39:08