# Accessing more than two coefficients in a plot

 Classic List Threaded
3 messages
Reply | Threaded
Open this post in threaded view
|

## Accessing more than two coefficients in a plot

 I've successfully plotted (in the plot and abline code below) a simple regression of Lambda1_2 on VV1_2. I then successfully regressed Lambda1_2 on VV1_2, VV1_22 and VV1_212 producing lm2.l. When I go to plot lm2.l using abline I get the warning: "1: In abline(lm2.l, col = "brown", lty = "dotted", lwd = 2) : only using the first two of 4 regression coefficients" Is there another function like abline that will produce a line using the constant and three coefficients from the lm2.l regression? lm.l <- lm(Lambda1_2 ~ VV1_2, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression lm2.l <- lm(Lambda1_2 ~ VV1_2 + VV1_22 + VV1_212, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression plot(VV1_2, Lambda1_2, ylim=yrange, tck=1, main="V(1) Parameters (V, V^2 & V^0.5)", xlab="VV1_2", ylab="Lambda & Beta1_2",pch=19,col="red") {abline(lm2.l, col="brown", lty="dotted", lwd=2) abline(wlm2.l, col="gold",lty="longdash", lwd=2) points(VV1_2, Beta1_2, pch=19, col="blue") abline(lm2.b, col="black",lty="dotted", lwd=2) abline(wlm2.b, col="blue", lty="longdash", lwd=2) legend("topright", inset=.05, title="Parameters", labels, lwd=2, lty=c(1, 1, 1, 1, 2), col=colors) }     ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Reply | Threaded
Open this post in threaded view
|

## Re: Accessing more than two coefficients in a plot

 Well, as a line in the plane is determined by 2 coefficients only, I'd guess that trying to find an R function that plots a line defined by 4 coefficients has about the same chance of success as finding a unicorn with 3 horns. You do understand that your linear model defines a hyperplane in your three covariates, do you not? Or do I misunderstand what you have requested? Cheers, Bert On Sun, Mar 25, 2012 at 2:32 PM, FJ M <[hidden email]> wrote: > > I've successfully plotted (in the plot and abline code below) a simple regression of Lambda1_2 on VV1_2. I then successfully regressed Lambda1_2 on VV1_2, VV1_22 and VV1_212 producing lm2.l. When I go to plot lm2.l using abline I get the warning: > > "1: In abline(lm2.l, col = "brown", lty = "dotted", lwd = 2) : only using the first two of 4 regression coefficients" > > Is there another function like abline that will produce a line using the constant and three coefficients from the lm2.l regression? > > > lm.l <- lm(Lambda1_2 ~ VV1_2, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression > > lm2.l <- lm(Lambda1_2 ~ VV1_2 + VV1_22 + VV1_212, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression > > plot(VV1_2, Lambda1_2, ylim=yrange, tck=1, main="V(1) Parameters (V, V^2 & V^0.5)", xlab="VV1_2", ylab="Lambda & Beta1_2",pch=19,col="red") > {abline(lm2.l, col="brown", lty="dotted", lwd=2) > abline(wlm2.l, col="gold",lty="longdash", lwd=2) > points(VV1_2, Beta1_2, pch=19, col="blue") > abline(lm2.b, col="black",lty="dotted", lwd=2) > abline(wlm2.b, col="blue", lty="longdash", lwd=2) > legend("topright", inset=.05, title="Parameters", > labels, lwd=2, lty=c(1, 1, 1, 1, 2), col=colors) > } > ______________________________________________ > [hidden email] mailing list > https://stat.ethz.ch/mailman/listinfo/r-help> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Reply | Threaded
Open this post in threaded view
|

## Re: Accessing more than two coefficients in a plot

 I needed to create my own forecast from the square root, linear and quadratic coefficients and then the abline() plot worked fine.   # Forecast l using non-linear regression coeffs - unweighted lm2.bforecast<- numeric(n) for (i in 1:n) { lm2.bforecast[i] <- lm2.b\$coeff["(Intercept)"]+lm2.b\$coeff["VV1_2"]*VV1_2[i]+lm2.b\$coeff["VV1_22"]*VV1_22[i]+lm2.b\$coeff["VV1_212"]*VV1_212[i] } lm2.bforecastline<-lm(lm2.bforecast ~ VV1_2, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted, non-linear regression forecast plot(VV1_2, Lambda1_2,  ylim=yrange, tck=1, main="Verizon V(1) Parameters (V, V^2 & V^0.5) Unweighted", xlab="VV1_2", ylab="Lambda1_2 & Beta1_2",pch=19,col="red") {points(VV1_2, lm2.lforecast, pch=19, col="brown") abline(lm2.lforecastline, col="brown", lty="longdash", lwd=2)  ...   > Date: Sun, 25 Mar 2012 15:36:20 -0700 > Subject: Re: [R] Accessing more than two coefficients in a plot > From: [hidden email] > To: [hidden email] > CC: [hidden email] > > Well, as a line in the plane is determined by 2 coefficients only, I'd > guess that trying to find an R function that plots a line defined by 4 > coefficients has about the same chance of success as finding a unicorn > with 3 horns. > > You do understand that your linear model defines a hyperplane in your > three covariates, do you not? Or do I misunderstand what you have > requested? > > Cheers, > Bert > > On Sun, Mar 25, 2012 at 2:32 PM, FJ M <[hidden email]> wrote: > > > > I've successfully plotted (in the plot and abline code below) a simple regression of Lambda1_2 on VV1_2. I then successfully regressed Lambda1_2 on VV1_2, VV1_22 and VV1_212 producing lm2.l. When I go to plot lm2.l using abline I get the warning: > > > > "1: In abline(lm2.l, col = "brown", lty = "dotted", lwd = 2) : only using the first two of 4 regression coefficients" > > > > Is there another function like abline that will produce a line using the constant and three coefficients from the lm2.l regression? > > > > > > lm.l <- lm(Lambda1_2 ~ VV1_2, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression > > > > lm2.l <- lm(Lambda1_2 ~ VV1_2 + VV1_22 + VV1_212, method = "qr", model = TRUE, x = FALSE, y = FALSE, qr = TRUE) # unweighted regression > > > > plot(VV1_2, Lambda1_2, ylim=yrange, tck=1, main="V(1) Parameters (V, V^2 & V^0.5)", xlab="VV1_2", ylab="Lambda & Beta1_2",pch=19,col="red") > > {abline(lm2.l, col="brown", lty="dotted", lwd=2) > > abline(wlm2.l, col="gold",lty="longdash", lwd=2) > > points(VV1_2, Beta1_2, pch=19, col="blue") > > abline(lm2.b, col="black",lty="dotted", lwd=2) > > abline(wlm2.b, col="blue", lty="longdash", lwd=2) > > legend("topright", inset=.05, title="Parameters", > > labels, lwd=2, lty=c(1, 1, 1, 1, 2), col=colors) > > } > > ______________________________________________ > > [hidden email] mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help> > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html> > and provide commented, minimal, self-contained, reproducible code. > > > > -- > > Bert Gunter > Genentech Nonclinical Biostatistics > > Internal Contact Info: > Phone: 467-7374 > Website: > http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm              [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.