forked from cecigongon/MatrizAgroecol-gica
-
Notifications
You must be signed in to change notification settings - Fork 0
/
analysis.R
28 lines (15 loc) · 969 Bytes
/
analysis.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
h0 <- c(5.5250,5.3500,5.2625,5.2250,4.9625,4.9875,4.7500,4.7875,4.6500,4.5500)
h1 <- c(5.2500,4.9750,4.7125,4.6000,4.3500,4.3500,4.0375,3.9875,3.5500,3.5250)
h2 <- c(4.9125,4.5625,4.2750,3.9375,3.6625,3.5000,3.2125,2.9875,2.7000,2.5375)
h3 <- c(4.8000,4.5500,4.0375,3.8000,3.4875,3.2125,2.7250,2.6625,2.5125,2.3375)
h4 <- c(4.7125,4.3375,3.9875,3.4875,3.1875,3.1500,2.6000,2.3500,2.1750,2.0375)
h0_30 <- c(6.7750,6.7750,6.7750,6.7500,6.7625,6.7250,6.7000,6.6875,6.6875,6.6500)
h1_30 <- c(6.7500,6.6625,6.5750,6.5375,6.5000,6.4625,6.4250,6.3375,6.2500,6.2250)
h2_30 <- c(6.6250,6.5125,6.4750,6.3500,6.2250,6.0500,6.0000,5.7625,5.6000,5.5375)
h3_30 <- c(6.6000,6.4750,6.4375,6.2125,6.0500,5.9250,5.6875,5.4000,5.3250,5.1125)
h4_30 <- c(6.6000,6.4750,6.4500,6.2250,6.0500,5.9250,5.7000,5.5125,5.3750,5.2625)
x <- c(0, .1, .2, .3, .4, .5, .6, .7, .8, .9)
fit1 <- lm(h0_30~poly(x, 1, raw=TRUE))
fit2 <- lm(h0_30~poly(x, 2, raw=TRUE))
summary(fit1)
summary(fit2)