我发布了一篇关于将x-y图控制为两条Normal曲线的帖子,后来发现我的事情太复杂了。我已经设法将其绘制为椭圆形,但这略微超过了估计误差;理想情况下可以绘制为菱形。
我要约会的代码是:
plot(c(-5,10), c(-5,5), xlab = expression(Age), ylab = expression(value), type="n")
draw.ellipse(Age, value, a=Age_error, b=value_error, col="grey70")
哪个情节:
有没有用菱形替换椭圆,其高度由2x value_error和宽度2x age_error控制?
我的数据框位于
之下 structure(list(Age = c(1L, 2L, 4L), value = c(3, -2, 0.01), Age_error = c(2,
1.4, 3), value_error = c(0.5, 1, 2.1)), .Names = c("Age", "value",
"Age_error", "value_error"), class = "data.frame", row.names = c(NA,
-3L))
非常感谢
答案 0 :(得分:1)
您可以使用TeachingDemos包中的my.symbols
和ms.polygon
函数绘制菱形:
library(TeachingDemos)
plot(c(-5,10), c(-5,5), xlab = expression(Age), ylab = expression(value),
type="n")
my.symbols( Age, value, ms.polygon, n=4, xsize=2*Age_error,
ysize=2*value_error, linesfun=polygon, col='grey' )
如果您不想填充菱形,请忽略linesfun
和col
。