WOLFRAM

Wolfram 语言与系统参考资料中心

狒狒下树的时间

在肯尼亚安博塞利保护区的一项研究对狒狒群体的进行了几天的观察，记录了半数群体从它们睡觉的树上下来的时间。某些天，观察员提前达到并准确记录。其他天中，他们在下树时间的中位数后达到，因此对当天的数据进行左截断.

    
Wolfram Language code: baboon = EventData[152];
对左截断的数据创建一个 Nelson–Aalen 模型并估计中位数：
Wolfram Language code: bModel = SurvivalModelFit[baboon, Method -> "NelsonAalen"]; bModel["MedianSurvival"]
计算下树时间的中位数大于9小时的概率：
Wolfram Language code: Probability[t > 9, tbaboon]
可视化累计下树概率的 95% 同时置信区间：
Wolfram Language code: Plot[{bModel["CDF"][t], bModel["LowerHallWellnerBand", "CDF"][t], bModel["UpperHallWellnerBand", "CDF"][t]}, {t, 6.5, 12}, PlotPoints -> 500, Exclusions -> None, AspectRatio -> 1, Filling -> {2 -> {3}}, PlotStyle -> {Blue, Purple, Purple}, PlotRange -> {{7, 10}, {0, 1}}]