Question

我是R的新手，无法找到这个（看似）简单问题的答案。我已经搜索了几天，并且阅读了几篇论文和帮助页面。

我已经能够绘制一条线（红色）。

red in my graph

我想绘制另一条适合后点的线条。我希望这条线看起来像这张图片中的黑线（Křivan和Priyadarshi，2015）。

然而，我无法绘制线条。

我尝试使用以下代码填充该行，但图表上没有显示任何内容：

我希望通过以下方式拟合的值：

Prey_isocline_x     <- c(8.2, 7.15, 7.65, 10.6, 7.947368421, 5.35,
                         6, 8.2, 7.473684211, 1.5, 1.3, 0.95, 1.85,
                         1.15, 0.6, 2.7, 1.3, 0.25, 0.25, 6.263157895,
                         4, 0.3, 5.1, 4.15, 1.15, 1.6, 1.6, 1.55)
Prey_isocline_y     <- c(0.45, 0.3, 0.2, 0.2, 0.105263158, 0.8, 0.5,
                         0.15, 0.052631579, 0.642857143, 1, 1, 1.15,
                         0.7, 0.55, 0.35, 0.8, 1.15, 1.55, 0.578947368,
                         0.5, 2.55, 0.15, 0.25, 0.45, 2.45, 2.45, 1.3)
Prey_isocline       <- data.frame(Prey_isocline_x, Prey_isocline_y)

Predator_isocline_x <- c(0.25, 0.15, 0.3, 0.7, 0.25, 0.25, 0.05, 0.5, 0.45,
                         0.5, 0.5, 0.15, 0.6, 1.4, 0.85, 0.15, 0.15, 0.6)
Predator_isocline_y <- c(2.35, 2.9, 3.6, 3.6, 2.35, 4.45, 1.45, 1.7, 1.65, 
                         1.7, 2.9, 1.8, 1.9, 2.35, 2.9, 2.8, 2.5, 3.05)
Predator_isocline   <- data.frame(Predator_isocline_x, Predator_isocline_y)

首次尝试绘图：

plot(Prey_isocline_x, Prey_isocline_y,
        axes = F,
        xlab= "",
        ylab= "",
        pch=1, col="black")
fit <- nls(Prey_isocline_y ~ SSlogis(Prey_isocline_x, Asym, xmid, scal), 
       data=Prey_isocline,
       trace = TRUE)
summary(fit)
curve(predict(fit, newdata = data.frame(Prey_isocline_y=x)), add=TRUE)

输出第一次尝试：

> par(new=T)
> plot(Prey_isocline_x, Prey_isocline_y,
+         axes = F,
+         xlab= "",
+         ylab= "",
+         pch=1, col="black")
> fit <- nls(Prey_isocline_y ~ SSlogis(Prey_isocline_x, Asym, xmid, scal), 
+            data=Prey_isocline,
+            trace = TRUE)
Error in nls(y ~ 1/(1 + exp((xmid - x)/scal)), data = xy, start = list(xmid = aux[1L],  : 
  step factor 0.000488281 reduced below 'minFactor' of 0.000976562
> summary(fit)
Error in summary(fit) : object 'fit' not found
> curve(predict(fit, newdata = data.frame(Prey_isocline_y=x)), add=TRUE)
Error in predict(fit, newdata = data.frame(Prey_isocline_y = x)) : 
  object 'fit' not found

第二次尝试：

model <- loess(formula=Prey_isocline_x~Prey_isocline_y, 
data=Predator_isocline)
abline(model, col="black")

第二次输出：

> model <- loess(formula=Prey_isocline_x~Prey_isocline_y, data=Predator_isocline)
> abline(model, col="black")

第三次尝试：

nls_fit <- nls(Prey_isocline_y ~ (b*Prey_isocline_x) - (b*Prey_isocline_x*Prey_isocline_x/K) -
              (Predator_isocline_y*(Prey_isocline_x^k/(x^k+C^k)*(l*x/(1+l*h*x)))),
               data = Prey_isocline,
               start = list(b = 2.2,
                            e = 1.5,
                            K = 30,
                            k = 20,
                            l = 0.1,
                            h = 0.25,
                            C = 1,
                            m = 1.0))
lines(Prey_isocline_x, predict(nls_fit), col = "green")

第三输出：

> nls_fit <- nls(Prey_isocline_y ~ (b*Prey_isocline_x) - (b*Prey_isocline_x*Prey_isocline_x/K) -
+               (Predator_isocline_y*(Prey_isocline_x^k/(x^k+C^k)*(l*x/(1+l*h*x)))),
+                data = Prey_isocline,
+                start = list(b = 2.2,
+                             e = 1.5,
+                             K = 30,
+                             k = 20,
+                             l = 0.1,
+                             h = 0.25,
+                             C = 1,
+                             m = 1.0))
Error in nlsModel(formula, mf, start, wts) : 
  singular gradient matrix at initial parameter estimates
In addition: There were 30 warnings (use warnings() to see them)
> lines(Prey_isocline_x, predict(nls_fit), col = "green")
Error in predict(nls_fit) : object 'nls_fit' not found

第四次尝试：

nls_fit <- nls(Prey_isocline_y ~ a + b * Prey_isocline_x^(-c), Prey_isocline,
               start = list(a = 80, b = 20, c = 0.2))
lines(Prey_isocline_x, predict(nls_fit), col = "green")

第四输出：

> nls_fit <- nls(Prey_isocline_y ~ a + b * Prey_isocline_x^(-c), Prey_isocline,
+                start = list(a = 80, b = 20, c = 0.2))
Error in nls(Prey_isocline_y ~ a + b * Prey_isocline_x^(-c), Prey_isocline,  : 
  step factor 0.000488281 reduced below 'minFactor' of 0.000976562
> lines(Prey_isocline_x, predict(nls_fit), col = "green")
Error in predict(nls_fit) : object 'nls_fit' not found

我完全迷失了，我希望有人可以帮助我。

Answer 1

以下是关于如何使用适合您的积分define('DRUPAL_ROOT', '/path_to_project'); require_once DRUPAL_ROOT . '/includes/bootstrap.inc';的情节的部分答案。

loess

基于Prey_isocline数据生成黄土模型：

# to prevent typing in messy codes, change "X_isocline_x" to "x" & "X_isocline_y" to "y"
names(Prey_isocline) <- c("x", "y")
names(Predator_isocline) <- c("x", "y")

为要绘制的黄土线创建新数据框：

model <- loess(y ~ x , Prey_isocline)

根据猎物等值线绘制黄土线：

new.prey <- data.frame(x=Prey_isocline$x)
new.prey$fit <- predict(model, new.prey)
new.prey <- new.prey[order(new.prey$x),]

重复捕食者的步骤

with(Prey_isocline, plot(x, y, ylim=c(0,5)))
with(new.prey, lines(x, fit))

为捕食者和黄土线添加点数：

model <- loess(y ~ x , Predator_isocline)
new.prd <- data.frame(x=Predator_isocline$x)
new.prd$fit <- predict(model, new.prd)
new.prd <- new.prd[order(new.prd$x),]

编辑：

如果两个数据帧合并，将更容易绘制。

with(Predator_isocline, points(x,y, col="red", pch=16))
with(new.prd, lines(x, fit))

如何在R中拟合非线性？

1 个答案: