R Notebook

Fuente: [https://docs.google.com/spreadsheets/d/1KxCgIUcBpiW3b8Pn-W3HqjEICs7e-jvNpF8y4vHSxxs/edit?usp=sharing]

Ver abajo la ED50 = 8.25 mg/kg

No se pudo encontrar diferencia entre la proporción de individuos mostrando retracción del ala de las cuatro distintas dósis. La dósis de 11mg/kg es la que mostró más individuos con mayor tiempo (hasta 20 min). Les hice un análisis binomial (regresión logística), tomando el umbral de 10 minutos, para ver si al menos la dósis 11mg/kg mostró diferencia con la dósis de 8mg/kg en cuanto a la proporción de individuos que lograron pasar los 10 minutos. Esta última es una forma más robusta de anaizarlo

Tienen un problema básico de diseño, ya que no tienen un grupo control (Sin Lidocaina). Al no tener un grupo sin Lidocaina, no pueden concluir si al menos la dósis máxima de Lidocaina tuvo algún efecto diferente. Existió un problema similar con el proyecto en ratas. El grupo control siempre es importante porque de otra forma no hay manera de saber si hay algun efecto de alguna dósis de lo que se está probando. No adjutaron los demás datos de las variables que en el protocolo dicen que analizarian (FC,FR, etc)

Reading Data

Plots

Call:
glm(formula = time.sens.min ~ dose, family = quasipoisson, data = lidoca)

Coefficients:
                Estimate Std. Error t value Pr(>|t|)    
(Intercept)    2.140e+00  1.782e-01  12.011 3.75e-14 ***
doseb.9mg.kg  -2.033e-10  2.520e-01   0.000   1.0000    
dosec.10mg.kg  1.112e-01  2.453e-01   0.454   0.6529    
dosed.11mg.kg  4.626e-01  2.275e-01   2.034   0.0494 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for quasipoisson family taken to be 2.698557)

    Null deviance: 114.573  on 39  degrees of freedom
Residual deviance:  98.546  on 36  degrees of freedom
AIC: NA

Number of Fisher Scoring iterations: 4

Passed 10 min

Call:
glm(formula = Y ~ dose, family = binomial, data = lidoca2)

Coefficients:
              Estimate Std. Error z value Pr(>|z|)  
(Intercept)    -2.1972     1.0541  -2.084   0.0371 *
doseb.9mg.kg    0.8109     1.3176   0.615   0.5383  
dosec.10mg.kg   1.3499     1.2599   1.071   0.2840  
dosed.11mg.kg   2.6027     1.2360   2.106   0.0352 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 6.6819e+00  on 3  degrees of freedom
Residual deviance: 6.6613e-16  on 0  degrees of freedom
AIC: 17.7

Number of Fisher Scoring iterations: 4

Pr. wing retraction

Proportion of wing retraction among groups

Regresión logística

Call:
glm(formula = wing ~ dose, family = binomial, data = wing2)

Coefficients:
                Estimate Std. Error z value Pr(>|z|)
(Intercept)   -4.701e-16  6.325e-01   0.000    1.000
doseb.9mg.kg   2.441e-16  8.944e-01   0.000    1.000
dosec.10mg.kg  4.055e-01  9.037e-01   0.449    0.654
dosed.11mg.kg  4.055e-01  9.037e-01   0.449    0.654

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 55.051  on 39  degrees of freedom
Residual deviance: 54.646  on 36  degrees of freedom
AIC: 62.646

Number of Fisher Scoring iterations: 4

ED50 (DE50)

(Intercept) 
       8.25 

No se pudo encontrar diferencia entre la proporción de individuos mostrando retracción del ala de las cuatro distintas dósis. La dósis de 11mg/kg es la que mostró más individuos con mayor tiempo (hasta 20 min). Les hice un análisis binomial (regresión logística), tomando el umbral de 10 minutos, para ver si al menos la dósis 11mg/kg mostró diferencia con la dósis de 8mg/kg en cuanto a la proporción de individuos que lograron pasar los 10 minutos. Esta última es una forma más robusta de anaizarlo

Tienen un problema básico de diseño, ya que no tienen un grupo control (Sin Lidocaina). Al no tener un grupo sin Lidocaina, no pueden concluir si al menos la dósis máxima de Lidocaina tuvo algún efecto diferente. Existió un problema similar con el proyecto en ratas. El grupo control siempre es importante porque de otra forma no hay manera de saber si hay algun efecto de alguna dósis de lo que se está probando.