data <- matrix(scan(),ncol=3,byrow=TRUE)
   1 250 16  1 325 18  1 475 32 
   2 250 19  2 475 46  2 550 45 
   3 325 26  3 400 39  3 550 61 
   4 400 21  4 475 35  4 550 55 
   5 325 19  5 475 47  5 550 48 
   6 250 20  6 400 33  6 475 31 
   7 250 13  7 325 13  7 400 34 
   8 250 21  8 400 30  8 550 52 
   9 250 24  9 325 10  9 550 50 
   10 325 24  10 400 31  10 475 37 

data

table92             <- data.frame(data)
colnames(table92)   <- c("run","pressure","conversion")
rm(data)

table92$run         <- as.factor(table92$run)  
table92$pressure    <- as.factor(table92$pressure) 


table92.lm          <- lm(conversion ~ run + pressure, data = table92)
anova(table92.lm)


 # check assumptions 

par(mfrow=c(2,2))  
plot(table92.lm)
par(mfrow=c(1,1))