0

基于我的data.frame的2个问题

structure(list(Collimator = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L), .Label = c("n", "y"), class = "factor"), angle = c(0L, 
15L, 30L, 45L, 60L, 75L, 90L, 105L, 120L, 135L, 150L, 165L, 180L, 
0L, 15L, 30L, 45L, 60L, 75L, 90L, 105L, 120L, 135L, 150L, 165L, 
180L), X1 = c(2099L, 11070L, 17273L, 21374L, 23555L, 23952L, 
23811L, 21908L, 19747L, 17561L, 12668L, 6008L, 362L, 53L, 21L, 
36L, 1418L, 6506L, 10922L, 12239L, 8727L, 4424L, 314L, 38L, 21L, 
50L), X2 = c(2126L, 10934L, 17361L, 21301L, 23101L, 23968L, 23923L, 
21940L, 19777L, 17458L, 12881L, 6051L, 323L, 40L, 34L, 46L, 1352L, 
6569L, 10880L, 12534L, 8956L, 4418L, 344L, 58L, 24L, 68L), X3 = c(2074L, 
11109L, 17377L, 21399L, 23159L, 23861L, 23739L, 21910L, 20088L, 
17445L, 12733L, 6046L, 317L, 45L, 26L, 46L, 1432L, 6495L, 10862L, 
12300L, 8720L, 4343L, 343L, 38L, 34L, 60L), average = c(2099.6666666667, 
11037.6666666667, 17337, 21358, 23271.6666666667, 23927, 23824.3333333333, 
21919.3333333333, 19870.6666666667, 17488, 12760.6666666667, 
6035, 334, 46, 27, 42.6666666667, 1400.6666666667, 6523.3333333333, 
10888, 12357.6666666667, 8801, 4395, 333.6666666667, 44.6666666667, 
26.3333333333, 59.3333333333)), .Names = c("Collimator", "angle", 
"X1", "X2", "X3", "average"), row.names = c(NA, -26L), class = "data.frame")

我希望绘制探测器计数与角度的关系图，无论是否安装准直器到设备上。我猜 geom_point 可能是总结数据的最佳方式

p <- ggplot(df, aes(x=angle,y=average,col=Collimator)) + geom_point() + geom_line()

与其在 y 轴上绘制平均计数，我更愿意重新调整数据的比例，使最大计数的角度对于准直器 Y 和 N 的值都为 1。我这样做的方式似乎很麻烦

range01 <- function(x){(x-min(x))/(max(x)-min(x))}
coly = subset(df,Collimator=='y')
coly$norm_count = range01(coly$average)
coln = subset(df,Collimator=='n')
coln$norm_count = range01(coln$average)
df = rbind(coln,coly)
p <- ggplot(df, aes(x=angle,y=norm_count,col=Collimator) + geom_point() + geom_line()

我确信这可以以更有效的方式完成，基于变量“准直器”将函数应用于 data.frame。我怎样才能做到这一点？

我还想为数据拟合一个函数，而不是使用 geom_line。我认为高斯函数在这种情况下可能有效，但不知道如何/是否可以在stat_smooth. 我也可以从这种拟合中提取蜂蜜酒/标准偏差吗？

4

2 回答 2

2

ggplot2与包plyr 齐头并进：

df <- ddply(df,.(Collimator),
            transform,
            norm_count1 = (average - min(average)) / (max(average) - min(average)) )
于 2013-07-12T15:56:23.977 回答
1

joran 的答案将最大值缩放为 1，将最小值缩放为 0；如果您只想缩放以使最大值为 1（并将 0 保留为 0），那就更简单了。

library("plyr")
df <- ddply(df, .(Collimator), transform,
            norm.average = average / max(average))

情节是

ggplot(df, aes(x=angle,y=norm.average,col=Collimator)) + 
  geom_point() + geom_line()

在此处输入图像描述

于 2013-07-12T17:48:01.373 回答