Farmacología: Universidad Nacional Autónoma de México Facultad de Química Departamento de Farmacia
Farmacología: Universidad Nacional Autónoma de México Facultad de Química Departamento de Farmacia
FACULTAD DE QUÍMICA
DEPARTAMENTO DE FARMACIA
FARMACOLOGÍA 1
DR. JOSÉ LUIS BALDERAS LÓPEZ
CÁLCULOS DE
FARMACOMETRÍA 1
UNIDAD 2.
Repetición 1 Repetición 2
C (mg/mL) ni ri ni ri
0 10 1 10 2
1 10 1 10 2
2 10 3 10 2
4 10 7 10 4
8 10 6 10 5
16 10 8 10 7
32 10 9 10 9
64 10 9 10 10
128 10 10 10 10
PRIMER PASO
1. IMPORTANTE: no sumar ni promediar los datos porque son categóricos nominales.
2. Calcule pi para cada concentración/dosis de cada repetición:
Repetición 1 Repetición 2
C (mg/mL)
ni ri pi ni ri pi
0 10 1 0.1 10 2 0.2
1 10 1 0.1 10 2 0.2
2 10 3 0.3 10 2 0.2
4 10 7 0.7 10 4 0.4
8 10 6 0.6 10 5 0.5
16 10 8 0.8 10 7 0.7
32 10 9 0.9 10 9 0.9
64 10 9 0.9 10 10 1
128 10 10 1 10 10 1
PRIMER PASO
3. En caso de que hubiera efecto en el control de una repetición, se corrigen todos los
datos de dicha repetición con la formula de Abbott.
PRIMER PASO
4. Calcular las probit (yi ) para cada proporción (pi ) con tablas o con formula.
Repetición 1 Repetición 2
C (mg/mL)
ni ri pi pi correg yi ni ri pi p correg yi
0 10 1 0.1 0 10 2 0.2 0
1 10 1 0.1 0 10 2 0.2 0
2 10 3 0.3 0.222 4.235 10 2 0.2 0
4 10 7 0.7 0.667 5.431 10 4 0.4 0.25 4.326
8 10 6 0.6 0.556 5.140 10 5 0.5 0.375 4.681
16 10 8 0.8 0.778 5.765 10 7 0.7 0.625 5.319
32 10 9 0.9 0.889 6.221 10 9 0.9 0.875 6.150
64 10 9 0.9 0.889 6.221 10 10 1 1
128 10 10 1 1 10 10 1 1
PRIMER PASO
Tabla de transformación de proporción (pi ) a unidades probit (yi )
y =DISTR.NORM.ESTAND.INV(proporción) + 5
PRIMER PASO
5. Se calcula el log10 de la dosis o concentración y se nombra a este valor como xi.
6. Se realiza la regresión lineal de xi vs yi utilizando todos los valores.
C (mg/mL) xi yi
2 0.301 4.235
4 0.602 5.431
4 0.602 4.326
8 0.903 5.140
8 0.903 4.681
16 1.204 5.765
16 1.204 5.319
32 1.505 6.221
32 1.505 6.150
64 1.806 6.221
PRIMER PASO
7. La ecuación de la recta de la regresión lineal es:
donde
despejando
PRIMER PASO
8. Se calcula la CL50 (0) mediante la siguiente formula:
( )
PRIMER PASO
9. Del problema:
.
( ) .
SEGUNDO PASO
1. A partir de los datos de regresión, se calcula la unidad probit asociada al logaritmo
de la dosis ( ).
Conc xi yi yi calc
Control
1 0.000 3.868
2 0.301 4.235 4.291
4 0.602 5.431 4.714
8 0.903 5.140 5.137
16 1.204 5.765 5.560
32 1.505 6.221 5.983
64 1.806 6.221 6.406
128 2.107 6.829
Control
1 0.000 3.868
2 0.301 4.291
4 0.602 4.326 4.714
8 0.903 4.681 5.137
16 1.204 5.319 5.560
32 1.505 6.150 5.983
DR. JOSÉ LUIS BALDERAS
64 LÓPEZ
1.806 6.406
Departamento de Farmacia 2.107
128 6.829
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
FACULTAD DE QUÍMICA
SEGUNDO PASO
2. Se transforma yi a su valor asociado en la distribución normal estándar (zi ) con la
siguiente formula:
-5
Conc xi yi yi calc zi
Control
1 0.000 3.868 -1.132
2 0.301 4.235 4.291 -0.709
4 0.602 5.431 4.714 -0.286
8 0.903 5.140 5.137 0.137
16 1.204 5.765 5.560 0.560
32 1.505 6.221 5.983 0.983
64 1.806 6.221 6.406 1.406
128 2.107 6.829 1.829
Control
1 0.000 3.868 -1.132
2 0.301 4.291 -0.709
4 0.602 4.326 4.714 -0.286
8 0.903 4.681 5.137 0.137
16 1.204 5.319 5.560 0.560
32 1.505 6.150 5.983 0.983
DR. JOSÉ LUIS BALDERAS
64 LÓPEZ
1.806 6.406 1.406
Departamento de Farmacia 2.107
128 6.829 1.829
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
FACULTAD DE QUÍMICA
SEGUNDO PASO
3. Se calcula la ordenada correspondiente al valor zi en la funcion normal estándar de
probabilidades con la formula:
SEGUNDO PASO
4. Se calcula la probabilidad acumulada (Pi ) de la distribución normal estándar, con
base al valor de zi :
SEGUNDO PASO
SEGUNDO PASO
4. Se calcula la probabilidad acumulada (Pi ) de la distribución normal estándar, con
base al valor de zi :
Pi =DISTR.NORM.ESTAND.N(zi , verdadero)
SEGUNDO PASO
5. Se calculan los factores de peso (wi ) con la formula:
SEGUNDO PASO
6. Se calculan los probits de trabajo (ni ) con la formula:
SEGUNDO PASO
7. Se calculan los valores: niwi, niwixi, niwixi2, niwini, niwinixi, niwi(xi-x)2 y sus sumatorias.
donde
niwi niwixi niwixi2 niwini niwinixi niwi(xi-)2
SEGUNDO PASO
8. Se calculan b1 y a1 con las siguientes formulas:
SEGUNDO PASO
Por lo tanto:
.
( ) .
ITERACIONES
Para la segunda iteración, se repite nuevamente el segundo paso del método,
calculando la unidad probit asociada al logaritmo de la dosis ( ) con:
ITERACIONES
Se continua con las iteraciones hasta:
( )
LÍMITES DE CONFIANZA
Se calcula
LIMITES DE CONFIANZA
Con los datos del problema, a la primera iteración:
a1 = 3.426
b1 = 1.798
CL50 (1) = 7.510
Sm2 = 0.00425
Sm = 0.0652 t(a=0.05,178) = 1.9734 t =INV.T.2C(probabilidad , grados de libertad)
minf = 0.728
msup = 1.023
LSC(1) = 10.546
LIC (1) = 5.349
DR. JOSÉ LUIS BALDERAS LÓPEZ
Departamento de Farmacia
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
FACULTAD DE QUÍMICA
Repetición 1 Repetición 2
C (mg/mL) ni ri ni ri
0 10 1 10 2
1 10 1 10 2
2 10 3 10 2
4 10 7 10 4
8 10 6 10 5
16 10 8 10 7
32 10 9 10 9
64 10 9 10 10
128 10 10 10 10
Repetición 1 Repetición 2
C (mg/mL)
ni ri pi ni ri pi
0 10 1 0.1 10 2 0.2
1 10 1 0.1 10 2 0.2
2 10 3 0.3 10 2 0.2
4 10 7 0.7 10 4 0.4
8 10 6 0.6 10 5 0.5
16 10 8 0.8 10 7 0.7
32 10 9 0.9 10 9 0.9
64 10 9 0.9 10 10 1
128 10 10 1 10 10 1
C (mg/mL) xi yi
1 0 3.355
2 0.301 4.235
2 0.301 3.355
4 0.602 5.431
4 0.602 4.326
8 0.903 5.140
8 0.903 4.681
16 1.204 5.765
16 1.204 5.319
32 1.505 6.221
32 1.505 6.150
64 1.806 6.221
64 1.806 6.645
128 2.107 6.645
donde y
despejando
.
.
donde:
. .
donde:
= El número total de unidades experimentales entre
pi = 0.07 y pi = 0.93
(∝ . , )
(∝ . , )
donde:
= número total de unidades experimentales utilizadas en el experimento.
(a=0.05, gl=178)
m
m
DR. JOSÉ LUIS BALDERAS LÓPEZ
Departamento de Farmacia
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
FACULTAD DE QUÍMICA
50 m
95 m
Repetición 1 Repetición 2
C (mg/mL) ni ri ni ri
0 10 1 10 2
1 10 1 10 2
2 10 3 10 2
4 10 7 10 4
8 10 6 10 5
16 10 8 10 7
32 10 9 10 9
64 10 9 10 10
128 10 10 10 10
Repetición 1 Repetición 2
C (mg/mL)
ni ri pi ni ri pi
0 10 1 0.1 10 2 0.2
1 10 1 0.1 10 2 0.2
2 10 3 0.3 10 2 0.2
4 10 7 0.7 10 4 0.4
8 10 6 0.6 10 5 0.5
16 10 8 0.8 10 7 0.7
32 10 9 0.9 10 9 0.9
64 10 9 0.9 10 10 1
128 10 10 1 10 10 1
log C pi
0.000 0
0.301 0.222
0.602 0.444
0.903 0.444
1.204 0.778
1.505 0.889
1.806 0.889
2.107 1
0.000 0
0.301 0
0.602 0.25
0.903 0.375
1.204 0.625
1.505 0.875
1.806 1
2.107 1
log C pi pi calc
0.000 0 0.005
0.301 0.222 0.161
0.602 0.444 0.316
0.903 0.444 0.472
1.204 0.778 0.627
1.505 0.889 0.783
1.806 0.889 0.938
2.107 1 1.094
0.000 0 0.005
0.301 0 0.161
0.602 0.25 0.316
0.903 0.375 0.472
1.204 0.625 0.627
1.505 0.875 0.783
1.806 1 0.938
2.107 1 1.094
0.1 0.032 0.035 0.038 0.041 0.044 0.047 0.049 0.052 0.055 0.057
0.2 0.060 0.062 0.065 0.067 0.070 0.072 0.074 0.076 0.078 0.081
0.3 0.083 0.084 0.086 0.088 0.090 0.092 0.093 0.094 0.096 0.098
0.4 0.099 0.100 0.101 0.102 0.103 0.104 0.104 0.104 0.104 0.105
0.5 - 0.895 0.896 0.896 0.896 0.897 0.897 0.898 0.899 0.900
0.6 0.901 0.902 0.904 0.905 0.907 0.908 0.910 0.912 0.914 0.916
0.7 0.917 0.919 0.922 0.924 0.926 0.928 0.930 0.933 0.935 0.938
0.8 0.940 0.943 0.945 0.948 0.951 0.953 0.956 0.959 0.962 0.965
0.9 0.968 0.971 0.974 0.977 0.980 0.984 0.987 0.990 0.993 0.997
84 m
16 m
CL50
m
m
MÉTODO DE SPEARMAN-KÄRBER
Es un método basado en la distribución de la tolerancia, donde la CL50 es el la media
de la distribución empírica de la tolerancia que, se asume simétrica.
Requisitos:
a) Debe haber concentraciones con p = 0 y p = 1, o asumir que debajo de xi será p = 0
y por arriba de xk será p = 1.
b) Si xi < xi+1 < … < xk, entonces pi < pi+1 < … < pk (“monotonized proportions”).
c) Se sugiere que los rangos logarítmicos de las concentraciones sean iguales a lo
largo de todo el experimento, aunque no es estrictamenmte necesario.
MÉTODO DE SPEARMAN-KÄRBER
a) Cuando los rangos de los logaritmos de las concentraciones son equivalentes
Donde:
xk = log[Conc] mínima a la cual pi = 1
pi = proporciones.
d = rango de los log[Conc] .
R = número de repeticiones
MÉTODO DE SPEARMAN-KÄRBER
b) Cuando los rangos de los logaritmos de las concentraciones NO son equivalentes
Donde:
xk = log[Conc] mínima a la cual pi = 1
xi = log[Conc].
pi = proporciones.
R = número de repeticiones
MÉTODO DE SPEARMAN-KÄRBER
C) Cálculo de la varianza y límites de confianza
Donde:
xk = log[Conc] mínima a la cual pi = 1
xi = log[Conc]. Algunos autores
ni = número de unidades experimentales en el grupo. sugieren cambiar el
pi = proporciones. valor de 4ni por ni-1
R = número de repeticiones
MÉTODO DE SPEARMAN-KÄRBER
C) Cálculo de la varianza y límites de confianza
( (∝/ ) ( ))
( (∝/ ) ( ))
Donde:
xk = log[Conc] mínima a la cual pi = 1
xi = log[Conc].
ni = número de unidades experimentales en el grupo.
pi = proporciones.
R = número de repeticiones
MÉTODO DE SPEARMAN-KÄRBER
Determine la CL50 y los límites de confianza del dicromato de
potasio en Artemia salina de acuerdo con los siguientes datos:
Repetición 1 Repetición 2
C (mg/mL) ni ri ni ri
0 10 1 10 2
1 10 1 10 2
2 10 3 10 2
4 10 7 10 4
8 10 6 10 5
16 10 8 10 7
32 10 9 10 9
64 10 9 10 10
128 10 10 10 10
MÉTODO DE SPEARMAN-KÄRBER
1. IMPORTANTE: no sumar ni promediar los datos porque son categóricos nominales.
2. Calcule pi para cada concentración/dosis de cada repetición:
Repetición 1 Repetición 2
C (mg/mL)
ni ri pi ni ri pi
0 10 1 0.1 10 2 0.2
1 10 1 0.1 10 2 0.2
2 10 3 0.3 10 2 0.2
4 10 7 0.7 10 4 0.4
8 10 6 0.6 10 5 0.5
16 10 8 0.8 10 7 0.7
32 10 9 0.9 10 9 0.9
64 10 9 0.9 10 10 1
128 10 10 1 10 10 1
MÉTODO DE SPEARMAN-KÄRBER
3. En caso de que hubiera efecto en el control de una repetición, se corrigen todos los
datos de dicha repetición con la formula de Abbott.
MÉTODO DE SPEARMAN-KÄRBER
4. Se aplica la corrección a las pi que no siguen un orden creciente, promediando las
pi y colocando el valor en ambas dosis. Además, se ajusta a una sola pi=0 y pi=1.
MÉTODO DE SPEARMAN-KÄRBER
5. Se calculan los valores (pi+1 - pi), (xi + xi+1) y el producto de ambos.
Repetición 1 Repetición 2
C (mg/mL) log C
p Correg SK pi+1 - pi xi + xi+1 producto p Correg SK pi+1 - pi xi + xi+1 producto
0
1 0 0 0.2222 0.3010 0.0669
2 0.3010 0.2222 0.4278 0.9031 0.3863 0 0.2500 0.9031 0.2258
4 0.6021 0.6500 0.0000 1.5051 0.0000 0.2500 0.1250 1.5051 0.1881
8 0.9031 0.6500 0.1278 2.1072 0.2693 0.3750 0.2500 2.1072 0.5268
16 1.2041 0.7778 0.1111 2.7093 0.3010 0.6250 0.2500 2.7093 0.6773
32 1.5051 0.8889 0.0000 3.3113 0.0000 0.8750 0.1250 3.3113 0.4139
64 1.8062 0.8889 0.1111 3.9134 0.4348 1
128 2.1072 1
MÉTODO DE SPEARMAN-KÄRBER
6. Se calcula el valor de CL50
MÉTODO DE SPEARMAN-KÄRBER
7. Para el cálculo de los límites de confianza, se calculan los valores de pi(1-pi)/4ni,
(xi+1 – xi)2 y el producto de ambos.
𝑝 (1 − 𝑝 )(𝑥 −𝑥 )
𝑆𝐸(log 𝐶𝐿 ) =
4𝑛
Repetición 1 Repetición 2
C (mg/mL) log C
ni p Correg SK pi(1-pi)/4ni (xi+1 – xi)2 producto ni p Correg SK pi(1-pi)/4ni (xi+1 – xi)2 producto
0
1 0 10 0.0000 0.0000 0.0906 0.0000
2 0.3010 10 0.2222 0.0043 0.0906 0.0004 10 0.0000 0.0000 0.0906 0.0000
4 0.6021 10 0.6500 0.0057 0.0906 0.0005 10 0.2500 0.0047 0.0906 0.0004
8 0.9031 10 0.6500 0.0057 0.0906 0.0005 10 0.3750 0.0059 0.0906 0.0005
16 1.2041 10 0.7778 0.0043 0.0906 0.0004 10 0.6250 0.0059 0.0906 0.0005
32 1.5051 10 0.8889 0.0025 0.0906 0.0002 10 0.8750 0.0027 0.0906 0.0002
64 1.8062 10 0.8889 0.0025 0.0906 0.0002 10 1
128 2.1072 10 1
MÉTODO DE SPEARMAN-KÄRBER
8. Se calcula el SE(log CL50) y los límites de confianza
( / ) .
( . ( . )( . ))
( . ( . )( . ))
RESUMEN
Valores calculados (mg/mL)
Método/Software
CL50 LSC LIC
Máxima verosimilitud 7.510 10.546 5.349
Miller-Tainter 7.963 12.174 3.752
Lichfield-Wilcoxon 9.237 13.299 6.415
Spearman-Kärber 7.457 9.930 5.600