BGFilterWeighted          package:beadarray          R Documentation

_W_e_i_g_h_t_e_d _B_a_c_k_g_r_o_u_n_d _F_i_l_t_e_r

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

     Finds local weighted means at each bead.

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

     BGFilterWeighted(E = NULL, neighbours, invasions = 20, weights = NULL)

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

       E: Error Image

neighbours: A Neighbours matrix. Required.

invasions: Integer - Number of invasions. This argument is passed to
          the function BGfilter.

 weights: Numerical vector - A vector of weights, from 0 to 1, to
          consider in the analysis. (see below.)

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

     This function finds the weighted mean of local bead intensities,
     using intensities from the given error image.

     To obtain our weighted mean for each bead, we use an invasion
     process. Links between beads are defined in the neighbours matrix.
     We define the local beads as those which can be reached in
     'invasions' steps from the first bead, and then collect their
     error values.

     We take a weighted mean of these error values, where the weights
     are calculated by taking the product of: a) 1/(the number of steps
     required to get to the bead from the central bead) b) (if
     supplied) the weights defined through the 'weights' parameter.

     This weighted mean is then assigned to the central bead.

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

     A vector - the weighted means. (NB: Whilst 'BGFilter' manipulates
     the error image and returns an updated error image, e.g.
     subtracting the local median, this function does not - it merely
     returns the local weighted means.)

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

     Jonathan Cairns

_S_e_e _A_l_s_o:

     'BGFilter'

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

     data(BLData)
     E <- generateE(BLData,1,method = "mean")
     neighbours <- generateNeighbours(BLData,1)
     ##bgf <- BGFilterWeighted(E, neighbours)

