Janeiro 2014 - R-br

Curso de estatistica com R gratuito em Stanford
by Wenceslau Teixeira 27 Jan '14

27 Jan '14

Excelente dica ! No curso tem também uma entrevista com John Chambers Um dos criadores da linguagem S e R https://www.youtube.com/watch?v=jk9S3RTAl38 E para os que não tiver tempo é possivel baixar o pdf do livro (legalmente) An introduction to statistical learning with applications in R (Gareth James, Daniela Witten, Trevor Hastie and Robert Tibshirani - Springer - 2013) http://www-bcf.usc.edu/~gareth/ISL/ Sds Wenceslau -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

Curso de estatistica com R gratuito em Stanford
by Fernando 25 Jan '14

25 Jan '14

Caros amigos Foi publicado no R-bloggers (http://www.r-bloggers.com/online-class-on-statistical-learning/) o anúncio de um curso on-line da universidade de Stanford sobre statistical learning (https://class.stanford.edu/courses/HumanitiesScience/StatLearning/Winter201…) com os autores do livro (http://www-bcf.usc.edu/~gareth/ISL/ <http://www-bcf.usc.edu>)). O curso iniciou dia 20/01 mas ainda pode ser feito inscrição. Basta cadastra (pelo que entendi). Abraços

2 1

Dúvida - elipse de confiança
by Mauro Sznelwar 24 Jan '14

24 Jan '14

Preciso traçar elipse de confiança no fim desta sequencia, queria saber como entro com os dados, que está no último comando abaixo: Data set:https://www.dropbox.com/s/405ggda0h6scnie/marajoara320amostras.txt marajoara<-read.table('marajoara320amostras.txt') =======NORMALIZAÇÃO======== marajoara.norm<-sapply(data.frame(marajoara),scale) =======CLUSTER========= win.graph() marajoara.norm.h<-hclust(dist(marajoara.norm)^2,method="ward") plot(marajoara.norm.h,main="marajoara320amostras") =======PCA======== win.graph() initial<-tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3),ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1<-kmeans(marajoara.norm,initial) marajoara.norm.pca<-princomp(marajoara.norm) marajoara.norm.px<-predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type="n",xlab="componente 1",ylab="componente 2") text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster),col=as.numeric(marajoara.norm.km1$cluster)) title("marajoara320amostras") # identificar amostras no gráfico de PCA identify(marajoara.norm.px[,1:2]) # Ellipse no gráfico de PCA (********DUVIDA*********) library(ellipse) plot(ellipse(marajoara.norm.px, centre = ???, t = ???), type = "???")

1 0

medida repetida?
by Marcelo Laia 20 Jan '14

20 Jan '14

Temos o seguinte: PlacaPetri Trat Time Mortos Vivos %Morto %Vivo 1 A 48 8 6 57,14 42,85 1 A 72 17 15 1 A 96 18 28 2 A 48 17 15 2 A 72 24 33 2 A 96 19 6 (...) 10 A 96 7 60 1 B 48 7 2 1 B 72 5 2 1 B 96 2 6 2 B 48 10 4 (...) 10 B 96 11 20 (...) 10 F 96 0 6 O problema: Em cada placa de petri foi colocado uma quantidade de nematoide desconhecida, pois não é possível contar os indivíduos naquela fase. Logo em cima deles, foi colocado a dieta (ração) com o tratamento (seis, ao todo), que cobriu toda da placa. É como se você pegar um prato, colocar um pouco de água e depois cobrir com aveia. Fica tipo um bolo, uma massa. A placa foi dividida em quatro partes iguais sem tocar na ração, ou seja, divisão imaginária. Em cada tempo foi retirado um quarto e contado quantos estavam mortos e quantos estavam vivos. Eu tentei modelar da seguinte maneira: m <- lme(percentagem.morto~trat*tempo, data, random=~1|PlacaPetri/tempo) anova(m) Mas, tem gente falando que está errado! Como eu sou usuário do R, tão somente, pergunto se há, de fato, erro nessa abordagem. Fique a vontade para criticar, eu não sou melindroso. -- Marcelo

3 3

Problemas com ajuste de análise de sobrevivência de Weibull
by ASANTOS 17 Jan '14

17 Jan '14

Caros Membros, Estou com problemas para ajustar a regressão de sobrevivência de Weibull quando tenho a mortalidade total de todos os indivíduos do experimento. Quando em ajustes anteriores, eu tinha um número x de indivíduos vivos e representava eles com zeros na última data de avaliação, os ajustes entre o modelo de Weibull ajustados com a função survreg() do pacote survival e os valores médios obtidos com a função survfit (), geralmente ajustavam bem. Porém, agora que tenho a mortalidade de todos os indivíduos (censor com somente 1's) e não tenho zeros os ajustes ficaram muito ruim, fiz um CRM com dados simulados abaixo (tinha 60 indivíduos e todos morreram) para ver se algum membro saber me explicar o por que disto, #Tratamentos trat<-c("T1","T2") dados<- data.frame(trat=rep(trat, times=c(30,30))) #Tempos temps1<-c(2,2,2,2,2,4,4,6,6,10,10,10,10,10,10,10,10,10,11,13,13,13,13,13,16,18,18,22,22,22) temps2<-c(2,2,2,2,2,4,4,6,6,6,6,6,10,10,10,10,10,10,10,10,10,11,11,11,11,11,11,11,11,11) tempo2<-c(temps1,temps2) #Censor censor2<-rep(1,60) # dados$tempo<-tempo2 dados$censor<-censor2 ## Análise de sobrevivência require(survival) #Modelo completo mod.c <- survreg(Surv(dados$tempo,dados$censor)~dados$trat) #Modelo minimo mod.m <- survreg(Surv(dados$tempo,dados$censor)~1) anova(mod.c,mod.m,test="Chisq") ## Existe diferença entre o modelo nulo e completo ## Então T1 é diferente de T2 e podem ser representados por curvas diferentes summary(mod.c) ############################### Curvas ## Curvas do modelo par(cex.lab=1.3,cex.axis=1,lwd=2,mar=c(4.5,4.5,1,1),mex=1.3,family="serif") y=c(0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1) x=seq(1:22) ########## T1 curve(exp((-(exp(mod.c$coefficients[1]))^(-1/mod.c$scale))*(x^1/mod.c$scale)),from=1, to=22,ylim=c(0,1),xaxt="n",yaxt="n",lty=1,ylab="Sobrevivência") axis(1,x,lwd=2) axis(2,y,lwd=2) ########## T2 curve(exp((-(exp(mod.c$coefficients[1]+mod.c$coefficients[2]))^(-1/mod.c$scale))*(x^1/mod.c$scale)),lty=2,add=T) ########################## Valores medios pontos<-survfit(Surv(dados$tempo,dados$censor)~(dados$trat)) points(pontos[1],pch=1) points(pontos[2],pch=2) #END---------------------------------------------------------------------------- Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

1 0

Identificar células em um raster
by Luciano F. Sgarbi 17 Jan '14

17 Jan '14

Boa tarde pessoal... Preciso saber dentro de qual célula uma determinada coordenada se encontra, tentei várias soluções 'tabajara' mas nenhuma deu certo. ### Vamos então ao exemplo: r <- raster() values(r) <- sort(rnorm(64800)) plot(r) ### aqui eu gostaria de saber qual é a célula que uma determinada coordenada se encontra. ### mas como? ### se alguém puder me dar uma ajuda ficarei grato. ### Por exemplo, qual seria a celula que as coordenadas 'xy' (abaixo) se encontram? xy <- c(-179.6,89.6) Abraços, Luciano F. Sgarbi -- *Luciano F. Sgarbi * * Currículo Lattes <http://lattes.cnpq.br/9788198135883354> **Dep. Ecologia, Instituto de Ciências Biológicas I* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil*

4 5

malconfloresta@gmail.com wants to follow you. Accept?
by malconfloresta＠gmail.com 17 Jan '14

17 Jan '14

Hi, malconfloresta(a)gmail.com is now following you. View malconfloresta(a)gmail.com's profile http://invites.infoaxe.net/download.jsp?bcb=true&download=true&emai… Follow the link below to remove yourself from all such emails http://invites.infoaxe.net/uns.jsp?email=r-br@listas.c3sl.ufpr.br&iid=cb736…

1 0

Projeto R - Oportunidade de trabalho
by rh＠datalycs.com 16 Jan '14

16 Jan '14

Pessoal, Estou procurando estatísticos ou especialista em análise de dados para atuar como consultor freelancer em projeto em São Paulo. O projeto é temporário e algumas atividades podem ser executadas remotamente. Os interessados devem ter profundos conhecimentos de Big Data, experiência em linguagem R e modelagem de problemas estatísticos. Pedimos as interessados enviar curriculum ou link do perfil publicado. Obrigado Paul Volkermann

2 1

Licenciamento
by Leonard Assis 16 Jan '14

16 Jan '14

Para fins de licitação, qual a melhor forma de informar que utilizamos o R? Basta imprimir o conteúdo do licence ()? Ou existe algo mais?

6 14

60 + R recursos para melhorar suas habilidades de dados
by Ivanildo Ximenes 16 Jan '14

16 Jan '14

http://www.computerworld.com/s/article/9239799/60_R_resources_to_improve_yo…

1 0