Walmes Marques Zeviani

"Experiência não é o que acontece com você, mas o que você fez com o que lhe aconteceu." (Aldous Huxley)

Engenheiro Agronômo pela Universidade Federal da Grande Dourados (2007) e Mestre em Estatística & Experimentação Agropecuária pela Universidade Federal de Lavras (2009). Doutorando do Programa de Pós-Graduação em Estatística & Experimentação Agropecuária - Professor Assistente I do Departamento de Estatística da Universidade Federal do Paraná (aprovado por Concurso e contratado em Julho 2010) onde atualmente leciona.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573                    VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br                      twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes    linux user number: 531218
==========================================================================

Avaliação de modelos de regressão não linear na cinética de liberação de potássio de resíduos orgânicos

Resumo: O potássio (K) é um elemento químico muito exigido pelas culturas, porém, sua ocorrência e disponibilidade em alguns solos brasileiros são inferiores à demanda pelas plantas. Seu fornecimento também pode ser feito a partir de fontes orgânicas de maneira eficiente, quando se conhece o padrão de liberação do nutriente. Modelos não lineares são adequados nessas situações, uma vez que estimam quantidades de interesse prático e apresentam boa qualidade de ajuste. Embora seu processo inferencial seja baseado em argumentos assintóticos, existem meios de se conhecer a intensidade da não linearidade. Neste trabalho avaliou-se a não linearidade, por meio das curvaturas de Bates & Watts, vício de Box e estudo das propriedades amostrais dos estimadores de mínimos quadrados por simulação, de dois modelos de regressão não linear. Tais modelos estimam o conteúdo de K prontamente liberável, o tempo de meia vida para a sua liberação e a taxa de liberação do K de lenta liberação. Os dados são provenientes do estudo, ao longo do tempo, da liberação de K de 11 resíduos orgânicos incubados em combinação com 3 tipos de solo. O modelo exponencial apresentou estimativas mais precisas para os parâmetros que o quociente e maior aproximação entre os resultados assintóticos e obtidos por simulação. O modelo exponencial foi mais adequado, em termos inferenciais e para aplicação prática, uma vez que, por todas as medidas, apresentou menor não linearidade.

Palavras-chaves: Curvatura de Bates e Watts. Vício de Box. Bootstrap. Tempo de meia vida.

Banca de Defesa:
Prof. Dr. Joel Augusto Muniz (Orientador) - UFLA
Prof. Dr. Carlos Alberto Silva (co-orientador) - UFLA
Profa. Dra. Taciana Villela Savian - UFLA
Prof. Dr. Luiz Alberto Beijo - UNIFAL-MG
Aprovada em 18 de setembro de 2009

O acesso ao LEG se dá pela escada/rampa localizada na face norte do bloco de Informática. Ao entrar pelo corredor do antigo salão de provas, procure por uma porta sinalizada do lado direito. Clique aqui para abrir a nossa localização em uma nova página com seta de indicação. <googlemap width="533px" height="400px" lat="-25.450603" lon="-49.231759" type="hybrid" zoom="19"> </googlemap>

• Site para produzir expressões matemáticas em latex (http://www.codecogs.com/latex/eqneditor.php).
  Permite a construção online de expressões matemáticas em código Latex, que podem ser copiadas e coladas no seu editor. O usuário tem opção de usar os menus com expressões pré-definidas, tais como matrizes, sistema de equações, derivadas, somatórios, simbologia de conjuntos, letras gregas, etc. Expressões comuns estão organizadas no botão "exemplos".
• Download do LaTable para confecção de tabelas para Latex (http://www.leg.ufpr.br/~walmes/docs/LaTable.exe).
  A versão disponível para Download tem compatibilidade com o wine e portanto pode ser usada por usuários Linux (usando o wine). Nos casos em as tabelas estão prontas em editores de texto (e.g. OpenOffice-Writer) ou em planilhas eletrônicas (e.g. OpenOffice-Calc) pode-se colar o conteúdo da tabela em um arquivo de extensão *txt e importar pelo LaTable. No LaTable é possível fazer a formatação de alinhamento de colunas e inserção linhas verticais e horizontais, usadas, em geral, para separar o cabeçalho da tabela do restante. Depois da edição é só copiar, em formato Latex, o conteúdo da tabela para área de transferência e colar no seu editor.
• Zotero: plugin do Firefox para extração de meta-dados de páginas de internet para criação de referências bibliográficas (http://www.zotero.org/).
  Zotero é capaz de extrair os metadados de sites da internet e produzir referências bibliográficas em diversos formatos. Sites Amazon.com (de venda de livros) e o ScienceDirect (para busca de artigos) são alguns dos compatíveis com o plugin. Basta fazer a busca por um livro/artigo no site que um ícone irá aparecer na sua barra de endereços. Ao clicar no ícone (que será um livro de capa azul no caso de livros) os metadados serão organizados pelo zotero. Para acessar os metadados, clica-se no nome "zotero" que fica no canto inferior direito do seu Firefox. Pode-se exportar as referências dos livros/artigos (e outros tipos de publicação) para diversos formatos, inclusive para o formato usado pelo bibtex, além de formatos para OpenOffice-Writer, html, etc. Mais detalhes consultar a documentação no site do plugin. Consulte também: tutorial em pdf, descrição, video.
• Site para criar senhas aleatórias, com diversas opções como tamanho e tipos de caracteres e criação de senhas pronunciáveis.
• Para converter as figuras em formato *pdf de um diretório para *jpg e fazer um zip delas no final (dica Linux). Salvar num arquivo com extensão *sh (e.g. converter.sh), dentro do diretório das figuras, o conteúdo abaixo (o tipo de aspas é importante)

for file in `ls *.pdf`
do
   convert $file `echo $file | sed 's/\.pdf$/\.jpg/'`
done
zip todojpg.zip *.jpg
rm *.jpg

$ chmod +x converter.sh
$ ./converter.sh

• Para fazer uma busca de uma expressão reugular dentro de uma árvore de diretórios apenas nos arquivos de determinada extensão. No caso, abaixo está ilustrada a procura dos arquivos de extensão R que possui a palavras 'try'

$ grep -r --include=*.R 'try' *

Trocar texto dentro de um arquivo. O texto 'walmes' será substituido por 'zeviani' dentro do arquivo experII.R

$ perl -i -pe 's/walmes/zeviani/;' experII.R 

Para redimensionar as fotos (*jpg, ou gráficos do R, rs) de um diretório, deixando as fotos com largura de 640 pixels e mantendo proporção original, e forçando a proporção 640x480.

$ mogrify -resize 640 *.jpg
$ mogrify -resize 640×480! *.jpg

Colocar alias no wiki. Logar na pata** via console

$ ssh -p <porta> <usuário>@<endereçomáquina>
$ sudo vim /etc/apache2/conf.d/alias (use esc : wq para sair e salvar)
$ sudo /etc/init.d/apache2 restart

• Para instalar pacotes do R no linux que possuem dependências de elementos do sistema operacional, geralmente pacotes gráficos (rgl, playwith, gWidgetsRGtk2), use o comando

$ sudo apt-get build-dep r-cran-nomedopacote

• para procurar e substituir no emacs dê alt+shift+5 (ou alt+%), abre o minibuffer, escreve nome de busca ENTER, escreve o nome de substituição ENTER, dê alt+% para avançar na tarefa, y para confirmar substituição, n para não substituir.
• para mudar o cursor de buffer no emacs use ctrl+c+o.
• para deixar o gedit com corretor ortográfico em português http://beaglenetwork.blogspot.com/2010/05/corretor-ortografico-em-portugues-para.html.
• para converter arquivos de texto da codoficação utf-8 para a iso-8859-15

#!/bin/bash 
# esse bash faz a conversão de arquivos R em codificação utf-8 para iso-8859-1
# $1 é para fornecer um nome para ir no nome arquivos, usar então o que está abaixo
$1
for a in *.R;
 do
  echo "Convertendo $a em utf-8 para ${a%%.R}-$1.R em iso-8859-15 ..."
   iconv -f utf-8 -t iso-8859-15 "${a}" -o "${a%%.R}-$1.R"
done

$ chmod +x utf2iso.sh
$ ./utf2iso.sh

• para montar pendrive com permissões de escrita

sudo mount -o remount,rw /media/HP\ v115w

1. Estatística Básica (ce003 - estatística II) para os cursos de Administração turmas N2 e N3 (2 semestre 2010)
   • Página da disciplina para as turmas N2 e N3;
   • Histórico das aulas para as tumas N2 e N3
   • Download da primeira avaliação de Estatística Básica com gabarito.
   • Download das notas e faltas referentes ao período até a primeira avaliação.
   • Download da lista com exercícios referentes ao conteúdo da segunda avaliação.
   • Download dos resultados dos exercícios referentes ao conteúdo da segunda avaliação.
   • Download das notas e faltas referentes ao período até a segunda avaliação.
   • Download da lista com exercícios referentes ao conteúdo da terceira avaliação.
   • Download das tabelas estatísticas (da normal, da t e outras).
   • Download das notas da terceira avaliação.
   • Download das notas do exame com as médias finais. Alunos que querem revisão de prova devem me procurar até as 19 horas de 20/12/2010.
2. Estatística Básica (ce003 - estatística II) para o curso de Engenharia Cívil turma B (1 semestre 2011)
   • Página da disciplina Turma B;
3. Estatística Computacional (ce223 - estatística II) para o curso de Estatística (1 semestre 2011)
   • Página da disciplina;

#==========================================================================================
# Estimação de parâmetros de modelos não lineares. Ajuste de modelos para a construção de
# intervalo hídrico ótimo
#==========================================================================================
 
#------------------------------------------------------------------------------------------
# pacotes necessários
 
install.packages("plyr", dep=TRUE, repos="http://cran-r.c3sl.ufpr.br/")
library(plyr)
 
#------------------------------------------------------------------------------------------
# entrada de dados
 
iho <- read.table("http://www.leg.ufpr.br/~walmes/docs/iho-data.txt", header=TRUE, sep="\t")
iho <- read.table("iho-data.txt", header=TRUE, sep="\t")
iho <- iho[,-c(1:4)]
names(iho) <- tolower(names(iho))
str(iho)
 
#------------------------------------------------------------------------------------------
# tratamento dos dados: usar valor médio das repetições para ajuste da curva
 
iho <- ddply(iho, .(tensao), "mean")
str(iho)
 
#------------------------------------------------------------------------------------------
# ajuste da curva característica segundo modelo theta = a*Psi^b*Ds^c
 
cra <- nls(theta~A*tensao^B*ds^C, data=iho,
           start=list(A=1, B=1, C=1))
summary(cra)
 
#------------------------------------------------------------------------------------------
# ajuste do modelo nulo (sem efeito de Ds e tensão)
 
nulo.cra <- nls(theta~mu, data=iho, start=list(mu=1))
 
#------------------------------------------------------------------------------------------
# testa a significância do modelo cra
 
an.cra <- anova(nulo.cra, cra)
an.cra
 
#------------------------------------------------------------------------------------------
# R^2 do modelo cra
 
R2.cra <- 1-deviance(cra)/deviance(nulo.cra)
R2.cra
 
#------------------------------------------------------------------------------------------
# verifica o ajuste observado vs predito (7 pontos apenas não dá pra fazer julgamento)
 
plot(iho$theta~fitted(cra))
 
#------------------------------------------------------------------------------------------
# ajuste do modelo de resistência a penetração
 
rp <- nls(rp~D*theta^E*ds^F, data=iho,
         start=list(D=1, E=1, F=1))
summary(rp)
 
#------------------------------------------------------------------------------------------
# ajuste do modelo nulo (sem efeito de theta e Ds)
 
nulo.rp <- nls(rp~mu, data=iho, start=list(mu=1))
 
#------------------------------------------------------------------------------------------
# testa a significância do modelo rp
 
an.rp <- anova(nulo.rp, rp)
an.rp
 
#------------------------------------------------------------------------------------------
# R^2 do modelo rp
 
R2.rp <- 1-deviance(rp)/deviance(nulo.rp)
R2.rp
 
#------------------------------------------------------------------------------------------
# verifica o ajuste observado vs predito (7 pontos apenas não dá pra fazer julgamento)
 
plot(iho$rp~fitted(rp))
 
#------------------------------------------------------------------------------------------
# tabela com as estimativas dos parâmetros e intervalos de confiânça
 
summary(cra)$coeff
summary(rp)$coeff
confint.default(cra)
confint.default(rp)
 
#------------------------------------------------------------------------------------------
# juntar tudo numa tabela
 
cbind(rbind(summary(cra)$coeff, summary(rp)$coeff),
      rbind(confint.default(cra), confint.default(rp)))
 
#------------------------------------------------------------------------------------------
# valor de F e R^2 retirar dos resultados de anova() e R^2
 
data.frame(Fval=c(an.cra$F[2], an.rp$F[2]),
           pval=c(an.cra$Pr[2], an.rp$Pr[2]),
           R2=c(R2.cra, R2.rp))
 
#------------------------------------------------------------------------------------------

<object width="636" height="516"><embed src="http://www.youtube.com/v/tDAB6BIz5nQ?version=3" type="application/x-shockwave-flash" allowscriptaccess="always" allowfullscreen="true" width="636" height="516"></embed></object>

Uma variação do problema dos aniversários http://pballew.blogspot.com/2011/01/getting-full-set.html

• Apostila do Curso (download);
• Dados usados durante Curso em planilha eletrônica (*.xls) (download);
• Dados usados durante Curso em arquivos de texto (*.txt) compactados (*.zip) (download);
• Questionário do curso (responder);

• https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
• http://r-br.2285057.n4.nabble.com/
• http://www.mail-archive.com/r-br@listas.c3sl.ufpr.br/