Dezembro 2011 - R-br

Re: [R-br] Ajuda para combinar Levels
by Ivan Bezerra Allaman 21 Dez '11

21 Dez '11

Segundo a análise combinatória, se a ordem dos fatores não importa, então você tem: combn(3,2) (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

Ajuda para combinar Levels
by Emmanuel Arnhold 21 Dez '11

21 Dez '11

Gostaria de concatenar combinações dos Levels de x (abaixo) com o character y (abaixo) para obter o “resultado” que apresento no final. Qual o melhor caminho para operacionalizar a obtenção do “resultado”? x=factor(c(1:3)) x y="vs" y resultado=c("1 vs 2","1 vs 3","2 vs 3") resultado

1 0

Transformação Box-Cox
by Alexandre Santos 21 Dez '11

21 Dez '11

Walmes, Descupe respondi meu email automaticamente e o envio da minha resposta foi para seu email pessoal ao invés de para o grupo. Continuando: Apesar da minha variável ser contínua, trata-se do peso de fezes de lagarta no tempo, sendo que para avaliarmos a velocidade com a qual uma determinada espécie de lagarta esta atacando um cultivo, colocamos vários quadrados de pano de 1m quadrado abaixo de uma área bastante atacada, porém nem todos os pontos amostrais possuem fezes em todas as datas de avaliação, apesar de no início do experimento (tempo1) todas as amostras possuem fezes. Meu problema é que se descartar os zeros produzo um desbalanceamento do meu DBC original e se usa-los vou ferir os pressupostos da distribuição normal. Por isso tentei a transformação Box-Cox, estava buscando uma abordagem inflacionadas com zeros, mais minha varíavel não tem distribuição de erros nem de Poisson, Binomial, Beta-binomial ... Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

1 0

Re: [R-br] Colocar legenda em histograma com curva de distribuição de weibull
by Ivan Bezerra Allaman 21 Dez '11

21 Dez '11

Utilize a library fdth para estatísticas descritivas. É uma excelente biblioteca. Baixe o pdf do library para maiores detalhes. library(fdth) x <- rnorm(1e5, mean=5, sd=1) plot(fdt(x, k=100), type='d', col=heat.colors(100)) curve(dnorm(x, mean=5, sd=1), col='darkgreen', add=TRUE, lwd=2) legend(par('usr')[1],par('usr')[4]-0.1,'Valores observados',pch=1,bty='n',col=heat.colors(100)) legend(par('usr')[1],par('usr')[4]-0.15,'Valores estimados',lty=1,bty='n',col='darkgreen') Como já exposto por Lucas, consulte o help da função legend para maiores detalhes! (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}

2 1

Colocar legenda em histograma com curva de distribuição de weibull
by Rodrigo Muniz 21 Dez '11

21 Dez '11

Caros amigo, Estou com os seguintes comandos: hist(Vento[[1]], seq(0, 18, 0.5), prob=TRUE, xlab="Velocidade (m/s)", ylab="Densidade", main="50 m") curve(dweibull(x, shape=2.223, scale=5.534), 0,18,add=T, col='red') Estes comando um gráfico com o histograma e uma linha representando a distribuição de weibull sobre o histograma, gostaria de saber como eu faço para adicionar uma legenda neste gráfico, representando os valores observados (histograma) e os valos estimado (Weibull). Desde já agradeço. 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)

2 2

Forçar a aparecer grau máximo em um modelo linear
by Henrique Ewbank 21 Dez '11

21 Dez '11

Bom dia a todos, Possuo dados gerados aleatoriamente e gostaria de modelá-los linearmente com diversos graus. A intenção é medir qual grau apresenta a melhor modelagem. Não trabalharei com previsão, mas sim buscarei as raízes reais desses polinômios. O problema ocorre quando é feita a modelagem para um polinômio de alto grau. Se tentar modelar com grau 20, por exemplo, não necessariamente a resposta terá grau 20, mas poderá apresentar qualquer grau inferior, como 19, por exemplo. É possível fixar que o maior grau exigido conste na resposta, ao modelar dados linearmente? Segue abaixo um exemplo de código. --- a<-NULL LT<-round(urtriang(10000,2.5,10,12.5)) DLT<-matrix(nrow=10000) for (i in 1:10000){ DLT[i]=sum(urtriang(LT[i],2.5,10,12.5)) } H<-hist(DLT[DLT>0],breaks=100,plot=FALSE) Dens<-H$density Int<-H$mids for (degreeTest in 2:20){ model<-lm(Dens~poly(Int,degreeTest,raw=TRUE)) Coeff1<-data.frame(model$coefficients[]) for (m in 1:nrow(Coeff1)) a <- c(a,Coeff1[m,1]) #end for m a <- c(a,rep(NA,21-m)) } b<-array(a,c(21,19),c("Coef","Grau")) ---- Os coeficientes de cada polinômio estão em *b*, sendo b[n,] os coeficientes de grau n-1. Exemplo: b[21,] = coeficientes de grau 20. Em *b* será possível observar que nem todos os graus (colunas) apresentarão coeficiente para seu maior expoente desejado. Obrigado, Henrique Ewbank Centro de Estudos em Logística COPPEAD

3 2

Dúvida no pacote ltm
by Felipe Buchbinder 21 Dez '11

21 Dez '11

Caros r-istas, Alguém utiliza o pacote ltm e pode me ajudar? Eu tenho usado a função grm para analisar respostas de 150 indivíduos a 4 questões e tenho recebido a seguinte mensagem de erro: Error: subscript out of bounds In addition: Warning messages: 1: glm.fit: algorithm did not converge 2: glm.fit: fitted probabilities numerically 0 or 1 occurred O que isso me diz? Como resolvo este problema? Obrigado, Felipe

3 4

Cartão de Natal No R [off-topic]
by Augusto Ribas 20 Dez '11

20 Dez '11

Ha, segue um cartão de natal feito com o R. Eu vi no Blog da revista Oikos: http://oikosjournal.wordpress.com/2011/12/19/merry-christmas/ Mas o código esta neste blog aqui: http://evol-eco.blogspot.com/2011/12/holiday-card.html Feliz natal a todos :) ###holiday caRd #install "rgl" package if necessary: install.packages("rgl") #load library library(rgl) ##Data r=1.0 h=5.0 m=h/r x=seq(0:99) y=x X=x*cos(y) Y=x*sin(y) Z=-1*(m*x) Z1=Z+rnorm(length(Z),Z,5) sub<-matrix(subset(c(X,Y,Z),Z==max(Z)),ncol=3) trunk<-matrix(c(0.54,0.84,-200,0.54,0.84,-650),byrow=TRUE,ncol=3) sno<-c(0,rep(15,9)) ##Plot card open3d(FOV=1) par3d(windowRect=c(100,100,600,600)) bg3d("aquamarine2") plot3d(X,Y,Z,col="darkgreen",type="l",box=FALSE,axes=FALSE,lwd=10,zlim=c(min(Z),0),xlab="",ylab="",zlab="",top=TRUE,border=TRUE) plot3d(X,Y,Z,col=c("yellow","red"),add=TRUE,type="s",radius=3) points3d(sub,col="orange",size=10) lines3d(trunk,col="brown",lwd=20) text3d(matrix(c(0.54,-100,300),ncol=3),text="Happy Holidays",font=5,cex=1.9,color="darkred",adj=c(0.5,1)) text3d(matrix(c(0.54,0.84,-650),ncol=3),text="and",font=5,cex=1.5,color="darkred",adj=c(0.5,1)) text3d(matrix(c(0.54,75,-780),ncol=3),text="a Happy New Year",font=5,cex=1.9,color="darkred",adj=c(0.5,1)) mat<-par3d("scale") ##Play card for(i in 1:10){ par3d(scale=mat,ignoreExtent=TRUE) points3d(matrix(c(sample(-160:160,sno[i]),sample(-160:160,sno[i]),sample(-500:150,sno[i])),ncol=3),col="white",size=3) play3d(spin3d(axis=c(0,0,1)),duration=1) if(i==10){points3d(sub,col="darkorange",size=17) plot3d(X,Y,Z,col=c("yellow","red"),type="s",radius=4,add=TRUE) } } text3d(matrix(c(0,140,-900),ncol=3),text="From the EEB & Flow",font=1,cex=1.0,color="darkgreen",adj=c(0.5,1)) ##The rgl device window allows you to use your mouse to manipulate the plot in 3-dimensions, so feel free to do so before closing it. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 1

(sem assunto)
by Edson Lira 20 Dez '11

20 Dez '11

Caros amigos, como faço para elimnar os NA´s? Quando li o banco usei a rotina bd[bd==""]<-NA. Não resolveu. > inap$ret_dias<-difftime(inap$dt_ret1, inap$dt_doa1, units = c("days")) > inap$ret_dias<-as.numeric(inap$ret_dias) > tapply(inap$ret_dias,inap$dt_doa_ano,summary) $`2010` Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 1 30 125 167 268 703 10318 $`2011` Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 1 4 35 75 107 23900 11468 []'s Edson Lira Estatístico Manaus-Amazonas

5 4

Função que transforma decimal em binario e vice-versa
by Claudio Carneiro 20 Dez '11

20 Dez '11

Olá todos, Estou a um tempo à procura, porém sem sucesso, por uma função, que já esteja disponivel no R, que transforme número decimal em binário, e uma que tranforme binário em decimal. Gostaria se vocês tem conhecimento de alguma? Grato desde já, Claudio Roberto Carneiro

3 2