Displaying 1 result from an estimated 1 matches for "meanbymsa".
2007 Feb 07
1
Problem with subsets and xyplot
...hivest)
cat <- CAT
cat[cat > 5] <- 6
msa <- as.numeric(MSA)
msa[msa == 7361] <- 7360
msa[msa == 7362] <- 7360
msa[msa == 7363] <- 7360
msa[msa == 5601] <- 5600
msa[msa == 5602] <- 5600
msa[msa == 6484] <- 6483
#### FIND MEANS FOR EACH MSA, FOR SUBSETTING LATER
meanbymsa <- aggregate(HIVEST, by = list(msa), FUN = mean, na.rm = T)
#### meanbymsa[,2] gives me the column I want; the 25%tile of this
column is about 3.1.
but when I try
plot1 <- xyplot(HIVEST~YEAR|as.factor(msa), pch = LETTERS[cat], subset
= (meanbymsa[,2] < 3.1))
plot1
I don't get wha...