# Copyright (c) 2014  Takafumi Kanamori [kanamori@is.nagoya-u.ac.jp]
# All rights reserved.  See the file COPYING for license terms.

source('functions.r')

## setup of sample distribution
ntrain <- 50   ## size of observations
dim <- 2       ## dim of x
contam_ratio <- 0.2           ## ratio of outliers: 1-c
in_mean  <- 0; in_var  <- 1;  ## parameter of target distribution
out_mean <- 3; out_var <- 1;  ## parameter of outlier distribution
## parameter of estimator

## observations
m <- sum(rbinom(ntrain,1,prob=contam_ratio)) ## num of outliers
n <- ntrain-m 
ix <- matrix(rnorm(n*dim,mean=in_mean,sd=sqrt(in_var)),  n,dim)
ox <- matrix(rnorm(m*dim,mean=out_mean,sd=sqrt(out_var)),m,dim) ## outliers
x <- rbind(ix,ox)  ## contaminated observations

## estimation
res <- est_UnnormalizedModel(x,gamma=0.5)
print(res)  ## print result