- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有这样的数据,可以在不同的日子给受试者打分(采用 POSIXct 格式)。
head(test)
Date Subject score
1 2012-08-10 Black6 0
2 2012-08-11 Black6 0
3 2012-08-12 Black6 0
4 2012-08-13 Black6 0
5 2012-08-14 Black6 0
6 2012-08-15 Black6 0
拟合黄土曲线很简单。
ggplot(test,aes(Date,score))+geom_smooth()+geom_point()
我想做的是拟合三阶多项式线。如果我输入以下内容,则会出现错误:
ggplot(test,aes(Date,score))+stat_smooth(method = "lm", formula = score ~ poly(Date, 3), size = 1)+geom_point()
Error in eval(expr, envir, enclos) : object 'score' not found
如果我在 stat_smooth() 中将 Date 指定为数字,我会得到同样的错误。有什么办法可以在 ggplot 中做到这一点吗?
这是数据:
test<-structure(list(Date = structure(c(1344556800, 1344643200, 1344729600,
1344816000, 1344902400, 1344988800, 1345075200, 1345161600, 1345248000,
1345334400, 1345420800, 1345507200, 1345593600, 1345680000, 1345766400,
1345852800, 1345939200, 1346025600, 1346112000, 1346198400, 1346284800,
1346371200, 1346457600, 1346544000, 1346630400, 1346716800, 1346803200,
1346889600, 1346976000, 1347062400, 1347148800, 1347235200, 1347321600,
1347408000, 1347494400, 1347580800, 1347667200, 1347753600, 1347840000,
1347926400, 1348012800, 1348099200, 1348185600, 1348272000, 1348358400,
1348444800, 1348531200, 1348617600, 1348704000, 1348790400, 1348876800,
1348963200, 1349049600, 1349136000, 1349222400, 1349308800, 1349395200
), class = c("POSIXct", "POSIXt"), tzone = "UTC"), Subject = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = "Black6", class = "factor"),
score = c(0, 0, 0, 0, 0, 0, 0, 0.25, 0.25, 0.25, 0.25, 0.25,
0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 1, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 1.5, 1.5, 1.25, 1.25, 1.25, 1,
1, 1, 1, 1, 1, 1, 1, 1.5, 1.5, 1.25, 1, 1, 1, 0.5, 0.5, 0.5,
0.25, 0.25)), .Names = c("Date", "Subject", "score"), row.names = c(NA,
57L), class = "data.frame")
最佳答案
stat_smooth
中的公式参数必须根据美学 y
和 x
指定,而不是映射到那些的原始变量美学。
ggplot(test,aes(Date,score)) +
stat_smooth(method = "lm", formula = y ~ poly(x, 3), size = 1) +
geom_point()
关于r - 以日期作为解释变量的ggplot多项式拟合,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/14779777/
gnuplot 中拟合函数的正确方法是什么 f(x)有下一个表格吗? f(x) = A*exp(x - B*f(x)) 我尝试使用以下方法将其拟合为任何其他函数: fit f(x) "data.txt
(1)首先要建立数据集 ? 1
测量显示一个信号,其形式类似于具有偏移量和因子的平方根函数。如何找到系数并在一个图中绘制原始数据和拟合曲线? require(ggplot2) require(nlmrt) # may be thi
我想将以下函数拟合到我的数据中: f(x) = Offset+Amplitudesin(FrequencyT+Phase), 或根据 Wikipedia : f(x) = C+alphasin(ome
我正在尝试使用与此工具相同的方法在 C# 中拟合 Akima 样条曲线:https://www.mycurvefit.com/share/4ab90a5f-af5e-435e-9ce4-652c95c
问题:开放层适合 map ,只有在添加特征之后(视觉),我该如何避免这种情况? 我在做这个 第 1 步 - 创建特征 var feature = new ol.Feature({...}); 第 2
我有一个数据变量,其中包含以下内容: [Object { score="2.8", word="Blue"}, Object { score="2.8", word="Red"}, Objec
我正在尝试用中等大小的 numpy float 组来填充森林 In [3]: data.shape Out[3]: (401125, 5) [...] forest = forest.fit(data
我想用洛伦兹函数拟合一些数据,但我发现当我使用不同数量级的参数时拟合会出现问题。 这是我的洛伦兹函数: function [ value ] = lorentz( x,x0,gamma,amp )
我有一些数据,我希望对其进行建模,以便能够在与数据相同的范围内获得相对准确的值。 为此,我使用 polyfit 来拟合 6 阶多项式,由于我的 x 轴值,它建议我将其居中并缩放以获得更准确的拟合。 但
我一直在寻找一种方法来使数据符合 beta 二项分布并估计 alpha 和 beta,类似于 VGAM 库中的 vglm 包的方式。我一直无法找到如何在 python 中执行此操作。有一个 scipy
我将 scipy.optimize.minimize ( https://docs.scipy.org/doc/scipy/reference/tutorial/optimize.html ) 函数与
在过去的几天里,我一直在尝试使用 python 绘制圆形数据,方法是构建一个范围从 0 到 2pi 的圆形直方图并拟合 Von Mises 分布。我真正想要实现的是: 具有拟合 Von-Mises 分
我有一个简单的循环,它在每次迭代中都会创建一个 LSTM(具有相同的参数)并将其拟合到相同的数据。问题是迭代过程中需要越来越多的时间。 batch_size = 10 optimizer = opti
我有一个 Python 系列,我想为其直方图拟合密度。问题:是否有一种巧妙的方法可以使用 np.histogram() 中的值来实现此结果? (请参阅下面的更新) 我目前的问题是,我执行的 kde 拟
我有一个简单的 keras 模型(正常套索线性模型),其中输入被移动到单个“神经元”Dense(1, kernel_regularizer=l1(fdr))(input_layer) 但是权重从这个模
我正在尝试解决 Boston Dataset 上的回归问题在random forest regressor的帮助下.我用的是GridSearchCV用于选择最佳超参数。 问题一 我是否应该将 Grid
使用以下函数,可以在输入点 P 上拟合三次样条: def plotCurve(P): pts = np.vstack([P, P[0]]) x, y = pts.T i = np.aran
我有 python 代码可以生成数字 x、y 和 z 的三元组列表。我想使用 scipy curve_fit 来拟合 z= f(x,y)。这是一些无效的代码 A = [(19,20,24), (10,
我正在尝试从 this answer 中复制代码,但是我在这样做时遇到了问题。我正在使用包 VGAM 中的gumbel 发行版和 fitdistrplus . 做的时候出现问题: fit = fi
我是一名优秀的程序员,十分优秀!