R中的框Cox变换，应用于列
R中的框Cox变换是一种常用的数据转换方法，用于处理偏态分布的数据。它通过对数据进行幂次变换，使得数据更加接近正态分布，从而满足许多统计模型的假设条件。
框Cox变换的一般形式为：
y' = (y^λ - 1) / λ
其中，y'是变换后的数据，y是原始数据，λ是变换参数。当λ为0时，变换公式变为对数变换，当λ为1时，变换公式变为平方根变换。
框Cox变换的优势在于可以改善数据的分布形态，使其更加接近正态分布。这有助于提高统计模型的准确性和可解释性。此外，框Cox变换还可以减小数据的离群值对统计分析的影响，提高模型的稳定性。
框Cox变换在许多领域都有广泛的应用，特别是在回归分析、生存分析和时间序列分析中常被使用。它可以用于处理各种类型的数据，包括连续型数据、计数型数据和百分比数据等。
对于R语言用户，可以使用R中的boxcox()函数来进行框Cox变换。该函数可以自动选择最佳的变换参数λ，并返回变换后的数据。
腾讯云提供了一系列与数据分析和机器学习相关的产品和服务，可以帮助用户进行数据处理和建模。例如，腾讯云提供的云服务器、云数据库、云存储等基础设施服务可以支持数据分析的计算和存储需求。此外，腾讯云还提供了人工智能平台、大数据平台和机器学习平台等高级服务，可以帮助用户进行数据挖掘、模型训练和预测分析。
更多关于腾讯云数据分析和机器学习相关产品的信息，可以参考以下链接：
请注意，以上答案仅供参考，具体的产品选择和使用需根据实际需求进行评估和决策。
相关·内容
2回答
是否有人在MATLAB中实现了Box Cox方法来转换多元数据？我已经找到了一个R ，但是它并不适用于MATLAB。谢谢。
浏览 1提问于2013-03-19得票数 1
这是我的输入数据： ? 这是对列r、f和m应用变换的所需输出，并将结果追加到原始数据 ? 代码如下： import pandas as pd
import numpy as np
from sklearn.preprocessing import StandardScaler
from sklearn.compose import ColumnTransformer
from sklearn.preprocessing import PowerTransformer
df = pd.DataFrame(np.random.randint(0,100,size
浏览 42提问于2019-02-08得票数 4
Python中有没有什么包可以进行数据转换:缩放、居中和Box-Cox转换，以消除数据的不对称性？在R中，这可以使用caret包来完成：
set.seed(1)
predictors = data.frame(x1 = rnorm(1000,
mean = 5,
sd = 2),
x2 = rexp(1000,
rate=10))
requ
浏览 1提问于2015-11-27得票数 6
回答已采纳
2回答
生成数据集的所有排列
、、、
我有一个数据框，看起来像这样：
df1 = pd.DataFrame({'Gene':['TP53', 'COX5', 'P16'], 'test':[1,3,0], 'Healthy':[0,0,2]})
Gene test Healthy
0 TP53 1 0
1 COX5 3 0
2 P16 0 2
我一直在尝试创建所有可能的值的排列。这个想法是将第一个基因"TP53“和它在”测试“列中的值映射到
浏览 16提问于2021-02-21得票数 0
回答已采纳
1回答
在R中运行Cox PH回归后，我需要在数据框中添加预测危险比的列。数据框是一个面板数据，其中numgvkey if公司标识符和age是时间标识符。您可以从以下链接下载日期的一小部分：
我做了以下几件事：
library(survival)
library(readstata13)
sme <- read.dta13("sme.dta")
reg<-coxph(Surv(age,EVENT2)~L1FETA+frailty(numgvkey), ties=c("efron"), data=sme)
summary(reg)
hr <- pred
浏览 0提问于2015-05-29得票数 5
5回答
反向Box-Cox变换
、、、、
我使用对连续变量执行。
from scipy.stats import boxcox
import numpy as np
y = np.random.random(100)
y_box, lambda_ = ss.boxcox(y + 1) # Add 1 to be able to transform 0 values
然后，我拟合一个统计模型来预测这个Box-Cox变换变量的值.模型预测在Box-Cox标度中，我想将它们转换为变量的原始标度。
from sklearn.ensemble import RandomForestRegressor
rf = RandomForestRegre
浏览 1提问于2014-10-15得票数 28
回答已采纳
我是R中的新手，我尝试在R中使用statistics ()函数，使用Nagelkerke R平方作为引导参数。我知道我需要一个函数来测量重新采样时原始数据的Nagelkerke R平方。然而，我不知道我应该把什么作为统计函数。
我知道Nagelkerke平方可以通过使用偏差和Null.deviance来计算，给定Logit回归。我编写了计算Nagelkerke R平方的函数。
NagR2 <- function(Objects){
n <- nrow(Objects)
reg <- glm(form,
family
浏览 70提问于2016-08-24得票数 1
1回答
我正在对线性回归模型执行Box转换，但在R中实现这个转换有一些困难。
我已经得到了Box-Cox的最优lambda值。我试图在R中执行以下操作
opt.lambda.test<- lapply(temp.list, function(x) with(x, x[which.max(y)]))
从Box-Cox中提取出最佳的lambda。temp.list是回归结果的列表。
opt.lambda.test
[[1]]
[1] 0.77
[[2]]
[1] 3.46
然后我把它转化成一个向量：
vec.lambda <- unlist(opt.lambda.test)
在本例中，我有一
浏览 1提问于2013-09-25得票数 2
回答已采纳
我正试着在画布‘DrawingContext’上用几个RenderTransforms绘制的Rect上做一些命中测试。我对RectangleGeometry应用了相同的旋转/缩放/平移变换，这样我就可以得到最终图形的边界。但是，RectangleGeometry上的.Bounds属性返回轴对齐的边界框，而我需要定向的边界框。
简而言之，我需要一些方法来找到实际转换后的矩形的角，这样我就可以对它执行命中/碰撞测试。
浏览 4提问于2011-01-13得票数 0
回答已采纳
我知道相机在0,0,0，我需要围绕它旋转世界，但我对平移和旋转的顺序感到困惑。
如果有一个理论上的x，y，z坐标系，相机在cx，cy，cz，它的方向是cox，coy，coz，我有一个立方体，它位于bx，by，bz方向是盒子，男孩，boz，那么需要什么系列的glTranslatef和glRotatef才能正确地旋转盒子，并在正确的位置远离相机？
以下是基本的操作，但我不知道如何放置它们，以及需要哪些其他操作才能使其按预期显示。
gl.glLoadIdentity();
// rotation and translation for cube
gl.glRotatef(box, 1,0,0);
浏览 0提问于2010-07-21得票数 1
回答已采纳
1回答
我使用DirectX11和DXUT显示点云，现在希望通过鼠标移动交互地选择点。特别地，当我拖动鼠标绘制矩形时，将选择该矩形中的点。我的方法是首先将鼠标位置转换为一个框，范围(-1，1)。然后利用mvp矩阵将点变换到投影空间。最后，将这两个点的x和y与方框进行比较。如果在框中，重点是selected.Here是我代码的某些片段。
CModelViewerCamera g_Camera;
XMMATRIX mView = g_Camera.GetViewMatrix();
XMMATRIX mProj = g_Camera.GetProjMatrix();
XMMATRIX mvp= mWor
浏览 9提问于2014-08-05得票数 0
回答已采纳
1回答
R中个体生存曲线的绘制
、、、、
在R中，我将我的数据库安排为一个计数过程，以应用扩展的Cox模型(具有时变协变量)：端点是事件的时间或审查的时间，而切点都是数据中的事件时间：
newdatabase <- survSplit(database,cut=eventTimes,
end=time_to_event_variable,
event=Status,start="start",id="newID")
object<-coxph(Surv(newdatabase$s
浏览 2提问于2014-07-24得票数 1
1回答
我在写一个简单的射线追踪器。当我遇到一个问题时，我正在研究我的相机矩阵。当我改变看起来向量，使它向下，图像，应该上升，实际上也下降了。同样地，如果我改变了看向量，使它朝上，图像也随之而去。据我所知，相机矩阵的工作是将世界坐标转换为眼睛坐标。由于我使用列矩阵，这是我的相机矩阵，$r$是正确的向量，$l$是lookAt向量，$u$是向上的向量
$$ \begin{b矩阵} r_x & u_x & -l_x &0\ r_y & u_y & -l_y &0\ r_z & u_z & -l_z &0\0&0&0&am
浏览 0提问于2018-02-17得票数 1
回答已采纳
1回答
我正在尝试建立R中具有时变协变量的Cox PH模型。我使用的是survival软件包中的coxph函数。虽然其中一个协变量的系数值太大，特别是2.5e+32，但在估计过程中没有任何问题。 ? 我猜不出这个问题的原因是什么，也不知道如何解决它。这个变量是非平稳的，并且违反了比例假设。这两个事实中的任何一个都可能导致这么大的系数值吗？
浏览 46提问于2021-10-26得票数 0
1回答
我试图在R中复制cox比例风险模型的估计，使用以下数据。
stata中的命令如下：
stset enddate2009, id(VPFid) fail(warends) origin(time startdate)
stcox HCTrebels o_rebstrength demdum independenceC transformC lnpop lngdppc africa diffreligion warage if keepobs==1, cluster(js_country)
Cox regression -- Breslow method for ties
No. of sub
浏览 3提问于2015-11-27得票数 5
回答已采纳
我需要将一个文件中的一个列条目与另一个包含列表(在shell中)的文件中的一个列条目相匹配。我使用的awk命令只匹配列表中的第一个单词，并不扫描列字段中的整个列表。
文件1如下所示：
chr1:725751 LOC100288069
rs3131980 LOC100288069
rs28830877 LINC01128
rs28873693 LINC01128
rs34221207 ATP4A
文件2如下所示：
Annotation Total Genes With Ann Your Genes With A
浏览 10提问于2016-08-23得票数 0
我有一个有许多列(字段)的表。在第一个字段中，我只需要保留唯一的值。在随后的列中，我需要计算第一列中存在的原始值数，但前提是给定列中的值大于0。
我已经成功地使用awk完成了部分任务，但我目前的尝试需要手动为表中的每一列创建一个数组，并手动为print命令键入每个数组。这不太可行。
任何帮助/建议(和解释一个潜在的解决方案如何工作)将不胜感激。
下面是输入表的一个子集(它已经在第1列上排序了)：
ATP6 93.883156 55.84006
COX1 230.708456 63.109
COX2 179.993226 74.224
浏览 4提问于2022-01-27得票数 1
回答已采纳
4回答
我需要将一些数据转换成“正常形状”，我读到了Box-Cox可以识别用于转换数据的指数。
因为我所理解的
car::boxCoxVariable(y)
用于线性模型中的响应变量，以及
MASS::boxcox(object)
为一个公式或合适的模型对象。因此，因为我的数据是dataframe的变量，所以我可以使用的唯一函数是：
car::powerTransform(dataframe$variable, family="bcPower")
对吗？还是我漏掉了什么？
第二个问题是在我获得
Estimated transformation parameters
dataframe$
浏览 11提问于2015-11-30得票数 40
我想用Cox比例风险模型的coxph R函数从惩罚的R包中得到对肺数据集的预测。
我有下面的示例代码。
library(survival)
library(pec)
library(penalized)
data("lung")
data <- lung
trainind <- sample(1:n,n*0.7)
testind <- (1:n)[-trainind]
frm <- as.formula(paste("Surv(time, status)~",paste(names(data[,-c(2,3)]), col
浏览 7提问于2017-09-11得票数 0
回答已采纳
我最近对空间点模式的研究使我开始研究LGCP过程。在我的例子中，我实际上有一系列的点模式，我想要适应一个单一的模型。我之前的一次询问让我用mppm来训练这样的模特(谢谢禤浩焯！)我的下一个问题是在mppm环境下使用这种Cox模型。
这是否可能将不均匀的LGCP Cox过程(或其他类型的Cox过程)与使用mppm的复制点模式相匹配？我看到了一些关于拟合Gibbs过程的信息，但实际上并不适用于Cox过程。似乎答案可能是“可能”通过创造性地使用“随机”论点。
举个例子，假设我用一个协变量X(这是一个im)来拟合一个使用点模式Y。对kppm的呼吁是：
myModel = kppm(Y ~ X,"
浏览 4提问于2021-03-31得票数 0
回答已采纳
1回答
我开始看openGL的教程，为了创建世界的实体，我使用了一个泛型类，它只保存了三件事:元素的位置，比例和旋转。有了它，任何时候都会渲染元素，创建一个模型矩阵，将这三个元素相乘。
我的问题是:我如何通过这三件事来获得这个元素的中心位置？
例如，我们在(0，0，0)位置有一个元素，其比例为(10，10，10)，旋转角度为(0，0，0)。我如何知道中心位置是什么？这取决于相机还是类似的东西？
编辑：
更多信息：
float ar = getAspectRatio();
cam->setOrthoParams(ar*-100.f, ar*100.f, -100.f, 100.f, -1
浏览 4提问于2014-02-15得票数 0
我的应用程序的MainActivity类由大约10个选项卡的ViewPager组成。在每个标签中，我使用android RecyclerView加载ListItem。当用户在ViewPager中从右向左滑动时出现问题。调试后，我发现这个问题出现在我的CustomListAdapter内部。
我发现如果我增加onCreateViewHolder内部布局充气的数量，滞后就会增加。
@Override
public RecyclerView.ViewHolder onCreateViewHolder(ViewGroup parent, int viewType) {
RecyclerView
浏览 0提问于2016-11-06得票数 0
我一直在编写一个游戏，只是为了在java上变得更好。在让球员轮换正常工作方面，我遇到了很多麻烦。我的第一个方法使用了以下内容
g2.setTransform(AffineTransform.getRotateInstance(radAngle,x_pos + (img.getWidth() / 2),y_pos+(img.getHeight() / 2)));
然而，这会导致所有图像随着播放器旋转，从而使射击和瞄准完全失灵。我正在研究，并看到有人使用这个代码来让他们的球员轮换。
Graphics2D g2 = (Graphics2D)g;
AffineTransf
浏览 1提问于2011-07-31得票数 2
1回答
我有一个包含许多列的数据集。我想对每一列进行Cox回归，所以我使用了"apply“。miRNA名称是预测变量，我想在分离的模型中进行Cox回归。这是代码：
mirna_names <-colnames(DB)[33:54]
cox_cont = sapply(mirna_names, function(x){
formula = as.formula(paste('Surv(years, AD)~', x))
coxFit = coxph(formula, data = DB)
summary(coxFit)$coefficients[,c(2,3,5
浏览 7提问于2022-06-01得票数 0
我想要创建一个虚拟的，其中包括NA在内的所有单元格变成1，而所有其他单元格-1。我试着做不同的事情，但都没有产生我想要的结果。
1)这只能导致NA值。
Cox.Reg$FirstTask <- ifelse(Cox.Reg$active_task_avg_depth == NA, 1, -1)
2)这只会导致新列中显示的-1值。
Cox.Reg$FirstTask <- ifelse(Cox.Reg$active_task_avg_depth == "NA", 1, -1)
浏览 1提问于2019-12-08得票数 0
回答已采纳
我正在建立一个Cox PH统计模型来预测乳腺癌患者复发的概率。我想为病人或医生创建一个在线界面，让他们能够输入相关的病人特征，并计算复发的概率。下面是一个完美的例子，尽管它适用于前列腺癌：
我的基本计划是用python创建工具，并根据用户的输入和我以前拟合的Cox PH模型用R计算概率。该工具需要下拉框和用户输入的数值。问题是，我从来没有用Python进行过任何web编程或GUI编程。我的经验是在Python的科学编程方面(例如Pylab等)。我的问题是：
我需要哪些与Python和R相关的软件包？从我做过的一些谷歌搜索来看，RPy和Tkinter似乎是明确的选择。
我如何存储统计模型
浏览 7提问于2012-07-11得票数 0
我正在尝试使用程序包MASS中的boxcox函数创建一个图。
但它正在创建一个rapache错误。
R代码：
<%
csvDF<- read.csv(GET$name1, header=TRUE)
a<-lm(csvDF[,GET$col_variable]~1)
require(MASS)
filename1 <- paste(tempfile(tmpdir='/var/www/images'), '.png', sep='')
png(filename1)
bx<-boxcox(a)
浏览 4提问于2011-04-15得票数 0
1回答
我想使用mlr在R中的右删失生存数据上运行xgboost，xgboost代码列出了一个目标函数存活:cox，它说：
生存: Cox :对右截尾生存时间数据的Cox回归(负值被认为是右截尾)。
我正在使用的Mlr 2只支持回归和分类学习者的xgboost。如果我尝试使用内置回归学习者的xgboost，它使用mse作为评估标准。因此，我尝试将度量更改为cindex，并得到了错误。
衡量标准: cindex索引
在乐趣中出错(X[i]，.)：测量cindex不支持任务类型regr！
于是我试着为xgboost写一个新的求生学习者，它只是回归学习者的一个副本，但随着"Regr“
浏览 2提问于2019-06-24得票数 5
回答已采纳
2回答
我正在使用R中的coxme函数{coxme}运行混合效果Cox模型，我想检查比例风险的假设。
我知道PH假设可以用cox.ph模型上的生存函数{ cox.zph }来验证。
但是，我找不到coxme模型的等价物。
2015年，类似的问题也在上发布，但没有得到答案。
我的问题是: 1)如何在混合效应cox模型coxme上检验PH假设? 2)如果coxme模型没有等价的cox.zph，那么在科学文章中发表运行混合效应coxme模型但在与coxme模型相同但没有随机效应的cox.ph模型上检验PH假设是否有效？
提前感谢你的回答。问候
浏览 10提问于2016-12-30得票数 0
1回答
我正在使用jQuery插件Datables：
下面是我的示例：
我的示例表中有几列，其中一列是"Name“。我添加了几个相同的名字，但一个以大写开头，另一个全部以小写开头。例如:阿什顿·考克斯阿什顿·考克斯。然而，我确实相信数据表中使用的默认排序是不区分大小写的，这意味着排序是“稳定的”(如果我错了，请纠正我)。
我的问题是，当我对它进行排序时，列变得不一致。如果它是升序的，大写字母将排在第一位，即第一行是Ashton cox，然后是ashton Cox。如果它是降序的，最后一行是ashton cox，倒数第二行是Ashton Cox。
因此，如果最后一行是相同的单词但大写，我如何编写
浏览 0提问于2014-08-19得票数 0
我想旋转一个JavaFX矩形。但是当我旋转它，试着得到它的X和Y时，这是和以前一样的X和Y。我使用以下代码：`
Rectangle r = new Rectangle(100,100,40,40);
r.getX();
r.getY();
r.setRotate(85);
r.getX();
r.getY(); `
getX()和getY()在rotate方法前后返回相同的结果。所以我的问题是，我可以做些什么来让左上角的X和Y位于他的新位置。下面是一个示例：
浏览 30提问于2018-01-12得票数 1
我很抱歉，这是一个非常普遍的问题，但我的谷歌技能完全失败了。
根据Cox PH测试，我想绘制每个效果大小(或至少0.1)的样本大小。一位同事(现在左)在R中产生了一个线性模型幂分析，作为Cox回归的占位符。有一些软件包(pwr)，但没有一个处理Cox回归。
我需要在R中这样做，这样我就可以将Cox分析的结果叠加在这个幂分析之上，来指示我的结果有多准确。
这是一个截图，我已经展示了什么，本质上是什么，我想重新创造，而不是通过眼睛弄虚作假。
浏览 1提问于2018-10-05得票数 0
回答已采纳
2回答
背景
我想绘制一个生存数据集的风险比率，包括它的置信区间。例如，我将从survival包中获取一个简化的数据集:冒号数据集。
library(survival)
library(tidyverse)
# Colon survival dataset
data <- colon %>%
filter(etype == 2) %>%
select(c(id, rx, status, time)) %>%
filter(rx == "Obs" | rx == "Lev+5FU") %>%
mutate(rx =
浏览 4提问于2021-03-02得票数 1
回答已采纳
我想通过计算方差通货膨胀因子(VIF)来评估cox比例风险模型中的多重共线性。像{car}这样的包中的vif-函数不接受coxph对象。
有没有办法计算R中cox模型的VIF？
浏览 12提问于2014-05-07得票数 3
回答已采纳
2回答
假设我有一个非正态分布的响应变量和一个解释性变量。让我们先创建这两个变量(以R编码)：
set.seed(12)
resp = (rnorm(120)+20)^3.79
expl = rep(c(1,2,3,4),30)
我运行了一个线性模型，我意识到残差不是正态分布的。(我知道运行夏皮罗可能不足以证明残差不是正态分布的，但这不是我问题的重点)
m1=lm(resp~expl)
shapiro.test(residuals(m1))
0.01794
因此，我想转换我的解释性变量(例如，寻找一个使用Box-Cox的转换)。
m2=lm(resp^(1/3.79)~expl)
shapiro.te
浏览 1提问于2013-08-29得票数 2
回答已采纳
我使用以下代码获得了一个名为bcvCindex$BootstrapCrossValCindexMat的元素列表： library(prodlim)
set.seed(13)
dat <- SimSurv(100)
library(survival)
library(randomForestSRC)
cox12 <- coxph(Surv(time,status)~X1+X2,data=dat,x=TRUE,y=TRUE)
cox1 <- coxph(Surv(time,status)~X1,data=dat,x=TRUE,y=TRUE)
cox2 <- coxph(S
浏览 9提问于2019-04-13得票数 0
回答已采纳
如果选择不可约多项式$m(x) = x^8+x^4+x^3+x+1$，甚至对于任何其他值，对于$01$将不存在乘法逆，因为$0000 0001$将完美地除以$m(x) = 100011011$，使余数为零。
那么，在AES S框中，01美元的价值如何映射为7c美元呢？
浏览 0提问于2018-03-10得票数 1
我有一个，它描述了三种类型的服务(酒吧、餐馆和快递)在十年时间内的死亡率。
数据包含三个变量:服务类型(1=saloon、2=restaurant和3=express)、年份(1到11的整数，其中11表示大于10年)和审查。
我有两个问题：
1)我已经拟合了Cox比例风险模型，但是有什么方法可以检验比例风险假设呢？也就是说，我们假设每个人的危险比率和基线危险与时间无关。
2)如何在R中拟合时间依赖的Cox模型？
下面是我的代码：
#Cox Proportional Hazards
cox <- coxph(Surv(Years, Censor) ~ data$`Service
浏览 4提问于2016-03-08得票数 1
我有图像，我希望它在布局中是全屏的。然后，我想将图标放在图像的特定位置；问题是当我在多种屏幕尺寸上运行应用程序时，图标并不位于同一位置。
唯一可以正常工作的是，我对xml布局文件中的图像使用了固定的宽度和高度。但这对我来说还不够，我希望它是全屏大小的，即使图像是拉伸的，也可以在所有屏幕尺寸和密度上使用。有人解决了这个问题吗？
实际上，我使用的是以下代码：
public class MainActivity extends Activity {
double [] X= {60,100,140};
double [] Y= {65,105,145};
//Parameter for icon
浏览 0提问于2013-05-16得票数 2
回答已采纳
从我看到的cox回归的例子来看：
，
数据还包括状态或事件列。
但是，如果我想对cox模型进行训练，然后使用经过训练的模型进行预测，我应该在测试数据中保留“状态”列吗？
此外，我们应该如何预测测试样本的可能存活时间？
浏览 0提问于2018-06-20得票数 0
回答已采纳
我想在R中使用BoxCoxTrans函数来解决偏斜问题。
但是，我有一个问题，不能得到结果作为数据帧。这是我的R代码。
df<-read.csv("dataSetNA1.csv",header=TRUE)
dd1<-apply(df[2:61],2,BoxCoxTrans) #Except independent variable that located first column, All variables are numeric variable.
dd1
$LT1Y_MXOD_AMT
Box-Cox Transformation
96249 data
浏览 0提问于2017-09-29得票数 0
1回答
我与R中的coxme包做了斗争，我想使用一个函数，比如survfit() --它通常用于coxph()模型--来绘制调整后的生存曲线，并在不同的参数值下找到中位生存期。
如果我使用coxph (无随机效应)来拟合模型，我可以做以下工作：
library(KMsurv)
data(burn)
my.surv <- with(burn, Surv(T1, D1))
cox_nr = coxph(my.surv ~ Z1 , data = burn)
survfit(cox_nr, newdata = data.frame(Z1 =1))
这提供了生存估计。但如果我把同样的模型和臀部搭配
浏览 2提问于2012-10-31得票数 10
2回答
如何从ITK配准中获得变换仿射？
、、、、
在三维核磁共振扫描A、B和C的情况下，我想在A上执行B的仿射(Co)配准，获取注册的仿射矩阵，并将其应用于C。
我的问题是，配准变换的仿射矩阵有错误的符号。也许是因为错误的定位？
TransformParameters包含12个值，其中前9个是按行大顺序排列的旋转矩阵，最后3个是平移值。
TransformParameters = [R1, R2, R3, R4, R5, R6, R7, R8, R9, Tx, Ty, Tz]
registration_affine = [[R1, R2, R3, Tx],
[R4, R5, R6, Ty],
浏览 67提问于2022-03-11得票数 2
回答已采纳
1回答
我已经为R中的一个物种计算了一个核密度估计(KDE)，为了简单起见，我在这里使用了来自spatstat的蚂蚁数据集(蚂蚁巢)，并且我希望有一个图例标度，每个测量单位的蚂蚁巢数(例如米)。如果我正确理解KDE图例显示0到1之间的概率密度，我如何在“真实世界”密度中转换这个概率，例如每平方米的点数(蚂蚁巢)？
这里的例子是：
require("spatstat")
data(ants)
dat <- ants
# estimate bandwith
h_cox <- bw.diggle(dat)
# calculate KDE
kd_cox <- densi
浏览 1提问于2020-05-11得票数 0
回答已采纳
1回答
数据解释和可视化
、、
我试着用母亲的年龄和婴儿的性别作为协变量来拟合考克斯比例风险模型。
我在用R语言表示数据时遇到了问题。
到目前为止，我将发布我的代码：
这些是我的包裹：
#Getting started:
# load libraries
pkgTest <- function(pkg){
new.pkg <- pkg[!(pkg %in% installed.packages()[, "Package"])]
if (length(new.pkg))
install.packages(new.pkg, dependencies = TRUE)
sapp
浏览 2提问于2022-04-04得票数 0
回答已采纳
在this tutorial之后，我正在尝试将R中的多变量Cox回归分析应用于我的数据集。特别是，我正在尝试应用以下函数coxph() install.packages(c("survival", "survminer"));
library("survival");
library("survminer");
data("lung");
res.cox <- coxph(Surv(time, status) ~ age + sex + ph.ecog, data = lung)
summary(res
浏览 49提问于2020-10-08得票数 1
回答已采纳
xgboost包启用了使用参数参数进行生存建模的方法: objective =“生存期:cox”和eval_metric = "cox-nloglik“。
结果模型的预测方法只输出风险分数(与type = "risk"在survival::coxph函数中的r值相同)。
如何使用xgboost来预测整个生存曲线？
浏览 0提问于2019-12-22得票数 4
回答已采纳
如果allowEditing为真，当用户从图库(或相机)加载照片时，他们可以裁剪图像以适合正方形。
是否可以更改该框的高度/宽度，这样它就可以不是正方形，而是一个矩形？
浏览 1提问于2016-06-15得票数 0
我进行了Cox回归分析，包括四个变量(性别、年龄和两个二元解释变量)，这些变量都与结果有显著的关联。我在R中使用了“生存”包中的coxph函数：
library(survival)
cox <- coxph(Surv(time, status_risk==1) ~ sex + age + stone_number +stone_size, data=cox_cut)
summary(cox1_3_cut)
Call:
coxph(formula = Surv(time, status_risk == 1) ~ sex + age +
stone_number + sto
浏览 2提问于2015-05-12得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
