- R-br - DINF UFPR

correlação de Spearman e respectiva linha no grafico
by Pedro Emmanuel Alvarenga Americano do Brasil 28 Nov '11

28 Nov '11

Amigos de R, Andei procurando RSiteSearch funçoes sobre a correlação de spearman. Eu achei diversas funçoes que fazem o rho de spearman e algumas extensões e correções, mas não consegui achar nenhuma que colocasse a linha da correlação em um grafico, a semelhança da reta na correlação de pearson. Alguem conhece alguma função ou tem algum script que faz isso? v1 v2 1 26 0.000 2 73 0.000 3 76 143.072 4 72 172.587 5 59 123.325 6 78 142.367 7 73 190.804 8 77 14.497 9 63 120.499 10 72 165.570 11 73 152.220 12 75 170.482 13 73 137.436 14 71 169.079 15 72 170.482 16 65 160.654 17 76 0.000 18 28 0.000 19 61 107.763 20 66 8.643 21 68 133.911 22 65 113.427 23 27 0.000 24 71 104.929 25 73 116.257 26 68 0.000 27 75 142.367 28 63 0.000 29 77 7.168 30 66 127.561 31 66 0.000 32 60 144.480 33 74 136.731 34 65 32.527 35 74 145.184 36 65 112.012 37 67 145.184 38 81 103.511 39 82 3.448 40 72 130.384 41 66 244.553 42 73 0.000 43 65 0.000 44 57 0.000 45 63 105.637 46 60 98.546 47 76 137.436 48 68 1.558 49 26 0.000 50 72 0.000 51 64 30.378 52 76 0.000 53 25 3.072 54 58 0.000 55 28 0.000 56 69 72.592 57 26 0.000 58 27 2.318 59 26 0.790 60 63 0.000 61 71 121.912 62 60 41.808 63 33 0.000 64 31 0.000 65 71 0.000 66 24 0.000 67 69 0.000 68 69 22.468 69 35 0.000 70 65 33.243 71 78 7.906 72 31 0.000 73 74 0.000 74 71 4.943 75 65 38.244 76 34 3.448 77 19 11.944 78 35 19.577 79 26 0.000 80 35 197.415 81 13 0.000 82 35 0.000 83 35 0.000 84 28 0.000 85 17 30.378 86 37 0.000 library('Hmisc') spearman2(v1 ~ v2, data=z, p=2) Spearman rho^2 Response variable:v1 rho2 F df1 df2 P Adjusted rho2 n v2 0.181 9.16 2 83 0.0003 0.161 86 Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

4 4

Inversão de matriz diagonal gigante
by Marcelo Cardoso mello 28 Nov '11

28 Nov '11

Olá pessoal. Gostaria de inverter uma matriz diagonal com algum processo que use menos memória que o do exemplo abaixo. (tenho uma matriz enorme e não esta invertendo via solve()). d<-c(1,1,0.5,0.75,0.5,0.5) d<-diag(d) i<-solve(d) i diag(i) # resultado que eu quero, mas de uma forma a gastar menos memória para diagonais gigantes!!!!

3 2

Versão R-2.14.0 no Ubunto
by Marcelo Cardoso mello 27 Nov '11

27 Nov '11

Pessoal, executei no terminal: $ sudo apt-get install r-base $ sudo apt-get install r-recommendad $ sudo apt-get etc etc que me foi indicado por aqui Tudo ok. Mas o R instalado é da versão 2.10 Tudo bem até aí. O problema é que não instala pacotes recentes. usei: instal.packages("vegan") tudo ok ... instalou mas não instala o lme4 por exemplo. O problema deve ser da versão do R que tenho? Como faço para instalar a versão R-2.14.0 no Ubunto? De forma rápida e simples. Obrigado por qualquer dica.

6 12

Re: [R-br] valor maximo e minimo
by Ivan Bezerra Allaman 27 Nov '11

27 Nov '11

Gilson, Busque por modelos que podem lhe dar respostas mais interessantes. Veja, A.Det=c(-0.595750,0.050250,0.270000,0.711500,1.608750,3.353500,4.937652,4.573750,4.955000) luz=c(0,25,50,75,100,500,1000,1500,2000) ajuste <- nls((A.Det~(a + 0.595750) * (1 - exp (-b * luz/(a-0.595750)))+0.595750), start=list(a=500, b=5)) plot(luz,A.Det) curve(3.9669 + 0.595750 * (1 - exp (-0.0061 * x/(3.9669-0.595750)))+0.595750,add=TRUE,col='blue') curve(SSasympOff(x,4.8292,-5.7353,23.452),add=TRUE,col='red') legend('bottomright', xjust=1,bty="n", c('Alternativo','Gilson'), lwd=2,col=c('red','blue')) Sobre este modelo 'alternativo', consulte "Mixed Effects Models in S and S-Plus, Pinheiro & Bates (2000)" OBS: Estes modelos não tem ponto de inflexão, pois quanto x tende ao infinito, y tende a zero, veja: D(expression(theta1*(1-exp(-exp(theta2)*(x-theta3)))),'x') curve(4.8292 * (exp(-exp(-5.7353) * (x - 23.452)) * exp(-5.7353)),to=0,from=2000,ylab='y') (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

valor maximo e minimo
by Gilson Sanchez 26 Nov '11

26 Nov '11

Pessoal mais uma vez solicitando ajuda de vc's O seguinte estou fazendo ajuste de curva com uma exponencial, então quero saber se da pra saber o valor do ponto maximo e minimo dos valores ajustados tanto no eixo X quanto no eixo Y? espe1 A.Det=c(-0.595750,0.050250,0.270000,0.711500,1.608750,3.353500,4.937652,4.573750,4.955000) luz=c(0,25,50,75,100,500,1000,1500,2000) ajuste = nls((A.Det~(a + 0.595750) * (1 - exp (-b * luz/(a-0.595750)))+ 0.595750), start=list(a=500, b=5)) summary(ajuste) obrigado! -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

2 3

Distribuição densidade de probabilidade de Weibull
by Rodrigo Muniz 25 Nov '11

25 Nov '11

Olá Pessoal, Sou usuário recente do R e como tal ainda tenho muitas dúvidas básicas na utilização do programa. Tenho uma série de dados de velocidade de vento na qual quero ajustar os parâmetros c e K da distribuição de Wuibull, no entanto quando uso as funções abaixo, da um erro que x deve ser um vetor numérico, no entanto se eu colocar um valor numérico ele gera uma série de números pseudoaleatórios, O que eu quero abrir os meus dados dentro da função. Se alguém puder me ajudar eu agradeço. Vento50m <- dweibull(Vento, shape=2, scale = 1) fitdistr(Vento50m, densfun="weibull") ATT. -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 6

Ajuda para realizar post-hoc em Anova de Medidas repetidas.
by Augusto Ribas 25 Nov '11

25 Nov '11

Ola, pessoal eu estou apanhando aqui para conseguir realizar uma analise de contrastes pra poder dizer que os resultados do Dia 1 são iguais do Dia 2, ai so no Dia 3 começa a aumentar. Eu vi um exemplo usando a função gls do pacote nlme, mas não consegui reproduzir pra conseguir verificar o que falo acima. Se alguém puder indicar um tutorial ou algum texto não muito complicado sobre como avaliar contrastes para anova de medidas repetidas, muito me ajudaria. Segue um exemplo dos dados: #Criandos dados set.seed(1) Resposta<-c(rnorm(20,2),rnorm(10,5),rnorm(20,5),rnorm(10,8)) Perigo<-rep(c("Baixo","Alto"),each=30) Id<-c(rep(1:10,3),rep(11:20,3)) Dia<-rep(rep(1:3,each=10),2) dados<-data.frame(cbind(Perigo,Dia,Id,Resposta)) dados$Resposta<-as.numeric(levels(dados$Resposta))[dados$Resposta] #Olhando os dados com um grafico str(dados) library(lattice) xyplot(Resposta~Dia|Perigo,data=dados, pch=16, col="black", cex=1.3, strip=strip.custom(bg="white"),type=c("p","a"),scale=list(tck=c(1,0),alternating=1)) #Anova de medidas repetidas modelo<-aov(Resposta~Perigo*Dia+Error(Id),data=dados) summary(modelo) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 5

Combinações de padrão específico
by Leonard de Assis 25 Nov '11

25 Nov '11

Estou com a seguinte situação: a <- letters[1:4] expand.grid(a,a) <- não serve pra mim, pois retorna os pares (a,a),(b,b), (c,c), (d,d) como combinações t(combn(a,2)) <- não serve pois retorna (a,b) mas não retorna (b,a) Preciso de algo q retorne um conjunto de pares como o expand.grid faz, mas excluindo as combinações repetidas (a com a, b com b, etc) Atualmente faço assim: (a <- letters[1:4]) (b <- expand.grid(a1=a,a2=a)) (final <- subset(b, a1 != a2)) alguma idéia melhor? Comandos entre parentesis para ilustrar o resultado -- []s Leonard de Assis assis<dot> leonard<at> gmail<dot> com

2 1

Cada alteração em um vetor sum () em outro
by ASANTOS 25 Nov '11

25 Nov '11

Boa tarde pessoal, Estou com uma dúvida de programação, onde tenho dois vetores de dados onde gostaria de obter a soma do primeiro, a cada alteração do segundo, tentei fazer: cont<-rpois(1000,4)####Cria os dados parc<-rep(1:length(freq),each=2)###Cria um novo vetor dados<-cbind(freq, parc) head(dados) ###Aqui gostaria de fazer a soma de cont a cada alteração de parc, usando um loop, for(i in 1:length(dados)) ##poderia ate especificar algo como dados$parc[i]>parc[i]+1 ## mas o problema seria fazer a sum() dentro de dados$freq Alguém teria alguma sugestão, Obrigado -- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304

2 5

Problemas numericos
by Carolina Pena 24 Nov '11

24 Nov '11

Boa tarde, Estou lendo um arquivo "txt" no R que contém o seguinte número: 698.4. Porém ao utilizar o seguinte comando: formatC(bd$x, digits = 40, format = 'e') [1] "6.9839999999999997724181577396507236699108e+02" O número 698.4 aparece conforme mostrado acima. Eu gostaria que o R entendesse como sendo: 6.984000000000000000000000000000000000e+02 Isso está me causando problemas ao realizar transformações pois preciso de um grau muito alto de precisão. Alguém sabe por que isso acontece e se é possível resolver? Agradeço desde já, Carolina

3 2