草庐IT

关于 r：为什么我的 lm 模型不显示线性关系，但在 geom_smooth 中显示？

codeneng 2023-03-28 原文

Why does my lm model not show a linear relationship but does in geom_smooth?

我试图建立一个线性模型来解释粒子浓度和荧光之间的关系。由于某种原因，我无法让模型使用 lm 来拟合数据，但它确实在 ggplot geom_smooth 函数中工作。

下面是对数荧光和对数粒子浓度的图...

我用下面的代码做了一个模型

1
2
3
Calicurve.M1 <- lm(Fluorescence~Particle.conc,
                na.action = na.exclude,
                data = Calicurve)

但是，当我使用此模型来预测值并添加到我的绘图(在 ggplot2 中)时，它看起来不正确

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
    ### predict values and put into dataframe
pdat <- expand.grid(Particle.conc = c(5, 50, 500, 5000, 50000, 500000, 5000000,
                                      50000000, 500000000, 5000000000, 50000000000,
                                      500000000000))
pred <- predict(Calicurve.M1, newdata = pdat, na.rm=T,
                type ="response", se.fit = TRUE)
predframe <- data.frame (pdat, preds=pred$fit, se=pred$se.fit)
predframe$upperse <- (predframe$preds+predframe$se)
predframe$lowerse <- (predframe$preds-predframe$se)

### plot calibration curve ###    
plot <- ggplot(data=Calicurve, aes(x=Particle.conc, y=Fluorescence)) +
  geom_point()+
  scale_y_log10(name ="Fluorecence (AFU)",
                     limits = c(1,1200))+
  scale_x_log10(name ="Particle concentration (particles/mL)")+
  theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
        strip.text = element_text(face ="italic"),
        legend.position = c(0.6, 0.75), legend.justification = c(1, 0)) +
  geom_line(data= predframe, aes(x=Particle.conc,y=preds), linetype=1) +
  geom_line(data= predframe, aes(x=Particle.conc,y=upperse), linetype=2) +
  geom_line(data= predframe, aes(x=Particle.conc,y=lowerse), linetype=2)

最后，当我使用 #geom_smooth(method='lm')+ 在 geom_smooth 中使用线性模型时，它会按照我的预期绘制关系...

我想了解第一个模型为什么不起作用以及如何修复它，以便我可以使用该模型来预测荧光值。

按要求提供样本数据

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
structure(list(Particle.conc = c(50, 50, 50, 500, 500, 500, 5000,
5000, 5000, 50000, 50000, 50000, 5e+05, 5e+05, 5e+05, 5e+06,
5e+06, 5e+06, 5e+07, 5e+07, 5e+07, 5e+08, 5e+08, 5e+08, 5e+09,
5e+09, 5e+09, 5e+10, 5e+10, 5e+10, 5e+11, 5e+11, 5e+11), Fluorecence = c(2.649,
2.671, 2.502, 3.926, 3.965, 4, 6.674, 6.337, 6.56, 12.204, 12.168,
12.209, 24.91, 25.54, 25.384, 38.232, 37.845, 37.979, 80.547,
80.343, 79.891, 168.693, 168.008, 168.826, 349.318, 351.304,
355.288, 556.081, 555.348, 554.112, 1105.749, 1103.063, 1097.552
), Average.FL = c(2.607333333, NA, NA, 3.963666667, NA, NA, 6.523666667,
NA, NA, 12.19366667, NA, NA, 25.278, NA, NA, 38.01866667, NA,
NA, 80.26033333, NA, NA, 168.509, NA, NA, 351.97, NA, NA, 555.1803333,
NA, NA, 1102.121333, NA, NA), Fl.Bl = c(1.463, NA, NA, 2.819333333,
NA, NA, 5.379333333, NA, NA, 11.04933333, NA, NA, 24.13366667,
NA, NA, 36.87433333, NA, NA, 79.116, NA, NA, 167.3646667, NA,
NA, 350.8256667, NA, NA, 554.036, NA, NA, 1100.977, NA, NA)), .Names = c("Particle.conc",
"Fluorecence","Average.FL","Fl.Bl"), class ="data.frame", row.names = c(NA,
-33L))

  • 请问您能提供一些示例数据吗？例如将 dput(Calicurve) 的输出粘贴到您的问题中
  • 您好 - 已按要求更新了样本数据

在将变量传递给 lm.

之前，您必须对其进行转换

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
df = tibble(x=10^(1:100), y=10^((1:100)*2) + rnorm(100))
lm(log(y, 10)~log(x, 10), df) %>% summary()

# Call:
# lm(formula = log(y, 10) ~ log(x, 10), data = df)
#
# Residuals:
#        Min         1Q     Median         3Q        Max
# -6.000e-05 -3.698e-05 -1.480e-05  8.320e-06  1.486e-03
#
# Coefficients:
#              Estimate Std. Error   t value Pr(>|t|)    
# (Intercept) 6.291e-05  3.086e-05 2.038e+00   0.0442 *  
# log(x, 10)  2.000e+00  5.305e-07 3.770e+06   <2e-16 ***
# ---
# Signif. codes:  0 a€?***a€? 0.001 a€?**a€? 0.01 a€?*a€? 0.05 a€?.a€? 0.1 a€? a€? 1
#
# Residual standard error: 0.0001531 on 98 degrees of freedom
# Multiple R-squared:      1,   Adjusted R-squared:      1
# F-statistic: 1.421e+13 on 1 and 98 DF,  p-value: < 2.2e-16

然后要预测，你必须像这样反转变换：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
df %>%
    mutate(y_hat = 10^predict(lm(log(y, 10)~log(x,10))))

# # A tibble: 100 x 3
#              x        y   y_hat
#          <dbl>    <dbl>   <dbl>
#  1          10  1.00e 2 1.00e 2
#  2         100  1.00e 4 1.00e 4
#  3        1000  1.00e 6 1.00e 6
#  4       10000 10.00e 7 1.00e 8
#  5      100000  1.00e10 1.00e10
#  6     1000000 10.00e11 1.00e12
#  7    10000000  1.00e14 1.00e14
#  8   100000000  1.00e16 1.00e16
#  9  1000000000  1.00e18 1.00e18
# 10 10000000000  1.00e20 1.00e20
# # ... with 90 more rows
geom_smoothsmoothnbspbrsectionggplot2glmlmrregression

有关关于 r：为什么我的 lm 模型不显示线性关系，但在 geom_smooth 中显示？的更多相关文章

  1. ruby - 为什么我可以在 Ruby 中使用 Object#send 访问私有(private)/ protected 方法？ - 2

    类classAprivatedeffooputs:fooendpublicdefbarputs:barendprivatedefzimputs:zimendprotecteddefdibputs:dibendendA的实例a=A.new测试a.foorescueputs:faila.barrescueputs:faila.zimrescueputs:faila.dibrescueputs:faila.gazrescueputs:fail测试输出failbarfailfailfail.发送测试[:foo,:bar,:zim,:dib,:gaz].each{|m|a.send(m)resc

  2. ruby-on-rails - Rails - 子类化模型的设计模式是什么？ - 2

    我有一个模型:classItem项目有一个属性“商店”基于存储的值，我希望Item对象对特定方法具有不同的行为。Rails中是否有针对此的通用设计模式？如果方法中没有大的if-else语句，这是如何干净利落地完成的？ 最佳答案 通常通过Single-TableInheritance. 关于ruby-on-rails-Rails-子类化模型的设计模式是什么？，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.co

  3. ruby-on-rails - Rails 编辑表单不显示嵌套项 - 2

    我得到了一个包含嵌套链接的表单。编辑时链接字段为空的问题。这是我的表格:Editingkategori{:action=>'update',:id=>@konkurrancer.id})do|f|%>'Trackingurl',:style=>'width:500;'%>'Editkonkurrence'%>|我的konkurrencer模型:has_one:link我的链接模型:classLink我的konkurrancer编辑操作:defedit@konkurrancer=Konkurrancer.find(params[:id])@konkurrancer.link_attrib

  4. ruby - 什么是填充的 Base64 编码字符串以及如何在 ruby​​ 中生成它们？ - 2

    我正在使用的第三方API的文档状态:"[O]urAPIonlyacceptspaddedBase64encodedstrings."什么是“填充的Base64编码字符串”以及如何在Ruby中生成它们。下面的代码是我第一次尝试创建转换为Base64的JSON格式数据。xa=Base64.encode64(a.to_json) 最佳答案 他们说的padding其实就是Base64本身的一部分。它是末尾的“=”和“==”。Base64将3个字节的数据包编码为4个编码字符。所以如果你的输入数据有长度n和n%3=1=>"=="末尾用于填充n%

  5. ruby - 解析 RDFa、微数据等的最佳方式是什么，使用统一的模式/词汇(例如 schema.org)存储和显示信息 - 2

    我主要使用Ruby来执行此操作，但到目前为止我的攻击计划如下:使用gemsrdf、rdf-rdfa和rdf-microdata或mida来解析给定任何URI的数据。我认为最好映射到像schema.org这样的统一模式，例如使用这个yaml文件，它试图描述数据词汇表和opengraph到schema.org之间的转换:#SchemaXtoschema.orgconversion#data-vocabularyDV:name:namestreet-address:streetAddressregion:addressRegionlocality:addressLocalityphoto:i

  6. ruby - 为什么 4.1%2 使用 Ruby 返回 0.0999999999999996？但是 4.2%2==0.2 - 2

    为什么4.1%2返回0.0999999999999996？但是4.2%2==0.2。 最佳答案 参见此处:WhatEveryProgrammerShouldKnowAboutFloating-PointArithmetic实数是无限的。计算机使用的位数有限(今天是32位、64位)。因此计算机进行的浮点运算不能代表所有的实数。0.1是这些数字之一。请注意，这不是与Ruby相关的问题，而是与所有编程语言相关的问题，因为它来自计算机表示实数的方式。 关于ruby-为什么4.1%2使用Ruby返

  7. ruby - ruby 中的 TOPLEVEL_BINDING 是什么？ - 2

    它不等于主线程的binding，这个toplevel作用域是什么？此作用域与主线程中的binding有何不同？>ruby-e'putsTOPLEVEL_BINDING===binding'false 最佳答案 事实是，TOPLEVEL_BINDING始终引用Binding的预定义全局实例，而Kernel#binding创建的新实例>Binding每次封装当前执行上下文。在顶层，它们都包含相同的绑定(bind)，但它们不是同一个对象，您无法使用==或===测试它们的绑定(bind)相等性。putsTOPLEVEL_BINDINGput

  8. ruby - Infinity 和 NaN 的类型是什么？ - 2

    我可以得到Infinity和NaNn=9.0/0#=>Infinityn.class#=>Floatm=0/0.0#=>NaNm.class#=>Float但是当我想直接访问Infinity或NaN时:Infinity#=>uninitializedconstantInfinity(NameError)NaN#=>uninitializedconstantNaN(NameError)什么是Infinity和NaN？它们是对象、关键字还是其他东西？ 最佳答案 您看到打印为Infinity和NaN的只是Float类的两个特殊实例的字符串

  9. ruby-on-rails - 如果 Object::try 被发送到一个 nil 对象，为什么它会起作用？ - 2

    如果您尝试在Ruby中的nil对象上调用方法，则会出现NoMethodError异常并显示消息:"undefinedmethod‘...’fornil:NilClass"然而，有一个tryRails中的方法，如果它被发送到一个nil对象，它只返回nil:require'rubygems'require'active_support/all'nil.try(:nonexisting_method)#noNoMethodErrorexceptionanymore那么try如何在内部工作以防止该异常？ 最佳答案 像Ruby中的所有其他对象

  10. ruby-on-rails - 使用 Sublime Text 3 突出显示 HTML 背景语法中的 ERB？ - 2

    所以我在关注Railscast，我注意到在html.erb文件中，ruby代码有一个微弱的背景高亮效果，以区别于其他代码HTML文档。我知道Ryan使用TextMate。我正在使用SublimeText3。我怎样才能达到同样的效果？谢谢! 最佳答案 为SublimeText安装ERB包。假设您安装了SublimeText包管理器*，只需点击cmd+shift+P即可获得命令菜单，然后键入installpackage并选择PackageControl:InstallPackage获取包管理器菜单。在该菜单中，键入ERB并在看到包时选择

随机推荐