rms

0热度

1回答

我想连续和分位数模型在同一小区联合收割机进行比较和对比两种方法（xtile是返回位数为因素的函数）： q.s <- cph(inc ~ rcs(exposure,3), data=data) q.q <- cph(inc ~ xtile(exposure,3), data=data) p.s <- Predict(q.s, exposure, fun=exp) p.q <- Predi

0热度

2回答

如何获得乐天认证密钥？

我正在为Rakuten.com开发一个Web应用程序。我发现他们的文档在https://developers.rakuten.com/但使用Rakuten RMS API我需要一个身份验证密钥。是否有另一种方式使用此API没有身份验证密钥（如演示api）？谢谢。

0热度

1回答

获得每个频率的RMS

我有一个输入信号，我计算了它的FFT。之后，我只需要在频率带宽上计算其均方根值，而不是针对所有频谱。我使用Parseval定理求解了整个频谱的RMS计算，但是如何计算这种RMS“选择性”？我已经正确地计算了索引以获得三个感兴趣的频率（F0，FC，F1），但是当将RMS应用于该频带时，似乎Parseval的定理不是完整的。我收到一个独特的10 KHz频率，从FFT总频谱的RMS是正确的，但其RM

1热度

1回答

如何找到一个整数列表（multiset，size n），其集合的均方根是一个整数？

我已经找到this one 蛮力当然是可能的，但有没有其他方法？有没有办法找到所有的多人组？有没有办法找出在一定限度内存在多少组合？也许这个问题对于SO来说太猥琐了，如果是这样的话我会移动它。我在javascript中创建了my own version，生成一列数字的所有可能组合，然后检查整数有效值。这些都是集合，而不是多集合。

1热度

1回答

来自predict.lrm的预测概率的置信区间

我正试图确定R中来自二项logistic回归的预测概率的置信区间。该模型使用lrm（来自包rms）进行估计，以允许在调查时将聚类标准误差受访者（每个受访出现最多3次的数据）： library(rms) model1<-lrm(outcome~var1+var2+var3,data=mydata,x=T,y=T,se.fit=T) model.rob<-robcov(model1,cluster

0热度

1回答

在J2ME中，RMS（记录管理系统）是否永久存储记录？

我正在使用Netbeans IDE 6.9 vesion开发J2ME移动应用程序。我的应用程序正常工作意味着将记录插入到RMS中，我可以列出系统默认电话中的记录。我的问题是，每当我运行应用程序刷新记录和新鲜再次存储。这意味着记录暂时存储。我想永久保存记录，并且在堆栈溢出中阅读了不同的答案，但不清楚。请为我提供解决此问题的方法。

1热度

1回答

“rms”中的survplot（）：整个组显示的风险数量，而不是每个曲线

我使用3210绘制生存曲线。但是，当我使用n.risk = TRUE来绘制风险表的数量时，R给了我整个数据集的数量，而不是每条曲线，我不知道为什么。 # initialize survival commands in R survive <- Surv(dat$dx_lastcontact_death_months, dat$event) library(rms) ff <- cp

3热度

1回答

使用受限三次样条的新数据的预测值

我有一些数据是使用受限三次样条建模的。我在rms包中使用rcs转换函数来生成线性模型的转换变量。这是一个使用5节的例子。 library('rms') my_df <- data.frame( y = -4 * -100:100 + -1.5 * (-100:100)**2 + 3 * (-100:100)**3 + rnorm(201, 0, 1e5), x = -100:

1热度

1回答

RMS包GLM：如何添加指标公式

在我试图建立R. 利用RMS包装GLM（）一GLM模型，但我不知道在GLM公式指标的正确表达。允许我使用虹膜数据为例，在基础R glm功能，我的代码工作，它是象下面这样： model1 = glm(Sepal.Length~Sepal.Width + Petal.Length + ifelse(Sepal.Width==3,1,0),data=iris) 但是，如果我用同样的公式中Glm，

2热度

1回答

如何将几条校准曲线添加到同一图上

set.seed(888) y<-sample(c(0,1), 100, replace = TRUE) x1<-rnorm(100,2,3) x2<-rnorm(100,1,4) library(rms) f1<-lrm(y~x1,,x=TRUE,y=TRUE) f2<-lrm(y~x2,,x=TRUE,y=TRUE) plot(calibrate(f1)) plot(calib