parametric bootstrap.htm
来自「matlab bootstrap程序设计方法」· HTM 代码 · 共 85 行
HTM
85 行
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<!-- saved from url=(0066)http://www-stat.stanford.edu/~susan/courses/b494/index/node55.html -->
<!--Converted with LaTeX2HTML 2K.1beta (1.57)original version by: Nikos Drakos, CBLU, University of Leeds* revised and updated by: Marcus Hennecke, Ross Moore, Herb Swan* with significant contributions from: Jens Lippmann, Marek Rouchal, Martin Wilck and others --><HTML><HEAD><TITLE>Parametric bootstrap</TITLE>
<META content="Parametric bootstrap" name=description>
<META content=t name=keywords>
<META content=document name=resource-type>
<META content=global name=distribution>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.2523" name=GENERATOR>
<META http-equiv=Content-Style-Type content=text/css><LINK
href="Parametric bootstrap.files/t.css" rel=STYLESHEET><LINK href="node56.html"
rel=next><LINK href="node54.html" rel=previous><LINK href="node52.html"
rel=up><LINK href="node56.html" rel=next></HEAD>
<BODY><!--Navigation Panel--><A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node56.html"
name=tex2html923><IMG height=24 alt=next
src="Parametric bootstrap.files/next.png" width=37 align=bottom border=0></A> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node52.html"
name=tex2html919><IMG height=24 alt=up src="Parametric bootstrap.files/up.png"
width=26 align=bottom border=0></A> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node54.html"
name=tex2html913><IMG height=24 alt=previous
src="Parametric bootstrap.files/prev.png" width=63 align=bottom border=0></A> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node124.html"
name=tex2html921><IMG height=24 alt=index
src="Parametric bootstrap.files/index.png" width=43 align=bottom border=0></A>
<BR><B>Next:</B> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node56.html"
name=tex2html924>Smoothed Bootstrap</A> <B>Up:</B> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node52.html"
name=tex2html920>The bootstrap : the</A> <B>Previous:</B> <A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node54.html"
name=tex2html914>S examples for simple</A> <B><A
href="http://www-stat.stanford.edu/~susan/courses/b494/index/node124.html"
name=tex2html922>Index</A></B> <BR><BR><!--End of Navigation Panel-->
<H2><A name=SECTION00313000000000000000>Parametric bootstrap</A> </H2>If one has
a parametric distribution <IMG height=17 alt=$F_\lambda$
src="Parametric bootstrap.files/img305.png" width=16 align=middle border=0>
which is well defined except for the pamameter <IMG height=34 alt=$\lambda$
src="Parametric bootstrap.files/img306.png" width=26 align=bottom border=0>,
instead of drawing with replacement from the original sample, one can draw from
the distribution <!-- MATH $F_{\widehat{\lambda}}$ --><IMG height=23
alt=$F_{\widehat{\lambda}}$ src="Parametric bootstrap.files/img307.png" width=16
align=middle border=0>, where <!-- MATH $\widehat{\lambda}$ --><IMG height=34
alt=$\widehat{\lambda}$ src="Parametric bootstrap.files/img308.png" width=26
align=bottom border=0> is the estimation obtained from the original sample.
<P><PRE>#Parametric Bootstrap
"boot.par"<-
function(model.sim = law.sim, n = 15, nboot = 66)
{
#Parametric bootstrap of correlation coefficient
#Law sim contains bivariate normal data
theta <- rep(0, nboot)
for(b in (1:nboot))
theta[b] <- corr(((b - 1) * n + 1):(b * n), law.sim)
return(theta)
}
#Example of using the functions and plotting the results
n1 <- rnorm(1000)
n2 <- rnorm(1000)
law.sim <- cbind(n1, 0.777 * n1 + sqrt(1 - (0.777^2)) * n2)
rbootpar <- boot.par()
hist(rbootpar)
par(new=T)
plot(density(rbootpar, width = 0.15), xlab = "", ylab = "",
lab = c(0, 0, 0),type = "l")
title("Parametric Bootstrap - Law School corr.", cex = 1.9)
</PRE>
<P><BR>
<HR>
<ADDRESS>Susan Holmes 2002-01-12 </ADDRESS></BODY></HTML>
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?