曲線の傾きと切片を取得するためにブラッドフォード キャリブレーション カーブを実行する関数を R で作成し、ラボ データを分析しました。ここに私が書いたコードがあります：
CalibrationLine <- function() {
BSAC <- readline("BSA concentrations (ug/uL): ")
BSAC <- str_split(BSAC, " ")
BSAC <- unlist(BSAC)
BSAC <- as.numeric(BSAC)
Absorbance <- readline("Absorbance (595 nm): ")
Absorbance <- str_split(Absorbance, " ")
Absorbance <- unlist(Absorbance)
Absorbance <- as.numeric(Absorbance)
Data <- data.frame(BSAC, Absorbance)
Plot <- ggplot(Data, aes(x = BSAC, y = Absorbance))+
geom_point()+
stat_smooth(method = lm)
Plot
paste("Correlation coefficient:", cor(Data$Absorbance, Data$BSAC), sep = " ")
model <- lm(Absorbance ~ BSAC, data = Data)
model
}
BSA 濃度とその吸光度の値を加算すると、線形回帰モデルの係数 (勾配と切片) が得られます。ただし、散布図はプロットされず、貼り付け機能は貼り付けられず、変数はグローバル環境に保存されません。これは私のコンソールに表示されるものです:
> CalibrationLine()
BSA concentrations (ug/uL): 0.125 0.5 1 2 4 5 7.5 10 14 18 22
Absorbance (595 nm): 0.006 0.032 0.065 0.203 0.33 0.364 0.55 0.661 0.846 1.022 1.123
Call:
lm(formula = Absorbance ~ BSAC, data = Data)
Coefficients:
(Intercept) BSAC
0.0719484 0.0524287
