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> &nbsp; <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"&lt;-
function(model.sim = law.sim, n = 15, nboot = 66)
{
#Parametric bootstrap of correlation coefficient
#Law sim contains bivariate normal data

        theta &lt;- rep(0, nboot)
        for(b in (1:nboot))
                theta[b] &lt;- corr(((b - 1) * n + 1):(b * n), law.sim)
        return(theta)
}


#Example of using the functions and plotting the results

n1 &lt;- rnorm(1000)
n2 &lt;- rnorm(1000)
law.sim &lt;- cbind(n1, 0.777 * n1 + sqrt(1 - (0.777^2)) * n2)
rbootpar &lt;- 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 + -
显示快捷键?