prise6/aVirtualTwins
1
0
Fork 0
mirror of https://github.com/prise6/aVirtualTwins.git synced 2024-05-01 19:52:43 +02:00
Code Issues Releases Wiki Activity

Nothing important

Browse source
This commit is contained in:
prise6 2015-06-10 14:38:36 +02:00
parent 8033da44eb
commit 82c56a4688
3 changed files with 9 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
R/tools.R
View file

@ -8,7 +8,7 @@ vt.getQAOriginal <- function(response, trt, ahat){
sum(response*ahat*(1-trt))/sum(ahat*(1-trt)) -
(sum(response*trt)/sum(trt) -
sum(response*(1-trt))/sum(1-trt))
}
}
return(tmp)
}
@ -32,7 +32,6 @@ vt.getIncidence <- function(df){
return(table.res)
}
vt.rr.snd <- function(vt.difft, selected){
if(sum(selected) == 0){
return(0)
@ -40,4 +39,8 @@ vt.rr.snd <- function(vt.difft, selected){
return((sum(vt.difft$twin1*selected*vt.difft$vt.object$data[, 2])/sum(selected*vt.difft$vt.object$data[, 2]))
/(sum(vt.difft$twin1*selected*(1-vt.difft$vt.object$data[, 2]))/sum(selected*(1-vt.difft$vt.object$data[, 2]))))
}
}
}

1
R/tree.R
View file

@ -12,6 +12,7 @@
#' @field tree rpart Rpart object to construct the tree
#' @field Ahat vector Indicator of beglonging to Ahat
#'
#' @import methods
VT.tree <- setRefClass(
Class = "VT.tree",

2
R/tree.class.R
View file

@ -3,6 +3,8 @@
#' A reference class to compute subgroups by classifiation tree
#'
#' @include tree.R
#'
#' @import methods
VT.tree.class <- setRefClass(
Class = "VT.tree.class",