28098

summary(pressure) ## podsumowanie zawarto?ci zbioru danych help(pressure) ## wyświetla informacji/ na temat zbioru dancyh pressure$temperature = pressure$temperature + 273.15 pressure$pressure = pressure$pressure * .1333 summa ry (p ressu re)

pres = p ressu re$p ressu re

temp = p ressu re$temperature

rm(pressure) ## usunil/cie zbl/dnych zmiennych

par(mfrow=c(1,4))    # jeden wiersz na 4 rysunki

plot(pres - temp, main="Ci?nienie pary \n dla rtłci", xlab="Temperatura (w stopniach Kelvina)", ylab="Ci?nienie(kPascal-e)")

plot(pres ~ temp, main="Ci?nienie pary \n dla rtl»ci", xlab="Temperatura (w stopniach Kelvina)", ylab="Ci?nienie(kPascal-e)", log="y")

plot(pres ~ temp, main="Ci?nienie pary \n dla rtłci", xlab="Temperatura (w stopniach Kelvina)\ ylab="Ci?nienie(kPascal-e)", log="xy")

plot (pres - temp, main="Ci?nienie pary \n dla rtlici", xlab="Temperatura (w stopniach Kelvina)", ylab="Ci?nienie(kPascal-e)", log="x")

#    jeden wiersz na dwa rysunki

#    model wykładniczy

#    rysunek reszt

par(mfrow=c(1,2))

Im.outl = lm(log(pres) - temp) Im.outl

plot(Im.outl$fitted. Im.outl$resid)

lm.out2 = lm(log(pres) - log(temp)) # transformacja log dla x i y lm.out2

plot(Im.out2$fitted, Im.out2$resid) # rysunek reszt

lm.out3 = Im (pres - temp + I(tenp^A2) + I(temp^A3)) summary(Im.out3)

par(mfrow=c(1,1))

plot(Im.out3$fitted, Im.out3$resid) plot(pres-temp)

curve(*475.8+3.804 *x -.009912*x^/v2+.00000844*x*3,add=T) plot(pres ~ I(l/temp)) libra ry ("MASS ") par(mfrow=c(1,2))

boxcox(pres - I(l/temp)) # podstawowy wykres box-cox boxcox(pres - I(l/temp), lambda=seq(-.2,.2,.01))

plot(log(pres) - I(l/temp))

lm.out4 = lm(log(pres) - I(l/temp)) summary(Im.out4)