0

我不知道这是否应该是一个新项目。

我将代码修改为：

library(Rcpp)

rm(list = ls())

datafile <- data.frame(
  N = c(1.5, 2.6, 0.555555555555556, 0.535714285714286, 0.418604651162791, 0.557377049180328, 0.463157894736842, 0.762589928057554, 0.583673469387755, 0.528350515463918, 0.649241146711636, 0.534764826175869, 0.556295802798135, 0.250856164383562, 0.202258726899384, 0.351266723598064, 0.226669475458184, 0.1275974583548, 0.0906183368869936, 0.123027510124284, 0.119124595871674)
)
View(datafile)
datafile$a<-dlogis(datafile$N, location = 0, scale = 1, log = FALSE)

数据文件不会变成 2 列 N 行。它变成 1 列 N+1 行，其中 datafile$a 变成一个条目。

我想知道是否有人会知道这是我的 dlogis 代码还是 optim 代码的问题。
1）控制台说我正在运行 dlogis 正确。
控制台同时说当我用 dlogis 调用 optim 时我没有正确运行它。使用以下代码调用 optim 的正确方法是什么？

2）另外，当我调用 dlogis 时，我想找到参数位置和比例，以便可以最大限度地减少我输入的数据的错误。还有其他我应该注意的事情吗？

谢谢

library(Rcpp)

rm(list = ls())

datafile <- data.frame(
  N = c(1.5, 2.6, 0.555555555555556, 0.535714285714286, 0.418604651162791, 0.557377049180328, 0.463157894736842, 0.762589928057554, 0.583673469387755, 0.528350515463918, 0.649241146711636, 0.534764826175869, 0.556295802798135, 0.250856164383562, 0.202258726899384, 0.351266723598064, 0.226669475458184, 0.1275974583548, 0.0906183368869936, 0.123027510124284, 0.119124595871674)
)
View(datafile)
datafile$a<-dlogis(datafile$N, location = 0, scale = 1, log = FALSE)
View(datafile)

optim(c(datafile$N, 0.5, 0.1), dlogis)

控制台给了我这个错误消息：

> optim(c(datafile$N, 0.5, 0.1), dlogis)
Error in optim(c(datafile$N, 0.5, 0.1), dlogis) : 
  objective function in optim evaluates to length 23 not 1
4

0 回答 0