Análise de regressão simples no R

Existe uma grande quantidade de análises que podem ser feitas encima de uma análise de regressão, como por exemplo uma análise de resíduos. Porém meu intuito aqui é só mostrar a parte introdutória.

Para iniciar precisamos de duas variáveis, uma dependente ou resposta e a outra independente ou explicativa, que chamaremos de x e y e que serão gerados de forma aleatória. São gerados 20 valores atravéz da distribuição uniforme entre 1 e 10 para x e entre 10 e 20 para y, após isso esses valores são ordenados em ordem crecente pelo comando sort().

set.seed(1)
x=sort(runif(20,1,10))
set.seed(2)
y=sort(runif(20,10,20))

Note que usei o comando set.seed antes de cada comando para que os valores aleatórios sejam os mesmo em qualquer lugar em que eles forem gerados, ou seja, os resultados daqui serão iguais pra qualquer pessoa que queira reproduzir o exemplo. Anyway, seguimos começando a análise onde a primeira coisa que olhamos é a correlação das variáveis.

cor(x,y)
[1] 0.982508

Veja que é uma ótima correlação. Em seguida fazemos a análise de regressão com o comando lm() de onde construímos nossa reta y=a+bx. Note que o y~x é a explicação de como deve ser a equação, ou seja, x explica y, ou ainda, y=x.

regressao=lm(y~x) summary(regressao)

Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -0.9339 -0.4361 0.1189 0.2904 0.9063 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 8.14171 0.33723 24.14 3.65e-15 *** x 1.16125 0.05188 22.38 1.37e-14 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.5823 on 18 degrees of freedom Multiple R-squared: 0.9653, Adjusted R-squared: 0.9634 F-statistic: 501.1 on 1 and 18 DF, p-value: 1.367e-14

Com isso nossa reta fica y= 8.14171 + 1.16125 x onde tanto o coeficiente angular quanto o linear são segnificantes à 1%. Maravilha de regressão! E pra completar o grafico de dispersão com a reta de regressão.

plot(x,y)
abline(regressao)

Guia cantado dos episódios de TNG

Sabem o que eu mais gosto deste vídeo? É saber que eu já assisti todos os episódios de TNG, HAHAHAHAAHAA, ai ai ah esquece... sim TNG = Star Trek: The Next Generation.

Distribuição de probabildiade de GURPS

Uma das primeiras coisas que quis fazer quando comecei a estudar probabilidade era construir a distribuição de probabilidade dos jogos de RPG que conheço. Essa semana consegui a de GURPS. Para vocês que não jogam RPG entenderem o sistema funciona da seguinte forma: O personagem possui um nível de habilidade (NH) em certa pericia, por exemplo salto. Joga-se 3 dados de 6 faces (normalmente conhecido como d6) e compare-se os resultados com o NH, se o resultado for menor ou igual ao NH o personagem é bem sucedido. Os testes de GURPS sofrem a critica que um redutor no NH causa uma redução da chance de acerto diferente para diferentes NHs, por exemplo ao retirar 1 de um NH 12 você causará um redutor de 11,6% na chance de passar em um teste, já um -1 em um NH 17 causará apenas 1,4% de penalidade. Mas apesar disso ainda acho GURPS meu sistema favorito.

Consegui esses resultados usando os seguintes comandos no R

d1=d2=d3=c(1:6)

resultados=expand.grid(d1,d2,d3)

soma.resultados=resultados$Var1+resultados$Var2+resultados$Var3

table(soma.resultados)

Isso ai me deu a tabela de frequência, o resto eu fiz no Excel mesmo, :)

O meu maior desafio será montar a FP de Storyteller...

Motores FTL

Bom dia a todos, humanos ou não, meu nome é Alysson Lívio e estou começando um bendito de um blog. Alguns amigos começaram essa mania e agora estou embarcando nela porque sinto vontade de me sentir perto das pessoas que convivi por tanto tempo bem próximo e agora, após a formatura, quase não vejo mais. Adoro compartilhar minha vida, e meus amigos sabem disso.

Estou escolhendo um nome pra o blog que me identifico com ele. Motores mais rápidos que a luz é um tema recorrente a qualquer seriado/filme de ficção cientifica que sou tão aficcionado. Espero também que essa seja a velocidade com que consiga produzir textos, pois não pretendo tornar o blog um eterno "A imortalidade e os livros", postagem única do blog do meu irmão (hahahaha).

Dentre os temas que estou pretendendo abordar aqui estão alguns assuntos Nerds, mas vou me dedicar principalmente à abordagem da Estatística como apostilas, questionários, e soluções estatísticas em R. Quero também ter a oportunidade de colocar algumas fotos, minhas ou de amigos, da tão bonita cidade de João Pessoa que é uma coisa que eu tenho vontade a muito tempo de fazer.

Bom... Espero que alguns gostem e que postem!