我使用最小二乘拟合的splines::ns的结果是不同的,没有我能看到的押韵或理由,我想我已经把问题追溯到了ns函数本身。
我已将问题减少到以下几个方面:
require(splines)
N <- 0
set.seed(1)
for (i in 1:100) N <- N + identical(ns(1:10,3),ns(1:10,3))
N
我的结果平均约39,范围34-44左右,但我预计每次100。为什么ns的结果应该是随机的?如果我在两个地方都用bs代替ns,我就能得到100,正如预期的那样。我的set.seed(1)希望证明我得到的随机性不是R想要的。
在一个干净的会话
我按照命令运行:
sudo apt-get install r-base
我得到了:
Reading package lists... Done
Building dependency tree
Reading state information... Done
Some packages could not be installed. This may mean that you have
requested an impossible situation or if you are using the unstable
distribution that some requi