liarozole               package:MCPAN               R Documentation

_M_a_r_k_e_d _i_m_p_r_o_v_e_m_e_n_t _o_f _p_s_o_r_i_a_s_i_s _a_f_t_e_r _a_p_p_l_i_c_a_t_i_o_n _o_f _l_i_a_r_o_z_o_l_e

_D_e_s_c_r_i_p_t_i_o_n:

     In a placebo controlled clinical trial, patients with psoriasis
     were randomly assigned to a placebo group and three dose groups
     (50 mg, 75 mg, and 150 mg). Variable of primary interest was the
     proportion of patients with marked improvement of psoriasis. This
     data.frame mimics how raw data could have been represented in a
     larger data frame.

_U_s_a_g_e:

     data(liarozole)

_F_o_r_m_a_t:

     A data frame with 137 observations on the following 2 variables.

     '_I_m_p_r_o_v_e_d' a factor with levels 'n',  'y', for "no" and "yes"

     '_T_r_e_a_t_m_e_n_t' a factor with levels 'Dose150',  'Dose50',  'Dose75', 
          'Placebo'

_D_e_t_a_i_l_s:

     For illustrative purpose only. Number of successes recalculated
     from proportions presented in the publication, while the number of
     patients in group 'Dose50' was not exactly clear.

_S_o_u_r_c_e:

     Berth-Jones J, Todd G, Hutchinson PE, Thestrup-Pedersen K,
     Vanhoutte FP: Treatment of psoriasis with oral liarozole: a
     dose-ranging study. British Journal of Dermatology 2000; 143:
     1170-1176.

_E_x_a_m_p_l_e_s:

     data(liarozole)
     head(liarozole)
     # create a contingency table:

     table(liarozole)

     # the order of the groups is alpha-numeric,
     # and "y" for success is of higher order than
     # to change the order:

     liarozole$Treatment<-factor(liarozole$Treatment,
      levels=c("Placebo", "Dose50", "Dose75", "Dose150"))

     liarozole$Improved<-factor(liarozole$Improved,
      levels=c("y", "n"))

     tab<-table(liarozole)
     tab

     # Approximate simultaneous confidence intervals
     # for the differences  pDose-pPlacebo:

     LCI<-binomRDci(tab, type="Dunnett", 
     alternative="greater", method="ADD1")

     LCI

     plot(LCI, main="Proportion of patients
      with marked improvement")

     # Perform a test on increasing trend 
     # vs. the placebo group:

     Ltest<-binomRDtest(tab, type="Williams", 
     alternative="greater", method="ADD1")

     summary(Ltest)

