Displaying 20 results from an estimated 900 matches similar to: "MICE for Cox model"
2013 Feb 14
2
Plotting survival curves after multiple imputation
I am working with some survival data with missing values.
I am using the mice package to do multiple imputation.
I have found code in this thread which handles pooling of the MI results:
https://stat.ethz.ch/pipermail/r-help/2007-May/132180.html
Now I would like to plot a survival curve using the pooled results.
Here is a reproducible example:
require(survival)
require(mice)
set.seed(2)
dt
2012 Mar 30
3
pooling in MICE
Hi everyone,
Does anyone here has experience using MICE to impute missing value? I am
having problem to pool the imputed dataset for a MANOVA test, could you
give me some advice please?
Here is my code:
> library(mice)
>
2006 Mar 01
1
mice library / survival analysis
Hello folks,
I am a relatively new user of R and created multiply imputed data sets
with the 'mice' library. This library provides two functions for
complete-data analysis on multiply imputed data set objects (lm.mids and
glm.mids). I am trying to estimate a series of Cox PH regression models
and cannot figure out the best way to do this. Is it possible with the
mitools library?
2007 May 31
0
Using MIcombine for coxph fits
R-helpers:
I am using R 2.5 on Windows XP, packages all up to date. I have run
into an issue with the MIcombine function of the mitools package that I
hoped some of you might be able to help with. I will work through a
reproducible example to demonstrate the issue.
First, make a dataset from the pbc dataset in the survival package
---------------
# Make a dataset
library(survival)
d <-
2013 Jan 30
3
Mac v Windows Mystery
Dear All,
I'm trying to track down a problem with my quantreg package reported by a user doing censored quantile regression.
When he runs the test4.R file attached below (which reads the csv file also attached) on his windows machine he gets
an error like this:
> Error in dimnames(B) <- list(c("tau", dimnames(x)[[2]], "Qhat"), NULL) :
> length of
2010 Jun 08
2
how to ignore rows missing arguments of a function when creating a function?
Hi,
I am relatively new to R; when creating functions, I run into problems with
missing values. I would like my functions to ignore rows with missing values
for arguments of my function) in the analysis (as for example is the case in
STATA). Note that I don't want my function to drop rows if there are missing
arguments elsewhere in a row, ie for variables that are not arguments of my
2009 Jun 04
0
Small mystery : passing a "subset=" argument to lme|lm through "..."
Dear list,
I have problems involving passing a "subset=" argument through "...".
I'm trying to augment the set of defined analyses for mice (homonymous
package) with a call to lme. This package create multiple imputations of
missing data in a "mids" object, each completed data set may be obtained
through the complete(data, set) function.
> sessionInfo()
R
2012 Jul 17
1
Threshold Quantile Regression code CRASHES in R
I am working on a two stage threshold quantile regression model in R, and my aim is to estimate the threshold of the reduced-form equation (call it rhohat), and the threshold of the structural equation (call it qhat), in two stages. On the first stage, i estimate rhohat by quantile regression and obtain the fitted values. I use these fitted values to estimate qhat on the second stage. The code is
2006 Sep 03
2
Running cox models
Hi,
I'm reading van Belle et al "Biostatistics" and trying to run a cox test using
a dataset from:
http://faculty.washington.edu/~heagerty/Books/Biostatistics/chapter16.html
(Primary Biliary Cirrhosis data link at top of the page),
I'm using the following code:
--------------- start of code
library(survival)
liver <-
2012 Jul 25
2
reshape -> reshape 2: function cast changed?
Hi,
I used to use reshape and moved to reshape2 (R 2.15.1). Now I tried some of my older scripts and was surprised that my cast function wasn't working like before.
What I did/want to do:
1) Melt a dataframe based on a vector specifying column names as measure.vars. Thats working so far:
dfm <- melt(df, measure.vars=n, variable_name = "species", na.rm = FALSE)
2) Recast the
2012 Nov 21
3
Reentrant NMIs, MCEs and interrupt stack tables.
Hello,
While working on a fix for the rare-but-possible problem of reentrant
NMIs and MCEs, I have discovered that it is sadly possible to generate
fake NMIs and MCEs which will run the relevant handlers on the relevant
stacks, without invoking any of the other CPU logic for these special
interrupts.
A fake NMI can be generated by a processor in PIC mode as opposed to
Virtual wire mode, with a
2009 Apr 22
1
How do I loop through strings?
I read in a CSV file with Data <-
read.csv(file="FileName.csv",head=TRUE,sep=",")
the file containts strings in the 2nd and 3rd columns and each has about
1000 rows. I need to either loop through the strings there looking for
strings that would trigger other logic or remove those rows and put them in
another array and loop through them in the new array. I can't seem
2020 Feb 11
2
[RFC PATCH 00/62] Linux as SEV-ES Guest Support
On Tue, Feb 11, 2020 at 03:50:08PM +0100, Peter Zijlstra wrote:
> Oh gawd; so instead of improving the whole NMI situation, AMD went and
> made it worse still ?!?
Well, depends on how you want to see it. Under SEV-ES an IRET will not
re-open the NMI window, but the guest has to tell the hypervisor
explicitly when it is ready to receive new NMIs via the NMI_COMPLETE
message. NMIs stay
2020 Feb 11
2
[RFC PATCH 00/62] Linux as SEV-ES Guest Support
On Tue, Feb 11, 2020 at 03:50:08PM +0100, Peter Zijlstra wrote:
> Oh gawd; so instead of improving the whole NMI situation, AMD went and
> made it worse still ?!?
Well, depends on how you want to see it. Under SEV-ES an IRET will not
re-open the NMI window, but the guest has to tell the hypervisor
explicitly when it is ready to receive new NMIs via the NMI_COMPLETE
message. NMIs stay
2005 Jun 09
1
single assignment affecting multiple sub-structures (PR#7924)
I'm trying to create a language structure that is a call to a function
with a number of arguments that is only known at run time. I do this by
using repeated indices to expand out a call with a single argument.
However, when I change one of the arguments, all are changed.
I don't see the same behavior when I initially create a call with
multiple arguments.
Even more strangely,
2012 Nov 22
41
[PATCH V3] vmx/nmi: Do not use self_nmi() in VMEXIT handler
The self_nmi() code cause''s an NMI to be triggered by sending an APIC
message to the local processor. However, NMIs are blocked by the
VMEXIT, until the next iret or VMENTER.
Volume 3 Chapter 27 Section 1 of the Intel SDM states:
An NMI causes subsequent NMIs to be blocked, but only after the VM exit
completes.
As a result, as soon as the VMENTER happens, an immediate VMEXIT
happens
2010 Oct 15
5
[LLVMdev] How do I find all memory allocations in an llvm ir code file?
I tried to compile this snippet of C++ code:
void FuncTest() {
int* a = new int;
int* b = new int[2];
}
using:
clang test.cpp -S -emit-llvm -o - > test.llvm
and obtained this:
define void @_Z8FuncTestv() {
entry:
%a = alloca i32*, align 4
%b = alloca i32*, align 4
%call = call noalias i8* @_Znwj(i32 4)
%0 = bitcast i8* %call to i32*
store i32* %0, i32** %a, align 4
2013 Mar 11
2
[LLVMdev] How to detect all free() calls
Thanks. isFreeCall() works well but for
%call2 = call i32 bitcast (i32 (...)* @free to i32 (i8*)*)(i8* %call1) nounwind, !dbg !16
So I tried to figure out when the above instruction occurred.
When <stdlib.h> is included, free(buf2R1); turn into call void @free(i8* %call1) nounwind, !dbg !16
when I forget to include <stdlib.h>, free(buf2R1); turn into %call2 = call i32 bitcast (i32
2012 Dec 12
7
[PATCH V5] x86/kexec: Change NMI and MCE handling on kexec path
xen/arch/x86/crash.c | 116 ++++++++++++++++++++++++++++++++++-----
xen/arch/x86/machine_kexec.c | 19 ++++++
xen/arch/x86/x86_64/entry.S | 34 +++++++++++
xen/include/asm-x86/desc.h | 45 +++++++++++++++
xen/include/asm-x86/processor.h | 4 +
5 files changed, 203 insertions(+), 15 deletions(-)
Experimentally, certain crash kernels will triple fault very early
2020 Feb 11
1
[PATCH 62/62] x86/sev-es: Add NMI state tracking
On Tue, Feb 11, 2020 at 5:53 AM Joerg Roedel <joro at 8bytes.org> wrote:
>
> From: Joerg Roedel <jroedel at suse.de>
>
> Keep NMI state in SEV-ES code so the kernel can re-enable NMIs for the
> vCPU when it reaches IRET.
This patch is overcomplicated IMO. Just do the magic incantation in C
from do_nmi or from here:
/*
* For ease of testing, unmask