arithscan                package:qtl                R Documentation

_A_r_i_t_h_m_e_t_i_c _O_p_e_r_a_t_o_r_s _f_o_r _g_e_n_o_m_e _s_c_a_n _r_e_s_u_l_t_s

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

     These binary operators perform arithmetic on genome scan result
     objects.

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

     scan1+scan2
     scan1-scan2

_A_r_g_u_m_e_n_t_s:

scan1, scan2: Genome scan results on the same set of chromosomes and
          markers.

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

     This is used to calculate the sum or difference of LOD scores of
     two genome scan results. The basic arithmetic operators were
     overrided to work on objects of class 'scanone', 'scantwo', and
     'scanqtl'.

     It currently only works for 'scanone' results.

_V_a_l_u_e:

     The same data structure as the input objects, with LOD scores
     recalcuated.

_A_u_t_h_o_r(_s):

     Karl W Broman, kbroman@jhsph.edu

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

     data(fake.bc)
     fake.bc <- calc.genoprob(fake.bc)
     # genome scans without covariates
     out.nocovar <- scanone(fake.bc)
     # genome scans with covariates
     ac <- fake.bc$pheno[,c("sex","age")]
     ic <- fake.bc$pheno[,"sex"]
     out.covar <- scanone(fake.bc, pheno.col=1,
                          addcovar=ac, intcovar=ic)
     # plot the difference of with and without additive covariate
     plot(out.covar-out.nocovar)

