生物统计与数学模型第五周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种曝气策略。
作图法。散点图、折线图,同一坐标内