Junho 2015 - R-br - DINF UFPR

Sistema é computacionalmente singular - Matriz/Mahalanobis
by Andre Oliveira 26 Jun '15

26 Jun '15

Prezados boa noite,estou com problemas em uma matriz inversa! Alguma sugestão? Precisava muito desta matriz D2! require(biotools)setwd("/home/andre")dados <- read.csv("Multivariada.txt", sep="")Cov=cov(dados)D2=D2.dist(dados, Cov) # Não hpa inversa!Error in solve.default(cov, (x[i, ] - x[j, ])) : sistema é computacionalmente singular: condição recíproca número = 1.53112e-19 André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

Dúvida glm (family)
by Karen Castillioni 25 Jun '15

25 Jun '15

Pessoal, estou tentando analisar dados de cobertura vegetal (5%, 10%, 15%, 20%..até 100%) e dados de biomassa (em gramas por m2) usando glm, pois nenhum dos dados apresenta distribuição normal, nem mesmo com transformação. Porém, estou em dúvida sobre qual família usar para cada variável. Alguém pode ajudar? Agradeço a ajuda.

3 3

Ciclos Codependentes
by Rafael Costa 25 Jun '15

25 Jun '15

Prezados colegas, Gostaria de saber se alguém dispõe dos códigos para rodar o teste de ciclos codependentes apresentado no artigo de 1997 dos autores Farshid Vahid e Robert F. Engle intitulado "Codependent cycles". Desde já, agradeço a todos! -- *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em Economia - CAENUniversidade Federal do Ceará - UFC*

1 0

Mahalanobis para o dendograma!
by Andre Oliveira 25 Jun '15

25 Jun '15

Pessoal como fazer para construir um dendrograma com a matriz de distancia de Mahalanobis? Já que Mahalanobis não esta na função dist! Já dei uma pesquisada na lista mesmo e no google, mas não consegui! x <- matrix(rnorm(100*3), ncol = 3)D2 <- mahalanobis(x, colMeans(x), cov(x))hc <- hclust(D2,method="ward.D") # erro in if (is.na(n) || n > 65536L) stop("size cannot be NA nor exceed 65536") : valor ausente onde TRUE/FALSE necessário André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

Fatorial 2^4 com ponto central
by ana paula coelho madeira 24 Jun '15

24 Jun '15

Prezados colegas, boa tarde. Tenho um planejamento fatorial 2^4 com um ponto central ( 4 repetições), totalizando 20 experimentos. Alguém tem sugestão de como fazer a análise no R (obs: consegui fazer no Statistica, mas queria apresentar no R)? Desde já agradeço Att., Ana Paula

5 9

Quantificar itens
by jose gomes de souza filho Gomes 24 Jun '15

24 Jun '15

Seja o banco de dados: http://www.datafilehost.com/d/5a8cd3a4# Ao usar a função *table*( ) para quantificar os itens, tive a seguinte situação: - Não diferenciou o item 2.1 do 2.10, ou seja, eram iguais. Como que faço para solucionar? Att Gomes

3 2

Re: [R-br] Histograma em escala logarítmica
by Romero Luiz M. Sales Filho 23 Jun '15

23 Jun '15

Caros amigos, Desculpem a demora em responder, estive com sérios problemas na net aqui de casa. Agradeço muito pelas valiosas dicas... Estarei agora testando todas elas e espero que consiga representar esses dados da melhor forma possível.. Mas realmente como foi dito nos comentários, é praticamente impossível conseguir um gráfico "bonito" com tão poucos dados, e principalemente quando esses dados apresenta alguns valores extremos... Muito obrigado! Abraços, Romero. -- *Profº Romero Luiz Mendonça Sales Filho* *Estatístico / Mestre em Engª de ProduçãoUFRPE - UAGCel.: (81) 9667-4477* *romero.sfilho(a)gmail.com* <romero.sfilho(a)gmail.com> *romero_sfilho(a)yahoo.com.br* <romero_sfilho(a)yahoo.com.br>

2 1

Ler arquivo .SHP
by Andre Oliveira 21 Jun '15

21 Jun '15

Boa tarde, baixei o arquivo .zip e extraí, mas não consigo ler o mesmo. Sempre li arquivos assim .... require(rgdal) mt=readOGR("51MUE250GC_SIR.shp",layer="51MUE250GC_SIR",encoding = "latin1") Mas, hoje baixei o arquivos e veio assim com estes formatos ... dir() [1] "1200401.DBF" "1200401.PRJ" "1200401.SHP" "1200401.SHX" require(rgdal) > fauna=readOGR("1200401.SHP",layer="1200401.SHX") Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv, : Cannot open layer O que eu poderia fazer no código de leitura para resolver?obrigado

3 2

Teste de Kruskal Wallis
by ana paula coelho madeira 19 Jun '15

19 Jun '15

Prezados, bom dia. Preciso fazer uma comparação entre4 tratamentos via teste de Kruskal Wallis. Utilizei o comando: kruskal.test(Resposta~trat, dados) Kruskal-Wallis rank sum test data: Resposta by trat Kruskal-Wallis chi-squared = 14.602, df = 3, p-value = 0.00219 indicando haver diferenças significativas. Como fazer comparações múltiplas entre os grupos? Desde já agradeço. Att., Ana Paula

5 7

Salvar gráficos via Rscript no Terminal
by Wagner Wolff 19 Jun '15

19 Jun '15

Olá pessoal estou com um problema, quero salvar uma figura .png via comando no Terminal por meio do comando Rscript, mas estou obtendo o seguinte erro: Error in plot.new() : could not open file '/opt/lampp/htdocs/Rphp/temp.png' Calls: hist -> hist.default -> plot -> plot.histogram -> plot.new Sendo assim alguém poderia me ajudar a contornar este erro? Segue a rotina (my_rscript.R): *args <- commandArgs(TRUE) N <- as.double(args[1])x <- rnorm(N,0,1)png(filename="/opt/lampp/htdocs/Rphp/temp.png", width=500, height=500)hist(x, col="lightblue")dev.off()* e o Comando no Terminal: *Rscript my_rscript.R 100*

4 4