Dezembro 2013 - R-br

regressão linear polinomial: classificar em simples ou múltipla (dúvida conceitual)
by Jose Claudio Faria 12 Dez '13

12 Dez '13

Caros membros, Vou usar da lista para possivelmente redimir uma dúvida conceitual, pois não vejo melhor local para fazê-lo, dado ao perfil dos usuários. Muitos autores classificam os modelos de regressão linear assim: 1- Simples (y = Bo + B1X) 2- Polinomial (y = Bo + B1X + B2X^2 + + B3X^3...) 3- Múltipla (y = Bo + B1X + B2Y + B3XY...) Outros (dando ênfase na variável preditora) classificam as polinomiais (dependentes de apenas uma variável fixa ou aleatória) como simples e usam r² (minúsculo) para o coef. de determinação, independente do grau do polinômio. Já outros (dando ênfase no modelo) classificam as polinomiais (dependentes de apenas uma variável fixa ou aleatória) como múltipla, e usam R² (maiúsculo) para o coef. de determinação, desde que tenha mais de um grau e independente do número de graus do polinômio. Gostaria muito da opinião do pessoal mais especializado sobre esse assunto! Não seria uma caso como o da "a cerveja que desce redondo"? Ou seja se a ênfase está no substantivo (cerveja) deveria ser: desce redonda. Outros pensam que "redondo" é um advérbio e se refere ao verbo "descer", portanto, que está correta a concordância no masculino. http://salp-pmf.blogspot.com.br/2012/06/cerveja-que-desce-redondo-ou-redond… Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

5 8

Re: [R-br] regressão linear polinomial: classificar em simples ou múltipla (dúvida conceitual)
by Ivan Bezerra Allaman 11 Dez '13

11 Dez '13

Para mim é claro e evidente que a diferença entre uma regressão linear simples e uma múltipla é o número de variáveis preditoras (ou independentes) indiferentemente do grau do polinômio. Vejamos a seguinte situação: Y = a + bX Y = a + bX + cZ O que os diferencia? - Ambos são polinômios do primeiro grau? Sim - Ambos tem o mesmo número de variáveis preditoras? Não . Então isso os diferencia? Sim Vejamos outra situação: Y = a + bX + cX^2 Y = a + bX + cX^2 + dZ + dZ^2 O que os diferencia? - Ambos são polinômios do segundo grau? Sim - Ambos tem o mesmo número de variáveis preditoras? Não . Então isso os diferencia? Sim Lembrando láaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa da matemática básica, todo *monômio *(comumente conhecido como regressão linear simples) é um *polinômio*. Portanto, o que diferencia um polinômio simples de um múltiplo é o número de variáveis preditoras. (s,f,p)

1 0

Converter dados de contagem em binomial
by ASANTOS 11 Dez '13

11 Dez '13

Caros Membros, Estou tentando converter dados de contagem de insetos mortos em notação binomial e estou fazendo algo errado, pois se tenho: trat<-gl(5,10) ## 5 tratamentos com dez repetições tempo<- gl(10,50) ## 10 avaliações no tempo dados<-expand.grid(trat, tempo) colnames(dados)<-c("trat","tempo") y<-rpois(length(dados[,1]),10) ## Simula a variável resposta dados<-cbind(dados,y) ## Converter os dados de contagem em binomial, por ex. se tenho 5, devo converter em 1 1 1 1 1 dados.bin <- NULL for (i in 1:nrow(dados)) dados.bin <- c(dados.bin, rep(c(1,dados[i,3]))) Por que para cada valor de contagem na coluna dados[,3] não tenho a transformação que espero? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

4 7

Como fazer uma Análise Fatorial
by Jodavid Ferreira 09 Dez '13

09 Dez '13

Alguém pode me ajudar a como realizar uma análise Multivariada no R? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

Análise Discriminante e Análise Fatorial
by Jodavid Ferreira 09 Dez '13

09 Dez '13

Como realizar uma Analise Fatorial e Análise Discriminante no software R? Onde posso encontrar material, ou algum algoritmo, script? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

filliben
by Hélio Gallo Rocha 07 Dez '13

07 Dez '13

Caros listeiros, bom dia. Estou com uma Li um comunicado técnico da EMBRAPA utilizando o Teste de Filliben, onde analisa os resíduos ortonormais para validação modelos geoestatísticos. neste link: https://www.dropbox.com/s/e9ww5h5w9dz8wed/CT60_FILLIBEN.pdf Fiz uma busca na lista e não achei nenhuma referencia. Achei neste link: http://genepi.qimr.edu.au/staff/davidD/R/filliben.R esta função: filliben <- function(x) { y<-sort(x) n<-length(y) p<-(1:n-0.3175)/(n+0.365) p[1]<-1-0.5^(1/n) p[n]<-0.5^(1/n) p<-qnorm(p) cor(y,p) } Para obter o resultado do teste: filliben(rstandard(reglin_esc)) onde, segundo o link reglin_esc seria produto de uma regressão, pensei entre os dados originais e os estimados pela krigagem, até ai tudo bem. No Comunicado técnico ele calcula o resíduo ortonomal. Alguém poderia dar uma ajuda? Agradeço antecipadamente, -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 3

Extração de coordenadas e tempo em múltiplas imagens
by ASANTOS 05 Dez '13

05 Dez '13

Caros Membros, Estou precisando de ajuda para criar um laço para extrair as coordenadas e o tempo (no nome do arquivo) em múltiplas figuras. No meu exemplo abaixo e usando apenas a figura 1, tenho uma regiões de pixels de interesse em branco e o nome da figura Fig1-0_00_00.png, preciso das coordenadas x e y das regiões em branco (Valor digital = 0), representado pelo objeto fig.data.pix e juntei ao tempo que é 0:00:00, pois o nome da Figura é 0_00_00. Conforme CRM abaixo, fazendo figura por figura esta OK, mas gostaria de fazer um laço para extrair essas informações das três figuras abaixo e armazená-las em um único objeto, alguém teria alguma sugestão? #------------------------------------------------------------------------------- # Rotina para extração de coordenadas e tempo em imagens #------------------------------------------------------------------------------- # Abertura das imagens require(RCurl) #Figura 1 ---------------------------------------------------------------------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/8b99gn9ve1yo7is/Fig1-0_00_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig1-0_00_00.png", "Fig1-0_00_00.png") #Figura 2 ---------------------------------------------------------------------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/8b99gn9ve1yo7is/Fig2-0_20_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig2-0_20_00.png", "Fig2-0_20_00.png") #Figura 3 --------------------------------------------------------------------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/bo6gbt3ca6xwxm0/Fig3-0_40_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig3-0_40_00.png", "Fig3-0_40_00.png") # Extraido as coordenadas dos pixels para figura 1 ----------------------------- require(raster) r <- raster("Fig1-0_00_00.png") image(r) r2<-as.data.frame(r) # Cria as coordenadas ---------------------------------------------------------- grid0 <- expand.grid(x=seq(0,length.out=r@nrows), y=seq(0,length.out=r@ncols)) # Coordenadas + valor dos pixels fig.data<-cbind(grid0,r2) # Região de interesse - Número digital = 0 fig.data.pix<-fig.data[fig.data[,3]==0,] # Extrai o tempo do nome do arquivo tempo2<-r@data@names tempo1<-gsub("Fig1.", "", tempo2) tempo<-gsub("_", ":", tempo1) # Cria o arquivo final dados<-cbind(fig.data.pix[,1:2],rep(tempo)) head(dados) #------------------------------------------------------------------------------- Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

3 3

zip
by ana paula coelho madeira 05 Dez '13

05 Dez '13

Boa tarde pessoal, Estou tentando ajustar um modelo com excesso de zeros, usando a função zeroinfl() do pacote pscl. No entanto, quando rodo o programa aparece a seguinte mensagem de erro: require(pscl) zip <-zeroinfl(Y~T+C+C:T,dist=c("poisson"), link=c("logit")) Error: could not find function "zeroinfl" Alguém sabe o que pode estar ocorrendo? Obrigada Ana Paula

3 6

Espelho do CRAN no Nordeste do Brasil: novo link
by Jose Claudio Faria 05 Dez '13

05 Dez '13

Olá para todos, Informo que ocorreu uma pequena alteração no link do novo repositório experimental do CRAN no Nordeste do Brasil (UESC/Ilhéus/Bahia). O dados corretos estão abaixo: - URL: http://nbcgib.uesc.br/mirrors/cran/ - Country: Brazil - Name: Brazil (BA) - Cidade: Ilhéus - Host: NBCGIB - Núcleo de Biologia Computacional e Gestão de Informações Biotecnológicas. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Gráfico de pontos
by Alex 05 Dez '13

05 Dez '13

Ao grupo r-br, por favor, eu criei um gráfico de pontos onde o tamanho e cor de cada um é baseado no quartil ao qual ele pertence, o trecho do código está abaixo. Entretanto estou com dificuldades de criar a legenda. Alguém poderia me orientar? *attach(dadosT1)* *par(mfrow=c(1,1))* *ramp<-colorRamp(c("yellow", "red"))* *points(dadosT1,xlab="Coordenada X (coluna)",* *ylab="Coordenada Y (linha)",* *pt.divide="quartiles", col=rgb(ramp(seq(0,1, length=5)), max=255))* Desde já agradeço, Alex Santos

3 2