%COMANDOS EM R REFERENTES AO EXEMPLO RESTAURANTES

  restaurante = read.table("restaurante.txt",header=TRUE)
  attach(restaurante)

  require(robustbase)

  adjbox(faturamento, ylab="Faturamento")
  adjbox(publicidade, ylab="Gastos Publicidade")
  plot(publicidade, faturamento, ylab="Faturamento", xlab="Gastos Publicidade")
  lines(smooth.spline(publicidade, faturamento, df=3))

  ajuste1.restaurantes = lm(faturamento ~ publicidade)
  summary(ajuste1.restaurantes)

  fit.model = ajuste1.restaurantes
  source("diag_resid_norm")
  source("envel_norm")

%AJUSTE REGRESSÃO PONDERADA NO LM E GAMLSS
  peso = -9226.002 + 7781.626*publicidade
  z = faturamento/sqrt(peso)
  u1 = 1/sqrt(peso)
  u2 = publicidade/sqrt(peso)

  ajuste2.restaurantes = lm(z ~ u1 + u2 - 1)
  summary(ajuste2.restaurantes)

  fit.model = ajuste2.restaurantes
  source("diag_resid_norm")
  source("diag_cook_norm")
  source("envel_norm")

  ajuste3.restaurantes = lm(z ~ u1 + u2 - 1,subset=-1)
  summary(ajuste3.restaurantes)

  require(gamlss)

  ajuste4.restaurantes = gamlss(z ~u1 + u2 -1, family=NO)
  summary(ajuste4.restaurantes)
  plot(ajuste4.restaurantes)

%AJUSTE DUPLO MÉDIA E VARIÂNCIA NO DGLM

  require(dglm)
  ajuste5.restaurantes = dglm(faturamento ~ publicidade, ~ publicidade, family=gaussian)
  summary(ajuste5.restaurantes)


%AJUSTE DUPLO MÉDIA E DESVIO PADRÃO NO GAMLSS
  
  ajuste6.restaurantes = gamlss(faturamento ~publicidade, family=NO)
  summary(ajuste6.restaurantes)
  plot(ajuste6.restaurantes)

  ajuste7.restaurantes = gamlss(faturamento ~publicidade,~publicidade, family=NO)
  summary(ajuste7.restaurantes)
  plot(ajuste7.restaurantes)