- R-br - DINF UFPR

Re: [R-br] modelos mistos - funçao para nested ancova
by Ivan Bezerra Allaman 13 Abr '12

13 Abr '12

Sem a mensagem de erro fica difícil ajudar! Poste a mensagem de erro, ou melhor, nos forneça um CRM para que a ajuda seja mais efetiva. (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

plotar gráfico com valores gerados em loop
by Vitor Aguiar 13 Abr '12

13 Abr '12

Caros, estou postando abaixo um loop com o qual estou trabalhando. É longo, mas se vocês olharem para as linhas finais, verão que estou querendo fazer um plot com os resultados. Mas fazendo dessa maneira, enquanto o loop corre, o R vai produzindo um gráfico com só um ponto (referente ao cálculo de 1 corrida do loop) e substituindo-o por outro gráfico a cada corrida do loop. Gostaria de plotar todos os resultados do loop num gráfico só, ou seja, a cada corrida do loop se adicionaria um ponto a um mesmo gráfico. Como fazer? desde já agradeço, Vitor counter = 1 while (counter <= length(MarkerNames)) { matObs = as.matrix(ObsMatD[[counter]]) ExpGen = matrix(0, nrow = nrow(matObs), ncol = ncol(matObs), dimnames = list(rownames(matObs), colnames(matObs))) N = sum(matObs) for(i in 1:nrow(ExpGen)) { for(j in i:ncol(ExpGen)) { if(i==j) { ExpGen[i, j] = sum(matObs[i, ], matObs[ ,j])^2/(4*N) } else if(i!=j) { ExpGen[i, j] = sum(matObs[i, ], matObs[ ,i])*sum(matObs[j, ], matObs[ ,j])/(2*N) } } } print(names(ObsMatD[counter])) popNArm = FrmNA(pop[,counter:(counter + 1)]) Alleles = c(popNArm[,1], popNArm[,2]) AlleleFreq = table(Alleles)/length(Alleles) Hobs = 1 - sum(diag(matObs))/sum(matObs) Hexp = 1 - sum(AlleleFreq^2) Fcoef = (Hexp - Hobs) / Hexp names(Fcoef) = "F" print(Fcoef) NAratio = (1 - nrow(popNArm)/nrow(pop[,counter:(counter + 1)])) names(NAratio) = "NA ratio" print(NAratio) plot(NAratio ~ Fcoef, pch = 19) counter = counter + 1 }

2 1

Análise multivariada
by Flávia Souza 13 Abr '12

13 Abr '12

Boa tarde a todos. Alguém pode me ajudar com material/link sobre análise multivariada utilizando o "R"? Gostaria de saber principalmente sobre Análise de Componentes principais, NMDS e Analise de Correlação Canônica. Desde já agradeço a ajuda Obrigada Flávia -- MSc. Flávia Kelly S. Souza Universidade Federal do Amazonas - UFAM Faculdade de Ciências Agrárias - Depto. de Ciências Pesqueiras Manaus - Amazonas flavia.souza(a)piatam.org.br flakel.souza(a)gmail.com flaviasouza(a)ufam.edu.br skype: fla_kel

8 8

Como estimar os parametros do modelo CAR no R
by Edmundo Caetano 13 Abr '12

13 Abr '12

Prezados, Gostaria de saber como estimar os parâmetros do modelo CAR usando o R. Edmundo do R. R. Caetano

2 1

Problema em loop
by ASANTOS 12 Abr '12

12 Abr '12

Bom dia pessoal, Estou tentando dividir uma área em 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, e 2048, para tira o CV entre as parcelas e armazenar o resultado em write.table(cv2, file ="area1.mod", row.names=F,append=T). Bom, estou conseguindo calcular um a um, mais eu gostaria de otimizar e fazer um loop em um vetor que especifica o número de parcelas, como em: n.parc<-c(2,4,8,16,32,64,128,256,512,1024,2048)##Número de parcelas for(k in 1:length(n.parc)) mas não sei porque razão não estou conseguindo fazer esse loop, já tentei com n.parc e n.parc[k] e nada, segue a minha rotina completa: ####Divisão em 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048 parcelas n.parc<-c(2,4,8,16,32,64,128,256,512,1024,2048)##Número de parcelas for(k in 1:length(n.parc)){ at2<-at/k parc2<-round(sqrt(at2)) N_S=seq(p0[1], by=parc2, length.out=10)##Define o centro das parcelas L_O=seq(p0[2],by=parc2, length.out=10) grid0 <- expand.grid(N_S=seq(p0[1], by=parc2, length.out=2), L_O=seq(p0[2],by=parc2, length.out=2))#transformar em área essas medidas do centro de parcela grid=GridTopology(c(304404, 6468774), c(parc2,parc2), c(10,10)) grid=SpatialGrid(grid=grid) plot(grid) points(ninhos) lines(limite) ##Passando para objetos da classe Spatial #Pontos pontos <- SpatialPoints(cbind(ninhos$x,ninhos$y)) pontos<- SpatialPoints(list(ninhos$x,ninhos$y)) pontos<- SpatialPoints(data.frame(ninhos$x,ninhos$y)) ##Limite da área borda=SpatialPolygons(list(Polygons(list(Polygon(limite)),"limite"))) ##Plotando plot(borda)#borda points(pontos,col='grey')#Pontos lines(grid,col='blue')#Malha #remover o que esta fora do limite fullgrid(grid)=FALSE clip<-overlay(grid,borda)##Jogando o limite sobre a malha sel.grid <-grid[!is.na(clip),]##Retirando os NA's que correspondem a malha fora do limite points(sel.grid, col='blue')#ver o resultado final grafico class(sel.grid) ## SpatialPixels ap3 <- as.SpatialPolygons.SpatialPixels(sel.grid) table(over(pontos, ap3)) ##Frequencia em cada parcela z<-table(over(pontos, ap3)) #contar quantos pontos tem em cada parcela z<-as.data.frame(z) #formato somente com colunas validas cc<- SpatialPoints(cbind(grid0$N_S,grid0$L_O))# cc<- SpatialPoints(list(grid0$N_S,grid0$L_O)) cc<- SpatialPoints(data.frame(grid0$N_S,grid0$L_O)) clip2<-overlay(cc,ap3) sel.cc<-cc[!is.na(clip2),]#seleção de pontos que nao estao dentro da area #criar matriz nula coords<-NULL#matriz sem valor para criar os pontos coords<-sel.cc[1:length(sel.cc)] #numero para cada parcela, ate o final do comprimento do vetor coords<-as.data.frame(coords)#transformar em data frame Var1<-1:length(sel.cc)#criar novo objeto para comparação, referenciando a coordenada do centroide coords<-cbind(coords,Var1) #Var1= numero das parcelas myList <- list(coords, z)#tabela com qtdade dentro de coords e de z final<-merge(coords, z, by.x='Var1', by.y='Var1', all.x=TRUE, all.y=TRUE)#comando para mostrar os NA para krigagem ordinal(limitação do software, não mostra as parcelas que tem 0 ninhos) final<-final[,-1]#não precisa da 1a coluna cv2<-sqrt(var(final$Freq))/mean(final$Freq) write.table(cv2, file ="area1.mod", row.names=F,append=T) result=read.table("area1.mod",header=F) } ########################################################## Alguém poderia me dar uma help, Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 11

Erro em subset
by Tiago Vieira 12 Abr '12

12 Abr '12

Caros Colegas de R, Minha dúvida parece ser simples: Eu gero um data.frame, visualizo um valor em uma coluna e quando peço para filtrar todas as linhas em que esta coluna apresenta este valor é retornado que não há linhas. Tente dois métodos de subset e o erro persiste. Já bati muito a cabeça e não acho razões para este erro. Abaixo segue os comando e, desde já, agradeço a atenção. > fxy <- function(x,y) exp(y/x) > x <- seq(0.1, 0.5, l=101) > y <- seq(min(x^3), max(x^2), l=101) > z <- seq(0, max(fxy(x,y)), l=101) > xyz <- expand.grid(x=x,y=y,z=z) > head(xyz) x y z 1 0.100 0.001 0 2 0.104 0.001 0 3 0.108 0.001 0 4 0.112 0.001 0 5 0.116 0.001 0 6 0.120 0.001 0 > xyz[xyz$y==0.001,] [1] x y z <0 rows> (or 0-length row.names) > subset(xyz, y==0.001) [1] x y z <0 rows> (or 0-length row.names) Att, Tiago Vieira.

5 6

ReadImage
by Natália Julieta 12 Abr '12

12 Abr '12

Olá pessoa, estou com um problema na instalação do pacote ReadImage no R, ta dando incompatibilidade de versão, eu to usando o R 2.15. Gostaria de saber se alguém sabe um jeito de instalar esse pacote para essa versão. Obrigada, -- Natália Julieta Soares de Souza http://lattes.cnpq.br/4179815655707887 Graduando em Ciência da Computação Instituto de Computação - IC Universidade Federal de Alagoas – UFAL "Para usufruir da liberdade de expressão é necessário a liberdade para se informar" (Julian Assange) Linux User #475951

1 0

Algoritmo K-protótipos (cluster)
by bruno vidigal 12 Abr '12

12 Abr '12

Boa tarde, alguém sabe informar se está implementado no R o algoritmo k-prototipos proposto por HUANG (1997,1998) ? É um algoritmo que trabalha com variáveis categóricas e contínuas. Obrigado -- Bruno C. Vidigal

1 0

Duvida String
by Sérgio Henrique almeida da silva ju 11 Abr '12

11 Abr '12

Olá alguns dias atrás vocês me ajudaram na criação de variáveis (x1 e y1), onde a mesmas me mostram a posição inicial de um certo valor (A360) em uma outra variável: A variável x1 mostra que o A360 inicia-se na segunda linha da variável X e assim por diante. Só que agora eu queria criar uma outra variável x2 e y2 que extrai o valor de A360 das variáveis exemplo Se x1 = 2 eu quero uma substring de x começando da 2ª linha até a 4ª me dando dessa forma o valor A360 Se y1 = 7 eu quero uma substring de y começando da 7ª linha até a 10ª me dando dessa forma o valor A360 e assim por diante. Como faço? x y x1 y1 *A360 *A340*A360 2 7 *A330 *A360 0 2 *A340 *A330*A340*360 0 12 *A340*A360 *A340 7 0 *A330*A340*360 *A330 12 0 Obrigado -- 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) 94429486/78101651 id: 123*20942

2 11

Ajuda
by Gilbert Queiroz 10 Abr '12

10 Abr '12

Pessoal, Estou trabalhando com dados distribuídos no tempo e no espaço. Alguém pode me indicar um pacote no R que trabalhe com modelos espaço-temporais? E alguma literatura em português? Abs.

2 1

Trabalhando com datas
by Sérgio Henrique almeida da silva ju 10 Abr '12

10 Abr '12

Olá Tenho duas variáveis que são: Inicio do tratamento e fim do tratamento, como faço para criar uma terceira variável que é o tempo de tratamento em anos? Ex: x = Inicio do tratamento y = Fim do tratamento z = Tempo de tratamento em anos x y z 01/01/2010 10/12/2012 ? 04/01/2010 05/01/2012 ? 10/01/2011 06/02/2012 ? 01/02/2011 07/02/2012 ? 07/02/2011 12/12/2012 ? Se ao invés de ano eu queira calcular em dia ou em mês, como faço? Obrigado mais uma vez. -- 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) 94429486/78101651 id: 123*20942

4 4

Divisão de área irregular em subáreas de mesma área
by ASANTOS 10 Abr '12

10 Abr '12

Bom dia pessoal, Estou tentando dividir um área em subáreas com mesma área, sendo que tenho: require(sp) ##Cria o limite da área limx<-c(304404.0000,304431.0000,304500.1960,304583.0000,304671.0000,304590.0000,304404.0000,304404.0000) limy<-c(6468847.0000,6468801.0000,6468773.6057,6468881.0000,6469056.0000,6469157.0000,6468847.0000,6468847.0000) lim=cbind(limx,limy) lim<-as.matrix(lim) ##Limite da área limite=SpatialPolygons(list(Polygons(list(Polygon(lim)),"lim"))) ##Plotando plot(limite)#Limite Sobre a qual gostaria de criar um grid, subdividindo a área em 4, 8, 16 e 32 partes, sendo que as subdivisões tem que ter a mesma área. Para tanto, gostaria de saber se alguém conhece alguma função em algum pacote que faz isto. Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0

Mapas Município
by Vinicius Brito Rocha 10 Abr '12

10 Abr '12

Senhores, onde posso encontrar shapes dos mapas dos municípios do Estado do RJ? 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.

6 6

[OFF-TOPIC] Sobre Outras Listas
by Listeiro 037 10 Abr '12

10 Abr '12

Bom dia a todos. Considero esta pergunta 99% off-topic porque este 1% refere-se ao que mais é necessário e feito para se fazer para o uso do R. Alguém recomendaria alguma lista de discussão sobre Linux, desde que não seja "oficial" de distro e também alguma de programação e desenvolvimento própria neste ambiente? O que mais está na moda são fóruns (não!) e o Yahoo! Groups está sucateado, virou repositório de spams. A antiga "prog-br" da Conectiva está parada há meses. Agradeço se alguem tiver algumas opções. Até mais.

1 0

erro em ggmap()
by D'Apolito 09 Abr '12

09 Abr '12

Olá pessoal, Se alguém souber desse erro que vou relatar e puder contribuir com alguma solução, agradeço. Estou usando Windows 7, e o R é o 2.15, 32-bit. quando uso o seguinte código: library(ggmap) hdf <- ggmap() me aparece o erro "Erro em readPNG(destfile, native = native) : file is not in PNG format". Procurei na internet e mais gente já teve o problema, mas não achei nenhuma solução. Quando eu rodava a mesma linha de código na versão anterior do R, funcionava. Já tentei inutilmente reinstalar o pacote. Desde já agradeço Carlos

1 0

listar funções de determinados pacotes
by Lucas Barbosa 09 Abr '12

09 Abr '12

Boa tarde pessoal, Preciso de uma listagem de todas as funções de vários pacotes (parallel, snow, multicore, bigmemory, synchronicity...). Alguma sugestão de como posso conseguir? Obrigado! Lucas Cusinato

4 11

dúvida: Mapas Temáticos
by Jean Carlos Quaresma Mariano 09 Abr '12

09 Abr '12

Bom Dia. Trabalho na UNESP de Ilha Solteira na Área de Hidráulica e Irrigação, temos uma rede de estações agrometeorológicas que medimos vários parâmetros tipo Temperatura, Umidade, Precipitação etc.. estes dados estão disponíveis em nosso Portal do CLIMA http://clima.feis.unesp.br o próximo passo do nosso projeto seria de disponibilizar mapas de interpolação exatamente com estes http://www.agritempo.gov.br/Editor/UserFiles/Image/RSE/RSE.BOLETIM_TMAX_Prev .GIF mas com a nossa área de estudos que é o noroeste paulista, a questão é alguém já teve alguma experiência na confecção deste mapas, algum exemplo, material? Obrigado Att Jean Carlos Quaresma Mariano Analista de Sistemas / Área de Hidráulica e Irrigação Faculdade de Engenharia de Ilha Solteira UNESP - Universidade Estadual Paulista www.agr.feis.unesp.br/irrigacao.php ( (18) 3743-1959 Skype: jean_quaresma MSN: jeanquaresma(a)gmail.com Twitter: @jean_quaresma Descrição: Descrição: assinatura AVISO LEGAL: As informações contidas nesta mensagem e seus anexos são CONFIDENCIAIS, protegidas por lei e por direitos autorais. A divulgação, distribuição, reprodução ou qualquer forma de utilização do teor destas informações depende de autorização do emissor, sujeitando-se o infrator às sanções legais. Caso tenha recebido esta mensagem por engano, favor avisar imediatamente o emissor respondendo esta mensagem e apagando-a, assim como todas as suas cópias e anexos.

3 2

Publicações da Área
by Listeiro 037 09 Abr '12

09 Abr '12

Olá pessoal! Cheguei agora. Gostaria de saber se existem publicações sobre a R. Desde já agradeço.

2 1

dúvidas sobre tratamentos equidistantes
by Fernando Antonio de souza 09 Abr '12

09 Abr '12

Caros amigos Gostaria de esclarecimento sobre a seguinte dúvida. Há necessidade de os tratamentos serem equidistantes para comparação de tratamentos e/ou ajuste de modelo?E no caso de variáveis instáveis?

2 1

Re: [R-br] Tobit espacial
by Diego andre 09 Abr '12

09 Abr '12

Wagner Hugo Bonat estava querendo modelar os tres tipos de dependencia: no erro (modelo SEM), na variavel dependente (modelo SAR) e em ambos (SARMA), quanto a referencia, no Lesage. -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

2 1

Problema de Memória em função
by João José de Simoni Gouveia 08 Abr '12

08 Abr '12

Prezados, Escrevi uma simples e tupiniquim função no R (em anexo com o nome ld.R) para umas análises que preciso rodar. Nos meus arquivos de teste (que são subset dos arquivos que terei que utilizar) a função está funcionando direitinho. O problema é que quando tento utilizar a função nos arquivos de trabalho (são 26 arquivos, um para cada cromossomo e cada arquivo com centenas de milhares de linhas) o R simplesmente trava depois de algum tempo. A impressão que me dá é que "estoura" a memória do meu notebook (um core i3 com 4 de ram, rodando ubuntu) Gostaria de saber se alguém teria alguma sugestão para que eu pudesse fazer isto funcionar... Anexado ao email envio a função e três arquivos teste (teste_CHR1.haplo.LD, teste_CHR2.haplo.LD, teste_CHR3.haplo.LD) A função está comentada. para rodar a função basta dar o comando: > ld('teste',3) Att JJ -- ============================================ 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 ============================================

1 0

R 2.15.0
by Benilton Carvalho 08 Abr '12

08 Abr '12

Comecem a preparar seus upgrades... e reclamacoes (ja que ha' muitos que querem pq querer atualizar o R sem atualizar seus pacotes): https://stat.ethz.ch/pipermail/r-announce/2012/000551.html benilton

8 9

Notepad++ e NppToR
by Leandro Marino 07 Abr '12

07 Abr '12

*Caros,* * * *alguém sabe como fazer a configuração de forma a tudo funcionar através de um pendrive?* * * *Não encontrei nenhum tutorial na internet que me ajudasse nisto...* * * *Abraços.* *Leandro* * *Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

2 7

dúvida simples: dar nome a variáveis dentro de um loop
by Vitor Aguiar 06 Abr '12

06 Abr '12

Caros senhores, tenho um loop com o seguinte formato: While (condição) { for (i in 1:n) { calculo } X2sum = sum(X2) print(X2sum) dfreedom = function(x) (nrow(x)*(nrow(x) - 1))/2 df = dfreedom(X2) print(df) pvalue = 1 - pchisq(sum(X2), dfreedom(X2)) print(pvalue) NAnumber = nrow(pop[,counter:(counter + 1)]) - nrow(FrmNA(pop[,counter:(counter + 1)])) print(NAnumber) } Esses objetos no final (X2sum, df, pvalue, NAnumber) não são impressos no R com o nome deles, mas do seguinte modo: [1] resultado [1] resultado [1] resultado [1] resultado Como eu faço para imprimir com os nomes ao invés desses [1]? Obrigado, Vitor Aguiar

2 2

Re: [R-br] dúvida function fts to create functional time series objects
by Mauro Sznelwar 06 Abr '12

06 Abr '12

Já achei! Este arquivo OCTOPUS está disponível em algum lugar? Novamente, um outro exemplo: dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porAno = dcast(dadosOriginais, MONTH ~ YEAR, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porAno) y0 = as.matrix(porAno[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0, plot.type='functions', plotlegend=TRUE) altere o codigo acima para os seus casos especificos. benilton

1 0

dúvida function fts to create functional time series objects
by carlos pombo sonderblohm 06 Abr '12

06 Abr '12

Caros colegas, Estou a tentar fazer uns graficos de series temporais com a libraria rainbow, a qual precisa de ter os dados arranjados numa functional time series (FTS). Eu ja sei transformar os dados numa ts normal, mas agora gostaria de mostrar todo num só gráfico com esta libraria rainbow: 1.- o comportamento da serie temporal por meses (no eixo x) dos diferentes barcos a longo dos anos, 2.- o comportamento dos desembarques por portos de desembarque, 3.- simplesmente os desembarques mensais por ano para todos os barcos e portos juntos (diferentes cores para cada ano). Os dados que tenho, ficheiro OCTOPUS.csv, são desembarques em kilogramos (QVENDA) de polvo, em preço (QVALOR, €), porto de desembarque (IPORTO), codigo de cada embarcação (CodEmbarc), DAY, MONTH, YEAR. Realmente não consigo perceber a ficheiro de ajuda quando escrevo ?fts, alguma ajuda disponivel? data= http://www.datafilehost.com/download-5bbb581a.html obrigado pelo tempo disponivel, Carlos -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

3 5

parametros booleano e numerico
by Cleber N.Borges 04 Abr '12

04 Abr '12

olá a todos, alguem sabe o motivo do valor 1 não ter o mesmo efeito de TRUE (booleano) na função memory.size? fiz alguns testes e acabei com essa duvida. Grato Cleber --------------------------- > memory.size() [1] 23.86 > memory.size(max=TRUE) [1] 25.62 > memory.size(max=FALSE) [1] 23.86 > > memory.size(max=1) [1] 7935.37 Warning message: In memory.size(max) : cannot decrease memory limit: ignored > memory.size(max=0) [1] 7935.37 Warning message: In memory.size(max) : cannot decrease memory limit: ignored > > memory.size function (max = FALSE) round(.Internal(memory.size(max)), 2L) <bytecode: 0x00000000054d5738> <environment: namespace:utils> > > memory.size(max=(1==TRUE)) [1] 25.62 > memory.size(max=(0==TRUE)) [1] 23.92 > > 23==TRUE [1] FALSE > 1==TRUE [1] TRUE > 23==TRUE [1] FALSE > > 0==FALSE [1] TRUE > 2==FALSE [1] FALSE > > sessionInfo() R version 2.15.0 (2012-03-30) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base >

3 4

DUVIDA ANOVA
by Kalill V.Páscoa 04 Abr '12

04 Abr '12

Nos meus dados existem 13 tratamentos porém quando gero a ANOVA ela indica que existe apenas 1, onde estou errando? > dados$Trat [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [41] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 [81] 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 7 7 7 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 [121] 8 8 8 8 8 8 8 8 8 8 8 8 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 10 10 10 10 10 10 10 10 [161] 10 10 10 10 10 10 10 10 10 10 10 10 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 12 12 12 12 12 12 12 12 [201] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13> mDAP <- aov(DAP~Bloco+Trat, data=dados)> class(mDAP)[1] "aov" "lm" > anova(mDAP)Analysis of Variance Table Response: DAP Df Sum Sq Mean Sq F value Pr(>F) Bloco 1 0.0 0.022 0.0013 0.970882 Trat 1 114.0 113.992 6.7906 0.009754 ** Residuals 233 3911.3 16.787 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 -- *Kalill José Viana da Páscoa* *Engenheiro Florestal * *Mestrando em Engenharia Florestal Universidade Federal de Lavras* * * *Rua Geraldo dos Passos, 81 * *Água Limpa, **Lavras-MG* *CEP: 37200-000* *(35)3821-6833* *(35)9152-4828* *(35)9928-6574*

2 2

PROBLEMA NA ANOVA E NO SKOTT-KNOTT
by Kalill V.Páscoa 04 Abr '12

04 Abr '12

Boa noite, estou tentando fazer a Anova e o teste Sk porém no fim ele gera um erro, alguém faz ideia do que seja? #ANAVA# mDAP <- aov(DAP~Bloco+Trat, data=dados) class(mDAP) anova(mDAP) summary(mDAP) summary.lm(mDAP) #TESTE DE NORMALIDADE# shapiro.test(dados$DAP) shapiro.test(residuals(mDAP)) library(ScottKnott) #TESTE SCOTT-KNOTT# as.factor(dados$Bloco) as.factor(dados$Trat) sk.DAP=SK(x=dados, y=dados$DAP, model="dados$DAP~dados$Bloco+dados$Trat",which="dados$Trat",sig.level=0.05) summary(sk.DAP) Erro em SK.aov(x = av, which = which, id.trim = id.trim, sig.level = sig.level) : No factors in the fitted model!Além disso: Mensagens de aviso perdidas:1: In replications(paste("~", xx), data = mf) : non-factors ignored: dados$Bloco2: In replications(paste("~", xx), data = mf) : non-factors ignored: dados$Trat -- *Kalill José Viana da Páscoa* *Engenheiro Florestal * *Mestrando em Engenharia Florestal Universidade Federal de Lavras* * * *Rua Geraldo dos Passos, 81 * *Água Limpa, **Lavras-MG* *CEP: 37200-000* *(35)3821-6833* *(35)9152-4828* *(35)9928-6574*

1 0

(sem assunto)
by Estevão Prado 04 Abr '12

04 Abr '12

Olá pessoal, Preciso carregar o mapa de uma cidade no R, já baixei os arquivos (.dbf, .shp...), porém quando vou carregar o arquivo ele dá erro. require(maptools) getwd() setwd("/home/est/Área de Trabalho/maps") city <- readShapePoly("maps/4105805") # 4105805 é código da cidade e nome dos arquivos (.dbf, .shp...) Dei uma olhada no site: http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:seufop11, mas não consegui obter êxito. Se alguém puder me ajudar desde já agradeço.

5 6

Re: [R-br] abrir mapa shp
by Estevão Prado 04 Abr '12

04 Abr '12

Alexandro, Com o arquivo que você me enviou consegui visualizar. Obrigado. Date: Tue, 3 Apr 2012 17:33:49 -0700 From: vl.alexandro(a)yahoo.com.br Subject: Re: [R-br] abrir mapa shp To: estevao.prado(a)hotmail.com Segue anexo o shape file de spse vc abri-lo isto significa que o seu shape file está ... vamos testar esta hipotése!tudo bem! De: Estevão Prado <estevao.prado(a)hotmail.com> Para: Lista R-br <r-br(a)listas.c3sl.ufpr.br> Enviadas: Terça-feira, 3 de Abril de 2012 20:36 Assunto: Re: [R-br] abrir mapa shp Alexandro, Tentei o que você me sugeriu, mas o mesmo erro ainda esta ocorrendo. (Erro em getinfo.shape(filen) : Error opening SHP file). Mas obrigado mesmo assim. Date: Tue, 3 Apr 2012 15:23:57 -0700 From: vl.alexandro(a)yahoo.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] abrir mapa shp Tente o seguinte: library(maptools)cidade<- readShapePoly(file.choose()) #procure onde está o arquivo shp #espero ter ajudado!!!Alex De: Estevão Prado <estevao.prado(a)hotmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Terça-feira, 3 de Abril de 2012 18:49 Assunto: [R-br] (sem assunto) Olá pessoal, Preciso carregar o mapa de uma cidade no R, já baixei os arquivos (.dbf, .shp...), porém quando vou carregar o arquivo ele dá erro. require(maptools) getwd() setwd("/home/est/Área de Trabalho/maps") city <- readShapePoly("maps/4105805") # 4105805 é código da cidade e nome dos arquivos (.dbf, .shp...) Dei uma olhada no site: http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:seufop11, mas não consegui obter êxito. Se alguém puder me ajudar desde já agradeço. _______________________________________________ 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 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Re: [R-br] Tobit espacial
by Diego andre 03 Abr '12

03 Abr '12

Caro Leonard de Assis, nestes pacotes que você citou tem o tobit, mas não na sua forma espacial (SPATIAL TOBIT), que é o que eu estou procurando. Mais uma vez, obrigado pela disposição de ajudar. -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

3 2

Re: [R-br] Intervalo de Confiança da Predição
by Danilo Scorzoni Ré 03 Abr '12

03 Abr '12

Perfeito! Deu certo, mto obrigado Walmes! Danilo 2012/4/3 <r-br-request(a)listas.c3sl.ufpr.br> > 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. Matriz de proximidade (Gilbert Queiroz) > 2. Re: Matriz de proximidade (Paulo J Ribeiro Jr) > 3. Matriz de proximidade (Gilbert Queiroz) > 4. Intervalo de Confiança da Predição (Danilo Scorzoni Ré) > 5. Re: Intervalo de Confiança da Predição (Walmes Zeviani) > 6. Re: Problemas com fixed breaks (Mauro Sznelwar) > 7. Re: Tobit espacial (Wagner bonat) > 8. Re: Matriz de proximidade (Wagner bonat) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 2 Apr 2012 12:07:08 -0700 (PDT) > From: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Matriz de proximidade > Message-ID: > <1333393628.70157.YahooMailNeo(a)web160603.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Amigos, > É possível gerar a matriz de proximidade no R, sem ter que importar de um > outro software, como o GeoDA ou ArcGis? > Quais os comandos? > Abs. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/20f4ab01/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 02 Apr 2012 16:26:47 -0300 > From: Paulo J Ribeiro Jr <paulojus(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br, Gilbert Queiroz > <gilbert_queiroz(a)yahoo.com.br> > Subject: Re: [R-br] Matriz de proximidade > Message-ID: <1333394807.16481.2.camel@atikum> > Content-Type: text/plain; charset="UTF-8" > > de uma olhada no pacote spdep e em particular em > > poly2nb() > nb2listw() > > > > Em Seg, 2012-04-02 às 12:07 -0700, Gilbert Queiroz escreveu: > > Amigos, > > É possível gerar a matriz de proximidade no R, sem ter que importar de > > um outro software, como o GeoDA ou ArcGis? > > Quais os comandos? > > Abs. > > > > _______________________________________________ > > 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. > > -- > Paulo Justiniano Ribeiro Jr > LEG (Laboratorio de Estatistica e Geoinformacao) > Universidade Federal do Parana > Caixa Postal 19.081 > CEP 81.531-990 > Curitiba, PR - Brasil > Tel: (+55) 41 3361 3573 > VOIP: (+55) (41) (3361 3600) 1053 1066 > Fax: (+55) 41 3361 3141 > e-mail: paulojus AT ufpr br > http://www.leg.ufpr.br/~paulojus > > > > > ------------------------------ > > Message: 3 > Date: Mon, 2 Apr 2012 12:25:52 -0700 (PDT) > From: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Matriz de proximidade > Message-ID: > <1333394752.69772.YahooMailNeo(a)web160604.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Amigos, > É possível gerar a matriz de proximidade no R, sem ter que importar de um > outro software, como o GeoDA ou ArcGis? > Quais os comandos? > Abs. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/b9fdfefa/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 2 Apr 2012 17:37:35 -0300 > From: Danilo Scorzoni Ré <danilo.scorzoni(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Intervalo de Confiança da Predição > Message-ID: > <CA+bJVrd5TuHhwPPjEYWngLKG7b+VaxFj20MHge4MUrDky0js-w(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Prezados, > > Estou ajustando um conjunto de dados ao modelo não-linear de Gompertz (G ~ > a * exp(b * exp(c * Idade))), onde a, b e c são parâmetros. > Consegui realizar o ajuste pelo comando nls() sem problema algum. O que eu > gostaria de obter agora são os intervalos de confiança da predição. Sei que > nos modelos lineares, eu consigo obter através da função predict(), > passando a informação para o parâmetro interval="prediction", no entanto, > ao que me parece, essa opção ainda não está implementada para o comando > nls(), pois quando eu faço o comando > > ic = predict(model2, newdata=g, interval="prediction", level=0.95) > > Ele me retorna apenas a resposta média, e não os intervalos de confiança da > predição. Testei a mesma sequência de comandos com o lm() e funcionou > corretamente. > Alguém sabe se existe outra maneira de estimar esses intervalos? > > Segue o comando mínimo reproduzível: > > > -------------------------------------------------------------------------------------------------------------------------- > > # Código Mínimo Reproduzível > > x = seq(0,15,0.1) > erros = rnorm(151, 0, 3) > y = 22 * exp(-3 * exp(-0.3 * x)) + erros > plot(y ~ x) # Comportamento não-linear > > # Predição e Intervalos de Confiança no Modelo Linear > model.linear = lm(y ~ x + I(x^2)) > summary(model.linear) > xest = as.data.frame(seq(0,15,1)) > colnames(xest) = "x" > pred = as.data.frame(predict(model.linear, newdata=xest, > interval="prediction")) > lines(pred$fit ~ xest$x, lwd=2) > lines(pred$lwr ~ xest$x, lty=2, col=2) > lines(pred$upr ~ xest$x, lty=2, col=2) > > # Predição e Intervalos de Confiança no Modelo Não-Linear > model.nlinear = nls(y ~ a * exp(b * exp(c * x)), start=list(a=25, b=-4, > c=-0.3)) > summary(model1) > > plot(y ~ x) > xest = as.data.frame(seq(0,15,1)) > colnames(xest) = "x" > pred = predict(model.nlinear, newdata=xest, interval="prediction") > str(pred) > lines(pred ~ xest$x, lwd=2) > > -------------------------------------------------------------------------------------------------------------------------- > > Desde já agradeço o empenho! > > Abraços, > > -- > Danilo Scorzoni Ré > Engenheiro Florestal > Mestre em Ciência Florestal > (14) 8180-2494 > http://about.me/dscorzoni > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/98c19b32/atta… > > > > ------------------------------ > > Message: 5 > Date: Mon, 2 Apr 2012 19:07:26 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Intervalo de Confiança da Predição > Message-ID: > <CAFU=EkZ+3nrxZjTuOmXLbkzsKuQJXxyzBXiR9Oc6GB6C=8r2Mw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Veja se isso pode ser útil > > > http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-d… > > À disposiçã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 > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/20f5e4cd/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 2 Apr 2012 22:08:08 -0300 > From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Problemas com fixed breaks > Message-ID: <C90D1BEEB10647F8A610AB827EBF250D@ACER4520> > Content-Type: text/plain; charset="iso-8859-1" > > > Tentei rodar de novo, mas deu dois problemas: > > library(maptools) > Carregando pacotes exigidos: foreign > Carregando pacotes exigidos: lattice > Checking rgeos availability: FALSE > Note: when rgeos is not available, polygon geometry > computations in maptools depend on gpclib, > which has a restricted licence. It is disabled by default; > to enable gpclib, type gpclibPermit() > > e depois > > > mapa = readShapePoly("saudemapa.shp") > Erro em getinfo.shape(filen) : Error opening SHP file > > E ainda faltou disponibilizar o arquivo > logo_pb.png > > > Em 1 de abril de 2012 21:05, Mauro Sznelwar <sznelwar(a)uol.com.br> > escreveuo: > > Estava tentando rodar este código, e vi que falta este > arquivoc:/users/Julianna/Documents/SecSaude/saudemapa.shp > Daria para enviar?. > > Mensagem abaixo em nome da Julianna... (dados disponiveis via > datafilehost) > > > ==== > > > E quando fui conferir os breaks, os intervalos estavam certos,mas o > mapa resultante estava errado. > Porque quando eu não importo o csv e uso o dataframe resultante da > função, os breaks ficam ok(da maneira que eu especifico em FixedBreaks) e, > do contrario, não??? > > > segue o código: > > > ####################################################### > # CONSTRUINDO MAPAS > ####################################################### > ##Code (Comments are preceded by ##) > ## Load required packages > > > # Limpando memoria. > rm(list=ls(all=TRUE)) > > > library(maps) > library(sp) > library(maptools) > library(spdep) > library(classInt) > library(RColorBrewer) > library(shape) > library(SDMTools) > library(latticeExtra) > library(ReadImages) > library(png) > > > rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0, > cols=1, cex=0.6,...) > { > cols <- rep(c("white","black"),8) > radii <- rep(size/c(1,4,2,4),4) > x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] > y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] > for (i in 1:15) > { > x1 <- c(x[i],x[i+1],loc[1]) > y1 <- c(y[i],y[i+1],loc[2]) > polygon(x1,y1,col=cols[i]) > } > polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) > b <- c("O","S","L","N") > for (i in 0:3) > text((size+par("cxy")[1]-0.40)*cos(bearing+i*pi/2)+loc[1], > > (size+par("cxy")[2]-0.42)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) > } > > > ## Lendo dados > ## http://www.datafilehost.com/download-e5cb7e86.html > dados = read.table("dadosJulianna.csv", sep = ",", header = TRUE) > > > # Lendo mapa. > mapa = > readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp") > proj4string(mapa) <- CRS("+init=epsg:4291") > scaleXY(mapa, 1836656) > > > > > # arquivo png > #img <- readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # > seu arquivo png aqui > #set.seed(123) > #par(xpd = TRUE) > > > > > colors <- brewer.pal(7, "Reds") > brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", > intervalClosure="right", > fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, > 10000), main="Fixed") > brks<- brks$brks > > > ID = > match(substr(mapa$CODIGO,1,6),dados$t_indicadores_municipio.f_municipio) > dados = dados[ID,] > > > png(filename="c:/mif2006.png", height=750, width=1500) > > > > > plot(mapa, col = colors[findInterval(dados$ano_2006, > brks,all.inside=TRUE)], axes = TRUE, cex = 1.5, mar=5) > ##, xlim=c(-39.7, -34.7), ylim=c(-8.5, -6) > > > ##invisible(text(getSpPPolygonsLabptSlots(mapa), > labels=as.character(mapa$NOMEMUN_1), cex=0.65)) > > > title(main = "PORCENTAGEM DE Ã BITOS INVESTIGADOS\ntipo: mulheres em > idade fÃ©rtil", > font.main = 2, cex.main = 1.8) > > > text(-35.7, -6.2, "ANO: 2006", cex = 1.8) > > > legend(-39.6,-7.6, c("Sem informaÃ§Ã£o", "0", "]0,20]", "]20,40]", > "]40,60]", "]60,80]","Acima de 80"), > fill = colors, bty="n", x.intersp = 1.1, y.intersp = 1.1, cex > = 1.6) > > > rosa_dos_ventos(loc=c(-34.5,-8),size=0.17,cex=1.20, bearing=0, cols=1) > > > SpatialPolygonsRescale(layout.scale.bar(), offset = c(-35,-8.30), > scale = 1, fill = c("transparent", "black"), > plot.grid= FALSE) > > > text(-34.5, -8.35, "110.3 KM", cex= 1.3) > > > #rasterImage(img, xleft = -39.5, ybottom = -6.25, xright = -38.9 , > ytop = -6) > > > dev.off() > > > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/e7f0095e/atta… > > > > ------------------------------ > > Message: 7 > Date: Tue, 3 Apr 2012 11:42:26 -0300 > From: Wagner bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Tobit espacial > Message-ID: > <CANt=4MietnxFQRVKX2FWTXnvNFVRt4NXrjtnBMLkPFMDFKb4=w(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Que tipo de dependência espacial vc quer modelar ? Onde encontrou este > modelo ? Pode não ter programado, mas com certeza pode ser programado. > > > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120403/2b5e44e9/atta… > > > > ------------------------------ > > Message: 8 > Date: Tue, 3 Apr 2012 11:43:55 -0300 > From: Wagner bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Gilbert Queiroz > <gilbert_queiroz(a)yahoo.com.br> > Subject: Re: [R-br] Matriz de proximidade > Message-ID: > <CANt=4MhG4Gug6tLEJLX-7PKb+_VP0SsQsYDi4B6iP_SjPK0JiA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Sim, o pacote spdep tem várias funções para montar matrizes de vizinhança > seguinte diveros critérios e de acordo com o tipo de dados que vc lê. > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120403/94bcdfd9/atta… > > > > ------------------------------ > > _______________________________________________ > 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 14, assunto 3 > **************************************** > -- Danilo Scorzoni Ré Engenheiro Florestal Mestre em Ciência Florestal (14) 8180-2494 http://about.me/dscorzoni

1 0

Matriz de proximidade
by Gilbert Queiroz 03 Abr '12

03 Abr '12

Amigos, É possível gerar a matriz de proximidade no R, sem ter que importar de um outro software, como o GeoDA ou ArcGis? Quais os comandos? Abs.

2 1

Re: [R-br] Problemas com fixed breaks
by Julianna Trindade 02 Abr '12

02 Abr '12

Obs: as faixas de valores com "n" casas decimais foram colocadas em função do erro que mostrei no email anterior, no entanto, ao selecionar o dataframe brks, vcs verão que 2 valores das classes de intervalos vem com erro. 2012/3/30 Julianna Trindade <julianna(a)jubalitpb.com> > arquivos em anexo > > > > Em 29 de março de 2012 21:33, Benilton Carvalho < > beniltoncarvalho(a)gmail.com> escreveu: > > Julianna, >> >> vc precisa dar-nos algo que possamos reproduzir para que possamos ajudar >> (esse e' o motivo de vc nao ter obtido resposta ao enviar esta mesma msg >> dois dias atras). >> >> Assumindo que 'dados' nao tenha muito mais que 30-50 linhas, seria muito >> util que vc postasse o resultado de: >> >> dput(dados$ano_2006) >> >> e tbm nos dissesse qual eh o pacote que te da' a funcao classIntervals(), >> afinal nao e' algo do base. >> >> benilton >> >> 2012/3/29 Julianna Trindade <julianna(a)jubalitpb.com> >> >>> Boa Tarde, >>> >>> Gostaria de uma solução para o FixedBreaks obedecer os limites >>> superiores e inferiores das minhas classes de intervalos. No entanto, ele >>> deve associar algum método estatístico e burla minha configuração. >>> >>> Eis o código que eu utilizo: >>> >>> brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", >>> intervalClosure="right", fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000), >>> main="Fixed") >>> brks<- brks$brks >>> >>> Eis o breaks gerado pelo R: >>> >>> c(-11.1, -0.5, 16.665, 41.665, 58.335, 83.335, 110.1) >>> >>> Alguém pode me ajudar, por favor? >>> >>> >>> >>> >>> -- >>> Julianna Trindade >>> >>> "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu >>> poder." >>> >>> Salmos 21; 13 >>> >>> _______________________________________________ >>> 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. >> > > > > -- > Julianna Trindade > > "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." > > Salmos 21; 13 > -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

4 10

Intervalo de Confiança da Predição
by Danilo Scorzoni Ré 02 Abr '12

02 Abr '12

Prezados, Estou ajustando um conjunto de dados ao modelo não-linear de Gompertz (G ~ a * exp(b * exp(c * Idade))), onde a, b e c são parâmetros. Consegui realizar o ajuste pelo comando nls() sem problema algum. O que eu gostaria de obter agora são os intervalos de confiança da predição. Sei que nos modelos lineares, eu consigo obter através da função predict(), passando a informação para o parâmetro interval="prediction", no entanto, ao que me parece, essa opção ainda não está implementada para o comando nls(), pois quando eu faço o comando ic = predict(model2, newdata=g, interval="prediction", level=0.95) Ele me retorna apenas a resposta média, e não os intervalos de confiança da predição. Testei a mesma sequência de comandos com o lm() e funcionou corretamente. Alguém sabe se existe outra maneira de estimar esses intervalos? Segue o comando mínimo reproduzível: -------------------------------------------------------------------------------------------------------------------------- # Código Mínimo Reproduzível x = seq(0,15,0.1) erros = rnorm(151, 0, 3) y = 22 * exp(-3 * exp(-0.3 * x)) + erros plot(y ~ x) # Comportamento não-linear # Predição e Intervalos de Confiança no Modelo Linear model.linear = lm(y ~ x + I(x^2)) summary(model.linear) xest = as.data.frame(seq(0,15,1)) colnames(xest) = "x" pred = as.data.frame(predict(model.linear, newdata=xest, interval="prediction")) lines(pred$fit ~ xest$x, lwd=2) lines(pred$lwr ~ xest$x, lty=2, col=2) lines(pred$upr ~ xest$x, lty=2, col=2) # Predição e Intervalos de Confiança no Modelo Não-Linear model.nlinear = nls(y ~ a * exp(b * exp(c * x)), start=list(a=25, b=-4, c=-0.3)) summary(model1) plot(y ~ x) xest = as.data.frame(seq(0,15,1)) colnames(xest) = "x" pred = predict(model.nlinear, newdata=xest, interval="prediction") str(pred) lines(pred ~ xest$x, lwd=2) -------------------------------------------------------------------------------------------------------------------------- Desde já agradeço o empenho! Abraços, -- Danilo Scorzoni Ré Engenheiro Florestal Mestre em Ciência Florestal (14) 8180-2494 http://about.me/dscorzoni

2 1

Matriz de proximidade
by Gilbert Queiroz 02 Abr '12

02 Abr '12

Amigos, É possível gerar a matriz de proximidade no R, sem ter que importar de um outro software, como o GeoDA ou ArcGis? Quais os comandos? Abs.

2 1

[dúvida] colocar única legenda em 4 plots no mesmo frame
by Marcelo Claro de Souza 31 Mar '12

31 Mar '12

Plotei quatro figuras em um mesmo frame utilizando: layout(matrix(1:4,nrow=2,ncol=2)) plot(x1~y1,data=dados) plot(x2~y2,data=dados) plot(x3~y3,data=dados) plot(x4~y4,data=dados) Eu gostaria de colocar uma única legenda nos eixos x e y para representar os eixos das 4 figuras ao mesmo tempo, sem a necessidade de chamar uma legenda para cada figura evitando redundâncias. Alguém teria alguma sugestão? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

Somatório com matrizes
by Guilherme Biz 31 Mar '12

31 Mar '12

Bom dia Senhores, exite alguma função pronta para fazer somatório de matrizes? Estou tentando somar, dentro de um loop, muitas matrizes e não encontrei uma função para este somatório. Vou ter que fazer outro loop para estas somas? Ex Tenho as matrizes A1,......A100. e quero fazer A1+A2+....+A100. Att, Guilherme Biz

3 3

Re: [R-br] Tobit espacial
by Diego andre 31 Mar '12

31 Mar '12

Obrigado Leonard de Assis pela resposta, mas infelizmente não tem *tobit espacial* neste pacote... -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

2 1

duvida integers e as.double
by Vitor Aguiar 30 Mar '12

30 Mar '12

Olá Benilton. Por favor, gostaria de tirar uma dúvida a respeito de um código que você me passou há uns dias. É esse código abaixo. grupos <- unique(gsub("(.*)\\.\\d{1}$", "\\1", names(Pop))) tabelaPorGrupo <- function(grp, dat){ cols <- grep(grp, names(dat)) Nomes <- levels(factor(c(dat[,cols[1]], dat[,cols[2]]))) table(factor(Pop[,cols[1]], levels=Nomes), factor(Pop[,cols[2]], levels=Nomes)) } tabelas <- lapply(grupos, tabelaPorGrupo, Pop) names(tabelas) <- grupos tabelas post: http://r-br.2285057.n4.nabble.com/R-br-problema-com-for-loop-td4465197.html Tenho que fazer muitos cálculos com os valores gerados por esse código e, algumas vezes, quando estou trabalhando com valores altos, há warning messages dizendo "NAs produced by integer overflow". Um outro programador me disse que eu poderia resolver isso estocando meus valores como doubles e não como integers, e que eu teria que fazer alguma modificação nesse código para esse fim. Você concorda e poderia me ajudar com isso? obrigado, Vitor Aguiar

2 5

Ajuda - Arima Sazonal
by Vinicius Macedo Magalhães 30 Mar '12

30 Mar '12

Fala, galera. Primeira vez usando este fórum. Tenho tido problemas para rodar meu código para previsão sazonal via ARIMA. Abaixo segue parte do código. >s.ts Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2004 118 143 169 158 143 135 135 140 135 125 120 120 2005 143 158 180 180 150 150 153 148 150 145 145 143 2006 165 180 190 168 163 163 168 175 180 160 155 150 2007 152 163 185 163 158 175 175 178 168 168 160 160 2008 155 170 195 195 195 185 185 165 165 153 154 159 ##### Início de ARIMA ##### linhaX <- 3 dados <- matrix() nome <- names(db) i=4 get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) { best.aic <- 1e8 n <- length(x.ts) for (p in 0:maxord[1]) for (d in 0:maxord[2]) for (q in 0:maxord[3]) for (P in 0:maxord[4]) for (D in 0:maxord[5]) for (Q in 0:maxord[6]) { fit <- arima(x.ts, order=c(p,d,q),seasonal=list(order=c(P,D,Q),period=12), method = "CSS") fit.aic <- -2 * fit$loglik + (log(n) + 1) * length(fit$coef) if (fit.aic < best.aic) { best.aic <- fit.aic best.fit <- fit best.model <- c(p,d,q,P,D,Q) } } list(best.aic,best.fit,best.model) } s.ts <- ts(db[,i],start=2004, fre=12) best.arima.ss <- get.best.arima(log(s.ts),maxord = c(2,2,2,2,2,2)) best.fit.ss <- best.arima.ss[[2]] best.arima.ss[[3]] pdf(paste("AR.", i-4, ".pdf", sep=""),width=7,height=5) ts.plot( cbind( window(s.ts, start=2004),exp(predict(best.fit.ss,24)$pred)),col=c("black","green")) #par(new=TRUE) #ts.plot(s.ts,best.fit.ss,col=c("black","red")) legend("bottomright",c("Real","Previsão_Passado","Previsão_Futuro"),col=1:3,lty=1,cex=0.7,merge = TRUE,bg = 'gray90') dev.off() report<-function(Actual,Forecast,Net,R2) { result <- new.env() if (R2==0) R2<-1-(sum((Forecast-Actual)^2)/sum((Actual-mean(Actual))^2)) result$Rsquare<-R2 result$MAPE <- 100*(sum(abs((Actual-Forecast))/Actual))/nrow(Forecast) ## Mean Absolute Percent Error result$CFE <- sum((Actual-Forecast)) ## Cumulative Forecast Error result$ME <- result$CFE/nrow(Forecast) ## Mean Error result$MSE <- sum((Actual-Forecast)^2)/nrow(Forecast) ## Mean Squared Error result$RMSE <- sqrt(result$MSE) ## Root Mean Squared Error result$MAD <- (sum(abs(Actual-Forecast))/nrow(Forecast)) ## Mean Absolute Deviation result$SUMMARY <- summary(Net) return(result) } ## Fim da função report pred.past<- db[,colY] previsao <- predict(best.fit.ss) forecast <- as.matrix(pred.past[1:length(pred.past)]) reportAR <- report(db[,colY],forecast,best.fit.ss,0) dados <- rbind(dados,"Produto:"=nome[colY],R2=reportAR$Rsquare,MAPE=reportAR$MAPE, CFE=reportAR$CFE,ME=reportAR$ME,MSE=reportAR$MSE,RMSE=reportAR$RMSE,MAD=reportAR$MAD) write.table(dados,sep=",",file="ReportAR.csv",) } ## Fim do for return(dados) ##### Fim de ARIMA ##### -- Att., Vinicius Macedo Magalhães (21) 9584-1533

2 1

loops no R
by Humberto Hazin 30 Mar '12

30 Mar '12

Prezados, Venho mais uma vez pedir ajuda a voces para resolver um problema! baixei uma série de images de satelite, por ano e mes (2004 a 2010) totalizando 163 arquivos ate ai tudo bem! Por termos de comodidade gostaria de criar um loop para nao ter que digitar essas 163 linhas ou mais para extrair latitude, longitude e SST. conforme o codigo abaixo: hdf4to5("S20040012004031.hdf");m0401 <- satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20, -50), itype="oceancolor") m0401<-(simplifySat(m0401, extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1; m0401$year<-2004 hdf4to5("S20040322004060.hdf");m0402 <- satinGet("S20040322004060.hdf5",lats=c(10, -30), lons=c(-20, -50), itype="oceancolor") m0402<-(simplifySat(m0402, extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2; m0402$year<-2004 Alguem poderia me ajudar! Nao coloquei o arquivo pois é muito grande quase 36gigas porém caso precisem posso colocar o parte do codigo para baixarem um ou mais arquivos e testarem Humberto

3 9

Tobit espacial
by Diego andre 30 Mar '12

30 Mar '12

Olá pessoal, alguém sabe como eu posso estimar um tobit espacial com o R? Obrigado. -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

2 1

BRugs_windows 7
by Daniele Portela de oliveira 29 Mar '12

29 Mar '12

Caros colegas eu troquei de computador e o novo veio com windows 7 64 bits. Depois disso não consegui mas usar o pacote BRugs, sempre dá erro quando peço a library. Alguém já teve problema semelhante e poderia me ajudar a solucioná-lo? Sem mai, agradeço > >Daniele Portela de Oliveira > >Doutoranda em Zootecnia >UNESP - Jaboticabal-SP >Cel: 16 8221 6815 > >

2 1

Problemas com fixed breaks
by Julianna Trindade 29 Mar '12

29 Mar '12

Boa Tarde, Gostaria de uma solução para o FixedBreaks obedecer os limites superiores e inferiores das minhas classes de intervalos. No entanto, ele deve associar algum método estatístico e burla minha configuração. Eis o código que eu utilizo: brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", intervalClosure="right", fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000), main="Fixed") brks<- brks$brks Eis o breaks gerado pelo R: c(-11.1, -0.5, 16.665, 41.665, 58.335, 83.335, 110.1) Alguém pode me ajudar, por favor? -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

2 1

Problemas para modificar vetor
by Augusto Ribas 29 Mar '12

29 Mar '12

Ola Pessoal, venho aqui pedir uma ajudinha. Eu tenho um conjunto de dados que são varios nomes de especies. Ai eu resumo eles em siglas para fazer um grafico. Com uma dica anterior do Benilton e do Paulo eu aprendi a resumir os nomes em 3 letras, a primeiro do genero e as duas primeira da especie. Meu problema é o seguinte eu tenho especies como: [1] "Aplectana crossodactyli" "Aplectana crucifer" Que as duas viram "Acr" O que eu gostaria de fazer é qd isso acontece, que são poucos casos dar o nome de Acr1 e Acr2, somente nesses casos. Seguem embaixo como estão os dados e como eu estava fazendo as coisas #dados nomes.lin<-c("Anuracanthorhynchus tritaxisentis", "Aplectana crossodactyli", "Aplectana crucifer", "Aplectana delirae", "Aplectana lopesi", "Aplectana membranosa", "Aplectana meridionalis", "Aplectana micropenis", "Aplectana pintoi", "Aplectana vellardi", "Capillaria recondita", "Catadiscus cohni", "Catadiscus cordeiroi", "Catadiscus inopinatus", "Catadiscus marinholutzi", "Catadiscus mirandai", "Catadiscus propinquus", "Catadiscus pygmaeus", "Catadiscus uruguayensis", "Centrorhynchus giganteus", "Centrorhynchus tumidulus", "cestoda não-identificado", "Choledocystus elegans", "Cosmocerca brasiliense", "Cosmocerca cruzi", "Cosmocerca freitasi", "Cosmocerca ornata", "Cosmocerca parva", "Cosmocerca podicipinus", "Cosmocerca rara", "Cosmocerca travassosi", "Cruzia tentaculata", "Cylindrotaenia americana", "Diplostomum compactum", "Falcaustra mascula", "Glypthelmins linguatula", "Glypthelmins palmipedis", "Glypthelmins parva", "Glypthelmins proximus", "Glypthelmins pseudium", "Glypthelmins repandum", "Glypthelmins simulans", "Glypthelmins vetelliniphilum", "Gorgoderina carioca", "Gorgoderina cedroi", "Gorgoderina cryptorchis", "Gorgoderina diaster", "Gorgoderina festoni", "Gorgoderina pigulevskyi", "Gorgodernina parvicava", "Gyrinicola chabaudi", "Haematoloechus freitasi", "Haematoloechus fuellborni", "Haematoloechus legrandi", "Haematoloechus lutzi", "Haematoloechus ozoroi", "Halipegus dubius", "Infindum infindum", "Maicuru solitarium", "Mesocoelium monas", "Mesocoelium travassosi", "Monostoma sulcatum", "Neohaematoloechus iturbei", "Neohaematoloechus neivai", "Ochoterenella convoluta", "Ochoterenella digiticauda", "Opistogliphe amplicavus", "Oswaldocruzia belenensis", "Oswaldocruzia lopesi", "Oswaldocruzia mazzai", "Oswaldocruzia proencai", "Oswaldocruzia subauricularis", "Oswaldocruzia vaucheri", "Oxyascaris oxyascaris", "Oxysomatium baylisi", "Paraoxyascaris travassosi", "Parapharyngodon alvarengai", "Physaloptera retusa", "Physalopteroides venancioi", "Plagiorchis lenti", "Plagiorchis rangeli", "Polystoma lopezromani", "Pseudoacanthocephalus lutzi", "Pteroxyascaris caudacutus", "Pteroxyascaris similis", "Raillietnema minor", "Raillietnema simples", "Raillietnema spectans", "Rhabdias androgyna", "Rhabdias elegans", "Rhabdias fuelleborni", "Rhabdias hermafrodita", "Rhabdias sphaerocephala", "Rudolphitrema rudolphi", "Schrankiana formosula", "Schrankiana freitasi", "Schrankiana fuscus", "Schrankiana inconspicata", "Schrankiana larvata", "Schrankiana schranki", "Schrankianella brasili", "Schulzia subventricosa", "Schulzia travassosi", "Strongyloides carinii", "Strongyloides pereirai", "Subulascaris falcaustriformis", "Thelandros owaldocruzia", "Travtrema stenocotyle") #fazendo siglas nomes.lin<-rownames(as.data.frame(matriz.final2)) for(i in 1:length(rownames(matriz.final2))) { nomes.lin[i]<-paste(substring(strsplit(nomes.lin[i]," ", fixed="TRUE")[[1]][1], 1,1), substring(strsplit(nomes.lin[i]," ", fixed="TRUE")[[1]][2], 1,2), sep="") } #ai o problema é assim: #no caso do Acr length(nomes.lin[which(nomes.lin==unique(nomes.lin)[2])])>1 #ele ocorre mais de uma vez, então qd isso for verdadeiro #tinha que virar Acr1 e Acr2 #então esses objetos which(nomes.lin==unique(nomes.lin)[2]) #que são iguais tinha que virar o Acr1 e Acr2 Eu tava tentando fazer algo tipo um loop pra testar cada situação unica, igual acima, se ela fosse verdadeira, ai eu tentei fazer um outro loop dentro pra ver qts casos era e ir adicionando 1 2 3 do lado da sigla Mas entrou em loop dentro de loop e me embananei todo :( Alguém poderia dar uma mão? -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Versão do R 2.14.2 mac os travando
by Edson Lira 29 Mar '12

29 Mar '12

Estou rodando um banco do excell xlsx com 1,6 mb. Consigo ler o arquivo, porém, quando uso a rotina names(), o computador está travando. Achei estranho, visto que minha máquina é um mac 2.7 ghz de clok com 4 gb de memória sendo um i7. Alguém tem alguma sugestão? Estou pensando em baixar novamente o R e instalar novamente. Ah! não desisntalei a versão anterior, pode ser isso? []'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Comparar modelos com gnls
by Jonas Brito 29 Mar '12

29 Mar '12

Quero comparar modelos para duas espécies de peixes. Consigo fazer usando nls mas não gnls. Alguém tem alguma idéia de como fazer? dados<-structure(list(Peixes = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"), Tempo = c(0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L, 270L, 270L, 0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L, 270L, 270L), Peso = c(6.09, 6.09, 6.09, 65.73, 52, 53.13, 102.88, 115.8, 129.31, 161.06, 174.94, 213.63, 233.67, 291.93, 301.33, 258.47, 326.84, 337.04, 8.92, 8.92, 8.92, 78.8, 72.07, 85.8, 167.82, 146.13, 179.31, 255.75, 230.73, 207.06, 320.67, 280.2, 243.93, 376.55, 441.96, 349.84)), .Names = c("Peixes", "Tempo", "Peso"), row.names = c(NA, -36L), class = "data.frame") library(nlme) m1 <-gnls(Peso~a*(1-b*(exp(-c*Tempo)))^3, weights=varPower(), data=dados, start=c(a=350,b=1, c=0.01)) m2<- nls(Peso~a*(1-b*(exp(-c*Tempo)))^3,data=dados, start=c(a=300,b=10, c=0.05)) m3<-gnls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,, weights=varPower(), data=dados, start=list(a=c(300,300),b=c(10,10), c=c(0.05,0.05))) m4<-nls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,data=dados, start=list(a=c(300,300),b=c(10,10), c=c(0.05,0.05))) anova(m1,m3) anova(m2,m4)

3 2

ajuda contar subperiodos no R
by Lucas Fernandes de Souza 29 Mar '12

29 Mar '12

Boa tarde. Tenho a seguinte planilha, onde as colunas são anos. Cada linha, represente um determinado dia do ciclo da cultura do milho. Os valores "0" representam que nesse dia, não houve chuva. Os valores "1" representam os dias que choveram. Pois bem, preciso contar os períodos de chuva e seca em cada ano. Porém, não é o perído total, e sim os subperíodos. Exemplo: No ano de 1979, tenho 1 dia de chuva (1° subperiodo), depois tenho 6 dias de seca (2° subperiodo), depois 1 dia de chuva (3 subperiodo), depois 2 dias de seca (4 subperiodo) e assim vai..... Existe um jeito de fazer isso no R automaticamente? Tenho vários anos para fazer isso..se for ficar contando, vou demorar muito... Desde já agradeço a ajuda. obrigado Exemplo: 1979 1980 1 0 0 0 0 1 0 1 0 1 0 0 0 0 1 1 0 0 0 0 1 1 0 1 1 1 0 1 1 0 1 0 1 0 1 0 1 1 0 1 0 1 1 1 1 0 1 0 1 0 1 0 1 0 -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

4 5