R muscle 氯化鈣對大鼠心髒肌肉收縮的影響

說明

本實驗的目的是評估溶液中的鈣對大鼠心肌收縮的影響。分離21隻大鼠心髒的左心耳，多次電刺激constant-length組織條並將其浸入不同濃度的氯化鈣溶液中，然後準確測量條帶的縮短作為響應。

用法

muscle

格式

該 DataFrame 包含以下列：

Strip

使用哪種心肌條？

Conc

氯化鈣溶液的濃度，為 2.2 mM 的倍數。

Length

帶材長度(縮短)的變化(據稱)以毫米為單位。

例子

## IGNORE_RDIFF_BEGIN
A <- model.matrix(~ Strip - 1, data=muscle)
rats.nls1 <- nls(log(Length) ~ cbind(A, rho^Conc),
   data = muscle, start = c(rho=0.1), algorithm="plinear")
(B <- coef(rats.nls1))

st <- list(alpha = B[2:22], beta = B[23], rho = B[1])
(rats.nls2 <- nls(log(Length) ~ alpha[Strip] + beta*rho^Conc,
                  data = muscle, start = st))
## IGNORE_RDIFF_END

Muscle <- with(muscle, {
Muscle <- expand.grid(Conc = sort(unique(Conc)), Strip = levels(Strip))
Muscle$Yhat <- predict(rats.nls2, Muscle)
Muscle <- cbind(Muscle, logLength = rep(as.numeric(NA), 126))
ind <- match(paste(Strip, Conc),
            paste(Muscle$Strip, Muscle$Conc))
Muscle$logLength[ind] <- log(Length)
Muscle})

lattice::xyplot(Yhat ~ Conc | Strip, Muscle, as.table = TRUE,
   ylim = range(c(Muscle$Yhat, Muscle$logLength), na.rm = TRUE),
   subscripts = TRUE, xlab = "Calcium Chloride concentration (mM)",
   ylab = "log(Length in mm)", panel =
   function(x, y, subscripts, ...) {
      panel.xyplot(x, Muscle$logLength[subscripts], ...)
      llines(spline(x, y))
   })

來源

Linder, A.、Chakravarti, I. M. 和 Vuagnat, P. (1964) 用不同漸近線擬合漸近回歸曲線。在對統計的貢獻中。在 P. C. Mahalanobis 教授 70 歲生日之際贈送給他，編輯。 C. R. Rao，第 221-228 頁。牛津：佩加蒙出版社。

參考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth Edition. Springer.