問題:一批植物材料,設定了不同梯度的溫度,每個處理重複3次,並測定了相對電導率,想計算各種基因型的半致死溫度。
方法:採用文獻《應用logistic方程確定植物組織低溫半致死溫度的研究》的方法進行計算。
前言:很遺憾,r語言自帶的nls函式(nls(ddl ~ sslogis(temp, a, b, c))無法執行,主要緣故是資料問題。因此,根據查詢的文獻,從手動開始計算,到程式設計一次性完成。
文獻涉及的logistic方程如下:
y =k
/(1+
a∗e(
−b∗x
))y=k/(1+a*e^)
y=k/(1
+a∗e
(−b∗
x))問題改為涉及對引數k,a
和bk,a和b
k,a和
b的求解,半致死溫度的求解公式如下:
l t50
=lna
/blt_=ln/b
lt50=
lna/
b關鍵步:參考文獻需要對原始的電導率作如下轉換:
y ′=
(k−y
)/y=
a∗e(
−b∗x
)y'=(k-y)/y=a*e^
y′=(k−
y)/y
=a∗e
(−b∗
x)上述方程進行對數轉換:
l ny
′=ln
a−b∗
xln=ln-b*x
lny′=l
na−b
∗x進一步lny
′=ln
(k−y
)/yln=ln
lny′=l
n(k−
y)/y
式中,y
yy是原始的相對電導率(去掉%的值),k
=100
k=100
k=100。
通過l ny
′lnlny′
即可對x
xx做普通的線性回歸,獲得lna
和bln和b
lna和
b的值。
最後要求解k
kk,採用等距離的3個y
yy值求解,方法如下:
k =(
y22∗
(y1+
y3)−
2∗y1
∗y2∗
y3)/
(y22
−y1∗
y3)k=(y_^2*(y_+y_)-2*y_*y_*y_)/(y_^2-y_*y_)
k=(y22
∗(y
1+y
3)−
2∗y1
∗y2
∗y3
)/(
y22
−y1
∗y3
)通過上述的步驟,即可獲得logistic方程的各項引數值。之後就是利用r語言程式設計的事了,以及如何實現批量分析的目的。
對原始的電導率資料做所需的對數轉換和均值計算:
df2<
-readxl:
:read_excel
('抗逆性.xls'
,sheet=
1,skip=1)
df2$y<
-log((
100-df2$ddl)
/df2$ddl)
df2a<
-df2 %
>
%select
(clone,y,ddl,temp)
%>
%group_by
(clone,temp)
%>
%summarise
(dm=
mean
(ddl)
,ym=
mean
(y))
lgssf <
-function
(dat,mod=ym~temp,x=temp,y0=dm,tn=3:
5,ts=3)
res2a <
- plyr:
:dlply
(df2a,
"clone"
,lgssf)
>
names
(res2a) # 各植物材料的基因型**[1
]"1"
"2""3"
"4""5"
"6""7"
"8""9"
# 第乙個基因型的結果
> res2a$`1`[c
('equat'
,'lt50'
,'fit')]
$`equat`[1
]"y=21.713/(1+1.148*exp^(-0.127*x))"
$lt50
lt50
-1.089
$fit
$fit$`r2`[1
]0.8356973
$fit$fv
value numdf dendf
20.34531
1.00000
4.00000
$fit$pv
value
0.01073532
圖形檢視res2a$`1`$fig
:
到此,所有求解和批量運算都結束了。
水質分析 電導率
1 電導是電阻的倒數,指在外電場作用下電解質溶液中的正負離子以相反的方向的移動的能力。2 電解質中電導測量,用兩個電極插入溶液中,測定兩極間的電阻值的倒數。1 電導率是指距離1cm,截面積為1cm2的兩電極間溶液所測得的電導值。2 水中的離解物質少時,其電導率值與溶液中的離子含量成比例變化,所以,測...
電導率測量儀表1066 C AN 60
電導率測量儀表1066 c an 60 外 殼 材質為聚碳酸脂,防護等級為nema 4x ip66 尺 寸 外形尺寸155 155 131公釐 6.10 6.10 5.15英吋 儀錶盤開孔尺寸1 2 din標準,139 139公釐 5.45 5.45英吋 進 線 孔 6個,英吋或pg13.5穿線管接...
常數乘以無窮大等於多少 電導率常數 水質分析儀
電導率是在整個工業過程控制中廣泛使用的水質的重要度量。電導率測量的應用包括鍋爐水處理,冷卻塔水處理和反滲透監測。有幾種不同型別的感測器可用於測量電導率。選擇適合您的應用的感測器將提高準確性並確保裝置的使用壽命。選擇電導率感測器之前要了解的乙個重要因素是如何選擇電導率池常數。電導率感測器的工作原理 可...