-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathperf_reg.r
72 lines (68 loc) · 2.32 KB
/
perf_reg.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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
observed <- c(0.22, 0.83, -0.12, 0.89, -0.23, -1.30, -0.15, -1.4,
0.62, 0.99, -0.18, 0.32, 0.34, -0.30, 0.04, -0.87,
0.55, -1.30, -1.15, 0.20)
predicted <- c(0.24, 0.78, -0.66, 0.53, 0.70, -0.75, -0.41, -0.43,
0.49, 0.79, -1.19, 0.06, 0.75, -0.07, 0.43, -0.42,
-0.25, -0.64, -1.26, -0.07)
residualValues <- observed - predicted
# summary(residualValues)
# observed values vs pred values
axisRange <- extendrange(c(observed, predicted))
plot(observed, predicted, ylim = axisRange, xlim = axisRange)
abline(0, 1, col="darkgrey", lty=2)
# residual vs predicted
plot(predicted, observed, ylab = "residual")
abline(h = 0, col="darkgrey", lty=2)
# Calculate R square and RMSE with caret package
library("caret")
R2(predicted, observed)
RMSE(predicted, observed)
# correlation
cor(predicted, observed)
cor(predicted, observed, method = "spearman")
library(lubridate)
library(reshape2)
body <- read.table("http://www.amstat.org/publications/jse/datasets/body.dat.txt")
boxplot(body)
keep.par <- par() # stored the original settings >par(keep.par)
par(mar=c(5, 4, 4, 2))
boxplot(body, las=3)
x <- 1:10
set.seed(23)
y <- x + rnorm(10)
plot(x, y)
plot(x, y, col=x)
plot(x, y, col=x, cex=x)
breaks <- seq(min(body$V22), max(body$V22), 5)
v22_group <- cut((body$V22), breaks)
body$group <- v22_group
plot(body$V11, body$V20, pch=body$V25, col=body$group)
＃the normal distribution using R
x <- seq(-4, 4, 0.01)
plot(x, dnorm(x), type="l")
plot(x, dnorm(x))
plot(x, pnorm(x))
qnorm(c(0.025, 0.975))
#one tailed test using Normal Distr
bt <- seq(60, 120, 1)
plot(bt, dnorm(bt, 90, 10), type="l")
pnorm(72, 90, 10)
abline(v=72)
cord.x <- c(60, seq(60, 72, 1), 72)
cord.y <- c(0, dnorm(seq(60, 72, 1), 90, 10), 0)
polygon(cord.x, cord.y, col="skyblue")
text(70, 0.005, "blue area = p = 0.0359")
#two-tailed test
bt <- seq(60, 120, 1)
plot(bt, dnorm(bt, 90, 10), type="l", xlim=c(60, 120), main="two-tailed test")
pnorm(72, 90, 10)
abline(v=72)
cord.x <- c(60,seq(60,72,1),72)
cord.y <- c(0,dnorm(seq(60, 72, 1), 90, 10),0)
polygon(cord.x,cord.y,col='skyblue')
cord.x1 <- c(108,seq(108,120,1),120)
cord.y1 <- c(0,dnorm(seq(108, 120, 1), 90, 10),0)
polygon(cord.x1,cord.y1,col='skyblue')
text(65, 0.005, round(pnorm(72, 90, 10), 3))
text(115, 0.005, round(pnorm(72, 90, 10), 3))
text(75, 0.02, "p = 0.072 " )