3

Estou tentando reproduzir um gráfico de animação de acordo com a evolução do parâmetro lambda. Abaixo está o código e a imagem.

Nesta parte são as funções

library(plotly)
library(splines)
bspline <- function(x,xl,xr,ndx,bdeg){
  dx <- (xr-xl)/ndx
  knots <- seq(xl-bdeg*dx,xr+bdeg*dx,by=dx)
  B <- spline.des(knots,x,bdeg+1,0*x,outer.ok=TRUE)$design
  output <- list(knots=knots,B=B)
  return(output)
}
tpoly <- function(x,t,p){
  B = NULL
  for (i in 1:length(t)){
    B <- cbind(B,(x-t[i])^p*(x>t[i]))
  }
  B
}


# P-spline regression fit with basis and coefficients
f <- sin(3*pi*x)
x <- seq(0,1,l=200)
B = list()
xl <- min(x)+0.01
xr <- max(x)+0.01
ndx = 15
pord <- 2
lam <- c(1e-3,1,100,1000)

Gráficos

graf = for(i in 1:4){
  plot(x,y,col="grey",pch=19,cex=.55); title(bquote(lambda == .(lam[i])))
  Basis <- bspline(x,0,1,ndx=ndx,bdeg=3)
  B <- Basis$B
  knots <- Basis$knots
  D <- diff(diag(ncol(B)),differences=pord)
  P <- lam[i]*t(D)%*%D
  theta <- solve(t(B)%*%B + P,t(B)%*%y)
  fit <- B%*%theta
  lines(x,fit,col=i,lwd=3,lty=1)
  matlines(x,B%*%diag(c(theta)),col=i,lty=2)
  points(knots,0*knots,col=i,cex=1.5,lwd=2.5,pch=15)
  points(knots[-c(1:2,length(knots)-1,length(knots))],theta,col=i,cex=1.5,lwd=2.5)
}
#####################################################
######################## Animation ##################
#####################################################
x11()
par(mfrow=c(2,2))
saveGIF(graf,interval=<0.05>,movie.name=<GRAFico>,outdir=<C:\Users\breni>)


O que estou tentando fazer é colocar um controlador para a variação do parâmetro Lambda e este gráfico se torna uma animação. Figure 1: enter image description here

1 Resposta 1

5

Se você colocar o código para produzir a imagem dentro da função saveGIF() funciona:

lam <- seq(1e-3, 1000, length.out = 100)
animation::saveGIF(
  {
    for(i in 1:length(lam)){
      plot(x, y, col = "grey", pch = 19, cex = .55, ylim = c(-1.5, 1.5))
      title(bquote(lambda == .(lam[i])))
      Basis <- bspline(x, 0, 1, ndx = ndx, bdeg = 3)
      B <- Basis$B
      knots <- Basis$knots
      D <- diff(diag(ncol(B)), differences = pord)
      P <- lam[i] * t(D) %*% D
      theta <- solve(t(B) %*% B + P, t(B) %*% y)
      fit <- B %*% theta
      lines(x, fit, col = i, lwd = 2, lty = 1)
      matlines(x, B %*% diag(c(theta)), col = i, lty = 2)
      points(knots, 0 * knots, col = i, cex = 1, lwd = 1.5, pch = 15)
      points(knots[-c(1:2, length(knots) - 1, length(knots))], theta, col = i, cex = 1, lwd = 1.5)
    }
  }, interval = 0.1, movie.name = "test.gif"
)

inserir a descrição da imagem aqui

4
  • 1
    Obrigado Willian!
    – user55546
    Commented 28/10/2019 às 13:38
  • Willian quem você utilizou como y?
    – user55546
    Commented 28/10/2019 às 13:45
  • Utilizei f como y pra criar o gif. Commented 28/10/2019 às 19:46
  • 1
    Ok!! Valeu pela ajuda!!
    – user55546
    Commented 29/10/2019 às 0:10

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .