- R-br - DINF UFPR

Teste Comparação - Pacote Agricolae
by Evelyn Alecrim 19 Jan '12

19 Jan '12

Olá a todos, Alguém sabe me informar qual o nome do teste de comparação feito após o Kruskal Wallis no pacote Agricolae? Obrigada! -- Evelyn da Fonseca Alecrim Bragion Bióloga/Universidade Federal de Lavras - UFLA Mestranda em Ecologia Aplicada - UFLA Lavras - MG/Brasil

1 0

Distribuição de pareto
by Vinicius Brito Rocha 19 Jan '12

19 Jan '12

Senhores, estou querendo gerar números a partir da distribuição de Pareto. Baixei o pacote VGAM. Meu problea é gerar dados de acordo com a proposta de que 80% da riqueza está na mão de 20% da sociciedade (no brasil é 90% em 2% rs) Alguem tem alguma noção de qual é a paremetrização para isso? No Pacote VGAM tem varias opões da distribuição de Pareto. 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.

1 0

Re: [R-br] Dúvida teste
by Ivan Bezerra Allaman 18 Jan '12

18 Jan '12

Bom dia Peterson! Sinceramente eu não entendi o seu problema. Ficou meio confuso. Você tem 4 linhas certo? Deu a entender que: Linha 1: Época 1 Linha 2: Época 2 Linha 3: Época 3 Linha 4: Época 4 É isso? Se for isso, cada planta constitui uma unidade experimental. Em outro momento, deu entender: Linha 1,2,3,4: Época 1 Linha 1,2,3,4: Época 2 Linha 1,2,3,4: Época 3 Linha 1,2,3,4: Época 4 Nesse caso, cada linha seria uma repetição. Talvez outros colegas tenham entedido, mas para agilizar a ajuda, se você "desenhar" fica melhor para ajudá-lo. Abraço! (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 2

Acessar funções de pacotes!
by Ivan Bezerra Allaman 18 Jan '12

18 Jan '12

Bom dia senhores! Gostaria de saber como faço para acessar as funções dos pacotes do R. Em um tempo remoto, bastava acessar a pastar "R", que lá estava todas as funções. Nos dias atuais, aparece aquirvos com extensões '.rdb' e '.rdx' que eu não consigo abrir. Será que as funções estão nestes arquivos? Caso sim, qual programa uso para abrí-los? Utilizo o ubuntu 10.04 lts. Desde já, agradeço a atenção de todos! (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 5

Dúvida teste
by Peterson Gardin 18 Jan '12

18 Jan '12

Pessoal Eu tenho um parreiral de uvas cabernet sauvignon. Nesse, tenho 4 linhas com 60 plantas cada. Fiz 4 tratamentos, que correspondem a época de colocação de uma cobertura plástica. Dependendo da época que coloquei a cobertura, ao final do ciclo, ocorrem diferenças de ataques de doenças, durante a maturação das uvas. Não existem repetições destas coberturas plásticas, pois cada linha custa mil reais, portanto, cada linha é um tratamento (ou seja, época de colocação da lona), pois se eu fizesse 4 repetições o experimento me comeria R$ 16.000,00, só de cobertura. Penso que 60 plantas sejam suficientes para realizar uma amostragem. Pergunto: Posso realizar uma amostragem aleatória em cada tratamento e fazer um teste de média como se fosse um DIC? Como faço se quiser avaliar uma doença que aumenta conforme transcorre o tempo de maturação, por exemplo, quanto mais madura a uva mais ataca. Coloco no modelo a variável tempo? Tenho dúvidas de analisar variáveis que evoluem com o tempo. Se tiverem alguma dica de bibliografia, para que eu veja as análises, tente fazer e retorne com as dúvidas eu agradeço. Obrigado. João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054

1 0

Problema na instalação do pacote geepck
by Olympio Neto 18 Jan '12

18 Jan '12

Pessoal, bom dia. Ontem tentei instalar o pacote 'geepack' no R.2.14.1 rodando no Ubuntu 11.10, mas ao terminar de compilar o pacote dá um erro que impossibilita a finalização da instalação. Alguém já passou por isso? No Windows instala corretamente. Alguém teve o mesmo problema? Abs Olympio

2 1

Gerar gráficos diferentes e em 3D
by Mauro Sznelwar 18 Jan '12

18 Jan '12

Tenho esta situação: f=rmultinom(30,3,prob=c(0.5,0.3,0.8))*dmultinom(x=1:3,prob=c(0.5,0.3,0.8)) n=5 for(i in 1:n){ g=f contour(g) } #preciso que gere 5 gráficos diferentes (mas não parece que gera gráficos #diferentes) jpeg(filename = "MeusGraficos%03d.jpg", width = 480, height = 480) for (i in 1:5) contour(g) dev.off() E como adapto o caso da normal para o multinomial para fazer gráfico 3D, já tentei mas acusa medidas diferentes. z <- matrix(rnorm(50*10),50,10) image(z) persp(z)

1 0

Dúvida equação GLM e Anova
by Caroline Figueiredo 17 Jan '12

17 Jan '12

Olá pessoal, Estou com uma dúvida, que na verdade deve ser bem ridícula, mas realmente não sei a resposta. A questão é a seguinte: Tenho dados binários como variável resposta ( 1 – vivo; 0 - morto), e variáveis contínuas como explicativas, e vi no The R book do Crawley que poderia usar GLM para construir um modelo, tendo binomial como família. Usei o seguinte código (os valores são fictícios, pois li meus dados direto da tabela, pois são muitas linhas). status<-c(0,1,0,1,1,1,1,1,0,1) altura<-c(2.2,1.3,4.5,6.7,1.3,4.5,1.2,2.0,7.8,1.2) newy<-c(5.5,1.2,7.8,1.4,1.3,2.3,3.2,1.3,6.7,1.4) m1<-glm(status~altura+newy+altura*newy, binomial) E o summary deu isso: Call: glm(formula = status ~ altura + newy + altura * newy, family = binomial) Deviance Residuals: Min 1Q Median 3Q Max -3.3122 0.1586 0.1969 0.2486 0.6394 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.24827 0.52797 0.470 0.63818 altura 0.80654 0.20006 4.032 5.54e-05 *** newy 0.23918 0.05101 4.689 2.74e-06 *** altura:newy -0.04578 0.01753 -2.611 0.00903 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 1705.0 on 7191 degrees of freedom Residual deviance: 1610.4 on 7188 degrees of freedom AIC: 1618.4 Number of Fisher Scoring iterations: 7 Minha dúvida é a seguinte: Qual é a equação ajustada para esse modelo? P. ex., se fosse uma regressão linear o modelo seria y~x, e a equação y=a + bx, mas e para esse modelo? A outra dúvida é se poderia fazer uma anova para comparar esse modelo com o modelo nulo, e qual teste poderia usar, pois no exemplo do livro ele usa o test= chi, porém está comparando dois modelos diferentes, e queria saber se posso usar para comparar com o modelo nulo também. Bem é isso pessoal, espero que possam me ajudar. Obrigada, Carol Figueiredo

4 3

Operações com datas
by ASANTOS 17 Jan '12

17 Jan '12

Boa tarde pessoal, Tenho dados de temperatura e precipitação diárias de 1994 a 2001 em uma determinada fazenda. Preciso retira as médias do dia 1º ao 15º dia e do 15º ao último dia de cada mês, porém estou tendo problemas em trabalhar com as datas, exemplificando através de uma rotina venho fazendo: ### Criando o objeto com datas e a variável climática x <- as.Date("1994-01-01") times <- 1:2865 datas <- x+times datas<-as.data.frame(datas)##Cria o vetor data de 01/01/1994 a 5/11/2001 variavel<-rnorm(2865)## Simula a variável medida no dia dados<-cbind(datas,variavel)##Cria o objeto ## Cria um vetor temporal a cada 15 dias e adiciona aos dados t<-1:191 for(i in 1:length(t)){xv<-sort(rep(t,15))} xv<-as.data.frame(xv) xv<-xv[1:2865,] tp.1994_2001<-cbind(dados,xv) head(tp.1994_2001) ## Calcula as médias a cada 15 dias md<-1:2865 res<-NULL for(k in 1:length(md)){Tm<-mean(tp.1994_2001[md[k],2]) res<-rbind(res,c(md[k],Tm)) } Porém minha abordagem esta incorreta, quando na verdade eu desejo a média do 1° ao 15º dia e do 15º ao último dia do mês, e não a cada 15 dias, peguei alguns materiais de séries temporais usando o R, mas todos tratam diretamente das análises e não da manipulação de dados. Alguém saberia que abordagem eu poderia empregar para obter os resultados corretos ou alguma literatura que envolva a manipulação de datas, 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

4 5

Como eu retiro so os 3 primeiro caracteres dos objetos de um Vetor?
by Augusto Ribas 17 Jan '12

17 Jan '12

Ola: #Vamos supor que eu tenha um vetor: teste<-c("Stomoxys.calcitrans","Musca.domestica","Syritta.flaviventris") Eu gostaria de estrair so as 3 primeiras letras, pra ter um resultado [1] "Sto" "Mus" "Syr" Meu interesse é em usar esse segundo vetor extraido para plotar um grafico, mas os nomes de especies inteiro ficam muito grandes e fica muito poluido, mas se der pra extrair os nomes assim fica mais facil fazer o grafico ja que são muitos nomes e além do trabalho é perigoso errar ordem de especies e tal ou escrever errado. Então se alguem souber como fazer. Outra coisa não sei se é muito complicado, mas todo nome de especie é o "Genero.especie" Daria pra mim tirar a primeira letra do Genero, e a primeira e segunda letra da especies, pra ter um resultado do tipo: [1] "Sca" "Mdo" "Sfl" E por ultimo, se o vetor for mistura especies e familias tipo: teste2<-c("Stomoxys.calcitrans","Musca.domestica","Syritta.flaviventris","Clusiidae") Neste caso o ultimo termo é uma familia, então ele é uma palavra so, e nao ta no esquema "Genero.especie" Daria pra nele ficar so as 3 primeiras letra, e chegar a um resultado final assim: [1] "Sca" "Mdo" "Sfl" "Clu" PS: "Genero.especie" Ta com um ponto separando os 2 mesmo. Não é espaço não. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 6