- R-br - DINF UFPR

Como colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R
by Ederson Civardi 30 Jan '14

30 Jan '14

Obrigado Fernando Toledo por responder. Porém o comando que me passou ainda não funciona. Da erro. Erro: símbolo inesperado in: "persp(xr,yr,z, main = "Banco", theta = 30, phi = 15, expand = 0.7, col = "darkorchid4" , xlab="Temperatura (ºC)", ylab="PMF (h)", zlab= expression ("Lesão (*mm^2*)",scale=T,ticktype = "detai persp" -- *_________________________________________* *Ederson A. Civardi* *Doutorando Produção Vegetal - UFG - PPGA - Goiânia* Cel. (Claro): (62) 9237-7817 ________________________________________________

3 2

Como colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R
by Ederson Civardi 30 Jan '14

30 Jan '14

Bom dia! Gostaria de uma ajuda básica para a formatação de gráficos no R. Preciso colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R, porém não estou conseguindo fazer isso. Tentei vários comandos, porém sempre aparece como por ex. (M^2). O comandos para o gráfico em 3D é: persp(xr,yr,z, main = "Banco", theta = 30, phi = 15, expand = 0.7, col = "darkorchid4" , xlab="Temperatura (ºC)", ylab="PMF (h)", zlab= expression ("Lesão (mm^2)"),scale=T,ticktype = "detailed") Outra dificuldade é inserir algumas linhas nos eixos x, y, z no gráfico que correspondam as escalas de cada eixo, para facilitar a visualização para que for ver. Alguém pode me ajudar? Desde já agradeço. -- *_________________________________________* *Ederson A. Civardi* *Doutorando Produção Vegetal - UFG - PPGA - Goiânia* Cel. (Claro): (62) 9237-7817 ________________________________________________

2 1

grafico de barras
by Natalia Martins 28 Jan '14

28 Jan '14

Boa tarde prezados, estou tentando fazer um gráfico de barras considerando 3 fatores, como p gráfico que o comando "bargraph.CI" permite fazer com dois fatores. Mas, não estou conseguindo. Alguém já fez este gráfico ou poderia me ajudar? Obrigada -- Natália Martins

3 8

Fábio, let's connect on LinkedIn
by Tiago Oliveira via LinkedIn 28 Jan '14

28 Jan '14

LinkedIn ------------ Tiago Oliveira requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Tiago Oliveira http://www.linkedin.com/e/-33g73l-hqz0ho8v-60/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Tiago Oliveira http://www.linkedin.com/e/-33g73l-hqz0ho8v-60/rso/295376692/x7Tr/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hqz0ho8v-60/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Re: [R-br] Dúvida - elipse de confiança
by Marcos Vital 27 Jan '14

27 Jan '14

Olá, Mauro Não sei se resolve no seu caso, mas dê uma olhada na função ordiellipse() do pacote vegan. Ela funciona bem para as ordenações feitas no próprio pacote, mas não sei se é compatível com outras ordenações, como no seu caso. Mas uma possibilidade é fazer a pca pelo vegan também, usando a função rda() com o argumento scale=T. Abraços Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

3 2

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