- R-br - DINF UFPR

None
by None 19 Abr '13

19 Abr '13

1 0

Transformação de dados (Box-Cox, Pot. est. variância de Tukey) e Parcela perdida
by Ederson Civardi 19 Abr '13

19 Abr '13

Olá, como vai? Alguém pode me ajudar a tirar umas dúvidas a respeito de quantas vezes posso transformar os dados. Já ouvi as duas situações, pode somente transformar uma vez seus dados e outros falam que posso transformar os dados quantas vezes forem necessárias para obter as pressuposições básicas. O que é o correto? Quantas transformações posso fazer nos mesmos dados? Outra dúvida é sobre parcela perdida, como deve usar ao digitar os dados para que o R entenda? Em branco, ponto, NA, Na ou na? Agradeço quem puder ajudar. Ederson

3 2

Re: [R-br] Pacote IBGEPesq e dicionarioIBGE
by Roney Fraga Souza 19 Abr '13

19 Abr '13

Olá Luciane, A vantagem de utilizar o pacote 'SAScii' é o fato dele importar os dados com o dicionário do SAS. Att Roney ## ---------------------------------------- # instalando o pacote 'SAScii' elaborado por Anthony Joseph Damico install.packages('SAScii') # carregando o pacote library(SAScii) # é necessário eliminar os comentários do dicionário, os /* … */, para realizar a leitura correta # lendo o dicionário da PNAD 2011 pessoas parse.SAScii( "2011_DOM_input.sas") # lendo o arquivo txt dom2011 <- read.SAScii('2011_DOM.TXT', '2011_DOM_input.sas', beginline=11, buffersize=1000) ## ---------------------------------------- # a função read.SAScii.sqlite() é muito útil para carregar os dados direto para SQLite source_url( "https://raw.github.com/ajdamico/usgsd/master/SQLite/read.SAScii.sqlite.R" ) drv <- dbDriver("SQLite") con <- dbConnect(drv, dbname='pnad') # transformar um arquivo do IBGE direto em uma base SQLite read.SAScii.sqlite(fn = "2011_PES.TXT", sas_ri = "2011_PES_input.sas", tl = TRUE, tablename = "pes2011", db = con)

3 2

raster e shape file
by Hélio Gallo Rocha 18 Abr '13

18 Abr '13

Caros, Também estou as voltas com resultados de krigagen. Estou acompanhando a msg do Alexandre Santos. Já pesquisei bastante, mas não achei solução para transformar o resultado de uma krigagem em um raster ou num shape file. Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

5 6

Espaços em branco
by Fátima Lima Paula 18 Abr '13

18 Abr '13

Pessoal, estou com um banco de 273 mil linhas. Envio o banco em .dbf para o excel, não mexo nele em nada. Salvo em xlsx e depois em csv. Quando eu leio o banco em csv (read.csv2) ele vem todo aparentemente bonitinho, mas a partir de uma variável ele zebra. Acho que o problema é no excel. Verifiquei que na coluna de diagobito e mae os espaços em branco não estão sendo reconhecidos como espaço em branco, no excel. Já alterei, no R, essas variáveis para factor, character, mas o problema persiste. Não tem como enviar o banco para vocês. Só queria saber se alguém já passou por isso e pode me ajudar. Vejam como fica o summary da variável: diagobito :233299 J960 :7018 etc................... Obrigada Fátima

6 13

Controlando ordem dos tratamentos em barplot2() do pacote gplots
by ASANTOS 18 Abr '13

18 Abr '13

Boa noite pessoal, Tenho dois tratamentos 1 casal e 10 casais mas na hora de plotar o barplot com a função barplot2() do pacote gplots, sempre aparece primeiro a barra 10 casais antes do 1 casal, tentei modificar inclusive a ordem do banco de dados e não dei certo, já li o help(barplot2) mas não encontrei nada que me permita controlar isto, alguém tem uma dica, Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

4 5

Espaço em branco
by Fátima Lima Paula 17 Abr '13

17 Abr '13

Pela vigésima vez. Agora deve ir. Tem uma coluna e 2000 linhas.

1 0

Pacote IBGEPesq
by Luciane Maria Pilotto 17 Abr '13

17 Abr '13

Olá pessoal, Gostaria de ajuda para instalar e carregar o pacote IBGEPesq no Wind XP para trabalhar com os micro dados das PNADs, está dando o erro descrito abaixo. Encontrei um debate com a mesma dúvida, porém, não vi a solução. Sou novata no R e nesta lista também. library(IBGEPesq) Erro em library(IBGEPesq) : package ‘IBGEPesq’ does not have a NAMESPACE and should be re-installed Att., Luciane

3 4

Matrix Hessiana - duvida
by Carlos Pereira 17 Abr '13

17 Abr '13

Oi pessoal to com um problema para calcular a matriz hessiana do seguinte modelo modelo = x1/((1+((x1/x2)**(x3*t) - 1))) sendo quet = rnorm(10,5,3)x1=8500x2=0.0156x3=0.26 Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

4 4

Gráfico do tipo stacked bar plot modificado
by ASANTOS 16 Abr '13

16 Abr '13

Boa tarde pessoal, Gostaria de fazer um gráfico do tipo stacked bar plot horizontal, mas gostaria que as cores das barras ficassem alinhadas pelo níveis, como por exemplo: counts <- table(mtcars$vs, mtcars$gear) barplot(counts, main="Car Distribution by Gears and VS", horiz=TRUE,xlab="Number of Gears", col=c("white","black"), legend = rownames(counts)) abline(v=7,col="red") # Que agora eu gostaria que as cores branca e preto ficassem alinhados por um valor arbitrário, que no meu caso seria 7 por exemplo, alguém sabe se existe um nome para este tipo de barra e se alguém do grupo já fez algo semelhante? Uma vez que help(barplot) não ajudou muito, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

2 2

ler datas no arquivo
by Helena Turon 16 Abr '13

16 Abr '13

Tenho um vetor de datas e horário idealizado, sem falhas. Preciso comparar este vetor com os dados observados, para contar "NAs". As datas observadas estão sendo reconhecidas pelo R como datas (Mas, só individualmente). dia(2), mês(2) e ano(4), tenho outro vetor com o horário das observações: 0, 12 18. Preciso juntar a data com a hora da observação. Quando colocadas num vetor numérico as datas ficam números (ex: -3271). Quando são colocadas em vetor lista, não podem ser transformadas para data, para serem comparadas. Assim, não consigo nem colocar horário, nem comparar. Helena

2 1

Simulação de dados
by Thiago de Andrade 15 Abr '13

15 Abr '13

Amigos, Gostaria de saber se vocês têm algum material pra me indicar sobre procedimentos de simulação de dados. Trabalho com psicometria e estou precisando fazer simulações de dados em que eu manipule aspectos como dimensionalidade, comunalidades, saturações de variáveis em componentes etc. Agradeço desde já! -- Thiago de Andrade Mestrando em Psicologia Social, PPGP-UFPB Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM)

4 4

Dúvidas com pressuposições da anova
by Ederson Civardi 15 Abr '13

15 Abr '13

Meu experimento está em DBC em parcelas subdivididas, sendo o Fator A (principal), com 5 níveis e o Fator B (secundário) com 4 níveis. Com 4 blocos ou repetições. E foi realizado a campo por 3 anos consecutivos no mesmo local e nas mesmas parcelas. Fiz inúmeras avaliações e tenho problemas com diversas variáveis respostas, pois pela característica do experimento apresenta muitos zeros. Com isso não consigo obter normalidade dos dados, principalmente. Além de não conseguir adequar as outras pressuposições básicas (aditividade, normalidade, homocedasticidade e independência) mesmo utilizando transformações ótimas de Box-Cox e Potencia estabilizadora da variância de Tukey; ou ambas juntas ou várias em conjunto. O que posso fazer? Tem alguma sugestão? Algumas dessas variáveis resposta até consigo atender as pressuposições através dessas transformações ótimas, obtendo uma diferente para casa ano. Mas não consigo usar a mesma para os 3 anos, para efetuar análise conjunta. Já que preciso usar a mesma transformação para os 3 anos para fazer a análise conjunta. É isso mesmo? Envio em anexo a planilha com 2 variáveis resposta que tenho mais problemas. Uma outra dúvida que tenho para as pressuposições é que devo verifica-las nos resíduos, correto? Devo testar a normalidade de forma geral (todos os resíduos) ou ajustados ao modelo (por exemplo: modeloMB1 <- aov(BraqEST2 ~ FatA*FatB*SAFRA+(REP:SAFRA) + Error(REP/FatA/FatB/SAFRA), data=Exp sala r-br2013) ). É isso mesmo? Tenho muitas outras dúvidas, mas deixarei para faze-las posteriormente. Estou aprendendo a utilizar essa maravilhosa ferramenta que é o R. Se puderem me ajudar ficarei muito agradecido. Abraço, Eder

1 0

Erro no calculo da matriz distancia
by Bernardo Silva 13 Abr '13

13 Abr '13

Olá a todos Estou a tentar calcular a distancia euclidiana e da me o seguinte erro distancia<-dist(dadoss) #Euclidiana Error in dist(dadoss) : cannot allocate vector of length 569278153 alguem me pode ajudar? nao estou a entender o erro., a dimensao dos dados é dim(dadoss) [1] 33743 17 > Obrigada;)

5 4

lendo dados diretamente da web
by Antonio Silva 13 Abr '13

13 Abr '13

Olá a todos, Estou com problemas para recuperar dados da web. Por exemplo, eu utilizava o comando dat.desc1 <- read.delim("http://dl.dropbox.com/u/755659/descarga.csv ",sep=",") para baixar dados que tenho em uma pasta pública no dropbox. Depois de uns meses sem utilizar o recurso fui tentar e recebi a mensagem Erro em file(file, "rt") : não é possível abrir a conexão Outros exemplos que peguei na internet como read.table(" http://archive.ics.uci.edu/ml/machine-learning-databases/wine/wine.data ",sep=",") Funcionam normalmente. Alguém já teve este problema? Há como contornar? Obrigado, Antônio Olinto

4 7

login
by Jeuel Alves 13 Abr '13

13 Abr '13

Estou tentando logar. A princípio confirmei pelo meu celular, agora estou no computador tentando "logar". Vou verificar os cookies, mas acredito que seja o ip.

1 0

Procurando números no nome
by Fátima Lima Paula 13 Abr '13

13 Abr '13

Olá pessoal. Como faço para buscar números em uma coluna com nomes? Obrigada Fátima

4 10

Soma de quadrados tipo III
by Thiago de Paula Protásio 12 Abr '13

12 Abr '13

Pessoal; Novamente recorro a ajuda de vcs... Preciso analisar um experimento em esquema fatorial duplo em blocos casualizados com diferentes números de repetições. Mas estou tendo dificuldades em obter as médias ajustadas para cada fator, pois a interação foi não significativa. Por favor, vejam a rotina abaixo: anava1=aov(QEA~especie*adesivo+Bloco) Anova(anava1,type="III") # ajuste do modelo aditivo com teste F marginal m1 <- lm(QEA~especie*adesivo+Bloco) drop1(m1, test="F") # estimativas dos efeitos (solução) sob a restrição do R summary(m1) # obtenção das médias ajustadas dos níveis de tratamento require(contrast) lapply(levels(especie), function(i){ contrast(m1, type="average", list(especie=i, adesivo=levels(adesivo))) } ) Erro em generateData(fit = list(coefficients = c(9.49304920888927, -3.79433616881114, : not enough factors Meus sinceros agradecimentos. -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

6 6

transf Box-Cox no R
by Andre Andres 12 Abr '13

12 Abr '13

> Prezados, > > Fiz minha "inscrição" a pouco e não sei se fiz correto, mas vamos lá. > > Sou iniciante no R, e estou fazendo umas análises de regressões nonlinear , LL, testando doses de um herbicida no comprimento da raiz do arroz (cm). > > Como fiz dois experimentos, estou vendo se ambos são iguais para discutir em conjunto. Já vi que não as curvas fitted não são iguais, tanto utilizando os dados originais (cm) ou usando dados transf. em %. Ou seja, o p <0.05. > > Tentei fazer uma transformação com box-cox, mas não fui"eficiente" em usar o R. Poderia me dar uma dica como usar esta ferramenta? > > Tenho no file uma cultivar, 7 doses, root , divididos em dois estudos > O que fiz inicialmente foi testar a logistica LL.3: > cult<drm(root~dose, data=cultivar, fct=LL.3()) > > mas com os dados originais ou em % não consigo P>0.05. POr isto gostaria de testar box-cox. > > > obrigado..Andre

3 3

(sem assunto)
by Helena Turon 12 Abr '13

12 Abr '13

Estou requisitando a inscrição novamente, pois deu erro na anterior.

1 0

exportar list
by Natalia Martins 12 Abr '13

12 Abr '13

Redistas, bom dia! Como eu poderia exportar um arquivo list para txt? Att -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestra em Ciências: Area de concentração: Estatística e Experimentação Agronômica - ESALQ/ USP Doutoranda em Ciências: Area de concentração: Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

3 2

um vetor com datas
by Helena Turon 12 Abr '13

12 Abr '13

Prezados, Tenho um banco de dados diários. Consegui ver os dados faltantes (NA) nos dias em que: houve observação na estação, mas não houve registro para aquela variável. Mas, quando o dia (ou os dias) foi pulado não sei como fazer. As datas observadas já estão sendo reconhecidas pelo R como datas, dia(2), mês(2) e ano(4). Tentei fazer um vetor com as datas e comparar com as datas das observações, mas não consegui colocar datas num vetor. Atenciosamente, Helena

2 1

distância de Mahalanobis par dois bancos
by Edson Lira 11 Abr '13

11 Abr '13

Bom dia a todos, tenho dois bancos, cada um com três variáveis. O banco 1 tem 5832 obsevações e 3 variáveis.(R,G,B) O banco 2 tem 6456 obsevações e 3 variáveis.(R,G,B) Abaixo as 6 primeiras observações de cada banco BANCO 1 R G B 1 142 137 255 2 139 138 251 3 113 122 232 4 124 120 243 5 133 128 255 6 121 126 231 BANCO 2 R G B 1 144 188 255 2 163 141 255 3 146 152 255 4 156 145 250 5 130 146 241 6 188 143 241 Gostaria de calcular a distância de Mahalanobis para os dois bancos, não estou conseguindo implementar o exemplo que tem no R. ?mahalanobis Alguem podria me ajudar? [ ]' s. Edson Lira Estatístico Manaus-Amazonas

1 0

obter media dos coeficientes de modelos
by Maurício Lordêlo 11 Abr '13

11 Abr '13

Caros, No CMR abaixo crio uma matriz de delineamento que estah associada ao nro de tratamentos, nro de repeticoes de cada tratamento e o tempo que cada resposta (neste caso binaria) foi observada. Suponha que o vetor y seja gerado levando em consideracao que no tempo=1, a probabilidade eh 0.78; nos tempos>1, esta probabilidade eh 0.35. n.trat=3 # número de tratamentos n.rep=2 # número de repetições n.time=4 # número de tempos n.ind=n.trat*n.rep # número de indivíduos n=n.trat*n.rep*n.time # número total de registros id=rep(1:n.ind, each = n.time) # identificação do indivíduo time=rep(1:n.time, times = n.ind) # identificação do tempo trat=rep(factor(LETTERS[1:n.trat]), each = n.rep*n.time) # tratamentos y = matrix(data = NA, nrow = n, ncol = 1, byrow = FALSE,dimnames = NULL) data1=data.frame(id,time,trat,y) for(i in 1:n) { y[i]=ifelse(data1[i,2]<2,rbinom(1,1,.78),rbinom(1,1,.35)) } data2=data.frame(id,time,trat,y) data2 Neste caso tenho formado meu conjunto de dados (data2) já com os valores de y gerados. Agora vem a minha dificuldade. Preciso rodar um mesmo modelo repetidas vezes, como por exemplo: modelo1=geeglm(y~trat, id=id, data=data2, family=binomial) e depois obter uma media para cada um dos coeficientes (com seus erros-padrao) obtidos. Ou seja, terei outros conjuntos de dados gerados (com o mesmo nro de tratamentos, repeticoes e tempos) sendo que a unica diferenca estarah nos valores de y. Preciso obter a média dos coeficientes de cada um dos modelos rodados a partir destes diferentes y´s. Agradeco quem puder auxiliar. Mauricio

2 1

Modelo misto fatorial
by Eric Ferreira 10 Abr '13

10 Abr '13

Caros colegas, Alguém poderia me ajudar com relação à análise de um fatorial duplo, em que apenas 1 dos fatores é aleatório? Suponha dois fatores qualitativos (categóricos): A um fator aleatório e B um fator fixo, compondo um fatorial duplo em DIC. 1. Vocês sugerem o uso da função lme(), lmer() ou tanto faz? 2. Como seria a sintaxe? Eu tentei: lme(resposta ~ A*B, random= ~1|A) ...mas não obtive a resposta desejada (ele inclui o intercepto na anava e gera NaN para o p-valor de A) Tentei também: lmer(Syrliglukt ~ 1|b*Produkt) ... mas ele retorna " * not meaningful for factors". Então os fatores não podem ser categóricos?? Não entendi... Obrigado pela atenção, Eric. -- Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

2 1

erros
by alanarocha＠sapo.pt 10 Abr '13

10 Abr '13

Bom dia, alguém me sabe esplicar que erro de aviso é este? e se é possível corrigir. setwd("f:/dados panakeia e não panakeia") dados<-read.table("CustosTratamentos.csv",header=TRUE,sep=";",stringsAsFactors=FALSE,dec=",") Warning message: In read.table("CustosTratamentos.csv", header = TRUE, sep = ";", : incomplete final line found by readTableHeader on 'CustosTratamentos.csv' names(dados) [1] "PK..." dados [1] PK... <0 rows> (or 0-length row.names) e assim não me aparecem os dados. existe alguma esplicação de porque esses avisos aparecem se sim onde? Cumprimentos Ana

1 0

Re: [R-br] Soma de quadrados tipo III
by Thiago de Paula Protásio 09 Abr '13

09 Abr '13

Walmes; Muito obrigado pela ajuda! Essa função facilitou demais! Consegui implementá-la. Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

Re: [R-br] Soma de quadrados tipo III
by Thiago de Paula Protásio 09 Abr '13

09 Abr '13

Walmes; Obrigado pela dica, mas em qual pacote eu encontro essa função? -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

Curso
by Marcio Pires 09 Abr '13

09 Abr '13

Prezados Membros Gostaria de saber se existe um curso introdutório do R em Salvador-Ba, pois existe um interesse de alguns estudantes em apreender o software. Márcio

1 0

Obter valor médio dos pixels de uma imagem gerada por krigagem
by ASANTOS 09 Abr '13

09 Abr '13

Boa noite membros da lista, Gostaria de receber ideias para calcular o valor médio dos pixel em uma imagem gerada por krigagem em sub-regiões de interesse na imagem, por exemplo usando um exemplo do Prof. Paulo, onde ajusta-se um semivariograma e faz a krigagem ordinária dos dados: # require(geoR) plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2), xlab="Coord X", ylab="Coord Y") loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0, 1.1), ncol=2) text(loci, as.character(1:4), col="red") polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2) # Fitting models with nugget fixed to zero ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T) # defining the grid pred.grid <- expand.grid(seq(0,1, l=51), seq(0,1, l=51)) # kriging calculations kc <- krige.conv(s100, loc = pred.grid, krige = krige.control(obj.m = ml)) # displaying predicted values image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)), xlab="Coord X", ylab="Coord Y") # #Agora vou dividir esta área em 25 partes e lançar os centroides destas subáreas: #Centro das 25 parcela 0.2 x 0.2 p0<-c(0.1,0.1) # c0 <- seq(p0[1],p0[1]+0.9,by=0.2) d0 <- seq(p0[2],p0[2]+0.9,by=0.2) d1 <- expand.grid(x=c0, y=d0) points(d1[,1],d1[,2],col="red") ver<-c(0.0,0.2,0.4,0.6,0.8,1.0) hor<-c(0.0,0.2,0.4,0.6,0.8,1.0) abline(v=ver, col="red") abline(h=hor, col="red") # E aqui vem minha pergunta, como calcular o valor médio dos pixels contidos em cada subárea de 0.2 x 0.2? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

4 8

Rotação de matrizes com coordenadas geográficas em UTM
by ASANTOS 08 Abr '13

08 Abr '13

Bom dia pessoal, Estou cá de volta com o problema de rotacionar matrizes para rotacionar um grid criado com expand.grid(), o exemplo numérico que havia postado na lista esta OK, porém agora quando emprego coordenadas geográficas em UTM não estou conseguindo plotar, pois os valores no eixo x estão estranhos e o eixo y esta negativo, alguém saberia me dar alguma sugestão, segue CRM abaixo, obrigado # p0<-c(524237,7978856)### Primeiro ponto (Coordenada geográfica conhecida em UTM) points(p0[1],p0[2], col="red") # # c0 <- seq(p0[1],p0[1]+400,by=100)### Criando o grid d0 <- seq(p0[2],p0[2]+500,by=100) d1 <- expand.grid(x=c0, y=d0) plot(d1[,1],d1[,2]) # d2<-as.matrix(d1) ang<- 60 #### Angulo a ser rotacionado Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang), cos(ang))) ## Matriz de rotação drt<-d2 %*% Q # Multiplicação de matrizes points(drt[,1],drt[,2], col="red") # -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

3 6

Modelo não linear- correção
by Fernando A. Souza 08 Abr '13

08 Abr '13

Caros colegas, Estou rodando uma série de equações não lineares, utilizando o modelo de gompertz. A minha dúvida é a seguinte: Quais estimativas eu devo fornecer para avaliação do modelo ? Eu penso em fornecer o "erro padrão residual" e o "log-likelihood". Quais a melhores informações a serem fornecidas? Abraço a todos -- Fernando Antônio de Souza Zootecnista, Dsc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

2 3

Criando uma mascara com shapefile sobre raster
by ASANTOS 06 Abr '13

06 Abr '13

Boa tarde pessoal, Estou com um problema e não consigo criar uma mascara com um arquivo no formato shapefile sobre raster usando a função crop() do pacote raster, por algum motivo meu objeto contorno que é um SpatialLinesDataFrame não esta servido para oferecer o polígono limítrofe da área, alguém poderia me dar um help, segue CRM: require("raster") A <- stack(c("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif"))## Imagem plotRGB(A) ## Plota a imagem contorno<- readShapeLines("Catas_Altas.shp") ## Contorno da área lines(contorno, col="red")Plota o contorno C <- crop(A,contorno) ## Recorta o que esta contido no contorno na imagem Erro em .local(x, y, ...) : nenhum slot de nome "legend" para esse objeto de classe "RasterStack" -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

2 13

Expandir valor!
by Sérgio Henrique almeida da silva ju 06 Abr '13

06 Abr '13

Prezados Como faço pra expandir um valor em uma coluna? Por exemplo, no banco abaixo, ele só me dá o valor da especie para a primeira linha que é o chefe do domicilio, mas eu gostaria de expandir esse valor para todos do cod 5001102003001. cod especie 5001102003001 1 5001102003001 5001102003001 5001102003001 5001102003001 5006309003002 2 5006309003002 5006309003002 5006309003002 alguém tem alguma ideia? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 / 94429486

2 3

Modelo não linear
by Fernando A. Souza 06 Abr '13

06 Abr '13

Caros colegas, Estou rodando uma série de equações não lineares, utilizando o modelo de gompertz. A minha dúvida é a seguinte: Quais estimativas eu devo fornecer para avaliação do modelo, uma vez que o Eu penso em fornecer o "erro padrão residual" e o "log-likelihood". Quais a melhores informações a serem fornecidas? Abraço a todos -- Fernando Antônio de Souza Zootecnista, Dsc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

1 0

importar dados SRTM usando getData() do pacote raster
by Adriano S. Melo 06 Abr '13

06 Abr '13

Prezados colegas, Estou tendo dificuldade com isto: getData("SRTM", download=T, lat=-16.68, lon=-49.256) tentando a URL ' http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_27_16.ZIP' Content type 'application/x-zip-compressed' length 35447524 bytes (33.8 Mb) URL aberta ================================================== downloaded 33.8 Mb Erro em .SRTM(..., download = download, path = path) : file not found Já tentei especificar o argumento 'path' mas não funcionou. O problema acontece no Ubuntu12.04 (testei em 2 computadores). No windows 7 não tive problemas. Por outro lado, este funciona: getData('alt', country="BRA") Alguma sugestão? Abraços, Adriano S. Melo Ecologia - UFG

1 1

Utilizar o Hershey do R
by Jodavid Ferreira 05 Abr '13

05 Abr '13

Alguém saberia informar como posso utilizar o Hershey do R. Rodando o comando: demo(Hershey) estou precisando utilizar o \\#H0907 que se encontra na ultima tabela que se chama: "More Raw Hershey Escape Sequences" mas não sei como fazer! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br 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>

1 0

atualização de pacotes
by Rui Svensson 05 Abr '13

05 Abr '13

Não é muito inteligente.Quando vc aciona a atualização de pacotes, deveria estar subentendido que os pacotes já atualizados não precisam ser baixados de novo. Tal como as dependências.Isso já existe no update do antivirus: se está uptodate, não precisa atualizar.Deu pra entender ? RUI SVENSSON FONSECA analista de sistemas & analista de custos & especialista em planilhas financeiras analista de O&M & crítico estrutural compulsivoZEN

6 8

exportar um vetor do tipo list
by Aline Lipsky 04 Abr '13

04 Abr '13

Boa tarde, Eu tenho um vetor do tipo list com vetores de tamanhos diferentes. Preciso exportar do r no formato "csv" ou "txt", mas não estou conseguindo.Alguém sabe se tem como fazer essa exportação? Grata, Aline F. Lipsky

5 5

Boxplots
by Edson Lira 04 Abr '13

04 Abr '13

Caros amigos, tenho uma variiável chamada número de casos e uma outra (fator) chamada tipo de doação. Usando a rotina: boxplot(casos ~ doação,data=dc,xlab="Tipo de Doação", ylab="No. de casos",cex.axis=0.8,cex.lab=1) Tenho os dois boxplots (lado a lado) por tipo de doação (voluntário e depósito). A rotina points(mean(dc$casos),col="blue") me mostra o ponto da média no box somente em um box, como fazer para mostrar nos dois plots Edson Lira Estatístico Manaus-Amazonas

3 5

índices de linha/coluna numa variável ou função no R
by alanarocha＠sapo.pt 04 Abr '13

04 Abr '13

Boa tarde, tenho um ficheiro com umas linhas e colunas: estes são os meus dados: 2012 JAN FEV MAR ABR MAI JUN JUL AGO SET OUT NOV DEZ Custo Tratamento (Fármacos, Análises e MCDs) Não Panakeia 23,62 28,31 22,46 22,03 27,66 20,92 21,28 25,64 20,38 19,96 30,07 17,6 % de doentes com Hemoglobina >= 10 gr/dl e <= 13 gr/dl Não Panakeia 75,0% 75,8% 77,6% 77,4% 77,1% 75,0% 75,3% 77,5% 78,3% 76,8% 75,8% 76,8% estou a fazer o seguinte código: setwd("f:/dados panakeia e não panakeia") dados<-read.table("Dados.csv",header=TRUE,sep=";",stringsAsFactors=FALSE,dec=",") for(i in 1:5){ dados[i]<-as.numeric(dados[i]) for(j in 1:13) { dados[j]<-as.real(gsub('%','',dados[j])) } } Error: (list) object cannot be coerced to type 'double' colnames(dados) <- c(2012,"Jan:Dez") dados não conheço o erro. Ana

2 1

file.append
by Alexandre Ribeiro Leichsenring 04 Abr '13

04 Abr '13

Olá a todos, Estou tendo dificuldades com a função file.append. Quando faço a concatenação de dois arquivos RData usando "file.append", vejo que o tamanho do arquivo resultante (em bytes) é igual à soma dos tamanhos dos arquivos originais. Porém, quando carrego o arquivo concatenado, aparece um objeto contendo somente os elementos do primeiro arquivo: > rm(list=ls()) > # Criando arquivo 1 > (df <- data.frame(x = letters[1:5], y = 1:5)) x y 1 a 1 2 b 2 3 c 3 4 d 4 5 e 5 > save(df, file = "df1.RData") > rm(df) > # Criando arquivo 2 > (df <- data.frame(x = LETTERS[1:5], y = 10*(1:5))) x y 1 A 10 2 B 20 3 C 30 4 D 40 5 E 50 > save(df, file = "df2.RData") > rm(df) > # Criando uma cópia do arquivo 1 > file.copy(from = "df1.RData", to = "df.RData", overwrite = T) [1] TRUE > > # Concatenando df2 em df. > file.append("df.RData", "df2.RData") [1] TRUE > file.info("df1.RData") *size *isdir mode mtime ctime atime exe df1.RData *177 *FALSE 666 2013-04-03 17:19:39 2013-04-03 16:52:49 2013-04-03 16:52:49 no > file.info("df2.RData") *size *isdir mode mtime ctime atime exe df2.RData *197 *FALSE 666 2013-04-03 17:19:39 2013-04-03 16:52:50 2013-04-03 16:52:50 no > file.info("df.RData") *size *isdir mode mtime ctime atime exe df.RData *374 *FALSE 666 2013-04-03 17:19:39 2013-04-03 16:56:50 2013-04-03 16:56:50 no > load("df.RData") > df x y 1 a 1 2 b 2 3 c 3 4 d 4 5 e 5 Alguém sabe dizer o que está acontecendo? Obrigado, Alexandre

3 7

R 3.0.0
by Thiago V. dos Santos 04 Abr '13

04 Abr '13

Pessoal, Alguém já testando o R v.3? Algum problema de falta de pacote, ou incompatibilidade de pacotes? Estou curioso para testar mas ao mesmo tempo não quero comprometer trabalhos em andamento se alguma coisa parar de funcionar. Saudações, Thiago.

5 7

Mapa com lat e long
by Edson Lira 03 Abr '13

03 Abr '13

Caros amigos, tenho um banco de dados com 910 observações com 17 variáveis, sendo duas a latitude e a longitude do local da coleta. O link abaixo é uma amostra do banco com algumas observações. http://www.datafilehost.com/download-9219b795.html O link abaixo é um arquivo zipado com os dados geográficos de Manaus (bairros.shp, bairros.dbf e etc) http://www.datafilehost.com/download-87046b14.html Gostaria de plotar mapas temáticos no R, um deles seria o mapa com a variável avaliação. No mapa colocar o nome e as avaliações por bairro e depois outro por zona geográfica Como plotar mapas com variaveis qualitativas? Exemplo, sexo, escolaridade, etc. Dei uma olhada no material do Jonas, mas não consegui replicar para meu conjunto de dados. Um abraço. Edson Lira Estatístico Manaus-Amazonas

1 0

[DIVULGAÇÃO] Versão preliminar de função para ajuste de modelos não lineares via interface gráfica
by Walmes Zeviani 03 Abr '13

03 Abr '13

Saudações listeiros, Venho por meio deste anunciar a versão preliminar da função que desenvolvi como parte de rotinas relacionado à minha Tese. A função permite o usuário ajustar modelos de regressão não linear controlando a posição da função por meio de deslizadores que alteram os valores dos parâmetros até que se possa a partir de então ajustar o modelo pela nls(). A interface gráfica é baseada nas funcionalidades do pacote rpanel. Essas implementações dão possibilidade para o usuário melhor conhecer o modelo, as propriedades da função e o significado dos parâmetros. Suas aplicações didáticas são inúmeras. A função está disponível no link https://dl.dropbox.com/u/48140237/rp.nlregression.R e um CMR com 7 exemplos está disponível em https://dl.dropbox.com/u/48140237/rp.nlregression_examples.R Gostaria que os interessados corressem os exemplos, que são reproduzíveis, e apontassem sugestões. Para convidá-los à experiência, a função cria uma interface gráfica tal como mostra essa figura animada em http://r-br.2285057.n4.nabble.com/ Dentre os 7 exemplos tem-se curvas de crescimento-temperatura, de água do solo, de lactação, de produtividade-nutrição, produtividade-desfolha. A intenção seria submeter à função ao mantenedor do rpanel para que componha o conjunto de funções deste. Uma implementação anterior à essa eu disponibilizei em https://ridiculas.wordpress.com/2012/12/07/interface-grafica-para-ajuste-de… Certo de vossa colaboração, agradeço. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

3 3

problema TwitteR
by Vinicius Brito Rocha 03 Abr '13

03 Abr '13

Senhores, estou tentando usar o twitteR, especificamente a função searchTwitter e estou tendo esse erro Erro em twInterfaceObj$doAPICall(cmd, params, "GET", ...) : OAuth authentication is required with Twitter's API v1.1 alguém pode ajudar ? abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* 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.

2 6

Variante do método seleção.
by Jackeline Bonetti Campos 02 Abr '13

02 Abr '13

Olá, O método de ordenação por seleção normal (onde o menor elemento é encontrado e colocado na primeira posição e assim consecutivamente até que todo o vetor seja ordenado) é dado pelo algoritmo: selecao<-function(vet){ tam<-length(vet) for(i in 1:tam) { pm<-which.min(vet[i:tam]) pm<-pm+i-1 temp<-vet[i] vet[i]<-vet[pm] vet[pm]<-temp } return(vet)} Como posso alterar esse código para que a busca seja feita pelo maior elemento e este seja colocado na última posição? Pensei em algo assim: selecao<-function(vet){ tam<-length(vet) for(i in tam:1) { pm<-which.max(vet[tam:i]) pm<-pm+ ?? temp<-vet[i] vet[i]<-vet[pm] vet[pm]<-temp } return(vet)} Porém, não sei que atribuição dar ao novo pm para que a ultima colocação seja desconsiderada. Agradeço desde já! Att,Jackeline

3 3

Identificando pontos em um plot
by Sérgio Henrique almeida da silva ju 02 Abr '13

02 Abr '13

Olá Qual comando eu uso pra identificar o nome de um ponto em um plot, quando eu uso o identify ele me mostra a posição da variável, mas eu gostaria que ele me mostrasse a variável "NOME" e não a posição de x. Como posso fazer? Banco NOME x 1 Agreste Alagoano 9 2 Agreste Paraibano 9 3 Agreste Pernambucano 30 4 Agreste Potiguar 3 5 Agreste Sergipano 4 6 Ara\xe7atuba 1 7 Araraquara 5 8 Assis 2 9 Baixadas 0 10 Baixo Amazonas 11 plot(x) identify????? -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 / 94429486

3 3

valores em percentagem
by alanarocha＠sapo.pt 02 Abr '13

02 Abr '13

Boa tarde, eu tenho este codigo: for (i in 2:6) { # Executa o comando a seguir para as colunas 2 a 6 (ano 2008 a 2012) dados[,i] <- as.real(gsub('%','',dados[,i])) # Tira o símbolo de percentual e transforma texto em número se eu quisesse que os valores apareçam com o simbolo de percentagem faço: for (i in 2:6) { # Executa o comando a seguir para as colunas 2 a 6 (ano 2008 a 2012) dados[,i] <- as.real('',dados[,i])) # Tira o símbolo de percentual e transforma texto em número } } exemplo 30% Com cumprimentos Ana Rocha

2 1

parcela subdividida
by Juraci Mendes Moreira 02 Abr '13

02 Abr '13

ola pessoal como faço anava em parcela subdividida com bloco agradeço a atenção

3 2