Abstract: The aim of the present invention is to estimate EDR directions in a single index model composed of a large number of variables with simple calculations without using the inverse matrix of the variance-covariance matrix and principle component analysis. Data conversion means 21 receives, from an input device 3, data to be analyzed, the data composed of sets of response variables and explanatory variables, standardizes the explanatory variables, and sends them to slice average calculating means 22. The slice average calculating means 22 divides the data into two slices with reference to the median of the response variables to calculate the mean vector of the explanatory variables on a slice basis. The calculated mean vectors are sent to EDR direction calculating means 23. The EDR direction calculating means 23 calculates the difference between the mean vectors for respective slices to estimate an EDR direction.

Abstract: More precise correction of global and local distortions of microarray data and correction of measurement errors caused by a difference in sensitivity between fluorescent dyes. A data standardization unit for a first process inputs gene expression intensity data from an input device, standardizes the gene expression intensity data by using grid-by-grid order statistics on the assumption that most genes are in a non-expression state, and outputs the standardized gene expression intensity data. A spot-position-based correction unit for a second process estimates a distortion depending on a spot position on a grid by grid basis by a nonparametric smoothing method and outputs gene expression intensity data whose distortion depending on the spot position has been corrected.