日志

置顶 ·分享 签到积分

仅好友可见 2018-1-16 00:58

tab= ; f_06clim = @(x) sum((mod(x,tab(:,1))==0).*tab(:,2:3)); accum=411;y=344; y=y+1; accum=accum+sum(f_06clim(y)); , jinqian /3 +1 /4 +2 /7 +2 /8 +3 /9 +5 /30 +20 /50 +30 /100+50 gongxian /10 +1 /20 +2 /30 +2 /50 +3 /100+5 24..16 27..16 28..14? 30..31+3 32..15 35..12 36..18 40..15+3 42 ...

个人分类: zhiding|0 个评论

分享 Bootstrap 的置信区间

2018-2-12 04:38

Bootstrap 是一种非参数的统计推断方法。 通过电脑随机再抽样(resample)，可以得到一个统计量的直方图(histogram)如下。（盗图自 Efron 1981, http://www.jstor.org/stable/3314608 ） 直方图的百分位数可以很直观的给出一个置信区间，这是其中一种最简单的置信区间估计。（Basic percentile method） 独立同分布(iid ...

个人分类: matlab|833 次阅读|0 个评论

分享 相关系数 r 的误差估计-Bootstrap

2018-2-8 12:23

Bootstrap的好处是不需要太多假设。（不假设正态分布） 论文请看 http://www.jstor.org/stable/2685844 Efron, B., Gong, G. (1983). A leisurely look at the bootstrap, the jackknife, and cross-validation. The American Statistician , 37 (1), 36-48. Matlab代码： load lawdata rng default % For ...

个人分类: matlab|777 次阅读|0 个评论

分享 [WRF] 精简版的ncdump

2018-2-5 02:55

ncdump是NetCDF自带的显示nc文件的命令。由于WRF输出的文件(wrfout, met_em, ...)带有很多没用的attributes（属性），使得显示有点太长。 我会在我的~/.bashrc 加上以下三行，定义一个bash函数 ncdwrf。 ncdwrf(){ ncdump -h $1 |grep -Ev 'FieldType|MemoryOrder|stagger|coordinates|sr_x|sr_y' } s ...

个人分类: technical|149 次阅读|0 个评论

本页有 7 篇日志因作者的隐私设置或未通过审核而隐藏