diff options
Diffstat (limited to 'R')
| -rw-r--r-- | R/interpweights.R | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/R/interpweights.R b/R/interpweights.R index f4f6cb50..c8702cdc 100644 --- a/R/interpweights.R +++ b/R/interpweights.R @@ -65,11 +65,11 @@ obj <- function(epsilon, vecpv, prob, support, cte){ return( 1 - crossprod(newprob, vecpv) - cte)
}
-optimize <- function(min, max, vecpv, prob, support, cte){
+tweak <- function(min, max, vecpv, prob, support, cte){
mid <- (min + max)/2
objective <- obj(mid, vecpv, prob, support, cte)
- while( abs(objective)>1e-6){
- if(objective>0){
+ while( abs(objective) > 1e-6){
+ if(objective > 0){
min <- mid
}else{
max <- mid
@@ -82,7 +82,7 @@ optimize <- function(min, max, vecpv, prob, support, cte){ interpweightsadjust <- function(w, v1, v2, vecpv){
interpweightsadjust <- interpweights(w, v1, v2)
- epsilon <- optimize(-0.5, 0.5, vecpv, interpweightsadjust, v2, 1)
+ epsilon <- tweak(-0.5, 0.5, vecpv, interpweightsadjust, v2, 1)
return( adjust_weights(interpweightsadjust, v2, epsilon) )
}
|
