R/create_phiU.R

In TimothyHyndman/deconvolve: Deconvolution Tools for Measurement Error Problems

create_phiU <- function(errors, errortype, sd_U){

	if(errortype == 'laplace' & errors == "hom") {
		phiU <- function(tt) {
			1 / (1 + sd_U^2 * tt^2 / 2)
		}
	}

	if(errortype == 'normal' & errors == "hom") {
		phiU <- function(tt) {
			exp(-sd_U^2 * tt^2 / 2)
		}
	}

	if(errortype == 'laplace' & errors == "het") {
		phiU <- c()
		for (sigUk in sd_U){
			phiUk <- function(tt) {
				1 / (1 + sigUk^2 * tt^2 / 2)
			}
			phiU <- c(phiU, phiUk)
		}
	}

	if(errortype == 'normal' & errors == "het") {
		phiU <- c()
		for (sigUk in sd_U){
			phiUk <- function(tt) {
				exp(-sigUk^2 * tt^2 / 2)
			}
			phiU <- c(phiU, phiUk)
		}
	}

	phiU
}