demo/debug_identification.R

In djhshih/orient-bias-filter: Orientation Bias Filter

library(io);
library(devtools);

load_all("..");

eps <- logistic(-15);

###

d <- qread("sim-data/obrs_low-signal_low-damage_data.tsv");

theta <- 0;

phis <- seq(eps, 1 - eps, by=1e-2);
lps <- lapply(phis, function(phi) -lp_bases_given_params(d, theta, phi));

phi_hat <- estimate_phi(d, theta, 0.5);

plot(phis, lps)

-lp_bases_given_params(d, theta, logistic(-15))
-lp_bases_given_params(d, theta, 0.5)

theta_hat <- estimate_theta(d, 0.5, phi_hat);
theta_hat
estimate_theta(d, 0.5, logistic(-15))




###

d <- qread("sim-data/obrs_low-signal_high-damage_data.tsv");

phi <- 0.1;


thetas <- seq(eps, 1 - eps, by=1e-2);
lps <- lapply(thetas, function(theta) -lp_bases_given_params(d, theta, phi));

plot(thetas, lps)

theta_hat <- estimate_theta(d, 0.5, phi);

-lp_bases_given_params(d, 0.1, phi)
-lp_bases_given_params(d, 0.0554461, phi)
-lp_bases_given_params(d, logistic(-15), phi)
-lp_bases_given_params(d, theta_hat, phi)

####

d <- qread("sim-data/obrs_no-signal_high-damage_data.tsv");

phi <- 0.1;

eps <- 1e-6;

thetas <- seq(eps, 1 - eps, by=1e-2);
lps <- lapply(thetas, function(theta) -lp_bases_given_params(d, theta, phi));

plot(thetas, lps)

theta_hat <- estimate_theta(d, 0.5, phi);

-lp_bases_given_params(d, 0.0295593, phi)
-lp_bases_given_params(d, theta_hat, phi)