Last data update: 2014.03.03

R: Interval-Censored Local Linear Density Estimation
iclldeR Documentation

Interval-Censored Local Linear Density Estimation

Description

This is the local linear version of ickde.

Usage

icllde(I, h, f, m, n.iterations = 10, x1, xm, right.limit = 10000, kernel="gaussian")

Arguments

I

A matrix with two columns, consisting of left and right endpoints of the interval data

h

A scalar bandwidth

f

An initial estimate of the density at a sequence of grid points (optional)

m

The number of (equally-spaced) grid points at which the density is to be estimated

n.iterations

The maximum number of iterations allowed

x1

The minimum grid point (optional)

xm

The maximum grid point (optional)

right.limit

For right-censored data, the value to be used as an artificial right endpoint for the intervals

kernel

character argument indicated choice of kernel; current choices are "gaussian", "epanechnikov", "biweight"

Value

An object of class IC

Examples

     data(ICHemophiliac)
     estimate <- icllde(ICHemophiliac, m=200, h=.9)
     plot(estimate, type="l")

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(ICE)
Loading required package: KernSmooth
KernSmooth 2.23 loaded
Copyright M. P. Wand 1997-2009
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ICE/icllde.Rd_%03d_medium.png", width=480, height=480)
> ### Name: icllde
> ### Title: Interval-Censored Local Linear Density Estimation
> ### Aliases: icllde
> ### Keywords: models
> 
> ### ** Examples
> 
>      data(ICHemophiliac)
>      estimate <- icllde(ICHemophiliac, m=200, h=.9)
>      plot(estimate, type="l")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>