R/chi2.R

Defines functions chi2

Documented in chi2

chi2 <-
function(data,alp=0.5,del=0.05){
#Phase 1
    dat <- data
    p <- ncol(dat)-1         
    
    alpha <- alp; delta <- del; d <- 0.1 
    inconRate <- incon(dat)
    Disc <- chiM(dat,alpha=alpha)$Disc.data
          
       while(inconRate<delta){
           if(alpha<=0.2) break
           alpha <- alpha-d
           Disc <- chiM(Disc,alpha=alpha)$Disc.data
               inconRate <- incon(Disc) 
       }
       
#Phase 2
    options(digits=3)
    eps <- 0.01
    sig <- array(alpha,p)
    d <- 0.05
    cutp=list()
    for(i in 1:p){
        while(TRUE){
             val <- value(i,Disc,alpha=sig[i])
             Disc <- val$disc
             inconRate <- incon(Disc)#inconRate for data
             sig[i] <- sig[i]-d
             if(inconRate>delta || sig[i]<=eps) break 
        }
     cutp[[i]] <- val$cuts
    }
return(list(cutp=cutp,Disc.data=Disc))
}

Try the discretization package in your browser

Any scripts or data that you put into this service are public.

discretization documentation built on June 9, 2022, 5:09 p.m.