- R-br - DINF UFPR

[Dúvida] Configuração de computador
by Pedro Rafael 04 Jan '12

04 Jan '12

Pessoal um processador intel Xeon E5620 é melhor que o core i7 para trabalhar com contas e grandes base de dados? Quero trabalhar com o R e algumas outras linguagens de programação neste computador. Mais detalhes segue abaixo: Processador: - Quantidade: 1 - Fabricante: Intel - Modelo: Xeon Nehalem E5520 - Freqüência de operação: 2,4GHz~2,66GHz - Cache L2: 12MB (6MB 6MB), compartilhado Memória: - Capacidade: 8GB (2 x 4GB) - Expanção Máxima Opcional: 36GB (9 x 4GB) - Voltagem: 1.5V - Arquitetura: DDR3 - Freqüência: PC3-1333Mhz - Número de pinos: 240 Disco Rígido (HD): - Capacidade: 1TB - Interface: Serial ATA 3.0Gbps - Velocidade: 7200 RPM Unidade Óptica: - Leitor/Gravador DVD-RW Placa Mãe: - Marca: Intel - Modelo: S5500HCV Chipset: - Intel® 5500 chipset with Intel® I/O Controller Hub ICH10R Soluções Termicas: - Intel® Thermal Solutions STS100 series -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Adicionado índices do teste de médias no barplot
by ASANTOS 04 Jan '12

04 Jan '12

Boa tarde pessoal, Estou fazendo um barplot com a função barplot2 do pacote gplot, porém eu gostaria de colocar os índices do teste de média do Tukey sobre as barras do erro padrão, para isso fiz: ### require(gplot) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, plot.ci=T, ci.u=TAB+TAB2,ci.l=TAB-TAB2, ylab="Number of Bradysia difformis", legend.text=T) # Adicionando os meses mtext(1, at = colMeans(mp), text = rep(c("April", "May", "June", "July","August", "September","October"), 7),line = 1, cex = 0.75) Porém gostaria de saber se existe algum comando que consigo especificar um vetor contendo os índices e faze-los aparecer sobre as barras do erro padrão, ao invés de especificar índice por índice como text(). 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

5 5

Teste de Dunnett
by Adriano Carvalho Costa 04 Jan '12

04 Jan '12

Boa noite, Venho pedir ajuda de voces para perguntar se alguém conhece alguma rotina para realizar um teste de Dunnet. Pois estou com um experimento fatorial com um tratamento adicional e devo comparar o adicional versus os demais (um a um). Muito obrigado!

3 2

Calcular sobreposição de poligonos de um "SpatialPolygonsDataFrame" do pacote "sp"
by Augusto Ribas 03 Jan '12

03 Jan '12

Ola, primeiro feliz ano novo todos da lista. Eu gostaria de calcular a sobreposição entre poligonos, mas não consigo. Nas verdade eu tenho varios poligos gostaria de uma matriz de sobreposição, como a que demonstro ali no final do script. Nao achei como fazer isso ainda. Procurando por Overlap pelas listas de R não tenho tido muito sucesso então acho que o problema deva ser como estou procurando. Os poligonos são areas de vida de alguns bixinhos e a sobreposição da a informação se eles deixam outro individuos entrarem no territorio deles ou não, então é isso que estou tentando fazer e nao consigo. Segue um CMR pra vcs verem como são os dados e um exemplo de como gostaria do do resultado mais ou menos: #Pacote com dados de exemplo, ele necessita do pacote sp library(adehabitatHR) library(sp) #Dados de exemplo data(puechabonsp) rel<- puechabonsp$relocs #Estimando Area de Vida cp <- mcp(rel[,1]) #Olhando o os poligonos plot(cp) text(coordinates(cp)[,1],coordinates(cp)[,2],label=c("A","B","C","D")) #Como estao os dados class(cp) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" #O output que eu gostaria # A B C D #A 100 0 0 0 #B 0 100 0 0 #C 0 0 100 X #D 0 0 X 100 -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 2

qcc por estrato
by Edson Lira 02 Jan '12

02 Jan '12

Caros amigos, tenho essa estrutura de dados abaixo, aonde os exames são HIV1 HIV2. Exame Total Reativo 141 HIV 1 4270 14 142 HIV 2 4270 18 143 HIV 1 3558 7 144 HIV 2 3558 9 145 HIV 1 3530 4 146 HIV 2 3530 10 147 HIV 1 3241 6 148 HIV 2 3241 9 149 HIV 1 3388 2 150 HIV 2 3388 13 151 HIV 1 4249 1 152 HIV 2 4249 13 153 HIV 1 3763 11 154 HIV 2 3763 30 155 HIV 1 3821 6 156 HIV 2 3821 13 157 HIV 1 3972 0 158 HIV 2 3972 18 159 HIV 1 3227 3 160 HIV 2 3227 16 161 HIV 1 3576 4 162 HIV 2 3576 10 163 HIV 1 2941 4 164 HIV 2 2941 13 Com a rotina abaixo, require(qcc) qcc(h1_05$Reativo,type="np",nsigmas=2,sizes=h1_05$Total) Me dá o gráfico para os dois exames conjuntamente. Como fazer separando HIV1 e HIV2? Edson Lira Estatístico Manaus-Amazonas

2 2

Re: [R-br] qcc por estrato
by Ivan Bezerra Allaman 02 Jan '12

02 Jan '12

?subset (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

Dúvida: Erro em shapiro test
by Flávia Machado 31 Dez '11

31 Dez '11

Olá! Estou tendo problemas com o teste shapiro wilk. Inseri o arquivo txt (com pontos ao invés de vírgulas) dos meus dados com o comando "dados<-read.table(file.choose(), header=T). Os dados foram carregados e ao digitar dados todos os objetos apareceram normalmente, mas quando fui fazer o teste de shapiro apareceu " Erro em sort(x[complete.cases(x)]) : objeto 'FRUTO' não encontrado". Para o teste de shapiro, digitei: "shapiro.test(FRUTO)" Estou achando bem estranho porque a pouco tempo fiz o mesmo teste, da mesma forma, com outra série de dados e deu certo. Hoje usei estes mesmos dados apra testar e deu o mesmo erro. Vocês saberiam me dizer como solucionar esse problema? Obs.: sou bem iniciante no R, acho que deu para perceber. Atenciosamente, -- Flávia Machado

2 3

modelo multinomial logistico
by Luiz André Zardo 30 Dez '11

30 Dez '11

Prezados, estou usando a função multinom para ajustar um modelo logistico. Porém não sei como calcular os valores da estatistica da razão de verossimilhança, do pseudo-R2 e da estatistica de Wald. Alguém saberia me indicar como conseguir esses calculos??? ou indicar algum outro pacote que possa disponibilisar essas estatisticas?? Obrigado Luiz Zardo

2 1

aggregate
by Edson Lira 30 Dez '11

30 Dez '11

Caros amigos do R, na rotina abaixo gostaria de calcular além da média, o desvio padrão e o CV=(sd/mean). Como fazer estes cálculos? ## example with character variables and NAs testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = c(11,33,55,77,88,33,55,NA,44,55,77,99) ) by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if you want to treat NAs as a group fby1 <- factor(by1, exclude = "") fby2 <- factor(by2, exclude = "") aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") FELIZ 2012 a todos os amigos e familiares da lista. Edson Lira Estatístico Manaus-Amazonas

2 1

Retirar elemento de uma lista
by Rodrigo Coster 29 Dez '11

29 Dez '11

Povo, como faço para retirar apenas um elemento de uma lista? No caso, eu tenho essa lista: list(fruta='melancia',x=matrix(rpois(100,5),10),y=matrix(rpois(200,7),10)) Eu quero uma lista só com x e y (no caso a minha lista tem um numero variado de 'termos', não apenas 3). É o equivalente a índice negativo em vetor, mas não funciona em list :( []'s

3 3

Encontrar o maior valor por linha entre colunas
by Daniel Marcelino 29 Dez '11

29 Dez '11

Olá a todos. Eu tenho um pequeno problema que é justamente encontrar qual é o maior valor entre as colunar em cada linha. Como eu poderia obter uma nova coluna ou mesmo um objeto com o nome da coluna para cada valor encontrado. Alguém tem alguma sugestão? Daniel

7 10

[Dúvida] Contar NA em Subset
by Pedro Rafael 28 Dez '11

28 Dez '11

Pessoal como fazer um comando subset contar os elementos NA? Considere o data.frame abaixo: dados = data.frame(x = c(1,2,3,7,12), y = c(0,10,2,3,NA)) Com este data.frame chamado dados gostaria de filtrar com o subset os elementos de y que são diferente de 10. Para o meu exemplo NA também é diferente de 10 e gostaria de fosse contabilizado. Com o mando abaixo o resultado segue logo após. > subset(dados, y != 10) x y 1 1 0 3 3 2 4 7 3 Veja que o NA não foi contabilizado. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 6

Ajuda para plotar vários gráficos em uma única janela
by Maria Papa 27 Dez '11

27 Dez '11

Olá pessoal, boa tarde, Eu estou tentanto, de uma forma muito deselegante, plotar vários gráficos na mesma janela, mas somente o último gráfico é impresso. Por favor, alguém poderia me dizer o que está errado? par(mfrow = c(2,5)) for(i in length(temp)) { p = temp[,i] ht = hist(p, freq = F, col = "gray86", main = "", xlab = "Temperaturas", ylab = "Densidade") dn = density(p) lines(dn, col = 4, lwd = 2) } temp é uma matriz (9X900) com valores de temperatura. desde já agradeço a todos. Célia

2 1

Comparações Múltiplas em modelo linear com co-variável.
by Danilo Scorzoni Ré 26 Dez '11

26 Dez '11

Prezados, Estou realizando uma análise estatística através de um modelo de análise de covariância: y ~ sítio + bloco + tratamento + densidade + sitio:bloco + sitio:tratamento onde: sitio, bloco e tratametno são qualitativos densidade é quantitativo, foi inserido no modelo para reduzir a variação dos dados Os resultados da análise apontam efeitos significativos de tratamento e densidade. O desafio, agora, é fazer comparações múltiplas entre os tratamentos, no entanto, fixando-se densidades médias que variam entre os tratamentos, ou seja, descontando os efeitos de sítios e blocos que não são significativos. Tenho densidades médias para os quatro tratamentos: (400, 1200, 1150 e 900 para os tratamentos A, B, C e D), desta forma, utilizando os parâmetros estimados para cada tratamento + efeitos da densidade, gostaria de saber como eu poderia comparar as estimativas e verificar se os tratamentos diferem entre si. Eu poderia usar o tukey? Lembrando que a variável resposta é contínua... Caso não tenha sido claro, trago alguns resultados para exemplificar o problema. Desde já agradeço o empenho de todos e Feliz Natal! Att. -- Danilo Scorzoni Ré Engenheiro Florestal Mestrando em Ciência Florestal FCA / UNESP - Botucatu (14) 8180-2494

2 1

Transformação de variáveis.
by Juliana Freitas de Mello e Silva 25 Dez '11

25 Dez '11

Eu achei que esse tópico já existia, porém não o encontrei. Quero codificar observações de uma variável no meu banco de dados.

6 13

Função que faz o código voltar até determinada linha
by Felipe Buchbinder 23 Dez '11

23 Dez '11

Prezados r-istas, Existe em R alguma função que, caso uma determinada condição seja satisfeita, faça o código continuar a correr a partir de uma linha anterior? Por exemplo, algo do tipo: if(........) go to line 126 Sei que na programação da HP isso existe. Queria saber se existe também em R. Obrigado, Felipe Buchbinder

4 3

[Dúvida] Automatizar um modelo de previsão.
by Pedro Rafael 23 Dez '11

23 Dez '11

Pessoal a função auto.arima do pacote forecast "presta" realmente para ajustar modelos AR, MA, ARIMA, SARIMA? Existe algum pacote que faz de forma automática e eficiente do ponto de vista estatístico as análises de independências e de auto-correlações? Fico com o pé atrás de usar pacotes de prometem ajustar modelos auto-regressivos de forma automática. No meu caso tenho vários dados de mortalidade por doenças de aparelho circulatório, neoplasias, doenças do aparelho respiratório, etc. Todos os meses aqui na secretaria de saúde tenho mais informações e minha série vai aumentando. No meu computador tenho vários dados de milhares de variáveis distribuídos mensalmente e gostaria de criar algum programa para ajustar essas séries de forma automática. Não necessito do melhor modelo possível, apenas queria um bom modelo. É possível automatizar de tal forma a "dispensar" uma análise cuidadosa de um estatístico? O único problema é tempo para se ajustar série por série, ficar fazendo gráficos de auto correlações, verificando independências etc. Hoje possuo informações mensais desde 1999 a novembro de 2011, 143 pontos para cada uma das variáveis que preciso estimar. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 4

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

23 Dez '11

Boa tarde pessoal, Tenho uma variável resposta constituída de medidas de peso, porém alguns valores são representados por zeros, estava tentando fazer uma transformação de box-cox para tentar normalizar os dados, porém quando fiz: boxcox(peso ~ trat, data = dados, plotit = T) Erro em boxcox.default(peso ~ trat, data = dados, plotit = T) : response variable must be positive Apareceu o seguinte erro, então gostaria de saber se são os zeros que estão causando o erro? 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

4 6

duvida cat ou paste
by Vinicius Brito Rocha 22 Dez '11

22 Dez '11

Senhores, estou tentando transformar essas linhas em um único objeto para realizar uma contagem, ou seja, preciso que as 6 colunas tornem-se apenas 1. tentei o cat e o paste combinados com o apply, mas não tive sucesso dado original: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 4 5 30 33 41 52 [2,] 9 37 39 41 43 49 [3,] 10 11 29 30 36 47 [4,] 1 5 6 27 42 59 [5,] 1 2 6 16 19 46 objetivo: [1] "4 5 30 33 41 52" [2] "9 37 39 41 43 49" ... Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

4 4

Boas festas de fim de ano
by Rodrigo Muniz 22 Dez '11

22 Dez '11

Desejo a todos os companheiros de R, um feliz natal, feliz ano novo. Gostaria de agradeçer a todos pela ajuda dada, começei a alguns meses a ser usuário de R e pretendo continuar usando, mais não teria persistido se não fosse a ajuda de vcs. Mais uma vez agradeço a todos. 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)

1 0

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

22 Dez '11

Usando a imaginação! rep <- combn(3,2) results <- paste(rep,'vs',combn(3,2),sep='')[-c(1,3,5)] (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}

3 3

Pacote Experimental Designs
by Eric Ferreira 21 Dez '11

21 Dez '11

Prezados colegas, Com satisfação, informo que o pacote ExpDes (Experimental Designs), dedicado a análise de variância de experimentos de modelo fixo balanceados, após 2 anos de desenvolvimento e divulgação, se encontra no CRAN. Sugestões de melhorias e indicação de "bugs" são muito bem vindas! Atenciosamente, -- Prof Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

3 2

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 20 Dez '11

20 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

soma condicional
by Luiz André Zardo 19 Dez '11

19 Dez '11

Como posso obter uma soma condicional no R?

4 4

Fwd: Duvida para an[alise
by Adriano Carvalho Costa 19 Dez '11

19 Dez '11

---------- Mensagem encaminhada ---------- De: Adriano Carvalho Costa <acarvalhocosta(a)gmail.com> Data: 19 de dezembro de 2011 13:03 Assunto: Duvida para an[alise Para: r-br-request(a)listas.c3sl.ufpr.br Boa tarde, Estou com uma duvida para analisar este experimento. Alguem pode me ajudar. Penso da seguinte forma. Contrastes 1. Pectina e Amido vs Prop.Ca e CaCl 2. Pectina vs Amido 3. Pectina 1,5 vs Pectina 7,5 e Pectina 15 4. Pectina 7,5 vs Pectina 15 5. Amido 15 vs amido 45 e amido 75 5. Amido 45 vs amido 75 6. Prop. Ca vs CaCl 7. Prop. Ca 0 vc Prop. Ca 15 e Prop. Ca 30 8. Prop. Ca 15 vs Prop. Ca 30 9. CaCl 0 vc CaCl 15 e CaCl30 10.CaCl 15 vs CaCl Ca 30 tratamento Concentrações Pectina 1,5 Pectina 7,5 Pectina 15 amido 15 amido 45 amido 75 Prop. Ca 0 Prop. Ca 15 Prop. Ca 30 CaCl 0 CaCl 15 CaCl 30

2 5

soma acumulada
by Luiz André Zardo 19 Dez '11

19 Dez '11

Prezados colegas, Como faço para criar uma variável em um dataframe que seja soma acumulada de outra variável? i Xi Yi 1 a a 2 b a+b 3 c a+b+c Obrigado.

2 1

H0 é a hipótese nula?
by Pedro Emmanuel Alvarenga Americano do Brasil 19 Dez '11

19 Dez '11

Amigos de R, Um pequeno divertimento para desejar boas festas aos amigos virtuais que ajudam em horas oportunas! http://www.phdcomics.com/comics.php?n=1457 Abraço forte e que a força esteja com vocês, 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

1 0

Medidas repetidas na mesma unidade experimental!!
by Ivan Bezerra Allaman 19 Dez '11

19 Dez '11

Boa noite senhores! Há alguns dias estou estudando a metodologia de medidas repetidas em uma mesma unidade experimental ao longo do tempo, com o intuito de ajudar uma pesquisadora. Encontrei materiais, livros, artigos, etc que me dão segurança em trabalhar com tal medotologia. No entanto, há uma dúvida que me deixou num impasse, e preciso da ajuda dos colegas mais experientes que trabalham com estes tipos de dados. O objetivo da pesquisadora foi avaliar o uso de dois tratamentos sobre a redução de um tumor cancerígeno em cães do sexo feminino. Inúmeras variáveis (nominais, contínuas, etc.) foram avaliadas (secreção, leucócitos, etc.). As medidas foram tomadas em um mesmo animal ao longo de um período de 30 dias. Até aí tudo bem. Pois posso utilizar a abordagem de modelos mistos considerando o animal como de efeito aleatório, e no modelo incluir uma estrutura de correlação do tipo Arima, sarima, entre outros. Ainda estou avaliando qual melhor estrutura. Mas a dúvida crucial é que além de ter tomadas várias medidas num mesmo animal, este mesmo animal recebeu mais de um tratamento. Por exemplo: Do dia 1 ao dia 15 recebeu o tratamento A e do dia 16 ao dia 30 o tratamento B. Aí que estou encontrando problema, pois creio que há efeitos residuais do tratamento A sobre o B, confundindo os verdadeiros efeitos de cada tratamento. Uma saída, seria tratar os dois tratamentos como apenas um. Alguém já trabalhou com estrutura semelhante? Alguém poderia me recomendar uma literatura para que eu possa ganhar mais este desafio? Desde já, agradeço a atenção de todos! Abraços! (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}

8 14

Comportamento estranho nd rStudio com o lme4 + laticce
by Leonard de Assis 19 Dez '11

19 Dez '11

Bom dia, Estava rodando o código R presente no paper: http://www.jstatsoft.org/v20/i02/paper require(lme4) data("lq2002", package = "multilevel") wrk <- lq2002 for (i in 3:16) wrk[[i]] <- ordered(wrk[[i]]) for (i in 17:21) wrk[[i]] <- ordered(5 - wrk[[i]]) lql <- reshape(wrk, varying = list(names(lq2002)[3:21]), v.names = "fivelev", idvar = "subj", timevar = "item", drop = names(lq2002)[c(2, 22:27)], direction = "long") lql$itype <- with(lql, factor(ifelse(item < 12, "Leadership", ifelse(item < 15, "Task Sig.", "Hostility")))) for (i in c(1, 2, 4, 5)) lql[[i]] <- factor(lql[[i]]) lql$dichot <- factor(ifelse(lql$fivelev < 4, 0, 1)) str(lql) summary(lql) (fm1 <- lmer(dichot ~ 0 + itype + (1 | subj) + (1 | COMPID) + (1 | item), lql, binomial)) rr <- ranef(fm1, postVar = TRUE) str(rr) qq <- qqmath(rr) print(qq$subj) O executei no terminal e rodou perfeitamente (R 2.14 dev compilado dos fontes - linux ubuntu) Quando fui executar no RStudio (Mesmo computador, mesmo usuário), encontro a seguinte mens de erro: > print(qq$subj) Error in x$xscale.components(x$x.limits, at = x$x.scales$at, used.at = x$x.used.at, : could not find function "calculateAxisComponents" Isso é normal ou loucura do RStudio? O resultado esperado do comando é o gráfico na pág 12 do Paper (o da esquerda). Isso eu consegui no R (terminal) -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

1 0

[Dúvida] Combinações de strings de um vetor.
by Pedro Rafael 18 Dez '11

18 Dez '11

Pessoal existe algum objeto no R para realizar combinações 2 a 2, 3 a 3. etc? Exemplo: Considere o vetor: vetor <- c("a","b","c","d") Queria obter todas as combinações possíveis com os elementos do vetor, ou seja, uma saída semelhante a: [1] "ab" "ac" "ad" "bc" "bd", "dc" Até poderia pensar em programar isto do zero usando o R mas será que já existe um objeto para fazer isto pronto e eficiente? Obrigado à todos desde já. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 11

[Dúvida] Tratamento de Erros no R.
by Pedro Rafael 18 Dez '11

18 Dez '11

Pessoal como tratar erro no R? Veja a função abaixo: dado = "fa\xE7ile" cod = combn(iconvlist(),2) i <- NULL conver <- vector() for(i in 1:dim(cod)[2]){ iconv(dado, from = as.vector(cod[,i])[1], to = as.vector(cod[,i])[2]) } Fiz todas as combinações possíveis de formato de arquivo e gravei no objeto cod. Tenho um outro objeto chamado dado que está em uma codificação até então "estranha", a sting correta seria façile. Queria ver qual codificação correta para converter o objeto dado para o formato correto. Fiz um FOR que usa o comando ivonc para todas as combinações possíveis. Só que a condição para quando ocorre um primeiro erro, ou seja, quando um tipo de conversão não é adequado. Queria saber como desconsiderar estes erros e dizer que pode continuar após a ocorrência de um erro. Mesmo havendo outra forma de resolver meu problema de codificação gostaria de saber como tratar problemas desse tipo, por exemplo dizer que algo pode ser continuado após um erro ou parado quando um erro ocorre ou uma quantidade específica de erros ocorrerem. A outra pergunta diz respeito a como saber a codificação correta do objeto de forma mais eficiente possível. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 2

Usando Sweave e MikTex no Windows [Off-Topic]
by Augusto Ribas 17 Dez '11

17 Dez '11

Ola. Bem eu começei a usar o Sweave no R e é uma excelente ferramenta pra escrever. Mas eu estou apanhando aqui um pouco. Bem, eu estou criando um Arquivo .Rnw num diretorio que concentro todos os meus arquivos referentes a um determinado projeto. Ai uso a Função Sweave() nele para criar o arquivo .tex So que quando vou usar o arquivo .tex com o miktex ele não vai e pede um arquivo, o Sweave.sty. Este que esta em /share/texmf/tex/latex Ai eu tenho que colocar esse arquivo no diretorio senão não da certo compilar(o termo compilar aqui ta certo, ou como eu devo falar ao usar o latex num arquivo tex pra fazer um dvi ou pdf?) o arquivo tex num pdf. Dai a Minha duvida é: Tem como eu dizer no MikTex onde ta esse aquivo Sweave.sty pra ele usar direto o da instalação? E não precisar estar no diretorio que estou trabalhando? Bem é isso. Obrigado :) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 3

Função para retornar uma lista de combinações de colunas de uma matriz
by Emmanuel Arnhold 16 Dez '11

16 Dez '11

Gostaria de alguma dica para operacionalizar a CMR abaixo em uma função. O resultado final é uma lista com combinações das colunas de uma matriz. x=matrix(c(1,2,3,4,5,6,7,8,9), nc=3);x l1=c(x[,1], x[,2]) l2=c(x[,1], x[,3]) l3=c(x[,2], x[,3]) Lista=list(l1,l2,l3);Lista Obrigado. Emmanuel

3 4

Uso do ifelse
by Edson Lira 16 Dez '11

16 Dez '11

R-istas, estou com um problema no uso do ifelse (Elias, finalmente postei o CMR). Tenho dois problemas: 1. ao avaliar percentual de gordura em homens e mulheres, temos os seguintes parâmetros Para Homens 15|-20 = Limite normal 20|-25 = Moderadamente obeso 25|-30 = Excessivamente obeso 10|-14 = Gordura ideal 03 = Gordura essencial Para Mulheres 25|-30 = Limite normal 30|-35 = Moderadamente obeso 35|-40 = Excessivamente obeso 14|-18 = Gordura ideal 12 = Gordura essencial Como fazê-lo sem criar subcojunto de homens e mulheres usando o ifelse? No banco as variáveis são sexo e perc_g. 2. Foi medido a pressão sitólica e distólica de um grupo, tenho os seguintes parâmetros ( que me foram passados): sistólica diastólica categoria <130 <85 normal 130-139 85-89 normal limítrofe 140-159 90-99 estágio 1 160-179 100-109 estágio 2 >180 >110 estágio 3 >=210 >=120 estágio 4 <140 <90 hipert sist isolada Como fazê-lo usando o ifelse? No endereço abaixo postei uma parte do banco de dados e o scritp que tentei implementar, mas não conseGui rodá-lo. https://gist.github.com/1481221 OU git://gist.github.com/1481221.git []'s. Edson Lira Estatístico Manaus-Amazonas

5 8

ordenar banco de dados
by Luiz André Zardo 16 Dez '11

16 Dez '11

Prezados, Alguém poderia me indicar como fazer para obter ordenação de um banco de dados Por exemplo, Como ordenar um banco de dados por uma variável específica. Obrigado!!

5 5

escalas em gráfico xyplot (lattice)
by Gustavo Marcatti 15 Dez '11

15 Dez '11

Pessoal, estou fazendo um gráfico de altura de uma planta em função de doses de N com a função xyplot. As doses são: 0, 120, 240, 360, 480, 600. Gostaria que o eixo x do gráfico viesse com as doses originais e não em uma sequência fixa (0, 100, 200, 300, 400, 500, 600) como o padrão. Alguem poderia ajudar Gustavo Marcatti Eng Florestal UFV

5 7

Consumo de Memória
by Fernando Neto 15 Dez '11

15 Dez '11

Estou com um problema de mudança de variáveis (transformar de discreta em colunas binárias), mas o programa tá consumindo memoria demais. Estou no ubuntu 64 e o consumo é tanto que chega um momento que nem o mouse mexe mais. Existe algum jeito de solicitar ao R que chame o garbage collector? Muito grato, PS.: qnto isso, continuo aceitando ajudas no problema de transformacao das variáveis (ja tenho uma thread em meu nome com isso =T) -- ---------------------------------------------------------------- *Fernando Neto* Twitter: @fernandompneto Facebook: facebook.com/fernandompneto Tecnologia de Ponte http://tecnologiadeponte.blogspot.com ---------------------------------------------------------------- *fmpn2 @ CIn - UFPE* http://cin.ufpe.br/~fmpn2 - Engenharia da Computação - Turma 2009.2 - CIn, UFPE. - Monitor de Estatistica e Probabilidade Para Engenharia da Computacao<http://www.cin.ufpe.br/%7Eet586/> <http://cin.ufpe.br/%7Eif672cc> ---------------------------------------------------------------- Confidencialidade*: *A informação contida nesta mensagem de e-mail, incluindo quaisquer anexos, é confidencial e está reservada apenas à pessoa ou entidade para a qual foi endereçada. Se você não é o destinatário ou a pessoa responsável por encaminhar esta mensagem ao destinatário, você está, por meio desta, notificado que não deverá rever, retransmitir, imprimir, copiar, usar ou distribuir esta mensagem de e-mail ou quaisquer anexos. Caso você tenha recebido esta mensagem por engano, por favor, contate o remetente imediatamente e apague esta mensagem de seu computador ou de qualquer outro banco de dados. Muito obrigado. Confidentiality Notice: The information contained in this email message, including any attachment, is confidential and is intended only for the person or entity to which it is addressed. If you are neither the intended recipient nor the employee or agent responsible for delivering this message to the intended recipient, you are hereby notified that you may not review, retransmit, convert to hard copy, copy, use or distribute this email message or any attachments to it. If you have received this email in error, please contact the sender immediately and delete this message from any computer or other data bank. Thank you. --------

5 6

[Dúvida] Como usar o Crontab com o R.
by Pedro Rafael 15 Dez '11

15 Dez '11

Pessoal como chamar um script em R utilizando a cada uma hora, por exemplo, usando os programas um crontab no Linux. As características do meu sistema e da versão do R estão abaixo: R version 2.14.0 (2011-10-31) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C LC_TIME=pt_BR [4] LC_COLLATE=pt_BR.UTF-8 LC_MONETARY=pt_BR LC_MESSAGES=pt_BR.UTF-8 [7] LC_PAPER=C LC_NAME=C LC_ADDRESS=C [10] LC_TELEPHONE=C LC_MEASUREMENT=pt_BR LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_2.14.0 Eu consegui agendar tarefas para escrever um texto em um arquivo .txt. Não estou entendendo como chamar o R e executar um determinado script em uma data específica. Muito obrigado à todos desde já. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 6

Re: [R-br] split-plot
by Fabiana Schneck 15 Dez '11

15 Dez '11

Oi Walmes, São 10 blocos, 5 deles em remansos e outros cinco em corredeiras. Em cada bloco existem dois plots colocados lado a lado, sendo um com exclusão elétrica e outro sem exclusão (ou seja, controle). Dentro de cada um desses dois plots existem substratos de acrílico lisos e rugosos. Preciso avaliar o efeito do habitat (remansos x corredeiras), o efeito da exclusão (exclusão x controle) e o efeito do tipo de substrato (liso x rugoso) levando em conta o design em spli-plot. Ficou mais claro? Obrigada, Fabiana Em 13 de dezembro de 2011 17:38, <r-br-request(a)listas.c3sl.ufpr.br>escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. split-plot (Fabiana Schneck) > 2. Re: Transformar Variáveis discretas em Numeros Binários -> > Pre-processamento para Redes Neurais (Elias T. Krainski) > 3. Rotina order (Edson Lira) > 4. Re: Rotina order (Elias T. Krainski) > 5. Re: MQM (Elias T. Krainski) > 6. Re: Modificar valores de dataframe com base em outro > dataframe (Elias T. Krainski) > 7. Categorização de variáveis (Edson Lira) > 8. Re: Categorização de variáveis (Luis Iván Ortiz Valencia) > 9. Re: Categorização de variáveis (Cesar Rabak) > 10. Re: Categorização de variáveis (Elias T. Krainski) > 11. Re: Categorização de variáveis (Edson Lira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 13 Dec 2011 12:25:13 -0200 > From: Fabiana Schneck <fabiana.schneck(a)gmail.com> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] split-plot > Message-ID: > <CALjK-AcaG3AELQAOopMDZ1X7oFYOzNg4957pn3-3BG9J+t07oA(a)mail.gmail.com > > > Content-Type: text/plain; charset="windows-1252" > > Prezados colegas, > > tenho um experimento com o seguinte design split-plot: > > fator 'ambiente': 5 corredeiras e 5 remansos > em cada plot de 'ambiente' (num total de 10) há outro fator: 'peixes' com > dois tratamentos: exclusão e controle > ainda, dentro de cada sub-plot do fator 'peixes' existe o fator 'substrato' > com dois tratamentos: liso e rugoso. > > A variável resposta é quantidade de clorofila em cada unidade experimental > (cada substrato). > Um exemplo dos dados: > > > > clorofila > > ambiente > > peixes > > substrato > > bloco > > 0.114 > > 1 > > 1 > > 1 > > 1 > > 0.160 > > 1 > > 1 > > 2 > > 1 > > 0.126 > > 1 > > 2 > > 1 > > 1 > > 0.149 > > 1 > > 2 > > 2 > > 1 > > 0.068 > > 1 > > 1 > > 1 > > 2 > > 0.137 > > 1 > > 1 > > 2 > > 2 > > 0.080 > > 1 > > 2 > > 1 > > 2 > > 0.160 > > 1 > > 2 > > 2 > > 2 > .... > .... > > 0.080 > > 2 > > 1 > > 1 > > 10 > > 0.068 > > 2 > > 1 > > 2 > > 10 > > 0.114 > > 2 > > 2 > > 1 > > 10 > > 0.124 > > 2 > > 2 > > 2 > > 10 > > Considerei que os fatores 'peixes' e 'substrato' estão em blocos (10), > enquanto o fator 'ambiente' está acima dos blocos (5 blocos são no ambiente > corredeiras e 5 no ambiente remansos). > > Analise os dados utilizando dois modelos distintos: > > *aov (log(cloro)~(amb+peixes+substrato)^2+Error(bloco/peixes))* > * > * > > Error: bloco > > Df Sum Sq Mean Sq F value Pr(>F) > > amb 1 0.1162 0.11618 0.128 0.7297 > > Residuals 8 7.2599 0.90749 > > > > Error: bloco:peixes > > Df Sum Sq Mean Sq F value Pr(>F) > > peixes 1 2.96190 2.96190 14.6147 0.005069 ** > > amb:peixes 1 0.30121 0.30121 1.4862 0.257528 > > Residuals 8 1.62133 0.20267 > > --- > > Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 > > > > Error: Within > > Df Sum Sq Mean Sq F value Pr(>F) > > substrato 1 1.9924 1.99244 10.2620 0.005209 ** > > amb:substrato 1 0.0013 0.00134 0.0069 0.934840 > > peixes:substrato 1 0.5934 0.59335 3.0560 0.098473 . > > Residuals 17 3.3007 0.19416 > > > * > * > *lme(log(cloro)~(amb+peixes+**substrato)^2, random=~1|bloco/peixes))* > > Linear mixed-effects model fit by REML > Data: NULL > AIC BIC logLik > 85.58484 100.5499 -32.79242 > > Random effects: > Formula: ~1 | bloco > (Intercept) > StdDev: 0.4197679 > > Formula: ~1 | peixes %in% bloco > (Intercept) Residual > StdDev: 0.06522495 0.4406333 > > Fixed effects: log(cloro) ~ (amb + peixes + substrato)^2 > Value Std.Error DF t-value p-value > (Intercept) -2.4342471 0.2647064 17 -9.196027 0.0000 > amb2 -0.2929029 0.3611524 8 -0.811023 0.4408 > peixes2 0.6142687 0.2448449 8 2.508808 0.0364 > substrato2 0.6783940 0.2413448 17 2.810891 0.0120 > amb2:peixes2 0.3471056 0.2847218 8 1.219104 0.2575 > amb2:substrato2 0.0231237 0.2786810 17 0.082976 0.9348 > peixes2:substrato2 -0.4871762 0.2786810 17 -1.748150 0.0985 > Correlation: > (Intr) amb2 peixs2 sbstr2 amb2:p2 amb2:s2 > amb2 -0.682 > peixes2 -0.462 0.229 > substrato2 -0.456 0.223 0.329 > amb2:peixes2 0.269 -0.394 -0.581 0.000 > amb2:substrato2 0.263 -0.386 0.000 -0.577 0.000 > peixes2:substrato2 0.263 0.000 -0.569 -0.577 0.000 0.000 > > Standardized Within-Group Residuals: > Min Q1 Med Q3 Max > -2.42000225 -0.41789060 -0.02451742 0.44159262 1.87090777 > > Number of Observations: 40 > Number of Groups: > bloco peixes %in% bloco > 10 20 > > > > Tenho dúvidas se estou definindo de forma correta os modelos e também se > devo usar aov ou lme, uma vez que os resultados dos termos isolados diferem > entre as duas funções (porém os resultados são indênticos para as > interações). > > Desde já obrigada pela atenção, > > Fabiana > > -- > ----- > Fabiana Schneck > Programa de Pós-Graduação em Ecologia > Universidade Federal do Rio Grande do Sul > CP 15007 - CEP 91501-970 > Porto Alegre, RS, Brasil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/6ca4ee6e/atta… > > > > ------------------------------ > > Message: 2 > Date: Tue, 13 Dec 2011 06:30:03 -0800 (PST) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Transformar Variáveis discretas em Numeros > Binários -> Pre-processamento para Redes Neurais > Message-ID: > <1323786603.22667.YahooMailNeo(a)web114716.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Fernando, talvez eu não tenha entendido exatamente o que você quer. Veja > dois exemplos, considerando as colunas 2 e 3 de 'warpbreaks' (disponivel no > 'datasets', como o 'iris'). > > > ### Considerando que um 'factor' é um inteiro, apenas retiramos os > 'levels', e vc pode manipular > ### Se ha 3 ou mais classes, o resultado nao é binário, é inteiro de 1 ao > número de classes > sapply(warpbreaks[,2:3], unclass) > sapply(warpbreaks[,2:3], unclass)-1 > > > ### converte em colunas binárias, como dummies em qualquer modelo regressão > > ### se há mais de duas classes, é criado mais de uma coluna > model.matrix(~wool + tension, warpbreaks)[,-1] > > > Att. > > Elias T. Krainski > > > >________________________________ > > De: Fernando Neto <fernandoneto7(a)gmail.com> > >Para: r-br(a)listas.c3sl.ufpr.br > >Enviadas: Segunda-feira, 12 de Dezembro de 2011 21:32 > >Assunto: [R-br] Transformar Variáveis discretas em Numeros Binários -> > Pre-processamento para Redes Neurais > > > > > >Boa noite, R-anos! > >Por favor, estou precisando ajustar os arquivos de uma base de dados para > processamento de uma RNA. > >Estou usando a função decodeClassLabels(...) mas a operacao tá ficando > muito custosa. > >Por exemplo > >__ > >Masculino Programador > >Masculino Médico > >Feminino Programador > >__ > >tem que ser > >__ > >0 1 0 1 > >1 0 1 0 > >__ > > > > > >Só que eu trato cada variável de cada vez. E acho que a operacao de > concatenacao que eu uso: > >cbind(...) tá tornando a operacao custosa (além de que a memória passa de > 2GB). > > > > > >PS.: eu faço o seguinte: > > > > > >for (i in (VETOR COM AS COLUNAS A SEREM TRANSFORMADAS)){ > > colAux = decodeClassLabels(vector[,i]) > > vector = vector [,-i] > > vector = cbind(vector, colAux) > >} > > > > > >Alguém sugere algo? > > > > > >Muito grato!!! > >Fernando > > > > > >-- > >---------------------------------------------------------------- > > > >Fernando Neto > >Twitter: @fernandompneto > >Facebook: facebook.com/fernandompneto > > > >Tecnologia de Ponte > >http://tecnologiadeponte.blogspot.com > >---------------------------------------------------------------- > >fmpn2 @ CIn - UFPE > >http://cin.ufpe.br/~fmpn2 > > > >- Engenharia da Computação - Turma 2009.2 - CIn, UFPE. > >- Monitor de Estatistica e Probabilidade Para Engenharia da Computacao > > > > > >---------------------------------------------------------------- > >Confidencialidade*: *A informação contida nesta mensagem de e-mail, > >incluindo quaisquer anexos, é confidencial e está reservada apenas à > pessoa ou entidade para a qual foi endereçada. Se você não é o destinatário > ou a pessoa responsável por encaminhar esta mensagem ao destinatário, você > está, por meio desta, notificado que não deverá rever, retransmitir, > imprimir, copiar, usar ou distribuir esta mensagem de e-mail ou quaisquer > anexos. Caso você tenha recebido esta mensagem por engano, por favor, > contate o remetente imediatamente e apague esta mensagem de seu computador > ou de qualquer outro banco de dados. Muito obrigado. > >Confidentiality Notice: The information contained in this email > message, including any attachment, is confidential and is intended only for > the person or entity to which it is addressed. If you are neither the > intended recipient nor the employee or agent responsible for delivering > this message to the intended recipient, you are hereby notified that you > may not review, retransmit, convert to hard copy, copy, use or distribute > this email message or any attachments to it. If you have received this > email in error, please contact the sender immediately and delete this > message from any computer or other data bank. Thank you. > > > >-------- > > > >_______________________________________________ > >R-br mailing list > >R-br(a)listas.c3sl.ufpr.br > >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/19fb4b33/atta… > > > > ------------------------------ > > Message: 3 > Date: Tue, 13 Dec 2011 07:22:54 -0800 (PST) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Rotina order > Message-ID: > <1323789774.75434.YahooMailNeo(a)web160704.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia a todos! > > Com a rotina abaixo eu consigo ordenar a tabela. > > inap<-table(tip_doa$INAPTIDAO,tip_doa$mês) > inap2 = inap[order(rowSums(inap), decreasing=TRUE),] > ftable(addmargins(inap2, c(1,2), > FUN=list(Total = sum))) > > Que me dá a saida abaixo(as primeiras linhas) > > jan > fev mar abr mai jun jul ago set out nov Total > > ESTEVE EM ZONA DE MALARIA 465 > 443 387 362 343 347 343 385 550 397 409 4431 > HEMATOCRITO IRREGULAR 299 > 387 478 473 576 393 266 307 362 381 333 4255 > VARIAÇÃO DE PARCEIRO SEM PRESERVATIVO (MAIS DE 4 EM 12 MESES) 249 > 329 240 152 100 183 150 86 110 124 106 1829 > VARIAÇÃO DE PARCEIROS SEM PRESERVATIVO (MAIS DE 1 EM 6 MESES) 190 > 206 147 83 73 104 87 58 103 56 102 1209 > VARIAÇÃO DE PARCEIROS SEXUAIS SEM PRESERVATIVO 0 > 0 20 106 238 135 109 138 134 142 160 1182 > GRIPE 119 > 112 125 98 105 112 75 96 90 84 102 1118 > > Gostaria de acrescentar o ano e ordenar. Quando acrescento o ano me dá o > seguinte problema. > > > inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > > inap2 = inap[order(rowSums(inap),decreasing=TRUE),] > Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : > número incorreto de dimensões > > inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > > inap2 = inap[order(rowSums(inap),decreasing=TRUE),] > Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : > número incorreto de dimensões > > ftable(addmargins(inap2, c(1,2,3), > + FUN=list(Total = sum))) > Erro em newdimnames[[margin]] : índice fora de limites > > Qual o problema? Alguém pode me ajudar? > > []'s. > > > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/d6d0f94c/atta… > > > > ------------------------------ > > Message: 4 > Date: Tue, 13 Dec 2011 07:43:37 -0800 (PST) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Edson Lira > <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] Rotina order > Message-ID: > <1323791017.65689.YahooMailNeo(a)web114705.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Com > > inap <- table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > 'inap' deixa de ter duas dimensões. Portanto, rowSums() é inadequado. > > > Elias T. Krainski > > > >________________________________ > > De: Edson Lira <edinhoestat(a)yahoo.com.br> > >Para: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > >Enviadas: Terça-feira, 13 de Dezembro de 2011 13:22 > >Assunto: [R-br] Rotina order > > > > > >Bom dia a todos! > > > >Com a rotina abaixo eu consigo ordenar a tabela. > > > >inap<-table(tip_doa$INAPTIDAO,tip_doa$mês) > >inap2 = inap[order(rowSums(inap), decreasing=TRUE),] > >ftable(addmargins(inap2, c(1,2), > > FUN=list(Total = sum))) > > > >Que me dá a saida abaixo(as primeiras linhas) > > > > jan > fev mar abr mai jun jul ago set out nov Total > > > >ESTEVE EM ZONA DE MALARIA 465 > 443 387 362 343 347 343 385 550 397 409 4431 > >HEMATOCRITO IRREGULAR 299 > 387 478 473 576 393 266 307 362 381 333 4255 > >VARIAÇÃO DE PARCEIRO SEM PRESERVATIVO (MAIS DE 4 EM 12 MESES) 249 > 329 240 152 100 183 150 86 110 124 106 1829 > >VARIAÇÃO DE PARCEIROS SEM PRESERVATIVO (MAIS DE 1 EM 6 MESES) 190 > 206 147 83 73 104 87 58 103 56 102 1209 > >VARIAÇÃO DE PARCEIROS SEXUAIS SEM PRESERVATIVO 0 > 0 20 106 238 135 109 138 134 142 160 1182 > >GRIPE 119 > 112 125 98 105 112 75 96 90 84 102 1118 > > > >Gostaria de acrescentar o ano e ordenar. Quando acrescento o ano me dá o > seguinte problema. > > > >> > inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > >> inap2 = inap[order(rowSums(inap),decreasing=TRUE),] > >Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : > > número incorreto de dimensões > >> inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > >> inap2 = inap[order(rowSums(inap),decreasing=TRUE),] > >Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : > > número incorreto de dimensões > >> ftable(addmargins(inap2, c(1,2,3), > >+ FUN=list(Total = sum))) > >Erro em newdimnames[[margin]] : índice fora de limites > > > >Qual o problema? Alguém pode me ajudar? > > > > []'s. > > > > > >Edson Lira > >Estatístico > >Manaus-Amazonas > >_______________________________________________ > >R-br mailing list > >R-br(a)listas.c3sl.ufpr.br > >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/53dc517b/atta… > > > > ------------------------------ > > Message: 5 > Date: Tue, 13 Dec 2011 08:00:29 -0800 (PST) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] MQM > Message-ID: > <1323792029.83765.YahooMailNeo(a)web114716.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Se em > > A b = x > A tem dimensão 3x3, então vc quer encontrar 6 valore. Para isso, vc tem 3 > valores em b e 3 valores em x. Isso é possível? No usual, quando vc quer > saber os 3 valores de b e conhece A e x, vc tem bem mais valores conhecidos > que desconhecidos. > > > ### Exemplo > > ### considerando o contexto de LM > n <- 5 > X <- cbind(1, runif(5)) > b <- c(10, 5) > y <- rnorm(n, X%*%b, 0.5) > > > ### estima b > (XX <- crossprod(X)) > (Xy <- crossprod(X, y)) > (bhat <- solve(XX, Xy)) > > ### com bhat e Xy podemos obter XX? > ### considere que > ### XX %*% bhat = Xy > ### para compatibilidade de operações seja > ### XX %*% (bhat %*% bhat') = Xy %*% bhat' > > require(Matrix) > bb <- forceSymmetric(tcrossprod(bhat)) > solve(bb, tcrossprod(Xy, bhat)) ## nao é igual a XX > > pode ser que este 'sistema' possua mais de uma solução... > > > Att. > > Elias T. Krainski > > > >________________________________ > > De: FHRB Toledo <fernandohtoledo(a)gmail.com> > >Para: R-Br <r-br(a)listas.c3sl.ufpr.br> > >Enviadas: Domingo, 11 de Dezembro de 2011 12:07 > >Assunto: [R-br] MQM > > > > > >Senhores, > > > >Tenho uma dúvida que não necessáriamente se relaciona com o programa, mas > que em muito tem a ver com teoria estatística e algebra de matrizes, espero > que aqui ache pelo menos uma sugestão para a resposta ou uma recomendação > de literatura. > > > >O problema é dado o sistema de equações na forma matricial: Xb = y, pelo > MQM, dada a matriz X e o vetor y, obtenho o vetor de soluções por X^-1y, > respeitando todas as pressuposições. > > > >Se caso conheço ao invés de X e y, conheço os vetores b e y, existe > alguma forma de conhecer a matriz X, podendo apenas assegurar que ela é uma > matriz quadrada simétrica e positiva definida? > > > >att, > >FH > > > >_______________________________________________ > >R-br mailing list > >R-br(a)listas.c3sl.ufpr.br > >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/018859e7/atta… > > > > ------------------------------ > > Message: 6 > Date: Tue, 13 Dec 2011 08:06:38 -0800 (PST) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Modificar valores de dataframe com base em outro > dataframe > Message-ID: > <1323792398.50539.YahooMailNeo(a)web114717.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > d1 <- data.frame(snp1=c("A","A","A","A","T"), snp2=c("A","G","G","A","G"), > > snp3=c("C","T","C","C","T"), > snp4=c("T","G","G","T","T")) > > > d2 <- data.frame(snp=paste("snp",1:4,sep=""), anc=c("A","G","C","T")) > > sapply(1:ncol(d1), function(i) > 1-(as.character(d1[,i])==as.character(d2[i,2]))) > > > Elias T. Krainski > > > >________________________________ > > De: João José de Simoni Gouveia <joao.gouveia(a)univasf.edu.br> > >Para: r-br(a)listas.c3sl.ufpr.br > >Enviadas: Sábado, 10 de Dezembro de 2011 22:21 > >Assunto: [R-br] Modificar valores de dataframe com base em outro dataframe > > > >Prezados, > > > >Possuo dois dataframes (nomeados aqui de 1 e 2). Preciso modificar os > valores do dataframe1 com base no dataframe2. > >Mais precisamente: > >Se os valores da coluna snp1 (dataframe1) forem iguais ao valor da linha > snp1 (dataframe2), preciso converter estes valores para 0 (zero), caso > contrário o valor deve ser convertido para 1 (um). > >Alguém tem alguma dica? > >Att > >JJ > >Abaixo uma ilustração do que eu tenho > >DataFrame 1 --> > > snp1 snp2 snp3 snp4 > >1 A A C T > >2 A G T G > >3 A G C G > >4 A A C T > >5 T G T T > > > >DataFrame 2 --> > > snp anc > >1 snp1 A > >2 snp2 G > >3 snp3 C > >4 snp4 T > > > >E o que eu preciso ter: > >Dataframe final --> > > snp1 snp2 snp3 snp4 > >1 0 1 0 0 > >2 0 0 1 1 > >3 0 0 0 1 > >4 0 1 0 0 > >5 1 0 1 0 > > > >-- > >============================================ > >João José de Simoni Gouveia > > >Médico Veterinário, Msc. > >Doutorando em Zootecnia - PDIZ (UFC/UFPB/UFRPE) > > >Professor Assistente I - Colegiado Acadêmico de Zootecnia > >Universidade Federal do Vale do São Francisco (UNIVASF) > > >---------------------------------------------------------------------------------------- > >Endereço: UNIVASF (Campus Ciências Agrárias) > >Rodovia BR 407, Km 12, Lote 543 > >Projeto de Irrigação Senador Nilo Coelho, s/nº-C1 > >Petrolina, PE - CEP 56.300-000 > > >Telefone: (87)3986-3800/3804-3801 > >Skype: joao_jose_de_simoni_gouveia > >E-mail alternativo: jjsgouveia(a)gmail.com > > >============================================ > >_______________________________________________ > >R-br mailing list > >R-br(a)listas.c3sl.ufpr.br > >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/7ca47fe5/atta… > > > > ------------------------------ > > Message: 7 > Date: Tue, 13 Dec 2011 09:54:57 -0800 (PST) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Categorização de variáveis > Message-ID: > <1323798897.47220.YahooMailNeo(a)web160703.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde, caros r-istas, estou com alguns probleminhas na categorização > de variáveis. > > Tenho este banco. > pad pas sexo pg > 155 134 83 feminino 15 > 156 72 99 feminino 30 > 157 82 151 feminino 10 > 158 23 70 feminino 20 > 159 75 128 masculino 13 > 160 73 125 feminino 05 > 161 125 80 masculino 08 > 162 65 110 feminino 10 > 163 132 92 masculino 13 > 164 114 68 masculino 11 > 165 83 153 feminino 26 > 166 166 107 masculino 22 > 167 114 76 feminino 03 > > > > Preciso classificar da seguinte forma: > se pad<85 e pas<130 é normal > se 85<pad<89 e 130<pas<139 Normal limítrofe > e assim por diante. > > > > a outra classificação é a seguinte: > > Para homens o percentual de gordura > > entre 15 e 20 é normal > > entre 20 e 25 levemente obeso > entre 25 e 30 excessivamente obeso > entre 10 e 14 gordura ideal > 3 gordura essencial > > > > > > Para homens o percentual de gordura > > entre 25 e 30 é normal > > entre 30 e 35 levemente obeso > entre 35 e 40 excessivamente obeso > entre 14 e 18 gordura ideal > 12 gordura essencial > > > > > Como fazer estaS categorizações? > > []'s. > > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/ddf3f837/atta… > > > > ------------------------------ > > Message: 8 > Date: Tue, 13 Dec 2011 15:58:43 -0200 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] Categorização de variáveis > Message-ID: > <CAJ+S-MHuH4RMQ35u092HTUwcCfMP3nqkR+WdzyouauFD9Yhhng(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > cut() > > Em 13 de dezembro de 2011 15:54, Edson Lira <edinhoestat(a)yahoo.com.br > >escreveu: > > > Boa tarde, caros r-istas, estou com alguns probleminhas na categorização > > de variáveis. > > > > Tenho este banco. > > pad pas sexo pg > > 155 134 83 feminino 15 > > 156 72 99 feminino 30 > > 157 82 151 feminino 10 > > 158 23 70 feminino 20 > > 159 75 128 masculino 13 > > 160 73 125 feminino 05 > > 161 125 80 masculino 08 > > 162 65 110 feminino 10 > > 163 132 92 masculino 13 > > 164 114 68 masculino 11 > > 165 83 153 feminino 26 > > 166 166 107 masculino 22 > > 167 114 76 feminino 03 > > > > > > Preciso classificar da seguinte forma: > > se pad<85 e pas<130 é normal > > se 85<pad<89 e 130<pas<139 Normal limítrofe > > e assim por diante. > > > > > > a outra classificação é a seguinte: > > > > Para homens o percentual de gordura > > > > entre 15 e 20 é normal > > entre 20 e 25 levemente obeso > > entre 25 e 30 excessivamente obeso > > entre 10 e 14 gordura ideal > > 3 gordura essencial > > > > > > > > Para homens o percentual de gordura > > > > entre 25 e 30 é normal > > entre 30 e 35 levemente obeso > > entre 35 e 40 excessivamente obeso > > entre 14 e 18 gordura ideal > > 12 gordura essencial > > > > > > Como fazer estaS categorizações? > > > > []'s. > > > > Edson Lira > > Estatístico > > Manaus-Amazonas > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > > > > -- > Luis Iván Ortiz Valencia > Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ > Estatístico Msc. > Spatial Analyst Msc. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/2702a502/atta… > > > > ------------------------------ > > Message: 9 > Date: Tue, 13 Dec 2011 15:14:02 -0300 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Categorização de variáveis > Message-ID: <4EE795EA.6090407(a)gmail.com> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Edson, > > Use o pacote car e a função recode. > > Em 13/12/2011 14:54, Edson Lira escreveu: > > Boa tarde, caros r-istas, estou com alguns probleminhas na categorização > > de variáveis. > > > > Tenho este banco. > > pad pas sexo pg > > 155 134 83 feminino 15 > > 156 72 99 feminino 30 > > 157 82 151 feminino 10 > > 158 23 70 feminino 20 > > 159 75 128 masculino 13 > > 160 73 125 feminino 05 > > 161 125 80 masculino 08 > > 162 65 110 feminino 10 > > 163 132 92 masculino 13 > > 164 114 68 masculino 11 > > 165 83 153 feminino 26 > > 166 166 107 masculino 22 > > 167 114 76 feminino 03 > > > > > > Preciso classificar da seguinte forma: > > se pad<85 e pas<130 é normal > > se 85<pad<89 e 130<pas<139 Normal limítrofe > > e assim por diante. > > > > > > a outra classificação é a seguinte: > > > > Para homens o percentual de gordura > > > > entre 15 e 20 é normal > > entre 20 e 25 levemente obeso > > entre 25 e 30 excessivamente obeso > > entre 10 e 14 gordura ideal > > 3 gordura essencial > > > > > > > > Para homens o percentual de gordura > > > > entre 25 e 30 é normal > > entre 30 e 35 levemente obeso > > entre 35 e 40 excessivamente obeso > > entre 14 e 18 gordura ideal > > 12 gordura essencial > > > > > > Como fazer estaS categorizações? > > > > []'s. > > > > Edson Lira > > Estatístico > > Manaus-Amazonas > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > -- > Cesar Rabak > GNU/Linux User 52247. > Get counted: http://counter.li.org/ > > > ------------------------------ > > Message: 10 > Date: Tue, 13 Dec 2011 10:22:32 -0800 (PST) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Categorização de variáveis > Message-ID: > <1323800552.38193.YahooMailNeo(a)web114707.mail.gq1.yahoo.com> > Content-Type: text/plain; charset=iso-8859-1 > > > > se vc enviasse um CMR seria mais fácil... veja este exemplo > > n <- 10 > x1 <- rnorm(n) > > ### classificação considerando apenas uma variavel > f1 <- factor(findInterval(x1, c(-Inf, -0.5, 0.5, Inf), labels=c("Ruim", > "Médio", "Bom")) > > quando vc tem: > > se pad<85 e pas<130 é normal > >se 85<pad<89 e 130<pas<139 Normal limítrofe > o que deve ser se pad<85 mas pas>130 por exemplo? > > Att. > Elias T. Krainski > > > > >________________________________ > > De: Edson Lira <edinhoestat(a)yahoo.com.br> > >Para: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > >Enviadas: Terça-feira, 13 de Dezembro de 2011 15:54 > >Assunto: [R-br] Categorização de variáveis > > > > > >Boa tarde, caros r-istas, estou com alguns probleminhas na categorização > de variáveis. > > > > > >Tenho este banco. > > pad pas sexo pg > >155 134 83 feminino 15 > >156 72 99 feminino 30 > >157 82 151 feminino 10 > >158 23 70 feminino 20 > >159 75 128 masculino 13 > >160 73 125 feminino 05 > >161 125 80 masculino 08 > >162 65 110 feminino 10 > >163 132 92 masculino 13 > >164 114 68 masculino 11 > >165 83 153 feminino 26 > >166 166 107 masculino 22 > >167 114 76 feminino 03 > > > > > > > > > >Preciso classificar da seguinte forma: > >se pad<85 e pas<130 é normal > >se 85<pad<89 e 130<pas<139 Normal limítrofe > >e assim por diante. > > > > > > > > > >a outra classificação é a seguinte: > > > > > >Para homens o percentual de gordura > > > > > >entre 15 e 20 é normal > > > >entre 20 e 25 levemente obeso > >entre 25 e 30 excessivamente obeso > >entre 10 e 14 gordura ideal > > 3 gordura essencial > > > > > > > > > > > > > > > >Para homens o percentual de gordura > > > > > >entre 25 e 30 é normal > > > >entre 30 e 35 levemente obeso > >entre 35 e 40 excessivamente obeso > >entre 14 e 18 gordura ideal > > 12 gordura essencial > > > > > > > > > > > >Como fazer estaS categorizações? > > > > > >[]'s. > > > > > >Edson Lira > >Estatístico > >Manaus-Amazonas > >_______________________________________________ > >R-br mailing list > >R-br(a)listas.c3sl.ufpr.br > >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > > > ------------------------------ > > Message: 11 > Date: Tue, 13 Dec 2011 15:37:02 -0400 > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, "Elias T. > Krainski" <eliaskrainski(a)yahoo.com.br> > Subject: Re: [R-br] Categorização de variáveis > Message-ID: <63652F35-F09C-4C24-8D38-0F462C8DBBF2(a)yahoo.com.br> > Content-Type: text/plain; charset=utf-8 > > O problema é que você tem duas variáveis, pas é pressão sistólica > e pad é pressão diastólica. Se pas é menor que 130 e pad é menor que 180 > o paciente é considerado normal e assim por diante. > [].s > Edson Lira > Estatístico > Ma-Am > > Em 13/12/2011, às 14:22, "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > escreveu: > > > > > > > se vc enviasse um CMR seria mais fácil... veja este exemplo > > > > n <- 10 > > x1 <- rnorm(n) > > > > ### classificação considerando apenas uma variavel > > f1 <- factor(findInterval(x1, c(-Inf, -0.5, 0.5, Inf), labels=c("Ruim", > "Médio", "Bom")) > > > > quando vc tem: > > > > se pad<85 e pas<130 é normal > > >se 85<pad<89 e 130<pas<139 Normal limítrofe > > o que deve ser se pad<85 mas pas>130 por exemplo? > > > > Att. > > Elias T. Krainski > > > > > > > >> ________________________________ > >> De: Edson Lira <edinhoestat(a)yahoo.com.br> > >> Para: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > >> Enviadas: Terça-feira, 13 de Dezembro de 2011 15:54 > >> Assunto: [R-br] Categorização de variáveis > >> > >> > >> Boa tarde, caros r-istas, estou com alguns probleminhas na > categorização de variáveis. > >> > >> > >> Tenho este banco. > >> pad pas sexo pg > >> 155 134 83 feminino 15 > >> 156 72 99 feminino 30 > >> 157 82 151 feminino 10 > >> 158 23 70 feminino 20 > >> 159 75 128 masculino 13 > >> 160 73 125 feminino 05 > >> 161 125 80 masculino 08 > >> 162 65 110 feminino 10 > >> 163 132 92 masculino 13 > >> 164 114 68 masculino 11 > >> 165 83 153 feminino 26 > >> 166 166 107 masculino 22 > >> 167 114 76 feminino 03 > >> > >> > >> > >> > >> Preciso classificar da seguinte forma: > >> se pad<85 e pas<130 é normal > >> se 85<pad<89 e 130<pas<139 Normal limítrofe > >> e assim por diante. > >> > >> > >> > >> > >> a outra classificação é a seguinte: > >> > >> > >> Para homens o percentual de gordura > >> > >> > >> entre 15 e 20 é normal > >> > >> entre 20 e 25 levemente obeso > >> entre 25 e 30 excessivamente obeso > >> entre 10 e 14 gordura ideal > >> 3 gordura essencial > >> > >> > >> > >> > >> > >> > >> > >> Para homens o percentual de gordura > >> > >> > >> entre 25 e 30 é normal > >> > >> entre 30 e 35 levemente obeso > >> entre 35 e 40 excessivamente obeso > >> entre 14 e 18 gordura ideal > >> 12 gordura essencial > >> > >> > >> > >> > >> > >> Como fazer estaS categorizações? > >> > >> > >> []'s. > >> > >> > >> Edson Lira > >> Estatístico > >> Manaus-Amazonas > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > >> > >> > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 10, assunto 19 > ***************************************** > -- ----- Fabiana Schneck Programa de Pós-Graduação em Ecologia Universidade Federal do Rio Grande do Sul CP 15007 - CEP 91501-970 Porto Alegre, RS, Brasil

1 0

Transformar Variáveis discretas em Numeros Binários -> Pre-processamento para Redes Neurais
by Fernando Neto 14 Dez '11

14 Dez '11

Boa noite, R-anos! Por favor, estou precisando ajustar os arquivos de uma base de dados para processamento de uma RNA. Estou usando a função decodeClassLabels(...) mas a operacao tá ficando muito custosa. Por exemplo __ Masculino Programador Masculino Médico Feminino Programador __ tem que ser __ 0 1 0 1 1 0 1 0 __ Só que eu trato cada variável de cada vez. E acho que a operacao de concatenacao que eu uso: cbind(...) tá tornando a operacao custosa (além de que a memória passa de 2GB). PS.: eu faço o seguinte: *for (i in (VETOR COM AS COLUNAS A SEREM TRANSFORMADAS)){* * colAux = decodeClassLabels(vector[,i])* * vector = vector [,-i]* * vector = cbind(vector, colAux) * *} * Alguém sugere algo? Muito grato!!! Fernando -- ---------------------------------------------------------------- *Fernando Neto* Twitter: @fernandompneto Facebook: facebook.com/fernandompneto Tecnologia de Ponte http://tecnologiadeponte.blogspot.com ---------------------------------------------------------------- *fmpn2 @ CIn - UFPE* http://cin.ufpe.br/~fmpn2 - Engenharia da Computação - Turma 2009.2 - CIn, UFPE. - Monitor de Estatistica e Probabilidade Para Engenharia da Computacao<http://www.cin.ufpe.br/%7Eet586/> <http://cin.ufpe.br/%7Eif672cc> ---------------------------------------------------------------- Confidencialidade*: *A informação contida nesta mensagem de e-mail, incluindo quaisquer anexos, é confidencial e está reservada apenas à pessoa ou entidade para a qual foi endereçada. Se você não é o destinatário ou a pessoa responsável por encaminhar esta mensagem ao destinatário, você está, por meio desta, notificado que não deverá rever, retransmitir, imprimir, copiar, usar ou distribuir esta mensagem de e-mail ou quaisquer anexos. Caso você tenha recebido esta mensagem por engano, por favor, contate o remetente imediatamente e apague esta mensagem de seu computador ou de qualquer outro banco de dados. Muito obrigado. Confidentiality Notice: The information contained in this email message, including any attachment, is confidential and is intended only for the person or entity to which it is addressed. If you are neither the intended recipient nor the employee or agent responsible for delivering this message to the intended recipient, you are hereby notified that you may not review, retransmit, convert to hard copy, copy, use or distribute this email message or any attachments to it. If you have received this email in error, please contact the sender immediately and delete this message from any computer or other data bank. Thank you. --------

4 10

split-plot
by Fabiana Schneck 14 Dez '11

14 Dez '11

Prezados colegas, tenho um experimento com o seguinte design split-plot: fator 'ambiente': 5 corredeiras e 5 remansos em cada plot de 'ambiente' (num total de 10) há outro fator: 'peixes' com dois tratamentos: exclusão e controle ainda, dentro de cada sub-plot do fator 'peixes' existe o fator 'substrato' com dois tratamentos: liso e rugoso. A variável resposta é quantidade de clorofila em cada unidade experimental (cada substrato). Um exemplo dos dados: clorofila ambiente peixes substrato bloco 0.114 1 1 1 1 0.160 1 1 2 1 0.126 1 2 1 1 0.149 1 2 2 1 0.068 1 1 1 2 0.137 1 1 2 2 0.080 1 2 1 2 0.160 1 2 2 2 .... .... 0.080 2 1 1 10 0.068 2 1 2 10 0.114 2 2 1 10 0.124 2 2 2 10 Considerei que os fatores 'peixes' e 'substrato' estão em blocos (10), enquanto o fator 'ambiente' está acima dos blocos (5 blocos são no ambiente corredeiras e 5 no ambiente remansos). Analise os dados utilizando dois modelos distintos: *aov (log(cloro)~(amb+peixes+substrato)^2+Error(bloco/peixes))* * * Error: bloco Df Sum Sq Mean Sq F value Pr(>F) amb 1 0.1162 0.11618 0.128 0.7297 Residuals 8 7.2599 0.90749 Error: bloco:peixes Df Sum Sq Mean Sq F value Pr(>F) peixes 1 2.96190 2.96190 14.6147 0.005069 ** amb:peixes 1 0.30121 0.30121 1.4862 0.257528 Residuals 8 1.62133 0.20267 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) substrato 1 1.9924 1.99244 10.2620 0.005209 ** amb:substrato 1 0.0013 0.00134 0.0069 0.934840 peixes:substrato 1 0.5934 0.59335 3.0560 0.098473 . Residuals 17 3.3007 0.19416 * * *lme(log(cloro)~(amb+peixes+**substrato)^2, random=~1|bloco/peixes))* Linear mixed-effects model fit by REML Data: NULL AIC BIC logLik 85.58484 100.5499 -32.79242 Random effects: Formula: ~1 | bloco (Intercept) StdDev: 0.4197679 Formula: ~1 | peixes %in% bloco (Intercept) Residual StdDev: 0.06522495 0.4406333 Fixed effects: log(cloro) ~ (amb + peixes + substrato)^2 Value Std.Error DF t-value p-value (Intercept) -2.4342471 0.2647064 17 -9.196027 0.0000 amb2 -0.2929029 0.3611524 8 -0.811023 0.4408 peixes2 0.6142687 0.2448449 8 2.508808 0.0364 substrato2 0.6783940 0.2413448 17 2.810891 0.0120 amb2:peixes2 0.3471056 0.2847218 8 1.219104 0.2575 amb2:substrato2 0.0231237 0.2786810 17 0.082976 0.9348 peixes2:substrato2 -0.4871762 0.2786810 17 -1.748150 0.0985 Correlation: (Intr) amb2 peixs2 sbstr2 amb2:p2 amb2:s2 amb2 -0.682 peixes2 -0.462 0.229 substrato2 -0.456 0.223 0.329 amb2:peixes2 0.269 -0.394 -0.581 0.000 amb2:substrato2 0.263 -0.386 0.000 -0.577 0.000 peixes2:substrato2 0.263 0.000 -0.569 -0.577 0.000 0.000 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -2.42000225 -0.41789060 -0.02451742 0.44159262 1.87090777 Number of Observations: 40 Number of Groups: bloco peixes %in% bloco 10 20 Tenho dúvidas se estou definindo de forma correta os modelos e também se devo usar aov ou lme, uma vez que os resultados dos termos isolados diferem entre as duas funções (porém os resultados são indênticos para as interações). Desde já obrigada pela atenção, Fabiana -- ----- Fabiana Schneck Programa de Pós-Graduação em Ecologia Universidade Federal do Rio Grande do Sul CP 15007 - CEP 91501-970 Porto Alegre, RS, Brasil

2 1