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.