Establishing relationship between discrete vector and categorical vector I have a vector of positive integer counts taken from a time-series (365 observations/days).
I believe that there is variance in these counts due to the day of week they are observed and as I am looking for abnormal observations I'd like to remove this variance before placing a threshold on what is abnormal.
I have tried creating a vector of days (1-7) which goes from t1 to t365 and counting the activity on those days. I then divided each of the count observations by the activity for that day.
Is this the right approach and how do I measure the strength of the correlation?
(I have looked at point-biserial correlation but is that just for binary categories?)
Ultimately I would like to see how much of the variance in the discrete vector is explained by the day of the week.
 A: How about this.  Let's make up some data where there really is a day effect
## dates of each day of 2012
days2012 <- seq(as.Date("2012-01-01"), as.Date("2012-12-31"), by="day")

## make variables representing interesting temporal units
dta <- data.frame(date=days2012, day=weekdays(days2012), month=months(days2012))

## specify average rates that are higher on Wednesdays
dta$muY <- ifelse(dta$day=='Wednesday', 10, 5)

## generate random count data from average rates  
dta$Y <- rpois(rep(1, nrow(dta)), dta$muY)

To spot the fact that Wednesdays are special in this data set you can use a generalised linear model.  This one asks whether rates vary by month or day or both:
mod <- glm(Y ~ day + month, data=dta, family=poisson)
summary(mod)

which gives a nice big effect for the dummy variable that represents Wednesday. 
To avoid the dummy variable contrasts obscuring the interpretation you can look at the marginal effects
library(effects) 
plot(effect("day", mod))

which again gives a big spike for Wednesday while operating on the scale of the observed counts.

A: Look into generalized linear models. You have a binomial or possibly Poisson dependent. What you will need to do is convert your categorical IV into indicator vars. Put another way, you will need to do an ANOVA using a non-normal DV. This is quite easy these days. In SAS, you use GENMOD.
