小编Max*_*don的帖子

获得使用 rms lrm 对象的信心

我无法从 rms 包中获取逻辑回归以与 confint() 一起使用,这是一个示例:

library(rms)
data(mtcars)
dd <- datadist(mtcars)
options(datadist = "dd")
fit <- lrm(am ~ gear + mpg, data=mtcars)
confint(fit)
Run Code Online (Sandbox Code Playgroud)

这给出了错误:

错误:$ 运算符对于原子向量无效

Traceback() 给出:

4: profile.glm(object, which = parm, alpha = (1 - level)/4, trace = trace)
3: profile(object, which = parm, alpha = (1 - level)/4, trace = trace)
2: confint.glm(fit)
1: confint(fit)
Run Code Online (Sandbox Code Playgroud)

我猜想 lrm() 模型没有实现限制。

我的问题

有没有方便的替代方法?是否有为 rms 包创建其他替代方案?

r glm

5
推荐指数
2
解决办法
3377
查看次数

确定脚本是否在RStudio笔记本中运行

在我的htmlTable包中,我一直在使用base::interactive,以确定输出应该作为字符串输出还是输出到浏览器窗口(这里是代码).不幸的是,在RStudio笔记本中运行时interactive()返回TRUE,在笔记本中确定执行环境的等效方法是什么?

我试过了ls(),search()但无论环境如何,它们看起来都一样.我也有默认情况下为knitr文档knit_print.htmlTable调用的S3函数,但它似乎在笔记本环境中没有正确检测到它.

r rstudio knitr r-markdown

5
推荐指数
1
解决办法
138
查看次数

rcs 在 lm() 模型中生成错误的预测

我正在尝试复制这篇关于过度拟合的博客文章。我想探索样条与测试多项式的比较。

我的问题:使用 rcs() - 限制三次样条 - 来自 rms 包我在常规 lm() 中应用时得到非常奇怪的预测。ols() 工作正常,但我对这种奇怪的行为感到有些惊讶。有人可以向我解释发生了什么吗?

library(rms)
p4 <- poly(1:100, degree=4)
true4 <- p4 %*% c(1,2,-6,9)
days <- 1:70

noise4 <- true4 + rnorm(100, sd=.5)
reg.n4.4 <- lm(noise4[1:70] ~ poly(days, 4))
reg.n4.4ns <- lm(noise4[1:70] ~ ns(days,5))
reg.n4.4rcs <- lm(noise4[1:70] ~ rcs(days,5))
dd <- datadist(noise4[1:70], days)
options("datadist" = "dd")
reg.n4.4rcs_ols <- ols(noise4[1:70] ~ rcs(days,5))

plot(1:100, noise4)
nd <- data.frame(days=1:100)
lines(1:100, predict(reg.n4.4, newdata=nd), col="orange", lwd=3)
lines(1:100, predict(reg.n4.4ns, newdata=nd), col="red", lwd=3)
lines(1:100, …
Run Code Online (Sandbox Code Playgroud)

r linear-regression

4
推荐指数
1
解决办法
2408
查看次数

加载非 js 文件时清除 Node.js 的缓存,例如 graphql 语言定义

我想从我的graphql定义中清除节点的缓存,但在我更改语句所在的文件之前什么也不会发生require('my.graphql')。我认为这个问题与这个问题相关,我已经尝试这样做:

Object.keys(require.cache).forEach(function(key) {
  delete require.cache[key];
});
Run Code Online (Sandbox Code Playgroud)

但它不会重新加载代码。一定还有其他一些我错过的缓存正在发生。

由于该示例有点长,因此我分叉了此存储库并创建了一个graphql 分支,您可以安装并测试调整该src/schema.graphql文件:

git clone git@github.com:gforge/graphql-relay-mongodb-pagination.git
cd graphql-relay-mongodb-pagination
git checkout git checkout mongoose-w-gql-lang 
npm install
Run Code Online (Sandbox Code Playgroud)

require('my.graphql')是通过babel-plugin-inline-import介导的。

node.js babeljs graphql

4
推荐指数
1
解决办法
3800
查看次数

在异步函数中使用参数default中的'await'

我想使用异步函数作为默认参数,但我得到一个错误:'await' is not a valid idientifier name in an async function.这是语言限制还是我错过了什么?

背景

支持以非常简洁的方式使用函数来默认参数:

> a = ({ no }, b = (function(c) { return c -1; })(no)) => console.log(no,b);
[Function: a]
> a({ no: 2 })
2 1
Run Code Online (Sandbox Code Playgroud)

这可以允许从昂贵的数据库调用传递结果,同时保留获取参数的选项(如果为null).这要求我们可以在参数调用中进行等待:

> a = async ({ no }, b = await (async function(c) { return c -1; })(no)) => console.log(no,b);
a = async ({ no }, b = await (async function(c) { return c -1; })(no)) => console.log(no,b);
                       ^^^^^

SyntaxError: …
Run Code Online (Sandbox Code Playgroud)

javascript async-await

3
推荐指数
1
解决办法
395
查看次数

当一个变量被分解时,ggplot2中的Qplot会导致多个回归线

当我做一个简单的qplot()我想要一个回归线的情节.我如何告诉stat_smooth忽略这些因素?

这是我的示例代码:

library("ggplot2")
qplot(y=wt, x=mpg, size=cyl, col=factor(gear), data=mtcars) +
    stat_smooth(method=lm, formula=y~x)
Run Code Online (Sandbox Code Playgroud)

这给了这个图像:

一条考虑顺畅的线条

当我删除因子时,我得到了我想要的图形(虽然我无法删除真实数据集中的因子):

qplot(y=wt, x=mpg, size=cyl, col=gear, data=mtcars) +
    stat_smooth(method=lm, formula=y~x)
Run Code Online (Sandbox Code Playgroud)

一个合适的形象

r ggplot2

2
推荐指数
1
解决办法
1933
查看次数