R软件.docVIP

分析下列每行程序，每个命令，运行个命令的结果，并对结果进行解读。 对程序的注释要写在每行程序的后面，建议采用ppt展示注释后的程序，并将结果列出。 params - function(N, breaks, p=seq(0.001, 1, length=100)) { list(N=N, T=1/breaks, p=p, q=1-p) } #创建一个名为params的函数，该函数执行创建一个含有四个分量的列表。 pdfcomp - function(comp, params) { n - params$T p - params$p q-params$q y - round(comp/n) choose(n, comp)*p^comp*q^(n-comp) / (1 - q^n) }#创建一个名为pdfcomp的函数，该函数执行令“params”中的值y=对comp、 /n取整”，以及计算 # Expected num sherds (for a vessel) [=completeness] expcomp - function(params) { params$T*params$p/(1-params$q^params$T) } #定义一个名为“expcomp”的函数，实行表达式的值 # Variance of num sherds (for a vessel) varcomp - function(params) { n - params$T p - params$p q - params$q # From Johnson Kotz (n*p*q / (1 - q^n)) - (n^2*p^2*q^n / (1 - q^n)^2) # n^2 times Thomas Yees formula # n^2*((p*(1 + p*(n - 1)) / (n*(1 - q^n))) - (p^2 / (1 - q^n)^2)) } #定义一个名为“varcomp”的函数，并执行下列表达式 ,,， # Expected value of completeness (for a sample of vessels) expmeancomp - function(params) { expcomp(params) } #定义一个名为“expmeancomp”的函数，执行表达式 # Variance of completeness (for a sample of vessels) # Use the expected number of vessels in sample as denominator varmeancomp - function(params) { varcomp(params)/(numvess(params)) } numvess - function(params) { params$N*(1-params$q^params$T) } # ecomp - function(p, T, comp) { q - 1 - p T*p/(1 - q^T) - comp } #定义一个含参数p，T，comp的名为ecomp的函数执行 estN - function(comp, broke, n) { T - 1/broke n / (1 - (1 - uniroot(ecomp, c(0.00001, 1), T=T, comp=comp)$root)^T) } ##定义一个含参数comp，broken，n的名为estN的函数执行 , nvessscale - function(params, xlim, ylim, new=TRUE) { if (new) par(new=TRUE) plot(0:1, c(1, params$N), type=n, axes=!new, ann=FALSE, xlim=xlim, ylim=ylim) }#定义一个名为“nvessscale”的函数，在函数中永久的令new为真，绘制区间的图形，type=n表示不让点显示，但设置坐