How to create a (total) Ghrelin calibration curve?

Version 11
by (unknown)
Version 12
by (unknown)

Deletions or items before changed

Additions or items after changed

1 == Computing Ghrelin calibration curve using R ==
2
3 The calibration curve can be created using R, and libracy "drc". Optionally, one can use library "sfmisc" for formatting of the labels on plot axis.
4 The R code is attached below.
5
6 The example assumes the data to be available in file "ghrelin_conc std_a std_b avg.csv"
7
8 The measured data:
9
10 || '''Ghrelin (ng/ml)''' || '''Standard a''' || '''Standard b''' ||
11 -
|| 1000000||-0.040596823||-0.052699697 ||
+
|| 1000000 ||-0.040596823||-0.052699697 ||
12 -
| |100000||0.136105144||0.119766263 ||
+
|| 100000 ||0.136105144||0.119766263 ||
13 -
||10000||0.61356354||0.606906959 ||
+
|| 10000 ||0.61356354||0.606906959 ||
14 -
||1000||0.846543873||0.839887292 ||
+
|| 1000 ||0.846543873||0.839887292 ||
15 -
||100||0.887693646||0.88345764 ||
+
|| 100 ||0.887693646||0.88345764 ||
16 ||0||0.896770802||0.896165658 ||
17
18 {{{
19
20 ##### Install libraries
21 install.packages("drc")
22 install.packages("sfsmisc")
23 require(drc)
24 library(sfsmisc)
25
26 ##### Read the data
27 hormone.data <- read.csv("ghrelin_conc std_a std_b avg.csv")
28 hormone.data <- hormone.data[,1:3]
29 colnames(hormone.data)[1:3] <- c("Concentration","Response_1", "Response_2")
30
31 ##### Reorganize the data
32 hormone.data <- reshape(hormone.data, varying=c("Response_1","Response_2"), direction="long", v.names=c("Response"))
33 hormone.data <- hormone.data[,c("Concentration", "Response")]
34
35 ##### Fitting the model (4-parameter log-logistic function)
36 hormone.data.model <- drm(Response ~ Concentration, data = hormone.data, fct = LL.4())
37 summary(hormone.data.model)
38
39 }}}
40
41 The calibration curve can be plotted using the commands below:
42
43 {{{
44 ##### Plotting a nice plot
45 par(pty="s", mar=c(5,5,1,1))
46 plot(hormone.data.model, type="confidence", cex.lab=2, axes=F, xlim=c(-10,10^6))
47 axis(side=1, at=hormone.data[1:6,1], labels=pretty10exp(hormone.data[1:6,1]), cex.axis=1.2)
48 axis(side=2, at=seq(0,1,0.2), labels=seq(0,1,0.2))
49 plot(hormone.data.model, type="all", add=T, pch=21, col="red", lwd=1, cex=2, bg="green")
50 }}}
51
52 -
+
[[Image(Ghrelin.png)]]
53
54
55 ##### Computing the concentration from the response, for instance for a response=0.1, and alpha=1-0.95
56 ED(hormone.data.model, respLev=0.1, interval="delta", type="absolute", level=0.95)