小编Pur*_*ina的帖子

R中的S形曲线拟合

我是 R 的初学者,我正在尝试将曲线拟合到(例如)可能如下所示的数据集上:

(x- value)  (y-value)
105 423
115 471
125 567
135 808
145 921.5
155 1040
Run Code Online (Sandbox Code Playgroud)

x 值代表刺激量,y 值代表运动反应(以 uV 为单位)。这些是 10 个科目的平均值,其中每个科目的 x 值都相同。

有人告诉我,这个数据集通常遵循 sigmoidal 拟合。我尝试用以下方法拟合它:

fit <- lm( y ~ poly(x, 3) )
Run Code Online (Sandbox Code Playgroud)

但我不确定这是否是执行此操作的合适方法:(

到目前为止,我的代码如下所示:

p <- ggplot (data, aes(x, y)) +
  geom_point(shape= 21, fill= "blue", colour= "black", size=2) +
  xlab("X value") + ylab("Y value") +
  geom_smooth(method= "lm", se= FALSE,  colour= "red", formula=y ~ poly(x, 3, raw=TRUE)) +
  geom_errorbar(aes(ymin=y-SE, ymax=y+SE), width=.9)+ 
  ggtitle ("Title")
p
Run Code Online (Sandbox Code Playgroud)

另外:一旦我拟合了曲线,我还想获得斜率(计算为曲线最陡点处的切线值) …

r curve-fitting ggplot2

5
推荐指数
1
解决办法
7324
查看次数

制作没有胡须的箱线图

我使用R中的box-plot函数为三个不同组的起始值制作了箱形图,如下所示:

boxplot(onset ~ group, data = pulse.dat, range = 0, col = "lightblue")
Run Code Online (Sandbox Code Playgroud)

但是,我想看看没有范围的数据看起来如何,所以我想创建一个没有胡须的盒子图.我也不介意任何类型的图表,只要它显示3组中的每一组的中位数,第25和第75四分位数.

有谁知道如何在R中做到这一点?

r boxplot

4
推荐指数
1
解决办法
6968
查看次数

标签 统计

r ×2

boxplot ×1

curve-fitting ×1

ggplot2 ×1