生物统计与数学模型第五周quiz（3-31）

1.1 计算总重：3.204+2.834+1.56=？g

答：两位有效小数，=7.60 g

1.2 计算C: 2.3850*6.46/2.224=?

答：相对舍入误差0.00005/2.3850=2.096E-5

0.005/6.46=0.00077=0.077%

0.0005/2.224=0.00022=0.022%

6.92765*0.077%=0.00536，即2位有效小数，=6.93 mmol

2.用循环语句求exp(2.5)，作散点图显示截断误差和n=5~15的关系

答：exp(2.5)=12.182494

以下是算出n=5~15时exp(2.5)的值的代码

S=[0];
for k=5:15;
 s1=1;for n=1:k,s1=s1+2.5^n/factorial(n);end,S(k,n)=s1;
end;
S

计算结果为11.67057,12.00966,12.13076,12.16860,12.17912,12.18174,12.18234,12.18246,12.18249,12.18249,12.18249

以下是计算截断误差&绘图代码

x=[5,6,7,8,9,10,11,12,13,14,15];
y=[11.67057,12.00966,12.13076,12.16860,12.17912,12.18174,12.18234,12.18246,12.18249,12.18249,12.18249];
z=[0];
for i=1:11;
  z(i)=12.182494-y(i);
end;
scatter(x,z)

3.产生一个n=60，μ=10，σ20=4.0的正态数组，求σ2，和σ20作0.05假设检验。

答：代码如下，不指定alpha就是默认的0.05，H为0代表接受假设。

pkg load statistics;x=normrnd(10,2.0,60,1);
S=var(x)
[H, PVAL, CI, STATS] = vartest (x,4.0)

计算结果为

S =  5.3944
H = 0
PVAL =  0.076847
CI =

   3.8758
   8.0245

STATS =

  scalar structure containing the fields:

    chisqstat =  79.567
    df =  59

关于vartest的matlab文档：

 -- Function File: [H, PVAL, CI, STATS] = vartest (X, Y)
 -- Function File: [H, PVAL, CI, STATS] = vartest (X, Y, NAME, VALUE)
     Perform a F-test for equal variances.

     If the second argument Y is a vector, a paired-t test of the
     hypothesis 'mean (X) = mean (Y)' is performed.

     The argument "alpha" can be used to specify the significance level
     of the test (the default value is 0.05). 

     If H is 0 the null hypothesis is accepted, if it is 1 the null
     hypothesis is rejected. 

     A 100(1-alpha)% confidence interval is returned in CI.
     
     STATS is a structure containing the value of the test statistic
     (TSTAT), the degrees of freedom (DF) and the sample standard
     deviation (SD).

每个数据的误差仅为舍入误差，即 E=0.5

绘图程序、及运行结果

N=[36;89;141;200;250];E=[0.5;0.5;0.5;0.5;0.5];errorbar(N,E)

1.举例练习各个分布函数语句（4种分布，5个函数）

2.举例练习正态分布的参数估计

1.经济竞争力。

2.生物质预处理的能耗。

3.发酵菌株的效率，高效率转化木质纤维素预处理液中所有的糖组分，即戊糖和己糖共利用。

4.发酵菌株对木质纤维素预处理产生的不同毒性物质的抑制物耐受性。

5-1.热化学转化中气化和催化的能耗成本。蒸汽爆破法的成本。

5-2.生物转化中纤维素酶的成本、效率，抑制物产量、副产物。

6.水解工艺步骤的能耗和可持续性。

7.纤维素、半纤维素糖的转化率和产量。

8.工艺过程放大。

9.预处理废料和废水的回收。

10.最低成本工艺的整合。

11.产业化初期所需的投资。

12.国与国之间发展政策、经费的差异。

问题稍难，可以参考Altered sterol composition renders yeast thermotolerant这篇Science上的论文，结合自己的思考给出

1.全基因组的测序。

2.全基因组基因表达。

3.代谢流分析。

**4.建立RNAi系统，从RNAi菌株库中筛选。**

> 5.先找到关键酶，再找到编码这些酶的基因（往往多个）。

### 关于🐴🐘编辑器

我阅读的文献：

Affiliation: \

DOI: 10.1038/srep25763

impact factor：4.011（2018）

publish date: 10 May 2016

本论文包含以下数据：

该研究运用的数据处理方法：

作图法。折线图，同一坐标内的六根折线，代表了6种曝气策略。

作图法。散点图、折线图，同一坐标内