r-mice

    1熱度

    1回答

    我有一個mids對象從mice創建。我想重新編碼一些推定變量並保留mids對象。我知道我可以用complete()將mids對象轉換爲「long」,但我想保留mids對象,因爲它有一些額外的用途。 以下是使用nhanes數據集的示例。運行mice()爲nhanes中的變量創建5個估算數據集。我正在關注hyp。 library(mice) names(nhanes) nhanes$hyp #[

    0熱度

    2回答

    我想在表格中計算一些缺失的數據,並在估算表上運行考克斯模型。 我可以得到插值運行在我的數據上,而cox模型運行在估算數據上,但我不明白如何從數據集中查看cox輸出,其中一些值被估算(即我特別需要輸出中的風險比和P值)。 的命令有: >library("mice") >Table <-read.table("TestTable",stringsAsFactors=TRUE,header=TRUE)

    1熱度

    1回答

    在小鼠中包提取完整的數據集可以使用complete()命令如下: install.packages("mice") library ("mice") imp1=mice(nhanes,10) fill1=complete(imp,1) fill2=complete(imp,2) fillall=complete(imp,"long") ,但可以有一個人告訴我如何提取阿梅利亞包完整的數

    0熱度

    1回答

    在使用MICE程序包對R進行插補之後,我想要生成應急表。適合顯示列表中的表,但如果我pool()他們,將引發以下錯誤:Error in pool(fit) : Object has no coef() method.我做錯了什麼? 這個簡單的例子重現錯誤: library("mice") imp <- mice(nhanes) fit <- with(imp, table(bmi, hyp)

    0熱度

    1回答

    所以我有一個包含很多缺失值的數據集。我想分離不同丟失模式的數據。我發現了包裝中的'老鼠',在總結缺失值模式時非常方便。但是,當我想選擇具有某種缺失模式的行時,選定行的數量遠少於缺少模式矩陣所建議的數量。 我的代碼如下。 要獲得丟失的模式: library(mice) # md.pattern returns a matrix, I convert the matrix into a data f

    1熱度

    1回答

    我希望能夠根據從字符串構造的公式,在mice()返回的mids對象上運行迴歸。對於一個普通的迴歸(沒有插補)這是很容易: ​​ 我覺得做對的推算數據集同樣的事情應該是這個樣子: imp <- mice(nhanes2,printFlag = FALSE) for (x in c('age','hyp','chl')) { regx <- with(data=imp,exp=lm(pas

    5熱度

    1回答

    我使用MICE開發了一個通過多重插補法建立的模型。我想使用這個模型來預測新觀察結果(不包含缺失數據)的響應,包括標準錯誤。將MICE中創建的模型對象傳遞給predict()函數不起作用 一個使用內置nhanes數據集的簡單示例。說我想開發一個邏輯迴歸模型的形式age==3 ~ bmi + hyp + chl,並使用此模型來預測,也就是說,概率(年齡= 3 | BMI = 20,HYP = 2,CH

    6熱度

    2回答

    我建立一個預測模型,現在用的是mice包在我的訓練集歸咎於港定居。因爲我需要重新使用相同的歸集方式爲我的測試集,我怎麼能重新應用到我的測試數據?在組合數據集 # generate example data set.seed(333) mydata <- data.frame(a = as.logical(rbinom(100, 1, 0.5)), b = as.logical

    1熱度

    1回答

    我有50個由鼠標創建的多重估算數據集,並且我在它們每個上運行Cox模型。我需要驗證一段時間內成比例危險的假設,爲此我將檢查cox.zph中的圖。 我想要做的是在同一圖上繪製來自每個數據集的樣條曲線。我知道這通常會在plot的電話後用lines()完成,但是當我用zph對象執行此操作時,R會連接每個點,這顯然不是我想要的。 我試圖繪製一次,然後在循環中繪製它。個人ZPH只是我所有的考克斯的名單分析

    1熱度

    1回答

    我已經在R中使用了鼠標包來計算我的數據中的某些缺失值,但並非針對所有變量。現在我想用原始數據中的列替換歸因數據中的列,如果列名相同的話。這裏是我的功能: replace_imp <- function(data,impdata) { for(i in 1:length(impdata)){ for(k in 1:length(data)){ if(colnames(