\documentclass[a4paper, oneside, 10pt]{article}
\usepackage[unicode]{hyperref}
\usepackage[utf8x]{inputenc}
\usepackage[english]{babel}
\usepackage{graphicx}
\graphicspath{{media/}}
\date{\today}
\title{}
\author{}
\begin{document}



\subsection{\texorpdfstring{Gisela Martini}{Gisela Martini}}
\label{sec:gisela_martini}


\includegraphics[keepaspectratio=true,width=0.8\textwidth]{bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/img-20150301-wa0009_-_copy}


Doutoranda em Ciências no Centro de Química e Meio Ambiente, IPEN - USP. Mestre pelo Instituto de Pesquisas Energéticas e Nucleares - Universidade de São Paulo (2013), atuou em projeto de pesquisa com ensaios de toxicidade e efeitos biológicos da radiação em organismos aquáticos. Graduada Ciências Biológicas com ênfase em Biologia Marinha pela Universidade Santa Cecília - Santos, SP (2010). Experiência nas áreas de: Gestão, Controle e Monitoramento em Ecotoxicologia. (\url{http://lattes.cnpq.br/0575674204746522})

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/ex_1_3.r}{ex\_1\_3.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/EXERCICIOS AULA 3 - NOTAR.txt}{exercicios\_aula\_3\_-\_notar.txt}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/4.2-4.5.r}{4.2-4.5.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/5.1-5.3.r}{5.1-5.3.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/106.2.r}{106.2.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/106.3.r}{106.3.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/exercicio_regressao_multipla.r}{exercicio\_regressao\_multipla.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/exercicio_8.r}{exercicio\_8.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/analise_exploratoria.r}{analise\_exploratoria.r}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/proposta_a.txt}{proposta\_a.txt}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/proposta_b.txt}{proposta\_b.txt}

\subsection{\texorpdfstring{PROPOSTA A}{PROPOSTA A}}
\label{sec:proposta_a}


Calculo do IQA e correlacao com ATIVIDADE ESTROGENICA em pontos de monitoramento de aguas superficiais do estado de Sao Paulo

O Indice de Qualidade das Aguas (IQA) incorpora nove variaveis consideradas relevantes para a avaliacao da qualidade das aguas, tendo como determinante principal a sua utilizacao para abastecimento publico. A proposta tem por objetivo criar uma funcao para o calculo do IQA utilizando como dados de entrada nove variaveis inseridas separadas por colunas em uma matriz, cujas linhas representam os valores amostrados (\textbf{locais de amostragem, certo?}). O IQA é calculado pelo produtorio ponderado da qualidade da agua correspondente as variaveis que integram o indice. A seguinte formula sera utilizada:   

IQA = prod(i=1,n)[qi\textasciicircum{}(wi)]

Onde: 
IQA: Indice de Qualidade das Aguas, um numero entre 0 e 100;   
qi: qualidade do i-esimo parametro, um numero entre 0 e 100, obtido da respectiva “curva media de variacao de qualidade”, em funcao de sua concentracao ou medida e,   
wi: peso correspondente ao i-esimo parametro, um numero entre 0 e 1, atribuido em funcao da sua importancia para a conformacao global de qualidade, sendo que:  

sum(i=1,n)[wi=1]

Em que: 
n: numero de variaveis que entram no calculo do IQA.
A qualidade da agua e determinada pelo resultado do calculo do IQA variando numa escala de 0 a 100 (CETESB, 2013).  
Disponível em: 
\url{http://www.cetesb.sp.gov.br/userfiles/file/agua/aguas-superficiais/apendices-relatorio-aguas-superficiais2013/Ap\%C3\%AAndice\%20C\%20\%20\%C3\%8Dndices\%20de\%20Qualidade\%20das\%20\%C3\%81guas.pdf}, acesso em 23 de fevereiro de 2015.

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/dados.csv}{dados.csv}

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/codigo_iqa.r}{codigo\_iqa.r}

Código da função:

\#Serão apresentados objetos para as variáveis utilizadas no cálculo do IQA, dentro de cada objeto estão os valores referentes ao ponto de amostragem

\#Lendo o arquivo com o conjunto de dados

read.table(,,dados.csv", header=T, sep=,,,", as.is=F)

\#Criando um data frame contendo as 9 variáveis utilizadas para o cálculo do índice de qualidade das águas - IQA

indice=read.table(,,dados.csv", header=TRUE, sep=,,;")
X=indice

\#Verificação dos dados de entrada

str(x)
dim(x)
D=dim(x)
head(x)

\#Criando os objetos para o cálculo do qi (qi é um valor calculado para fornecer a qualidade da variável isolada, esse número pode variar de 0-100)
\#qi.(cada variável)foi definida como sendo = 0 (zero) para que prossiga com o cáculo de cada parâmetro dentro da função

qi.ct=0
qi.ph=0
qi.dbo=0
qi.nt=0
qi.ft=0
qi.temp=0
qi.turb=0
qi.st=0
qi.od=0

\#Objetos contendo os dados dos pontos amostrados

x1=X[1,]
x2=X[2,]
x3=X[3,]
x4=X[4,]
x5=X[5,]
x6=X[6,]
x7=X[7,]
x8=X[8,]
x9=X[9,]
x10=X[10,]

\#Criando a função

IQA= function(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10)
\begin{verbatim}
{
\end{verbatim}


\#Definindo as variáveis de entrada

ct=(x1[1])
ph=(x1[2])
dbo=(x1[3])
nt=(x1[4])
ft=(x1[5])
temp=(x1[6])
turb=(x1[7])
st=(x1[8])
od=(x1[9])

\#Calculando os valores de qi x1

\#Calculo do qi para variavel ct - Coliformes termotolerantes

lg=log10(ct)

if(lg$\Leftarrow$1) 
\begin{verbatim}
                              { 
                              qi.ct=100-(33*lg) 
                              }
\end{verbatim}


if(lg\textgreater 1 \& lg$\Leftarrow$5) 
\begin{verbatim}
                              { 
                              qi.ct=100-(37.2*lg)+(3.60743*log10(ct^2)) 
                              }
\end{verbatim}


if(lg\textgreater 5) 
\begin{verbatim}
                              { 
                              qi.ct=3
                              }
\end{verbatim}


\#Calculo do qi para variavel ph - pH

if(ph$\Leftarrow$2) 
\begin{verbatim}
                              { 
                              qi.ph=2 
                              }
\end{verbatim}


if(ph\textgreater 2 \& ph$\Leftarrow$4) 
\begin{verbatim}
                              { 
                              qi.ph=13.6-(10.6*ph)+(2.4364*(ph^2)) 
                              }
\end{verbatim}


if(ph\textgreater 4 \& ph$\Leftarrow$6.2) 
\begin{verbatim}
                              { 
                              qi.ph=155.5-(77.36*ph)+(10.2481*(ph^2))
                              }
\end{verbatim}


if(ph\textgreater 6.2 \& ph$\Leftarrow$7) 
\begin{verbatim}
                              { 
                              qi.ph=-657.2+(197.38*ph)-(12.9167*(ph^2))
                              }
\end{verbatim}


if(ph\textgreater 7 \& ph$\Leftarrow$8) 
\begin{verbatim}
                              { 
                              qi.ph=-427.8+(142.05*ph)-(9.695*(ph^2))
                              }
\end{verbatim}


if(ph\textgreater 8 \& ph$\Leftarrow$8.5) 
\begin{verbatim}
                              { 
                              qi.ph=216-(16*ph)
                              }
\end{verbatim}


if(ph\textgreater 8.5 \& ph$\Leftarrow$9) 
\begin{verbatim}
                              { 
                              qi.ph=1415823*exp(-1.1507*ph)
                              }
\end{verbatim}


if(ph\textgreater 9 \& ph$\Leftarrow$10) 
\begin{verbatim}
                              { 
                              qi.ph=228-(27*ph)
                              }
\end{verbatim}


if(ph\textgreater 10 \& ph$\Leftarrow$12) 
\begin{verbatim}
                              { 
                              qi.ph=633-(106.5*ph)+(4.5*(ph^2)) 
                              }
\end{verbatim}


if(ph\textgreater 12 \& ph$\Leftarrow$14) 
\begin{verbatim}
                              { 
                              qi.ph=3
                              }
\end{verbatim}


\#Calculo do qi para variavel dbo - DBO

if(dbo$\Leftarrow$5) 
\begin{verbatim}
                              { 
                              qi.dbo=99.96*exp(-0.1232738*dbo) 
                              }
\end{verbatim}


if(dbo\textgreater 5 \& dbo$\Leftarrow$15) 
\begin{verbatim}
                              { 
                              qi.dbo=104.67-(31.5463*(log10(dbo)))
                              }
\end{verbatim}


if(dbo\textgreater 15 \& dbo$\Leftarrow$30) 
\begin{verbatim}
                              { 
                              qi.dbo=4394.91*(dbo^-1.99809)
                              }
\end{verbatim}


if(dbo\textgreater 30) 
\begin{verbatim}
                              { 
                              qi.dbo=2
                              }
\end{verbatim}


\#Calculo do qi para variavel nt - Nitrogenio total

if(nt$\Leftarrow$10) 
\begin{verbatim}
                              { 
                              qi.nt=100-(8.169*nt)+(0.3059*(nt^2))
                              }
\end{verbatim}


if(nt\textgreater 10 \& nt$\Leftarrow$60) 
\begin{verbatim}
                              { 
                              qi.nt=101.9-23.1023*log10(nt)
                              }
\end{verbatim}


if(nt\textgreater 60 \& nt$\Leftarrow$100) 
\begin{verbatim}
                              { 
                              qi.nt=159.3148*exp(-0.0512842*nt)
                              }
\end{verbatim}


if(nt\textgreater 100) 
\begin{verbatim}
                              { 
                              qi.nt=1
                              }
\end{verbatim}


\#Calculo do qi para variavel ft - Fósforo total

if(ft$\Leftarrow$1) 
\begin{verbatim}
                              { 
                              qi.ft=99*exp(-0.91629*ft)
                              }
\end{verbatim}


if(ft\textgreater 1 \& ft$\Leftarrow$5) 
\begin{verbatim}
                              { 
                              qi.ft=57.6-(20.178*ft)+(2.1326*(ft^2))
                              }
\end{verbatim}


if(ft\textgreater 5 \& ft$\Leftarrow$10) 
\begin{verbatim}
                              { 
                              qi.ft=19.8*exp(-0.13544*ft)
                              }
\end{verbatim}


if(ft\textgreater 10) 
\begin{verbatim}
                              { 
                              qi.ft=5
                              }
\end{verbatim}


temp=94 
\begin{verbatim}
                              { 
                              qi.temp=94
                              }
\end{verbatim}


\#Calculo do qi para variavel turb - Turbidez

if(turb$\Leftarrow$25) 
\begin{verbatim}
                              { 
                              qi.turb=100.17-(2.67*turb)+(0.03775*(turb^2))
                              }
\end{verbatim}


if(turb\textgreater 25 \& turb$\Leftarrow$100) 
\begin{verbatim}
                              { 
                              qi.turb=84.768*exp(-0.016206*turb)
                              }
\end{verbatim}


if(turb\textgreater 100) 
\begin{verbatim}
                              { 
                              qi.turb=5
                              }
\end{verbatim}


\#Calculo do qi para variavel st - Sólidos totais

if(st$\Leftarrow$150) 
\begin{verbatim}
                              { 
                              qi.st=79.75+(0.166*st)-(0.001088*(st^2))
                              }
\end{verbatim}


if(st\textgreater 150 \& ft$\Leftarrow$500) 
\begin{verbatim}
                              { 
                              qi.st=101.67-(0.13917*st)
                              }
\end{verbatim}


if(st\textgreater 500) 
\begin{verbatim}
                              { 
                              qi.st=32
                              }
\end{verbatim}


\#Calculo do qi para variavel od - oxigênio dissolvido

if(od$\Leftarrow$50) 
\begin{verbatim}
                              { 
                              qi.od=3+(0.34*od)+(0.008095*(od^2))+1.35252*(0.00001*(od^3))
                              }
\end{verbatim}


if(od\textgreater 50 \& od$\Leftarrow$85) 
\begin{verbatim}
                              { 
                              qi.od=3-(1.166*od)+(0.058*(od^2))-3.803435*(0.0001*(od^3))
                              }
\end{verbatim}


if(od\textgreater 85 \& od$\Leftarrow$100) 
\begin{verbatim}
                              { 
                              qi.od=3+(3.7745)*(od^0.704889)
                              }
\end{verbatim}


if(od\textgreater 100 \& od$\Leftarrow$140) 
\begin{verbatim}
                              { 
                              qi.od=3+2.9*(od)-(0.02496)*(od^2)+5.60919*(0.000018(od^3))
                              }
\end{verbatim}


if(od\textgreater 140) 
\begin{verbatim}
                              { 
                              qi.od=3+47
                              }
\end{verbatim}


\#Criando objeto com os valores de peso para cada parâmetro (wi) (o peso wi, é uma constante aplicada a cada variável)

wi=c(0.15,0.12,0.10,0.10,0.10,0.10,0.08,0.08,0.17)

\#Criando um objeto com os valores de qi calculados para cada parâmetro na função

qi=c(qi.ct,qi.ph, qi.dbo,qi.nt, qi.ft, qi.temp, qi.turb, qi.st, qi.od)

\#Calcular o qi\textasciicircum{}wi para cada variável 

iqa.ct= (qi.ct\textasciicircum{}wi[1])
iqa.ph= (qi.ph\textasciicircum{}wi[2])
iqa.dbo= (qi.dbo\textasciicircum{}wi[3])
iqa.nt= (qi.nt\textasciicircum{}wi[4])
iqa.ft= (qi.ft\textasciicircum{}wi[5])
iqa.temp= (qi.temp\textasciicircum{}wi[6])
iqa.turb= (qi.turb\textasciicircum{}wi[7])
iqa.st= (qi.st\textasciicircum{}wi[8])
iqa.od= (qi.od\textasciicircum{}wi[9])

\#Calculando o IQA para x1

IQA1 = iqa.ct*iqa.ph*iqa.dbo*iqa.nt*iqa.ft*iqa.temp*iqa.turb*iqa.st*iqa.od

(OBS: repetição para cada ponto de amostragem, vide arquivo: iqa.r)

\href{media/bie5782/01_curso_atual/alunos/trabalho_final/gisela.martini/help.r}{help.r}

\end{document}
