Regressão linear usando Shell Bash + AWK
Publicado por Perfil removido (última atualização em 27/06/2016)
[ Hits: 2.362 ]
Ouvi falar do AWK enquanto estudava Programação Shell Linux, Júlio César; e de maneira bem cômica, ele explica um pouco dessa "Linguagem/Ferramenta/Tudo de bom". kkk
Então vi grande potencial para resolver alguns problemas que gastaria muito tempo, as vezes.
Bem, para quem faz alguma Análise de Dados, e desses dados necessite fazer regressão linear com cálculo das incertezas, vi que o AWK é excelente para fazer tal coisa.
E resolvi compartilhar com a comunidade esse pequeno trabalho, fiquem a vontade para compartilhar, melhorar e contribuir para o maior uso dessa fantástica ferramenta.
Então, basta fazer:
./linear_simples.sh "dados.dat"
por exemplo, e pronto, um arquivo com tudo. Ou quase tudo. ;)
#/bin/bash #---Regressao linear simples------------# # usando awk e shell # Autor: Marcos E. O. Sales # email: [email protected] # # feito:03/05/2016 # ultima atualizacao:04/05/2016 #---Regressao linear simples------------# #---(Utilidade)--------------------------------# #Fazer regressao atraves simples atravez # da reta y = a + b*x # atraves de uma entrada de dados. #-------------------------------------------------# # versao: 1.1 #==============================# #-------Definindo funcoes------------------------------------------------------# Uso() { sed -n '/^# ultima/,/^# versao:/ {/^# ultima/b; /^# versao:/b; 1,$p;}' $(basename $0) echo echo "Uso: $(basename "$0") ARQ.dat" echo "Onde ARQ.dat e um arquivo com os valores: Dados x Dados Y" echo -e "\nPara retornar um arquivo formatado com os dados,e com os parametros analizados" exit 0 } # Regressao() { clear sed '1d' "$Dados" |\ awk ' BEGIN { printf "%8s %13s %13s \n", "Dados X", "Dados Y", "sigma Y" } #colocar primeira linha { printf "%8e %13e %8e\n",$1, $2, $3 #printa dados formatados } #operacoes { Somax = Somax + $1 ; Somay = Somay + $2 Mediax = Somax / NR Mediay = Somay / NR somaquadx+=$1*$1 somaquady+=$2*$2 somaxy+=$1*$2 } { cov = ($1-Mediax)*($2-Mediay) desvx = sqrt(($1-Mediax)^2) desvy = sqrt(($2-Mediay)^2) delta = (NR*somaquadx-Somax^2) } END{ a = ((somaquadx*Somay)-(Somax*somaxy))/(delta) b = ((NR*somaxy) - (Somax*Somay))/(delta) sigmay = sqrt((1/(NR-2))*(($2-a-b*$1)^2)) sigmaa = sigmay*sqrt(somaquadx/delta) sigmab = sigmay*sqrt(NR/delta) r = cov/(desvx*desvy) printf "\n%7s\n%8s %7e\n%8s %10e\n %15s \n %1s %8e \n %8s %8e \n %1s %8e\n %8s %8e\n %8s %8e \n" , "\nParametros:\n", "Media X:", Mediax , "Media Y:", Mediay, "Regressao: y = b*x + a", "a:", a, "sigma a:",sigmaa,"b:",b,"sigma b:",sigmab, "Correlacao",r printf "\n%7s %8e\n%7s %8e \n%7s %8e\n", "covariancia:",cov, "desvio x",desvx, "desvio y", desvy}' |tee resultado.txt } #----Definindo o case----------------------------------------------------------# case $# in 0) Uso ;; 1) Dados="$1"; Regressao ;; *) echo -e "\033[01;31mNumero de parametros errados.\033[m" Uso ;; esac #--------------------------------------------------------------------------------------# # Fim do User #--------------------------------------------------------------------------------------#
Script Para Formatação de Texto Com Cabeçalho, Rodapé e Campos de Dados Originados de Arquivo em AWK
Pesquisa de grupos do sistema por nome ou GID em AWK
Mostrar as partições existentes no HD
Mostrar as partições existentes no HD
Script Para Formatação de Texto Com Cabeçalho, Rodapé e Campos de Dados Originados de Arquivo em AWK
Nenhum comentário foi encontrado.
Instalação do Gentoo GNU/Linux em um SSD NVMe (UEFI/GPT)
Adicionar mensagens iniciais no terminal nas sessões de sudo
Desenhar com o mouse na tela do desktop
Gerencie suas contas financeiras pessoais com Terminal Finances
Usando "sscanf" para Dividir String com Delimitadores em C
Processo aleatório consumindo CPU por meio do usuário Postgres
Linux trava Inexplicavelm... (10)
Meu note esquenta muito com linux (3)
Como escrever um artigo científico? (2)
Teclado parando de funcionar (0)
Usar notify-send pelo Root enviando para todos usuarios do sistema da ... (0)