r语言如何绘制带标准误的线性关系图

Python024

r语言如何绘制带标准误的线性关系图,第1张

r语言绘制带标准误的线性关系

1.打开文件,输入几组具有线性关系的数据。

2.用鼠标选中这些数据,点击菜单栏中的“插入”选项

3.在插入菜单中,选择一种散点图。

4.右击图表中的散点,在其右键菜单中点击“添加趋势线”选项。

5.在出来的页面中,选择“线性”,勾选“显示公式”选项,关闭窗口。

6.完成以上设置后,即可用excel制作线性关系图图表。

可能稍微有点麻烦。。。比如在向量v=(v[1],v[2],...,v[n])的第i个元素前加一个元素x变成(v[1],v[2],...,v[i-1],x,v[i],v[i+1],...,v[n])可以这样:

v<-c(v[1:(i-1)],x,v[i:length(v)])

LagrangePolynomial <- function(x,y) {

  len = length(x)

  if(len != length(y))

    stop("length not equal!")

   

  if(len < 2)

    stop("dim size must more than 1")

   

  #pretreat data abd alloc memery

  xx <- paste("(","a -",x,")")

  m <- c(rep(0,len))

   

  #combin express

  for(i in 1:len) {

    td <- 1

    tm <- "1"

    for(j in 1:len) {

      if(i != j) {

        td <- td*(x[i] - x[j])

        tm <- paste(tm,"*",xx[j])

      }

    }

    tm <- paste(tm,"/",td)

    m[i]<-tm #m[i] <- parse(text=tm)

  }

   

  #combin the exrpession

  m <- paste(m,"*",y)

  r <- paste(m,collapse="+")

   

  #combin the function

  fbody <- paste("{ return(",r,")}")

  f <- function(a) {}

   

  #fill the function's body

  body(f) <- parse(text=fbody)

   

  return(f)

}

这是拉格朗日多项式插值算法  你参考下吧