Nothing
surface <- ice[ice$Location=='surface',c("Treatment","T1930")]
intra <- ice[ice$Location=='intramuscular',"T1930"]
newdata <- cbind(surface,intra)
names(newdata) <- c('Treatment','SurfTemp','IntraTemp')
anova(lm(SurfTemp - IntraTemp ~ Treatment, newdata))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.