小编Dan*_* M.的帖子

我ggplot用来绘制时间过程数据(随着时间的推移固定比例到屏幕上的不同对象)并想要使用功能区来显示SE,但是功能区本身在顶部和底部边缘都有线条,这使得读取图形a有点难.我无法弄清楚如何摆脱这些边缘线.这是我的情节代码:

ggplot(d, aes(Time, y, color = Object, fill = Object)) +
  stat_summary(fun.y = "mean", geom = "line", size = 2) +
  stat_summary(fun.data = "mean_se", geom = "ribbon", alpha = .3)

有什么建议？

这是一个最小的工作示例.我已将我的数据压缩为:

   Time Object          y      lower     upper
 1 1000      C 0.12453389 0.04510504 0.2039627
 2 1000      T 0.58826856 0.37615078 0.8003864
 3 1000      U 0.09437160 0.03278069 0.1559625
 4 1100      C 0.12140127 0.03943988 0.2033627
 5 1100      T 0.64560823 0.44898727 0.8422292
 6 1100      U 0.06725172 0.01584248 0.1186610

d <- structure(list(Time = c(1000L, …

我试图使用nlmer与SSfpl用四参数逻辑函数,以适应一些数据.使用以下方法,我可以很好地适应整体数据:

nm.fpl <- nlmer(meanFix ~ SSfpl(Time, A, B, xmid, scal) ~ (scal | Subject), 
          data = dataSubset, start = c(A = 0.2, B = 0.7, xmid = 600, scal = 100))

现在我想添加一个固定的效果Condition,它有2个主题内的水平.我想评估两个条件是否在4个参数(A,B,xmid,scal)中的任何一个方面有所不同,但我不知道如何在这个公式中指定它.我可以将模型分别拟合到两个子集(条件A和条件B),然后比较参数,但这似乎不是正确的方法.