Argthtjtr

I'm having problems with the following code, I do not know how to introduce this instruction below inside a cat instruction which is in the function.

df <- data.frame(x,y)

plot(y,x)

abline(lm(y ~ x))

R marks an error saying that a comma is needed between plot(x,y) and abline(lm(y~ x)) but when adding it, the display is

Warning messages:

1: In doTryCatch(return(expr), name, parentenv, handler) :

  display list redraw incomplete

2: In doTryCatch(return(expr), name, parentenv, handler) :

  invalid graphics state

3: In doTryCatch(return(expr), name, parentenv, handler) :

  invalid graphics state

What can I do to fix the problem?

Please someone help me, thank you in advance.

Summary<-function(x,y,print=TRUE)

{

  p<-2

  n<-length(x)





  x<-matrix(c(rep(1,n),x),n,p)

  bg<-solve(t(x)%*%x,t(x)%*%y)

  invx<-solve(t(x)%*%x)

  xty<-t(x)%*%y

  e<-y-x%*%bg

  SCT<-sum(y^2)-n*(mean(y)^2)

  SCE<-sum(e*e)

  SCRegression<-SCT-SCE



  sigma<-SCE/(n-p)



  varbeta<-sigma*invx





  seb0 <- sqrt(varbeta[1,1])

  seb1 <- sqrt(varbeta[2,2])



  alfa <- 0.05



  valuet <- qt(1-alfa/2,n-p) 



  valuechi1 <- qchisq(1-alfa/2,n-p)

  valuechi2 <- qchisq(alfa/2,n-p)







  bg[1]-valuet*seb0

  bg[1]+valuet*seb0







  bg[2]-valuet*seb1

  bg[2]+valuet*seb1





  tvalueb1 <- bg[2]/seb1

  tvalueb1

  tvalueb0 <- bg[1]/seb0

  tvalueb0





  valuepb0 <- 2*pt(1-tvalueb0,n-p)

  valuepb0



  valuepb1 <- 2*pt(tvalueb1,n-p)

  valuepb1





  SCE/valuechi1

  SCE/valuechi2





  x0 <- mean(x) 

  yhat <- bg[1]+bg[2]*x0

  x0



  varmu <- varbeta[1,1]+x0^2*varbeta[2,2]+2*x0*varbeta[1,2]





  semu <- sqrt(varmu)

  semu



  yhat - valuet*semu

  yhat + valuet*semu



  Rsq <- SCRegression/SCT

  Rsq



  yhatt <- bg[1]+bg[2]*x

  yhatt



  ee <- y-yhatt



  df <- data.frame(x,y)



  results <- list(TSS=SCT, ESS=SCE, p=p, x=x, y=y)



  if (SCE == 0) warning("Error square sum is zero", call.=FALSE)



  if (print) {

    cat("Results for the variables", "nt",

        deparse(match.call()$x), " and ", deparse(match.call()$y),

        "nn", sep="")

    cat("The total square sum is: ", SCT, "nn",

        "The error square sum is: ", SCE, "nn",

        "The errors graph to determine homogeneity or heterogeneity is: ",plot(x,ee) , "nn",

        "The linear model plot: ",plot(y,x) abline(lm(y ~ x)) , "nn",

        sep="")

    invisible(results)

  } else {

    results

  }      

}