220 likes | 327 Views
Transforming the data. Modified from: Gotelli and Allison 2004. Chapter 8; Sokal and Rohlf 2000 Chapter 13. What is a transformation?. It is a mathematical function that is applied to all the observations of a given variable
E N D
Transforming the data Modified from: Gotelli and Allison 2004. Chapter 8; Sokal and Rohlf 2000 Chapter 13
What is a transformation? It is a mathematical function that is applied to all the observations of a given variable • Y represents the original variable, Y* is the transformed variable, and f is a mathematical function that is applied to the data
Most are monotonic: • Monotonic functions do not change the rank order of the data, but they do change their relative spacing, and therefore affect the variance and shape of the probability distribution
There are two legitimate reasons to transform your data before analysis • The patterns in the transformed data may be easier to understand and communicate than patterns in the raw data. • They may be necessary so that the analysis is valid
They are often useful for converting curves into straight lines: The logarithmic function is very useful when two variables are related to each other by multiplicative or exponential functions
Statistics and transformation Data to be analyzed using analysis of variance must meet to assumptions: • The data must be homoscedastic: variances of treatment groups need to be approximately equal • The residuals, or deviations from the mean must be normal random variables
Lets look an example • A single variate of the simplest type of ANOVA (completely randomized, single classification) decomposes as follows: • In this model the components are additive with the error term εij distributed normally
However… • We might encounter a situation in which the components are multiplicative in effect, where • If we fitted a standard ANOVA model, the observed deviations from the group means would lack normality and homoscedasticity
The logarithmic transformation • We can correct this situation by transforming our model into logarithms Wherever the mean is positively correlated with the variance the logarithmic transformation is likely to remedy the situation and make the variance independent of the mean
We would obtain • Which is additive and homoscedastic
The square root transformation • It is used most frequently with count data. Such distributions are likely to be Poisson distributed rather than normally distributed. In the Poisson distribution the variance is the same as the mean. Transforming the variates to square roots generally makes the variances independents of the means for these type of data. When counts include zero values, it is desirable to code all variates by adding 0.5.
The box-cox transformation • Often one do not have a-priori reason for selecting a specific transformation. • Box and Cox (1964) developed a procedure for estimating the best transformation to normality within the family of power transformation
The box-cox transformation • The value of lambda which maximizes the log-likelihood function: yields the best transformation to normality within the family of transformations s2Tis the variance of the transformed values (based on v degrees of freedom). The second term involves the sum of the ln of untransformed values
box-cox in R (for a vector of data Y) >library(MASS) >lamb <- seq(0,2.5,0.5) >boxcox(Y_~1,lamb,plotit=T) >library(car) >transform_Y<-box.cox(Y,lamb) What do you conclude from this plot? Read more in Sokal and Rohlf 2000 page 417
The arcsine transformation • Also known as the angular transformation • It is especially appropriate to percentages
The arcsine transformation Transformed data It is appropriate only for data expressed as proportions Proportion original data
Since the transformations discussed are NON-LINEAR, confidence limits computed in the transformed scale and changed back to the original scale would be asymmetrical