- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试使用 R 中的 ggplot 在公式 y ~ a*x^b 中绘制幂回归线。
看起来我需要在 stat_smooth
下使用 method = "nls"但我收到一条警告消息,指出在评估模型时会产生缺失值或无穷大值。
如果我从我的数据中排除一些点,我可以让它工作,但我看不出这些点和那些有效的点之间有任何区别
这里发生了什么?
数据
library(dplyr)
library(ggplot2)
mydf <- structure(list(x = c(3.09090909090909, 1.09677419354839, 1.47826086956522,
2.125, 5.23076923076923, 2.42857142857143, 45.3333333333333,
2.3448275862069, 2.51851851851852, 45.3333333333333, 6.18181818181818,
1.78947368421053, 9.71428571428572, 13.6, 2.95652173913043, 4,
11.3333333333333, 1.51111111111111, 3.23809523809524, 6.8, 22.6666666666667,
5.66666666666667, 1.83783783783784, 2.61538461538462, 1.30769230769231,
1.7, 7.55555555555556, 1.03030303030303, 2, 1.88888888888889,
1.17241379310345, 1.13333333333333, 1.4468085106383, 8.5, 2.26666666666667,
1.01492537313433, 1.61904761904762, 3.57894736842105, 4.68965517241379,
2.19354838709677, 3.77777777777778, 17, 4.25, 1.07936507936508,
1.33333333333333, 1.74358974358974, 1.28301886792453, 2.06060606060606,
4.68965517241379, 1.94285714285714, 2.83333333333333, 1.04615384615385,
2.72, 1.58139534883721, 1.23636363636364, 1.41666666666667, 1.65853658536585,
1.15254237288136, 1.11475409836066, 1.38775510204082, 1.0625,
1.36, 3.4, 1.19298245614035, 1.21428571428571, 1.25925925925926,
1.54545454545455), y = c(193.571428571429, 289, 249, 215.285714285714,
179.857142857143, 204, 141.142857142857, 206.142857142857, 201.428571428571,
141.142857142857, 178.285714285714, 224.714285714286, 165.714285714286,
160, 198.428571428571, 184.857142857143, 164.285714285714, 245.857142857143,
192.714285714286, 172, 150.714285714286, 179.142857142857, 223.571428571429,
200.714285714286, 259.285714285714, 230, 171.714285714286, 315.714285714286,
216, 220.571428571429, 277, 285, 250, 170.285714285714, 209.142857142857,
335.714285714286, 238.285714285714, 188.571428571429, 181.142857142857,
210.857142857143, 185.714285714286, 151.142857142857, 181.428571428571,
295.857142857143, 257.142857142857, 225.857142857143, 261.714285714286,
215.714285714286, 181.142857142857, 217.857142857143, 199.714285714286,
313.571428571429, 200.428571428571, 242.857142857143, 269.714285714286,
250.142857142857, 234.857142857143, 282, 285.142857142857, 253.285714285714,
309.857142857143, 253.857142857143, 189.142857142857, 274.714285714286,
273.285714285714, 261.857142857143, 244.142857142857)), row.names = c(NA, -67), class = "data.frame")
mydf %>%
ggplot(aes(x = x, y = y))+
geom_point()+
stat_smooth(
method = "nls",
formula = "y~a*x^b",
method.args = list(start= c(a = 1,b=1)),
se=FALSE
)
# Warning message:
# Computation failed in `stat_smooth()`:
# Missing value or an infinity produced when evaluating the model
mydf_trunc <- mydf[c(1:51, 53:58, 60, 63, 67),]
mydf_trunc%>%
ggplot(aes(x = x, y = y))+
geom_point()+
stat_smooth(
method = "nls",
formula = "y~a*x^b",
method.args = list(start= c(a = 1,b=1)),
se=FALSE
)
bad_pts <- mydf[c(52,59, 61, 62, 64, 65, 66),]
bad_pts
# # A tibble: 7 x 2
# x y
# <dbl> <dbl>
# 1 1.05 314.
# 2 1.11 285.
# 3 1.06 310.
# 4 1.36 254.
# 5 1.19 275.
# 6 1.21 273.
# 7 1.26 262.
最佳答案
我怀疑 a) 有问题的点至少对于它们的 x 值是其他一些点的重复,因此在估计过程的某些部分为您提供无限斜率,并且 b) 应该允许渐近线高于零:
ggplot(data=mydf, aes(x = jitter(x), y = y))+
geom_point()+
geom_smooth(
method="nls",
formula= y~ a*x^(-b) +cc,
method.args = list(start= c(a = 1,b=2 ,cc=100)),
se=FALSE
)
关于r - 使用选择点绘制幂回归时 stat_smooth 中的错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/59329935/
是否有某种方法可以使用 JPA 或 Hibernate Crtiteria API 来表示这种 SQL?或者我应该将其作为 native 执行吗? SELECT A.X FROM (SELECT X,
在查询中, select id,name,feature,marks from (....) 我想删除其 id 在另一个 select 语句中存在的那些。 从 (...) 中选择 id 我是 sql
我想响应用户在 select 元素中选择一个项目。然而这个 jQuery: $('#platypusDropDown').select(function () { alert('You sel
这个问题在这里已经有了答案: SQL select only rows with max value on a column [duplicate] (27 个回答) 关闭8年前。 我正在学习 SQL
This question already has answers here: “Notice: Undefined variable”, “Notice: Undefined index”, and
我在 php 脚本中调用 SQL。有时“DE”中没有值,如果是这种情况我想从“EN”中获取值 应该是这样的,但不是这样的 IF (EXISTS (SELECT epf_application_deta
这可能是一个奇怪的问题,但不知道如何研究它。执行以下查询时: SELECT Foo.col1, Foo.col2, Foo.col3 FROM Foo INNER JOIN Bar ON
如何在使用 Camera.DestinationType.FILE_URI. 时在 phonegap camera API 中同时选择或拾取多个图像我能够一次只选择一张图像。我可以使用 this 在
这是一个纯粹的学术问题。这两个陈述实际上是否相同? IF EXISTS (SELECT TOP 1 1 FROM Table1) SELECT 1 ELSE SELECT 0 相对 IF EXIS
我使用 JSoup 来解析 HTML 响应。我有多个 Div 标签。我必须根据 ID 选择 Div 标签。 我的伪代码是这样的 Document divTag = Jsoup.connect(link
我正在处理一个具有多个选择框的表单。当用户从 selectbox1 中选择一个选项时,我需要 selectbox2 active 的另一个值。同样,当他选择 selectbox2 的另一个值时,我需要
Acme Inc. Christa Woods Charlotte Freeman Jeffrey Walton Ella Hubbard Se
我有一个login.html其中form定义如下: First Initial Plus Last Name : 我的do_authorize如下: "; pri
$.get( 'http://www.ufilme.ro/api/load/maron_online/470', function(data
我有一个下拉列表“磅”、“克”、“千克”和“盎司”。我想要这样一种情况,当我选择 gram 来执行一个函数时,当我在输入字段中输入一个值时,当我选择 pounds 时,我想要另一个函数来执行时我在输入
我有一个 GLSL 着色器,它从输入纹理的 channel 之一(例如 R)读取,然后写入输出纹理中的同一 channel 。该 channel 必须由用户选择。 我现在能想到的就是使用一个 int
我想根据下拉列表中的选定值生成输入文本框。 Options 2 3 4 5 就在这个选择框之后,一些输入字段应该按照选定的数字出现。 最佳答案 我建议您使用响应式(Reac
我是 SQL 新手,我想问一下如何根据首选项和分组选择条目。 +----------+----------+------+ | ENTRY_ID | ROUTE_ID | TYPE | +------
我有以下表结构: CREATE TABLE [dbo].[UTS_USERCLIENT_MAPPING_USER_LIST] ( [MAPPING_ID] [int] IDENTITY(1,1
我在移除不必要的床单时遇到了问题。我查看了不同的论坛并将不同的解决方案混合在一起。 此宏删除工作表(第一张工作表除外)。 Sub wrong() Dim sht As Object Applicati
我是一名优秀的程序员,十分优秀!