大家能不能帮忙画个图，实在不会了。

carnival

画图的代码不会写了。。。求大家帮忙画个图

setwd("c:/work")
data <- read.table ("Hollywood.txt", header = TRUE)
attach(data)
y<-data$X1
x1<-data$X2
x2<-data$X3
x3<-data$X4
fit1 <- lm (y~x1)
fit2 <- lm (y~x2)
fit3 <- lm (y~x3)
fit4 <- lm (y ~ x1 + x2 +x3)
fit5 <- lm (y ~ x1 + x2 )
summary (fit1)
summary (fit2)
summary (fit3)
summary (fit4)
summary (fit5)

Hollywood.txt里的内容是
X1 X2 X3 X4
85.1  8.5  5.1  4.7
106.3  12.9  5.8  8.8
50.2  5.2  2.1  15.1
130.6  10.7  8.4  12.2
54.8  3.1  2.9  10.6
30.3  3.5  1.2  3.5
79.4  9.2  3.7  9.7
91.0  9.0  7.6  5.9
135.4  15.1  7.7  20.8
89.3  10.2  4.5  7.9


多谢大家了！！

carnival

找了一个例子，但是不知道数据怎么换。。 大家帮帮忙吧 谢谢
x11(8,4)
op <- par (bg = "white")
par ( mar = c(5,4,4,2) - c (1,0,1,1) + 0.6 , mfrow = c (1,2) , cex = 0.8 ,
mex = 0.8 , pch = 16 , lab = c (5,5,7))

par (cex = 0.8)
plot (y~x1)
abline (crime4, col = "red")

小小七nx

你这是要分别画出每一张图吗？

carnival

小小七nx 发表于 2016-4-27 18:52
你这是要分别画出每一张图吗？

新手不太懂呢  算出来y ~ x1 + x2是正确的。能帮我画这个么

Tracy-hu

不知道你的意思是什么？

zhigaowu

setwd("c:/work")
data <- read.table ("Hollywood.txt", header = TRUE)
data
attach(data)
y<-data$X1
x1<-data$X2
x2<-data$X3
x3<-data$X4
fit1 <- lm (y~x1)
fit2 <- lm (y~x2)
fit3 <- lm (y~x3)
fit4 <- lm (y~x1 + x2 + x3)
fit5 <- lm (y~x1 + x2 )
summary (fit1)
summary (fit2)
summary (fit3)
summary (fit4)
summary (fit5)

#####################################################
#refrence: http://www.walkingrandomly.com/?p=5254
#只绘了前三项图
par(bg=colors()[1])
plot(x1, y)
new = data.frame(x1 = seq(min(x1),max(x1),len=2))
lines(new$x1,predict(fit1,newdata=new))
Sys.sleep(5)

par(bg=colors()[399])
plot(x2, y)
new = data.frame(x2 = seq(min(x2),max(x2),len=2))
lines(new$x2,predict(fit2,newdata=new))
Sys.sleep(3)

par(bg=colors()[652])
plot(x3, y)
new = data.frame(x3 = seq(min(x3),max(x3),len=2))
lines(new$x3,predict(fit3,newdata=new))
Sys.sleep(3)

hou0922

楼主问题表达不明确啊，画什么图都没说清楚

xydf_1992

你是要把回归结果做散点图么？把下面代码根据自己情况改改应该就行

plot(NDBI$V3, LST$V3, #X轴为NDBI，Y轴为LST数据，V3是因为前两行为空
     xlab = "NDBI",  #X轴显示NDBI
     ylab = "LST(摄氏度)", #Y轴显示LST摄氏度
     col="black",pch=1) #指定点颜色和样式，可以改的
abline(huigui) #绘制斜线

xydf_1992

你是要把回归结果做散点图么？把下面代码根据自己情况改改应该就行

plot(NDBI$V3, LST$V3, #X轴为NDBI，Y轴为LST数据，V3是因为前两行为空
     xlab = "NDBI",  #X轴显示NDBI
     ylab = "LST(摄氏度)", #Y轴显示LST摄氏度
     col="black",pch=1) #指定点颜色和样式，可以改的
abline(huigui) #绘制斜线