Efficacy_test {CUtools}R Documentation

Test to compare efficacy of two models for a percentage of misclassified events

Description

Test to compare the efficacy of two markers for paired or unpaired cases

Usage

Efficacy_test( paired, Prob1,Prob2,yt1,yt2,z)

Arguments

paired

if sample is paired 1 else 0

Prob1

A vector with the event probability values provided by the biomarker 1

yt1

A vector with the actual event values for the biomarker 1

Prob2

A vector with the event probability values provided by the biomarker 2

yt2

A vector with the actual event values for the biomarker 2

z

The misclassification rate at which the effectiveness of the marker will be estimated.

Details

Prob1 and Prob2 must be numeric vectors with values between 0 and 1, yt1 and yt2 numeric vectors with dichotomic values 0/1 and z a numeric value between 0 and 100. in a case of a paired comparison, yt1 and yt2 must be the same vector.

Value

The returned results of a test.

Efficacy_test

It gives the result of the comparison test between markers in terms of efficacy

Author(s)

Maria Escorihuela, Luis Mariano Esteban, Gerardo Sanz, Angel Borque

Examples

###We generate a marker to serve as an example
Prob1<-c(rnorm(10000,0.4,0.1),rnorm(10000,0.5,0.05))
Prob2<-c(rnorm(10000,0.4,0.1),rnorm(10000,0.5,0.05))
yt1<-rep(c(0,1),c(10000,10000))
yt2<-rep(c(0,1),c(10000,10000))
#We choose a rate of 10% for misclassified events.
##For a paired test

Efficacy_test(paired=1,Prob1,Prob2,yt1,yt2,z=10)

##For a unpaired test

Efficacy_test(paired=0,Prob1,Prob2,yt1,yt2,z=10)

[Package CUtools version 0.1.0 Index]