- R-br - DINF UFPR

Contrastes e função
by Leandro Marino 18 Ago '11

18 Ago '11

#Caros, existe alguma forma de melhorar este exemplo: vetor <- rep('0',108) z=1 for(i in 0:5){ for(j in 0:2){ for(k in 1:2){ for (l in 0:2){ vetor[z] <- paste(i,j,k,l,sep='') z = z + 1 } } } } vetor #Aproveitando, se quero usar um tapply extraindo média, #desvio-padrão e uma função personalizada, alguém saberia com fazer as #três funções de uma vez só? #Gostaria de executar o comando abaixo de uma única vez em um tapply. #Alguma sugestão?! set.seed(1) x <- rnorm(1000,500,100) y <- trunc(runif(1000,1,10)) fper <- function(x){mean(x-500)/100} dataframe <- t(rbind(tapply(x,y,mean),tapply(x,y,sd),tapply(x,y,fper))) colnames(dataframe) <- c('media','desvp','fper') dataframe Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 R (21) 2103-9600 R.:236 ( leandro(a)cesgranrio.org.br  Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

3 3

Material Novo no CRAN
by Marcos Silva 17 Ago '11

17 Ago '11

Só para avisar que existe um novo material sobre o R em portugues no CRAN. http://cran.r-project.org/doc/contrib/Provete-Estatistica_aplicada.pdf Valeu. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

2 1

Malha sobre poligono
by Eder David Borges da Silva 17 Ago '11

17 Ago '11

#Caros colegas, #minha duvida é uma implementação para gerar uma malha da seguinte forma: # 1) Imaginem que tenho um poligono como esse? x <- c(-0.512320711, -0.503592775, -0.669423554, -0.841072957, -0.459953097, -0.186477777, -0.221389520, -0.000281815, 0.319742495, 0.569943319, 0.503029145, 0.209188642, 0.022992680, -0.169021906, -0.241754703, -0.305759565, -0.352308556, -0.390129611, -0.422132042, -0.422132042,-0.512320711) y <- c( 0.73217213, 0.38635022, 0.08640264, -0.43585901, -0.80638249, -0.44997419, -0.03357638, 0.12169059,-0.13238264, 0.07228747, 0.48162768, 0.81333441 , 1.10622439, 1.13092596, 1.14504114, 1.13445475, 1.12739716, 1.12033957, 1.06035006, 1.00388934, 0.73217213) plot(x,y,type='l') # 2) dentro deste poligono quero fazer quadrados irregulares como este px1 <- c(-0.6131326, -0.4633484, -0.3762646, -0.5034070) py1 <- c(-0.6037305,-0.7381249, -0.5774359, -0.3992173) polygon(px1,py1) # 3) quero fazer isso com o mouser clicando 4 vezes e guardando as coordenadas e cada coordenada ganha um ID, ficando assim: data.frame(id=1:length(px1),px1,py1) # 4) este quatros pontos forman um elemento (quadrado irregular), onde o elemento guarda o ID os pontos clicados. elemento <- c(1,1,2,3,4)#Idelemento - 4 ids dos pontos # 5) Esse procedimento é repetido n vezes cobrindo toda a região com poligonos irregulares, como resultado final: # um data frame com id e x, y do ponto # um data frame com id do elemento # restrições para o segundo elemento eu devo clicar por exemplo no ponto 2 e 3 e fazer mais o 5 e 6 para fechar ele, desta forma, # os elementos tem de ficar plotados para aproveitar os pontos já feitos. #OBS: para os deterministicos de plantão isso é uma malha para modelos numericos de mecanica de fluidos. #OBS1: Sugestoes de softwares livres que faz isso são bem vindo. # Diante do problema alguem tem sugestoes de funçoes que posso utilizar, na verdade é um mix de locator + identifique, pois não # consigui fechar todo o raciocinio de como fazer isso. # Obrigado pessoal #Att

2 2

Modelo multionomial
by Rodrigo Pereira 17 Ago '11

17 Ago '11

Prezados colegas, bom dia. Estou tentando ajustar um modelo logístico multinomial no R, mas creio que não estou tendo sucesso. A variável resposta tem três níveis (0, 1 e 2) e com 3 variáveis resposta, sendo duas na forma dicotômica e uma contínua. Tentei ajustar um modelo utilizando a função 'glm', mas para a opção de família não há explícita a multinomial. Devo utilizar a binomial e automaticamente a função entenderá que deverá ocorrer uma modificação? Há algum pacote que forneça uma análise completa para uma regressão multinomial? Além disso, preciso de algumas outras estatísticas, tais como odds ratio ou risco relativo (se possível numa comparação duas a duas), intervalos de confiança para os parâmetros, avaliação de efeitos, etc. Desde já agradeço a atenção dispensada. Atenciosamente, Rodrigo Pereira

6 9

Ajuda para simulações
by sznelwar 17 Ago '11

17 Ago '11

4 5

[OFF-Topic] Mendelev
by FHRB Toledo 17 Ago '11

17 Ago '11

Pessoal da R-Br, olá Não sei se todos estão altamente "antenados", mas parece que estamos vivendo uma revolução... Surgem grandes oportundiades p/ o pessoal envolvido em pesquisa ao usar o Mendelev! Pois bem, estou cursando uma disciplina e o professor faz parte de algum "circulo maior" da organização do mendelev! Para incentivar o seu uso ele agendou um webnário sobre o ambiente... Seguem as informações: Atendendo ao foi solicitado, o pessoal do Mendeley vai apresentar um seminário em português especial para os alunos da disciplina e estou tomando a liberdade de extender p/ turma da lista! Data e horário: Terça-feira, dia 16 de Agosto, 11:00 da manhã (hora de Brasília). Local: internet (é um "webnário"). Vocês devem se cadastrar antes no link https://mendeley.webex.com/mendeley/onstage/g.php?t=a&d=663141002 e na hora do seminário basta entrar na mesma página. **IMPORTANTE:** uma vez que isso será feito porque nós pedimos, será bem chato se poucos participarem, então por favor não deixem de se inscrever e entrar na página no horário correto. Estou esperando pelo menos uns 30 participantes. Abraços, FH

9 14

Ajuda com modelo não linear
by Eduardo Lucas Terra Peixoto 17 Ago '11

17 Ago '11

Boa noite a todosEstou fazendo doutorado na UEL no curso de Ciência animal. Ainda sou novato no uso do R e também não tenho base estatística aprofundada em modelos não lineares.E preciso de um help, estou adaptando uma fórmula de modelo não linear do proc do SAS para o R, mas li o help, busquei dicas no site http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase30.html#x32-20100030 e fui seguindo o modelo explicado por esse autor, mas tem algumas coisas que não entendi e não estou conseguindo rodar o modelo.Primeiramente explicarei resumidamente o que se trata o modelo. Esse modelo é para poder determinar cinco outras variáveis, através da produção de gás pela fermentação de bacteriais do líquido ruminal in vitro, simulando assim o rumen do animal. As variáveis que tenho em mãos são o volume de gás produzido e o tempo da leitura, a partir deles preciso determinar:VcnF = volume de gás produzido pelo carboidratos não fibrosos (mL de gás)kdcnf = taxa de degradação dos carboidratos não fibrosos (%/hora)VcF = volume de gás produzido pelo carboidratos fibrosos (mL de gás)kdcf = taxa de degradação dos carboidratos fibrosos (%/hora)L= tempo de latência (h), este é o tempo que as bactérias gastam para começar a degradar o substrato.Mas ao rodar o modelo abaixo aparece o seguinte erro: fit30=nls(Y~ (VcnF/(1+exp(2-4*kdcnf*I(L-Tempo)))+ VcF/I(1+exp(2-4*kdcf*(L-Tempo)))), + data=cra30, + start=list(VcnF = 10,kdcnf = 0.01,L =0.1,VcF = 50,kdcf =0.01, Tempo=1), + alpha= (0.05))Error in qr.qty(QR, resid) : NA/NaN/Inf in foreign function call (arg 5)In addition: Warning message:In Ops.factor(lhs, rhs) : - not meaningful for factors Aproveitando o ensejo gostaria de perguntar como faço para limitar os parâmetros pois no proc do SAS está assim:VcnF = 10 to 200 by 50 kdcnf = 0.01 to 0.8 by 0.05 L = 0.1 to 5 by 0.2 VcF = 50 to 150 by 30 kdcf = 0.01 to 0.8 by 0.05; bounds L >=0; Caso alguém tenha materiais relativos a modelos não linear para que possa aprofundar melhor, também contribuirá bastante. Agradeço a todos, e espero que alguém possa me ajudar! Eduardo Lucas Terra PeixotoDiscente de doutorado - Ciência AnimalUniversidade Estadual de Londrina Cel. (43) 9653-5574

2 1

Fwd: [Computing] package parallel in R-2.14.0
by Leonard Assis 16 Ago '11

16 Ago '11

Boa notícia lmassis <at> yahoo <dot> com <dot> br assis.leonard <at> gmail <dot> com ---------- Forwarded message ---------- From: Jan de Leeuw <deleeuw(a)stat.ucla.edu> Date: Tue, Aug 16, 2011 at 1:40 PM Subject: [Computing] package parallel in R-2.14.0 To: statcompute <computing(a)stat.ucla.edu> parallel is a recommended package in R-2.14.0 -- I attach the vignette. === Jan de Leeuw; Distinguished Professor and Chair, UCLA Department of Statistics; Editor: Journal of Multivariate Analysis, Journal of Statistical Software; US mail: 8125 Math Sciences Bldg, Box 951554, Los Angeles, CA 90095-1554 phone (310)-825-9550; fax (310)-206-5658; email: deleeuw(a)stat.ucla.edu .mac: jdeleeuw ++++++ aim: deleeuwjan ++++++ skype: j_deleeuw homepages: http://gifi.stat.ucla.edu ++++++ http://www.cuddyvalley.org ------------------------------------------------------------------------------------------------- No matter where you go, there you are. --- Buckaroo Banzai http://gifi.stat.ucla.edu/sounds/nomatter.au ------------------------------------------------------------------------------------------------- _______________________________________________ Computing mailing list Computing(a)lists.stat.ucla.edu http://lists.stat.ucla.edu/mailman/listinfo/computing

1 0

dúvida em glm
by Vinicius Brito Rocha 16 Ago '11

16 Ago '11

Senhores, ajustei um glm com família binommial e link logit, onde minha variável resposta é a probabilidade do evento (chamadas telefônicas repetidas em 24hs / total de chamadas telefônicas), obtidos a aprtir de um aggregate em um banco de dados com informação de telefones. A minha varíavel independente é um fator (UF), Não estou conseguindo interpretar os coeficientes do meu modelo. tenho 27 coef, devido ao meu fator ser UF +1 intercepto. Alguém pode me dar uma força? 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.

3 6

Re: [R-br] Análise Envoltória de Dados (DEA)
by Roney Fraga Souza 16 Ago '11

16 Ago '11

Daniel, bom dia! Ainda não trabalhei com DEA SBM sem orientação, contudo estou aqui pensando se vc utilizar no comando dea() as opções ORIENTATION='in-out' com DIRECT diferente de NULL será que chega perto do seu objetivo? Outra dúvida, em um tópico anterior também sobre DEA, no link, http://r-br.2285057.n4.nabble.com/R-br-Modelo-linear-generalizado-misto-td3… vc postou que, Pacote Benchmarking: Vantagem - Te retorna os pesos(contribuição) de cada DMU onde a soma é igual a 1. Desvantagem - Não te fornece os slacks (quanto a DMU deve aumentar ou diminuir de produto ou insumo para se tornar eficiente) com relação a não retornar os slacks, quando utilizo em ccr.in <- dea(insumos, produtos, RTS='crs', ORIENTATION='in', SLACK=TRUE) o comando summary(ccr.in) já aparece se existe folga em x ou em y. E quando chamo names(ccr.in) cbind(ccr.in$sx,ccr.in$sy,ccr.in$slack) tenho as folgas de x e y. Minha pergunta é, minha análise está equivocada, ou realmente é possível obter as folgas pelo pacote Benchmarking? Atenciosamente Roney

1 0

[OFF-TOPIC] (webnário de) o que github?
by Walmes Zeviani 16 Ago '11

16 Ago '11

Pessoal, Aproveitando o embalo sobre soluções (web livres) para facilitar a nossa vida em colaboração com outras pessoas, gostaria de anunciar que vai haver um webnário de GitHub amanhã também ( https://github.com/blog/906-github-training-tomorrow) Gostaria de pedir para quem conhece o GitHub (Benilton Carvalho, Fernando Mayer, ...) para dar uma pequena noção para a gente da ferramenta. Inclusive, pelo que andei lendo, é possível hospedar dados lá e importá-los por link, ótimo para CMRs. Parece ser útil para quem tá desenvolvendo dissertação/tese por ser gerenciador de versões, tô certo? Grato. 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 ==========================================================================

4 3

Análise filogenética (pPCA)
by Marcelo Claro de Souza 16 Ago '11

16 Ago '11

Preciso rodar uma análise de componentes principais filogenética (ppca) utilizando o pacote adephylo. Gostaria de saber como montar uma árvore filogenica de Newik e salvar o resultado em formato phylo, para depois poder rodar a ppca. Possuo um conjunto de dados formado por 36 espécies e seis variáveis. Outra dúvida, o nome das espécies deve ser o nome das linhas ou devo manter uma coluna com o nome das espécies, ficando com 7 colunas ao invés de 6. Muito obrigado.

1 0

Re: [R-br] Análise Envoltória de Dados (DEA)
by Roney Fraga Souza 16 Ago '11

16 Ago '11

Daniel, boa noite! Não compreendi o seu problema postado na lista, possivelmente pq não utilizo o pacote DEA, onde posso baixar o pacote DEA para testes? Também utilizo o R para analisar eficiência por meio de análise envoltória de dados. Estou trabalhando com o pacote Benchmarking, que além do excelente help tem o livro do mesmo autores que criaram o pacote que é bastante útil: Bogetoft and Otto (2011), Benchmarking with DEA, SFA, and R. Como exemplo do uso do pacote Benchmarking segue medidas de eficiência para o modelos básicos de DEA utilizando os dados da Tabela 2.1 de Coll e Blasco (2006) Evaluación de la Eficiencia mediante el Análisis Envolvente de Datos. Vale ressaltar que no seu caso o arquivo é grande, e nesse exemplo o arquivo .csv tem poucas linhas, mas acredito que não terá problema em manipular suas matrizes ou data.frames utilizando o Benchmarking. Espero que o exemplo sejá útil, caso encontre algum erro pode corrigir e postar a versão corrigida. Att Roney #################### # TODO: Dados da Tabla 2.1 Valores observados concesionarios, do livro de Coll e Blasco (2006) # INPUTS # x1 = Número de empleados # x2 = Depreciación del Inmovilizado, como proxy del capital # OUTPUTS # y1 = Número de vehículos vendidos # y2 = Número de órdenes de trabajo recibidas en taller # # Coll e Blasco (2006) Evaluación de la Eficiencia mediante el Análisis Envolvente de Datos # Ferreira e Gomes (2009) Introdução à Análise Envoltória de Dados. Teoria, Modelos e Aplicações # Bogetoft, Otto (2011) Benchmark and frontier analysis using DEA and SFA and R # Author: roney ############################################################################### rm(list=ls(all=TRUE)) getwd() setwd("/Users/roney/Documents/Economia/R-workspace/DEA") dir() library(Benchmarking) help(package="Benchmarking") ?dea ## dea(X, Y, RTS="vrs", ORIENTATION="in", XREF=NULL, YREF=NULL, # FRONT.IDX=NULL, SLACK=FALSE, DUAL=FALSE, DIRECT=NULL, param=NULL, # TRANSPOSE=FALSE, FAST=FALSE, LP=FALSE, CONTROL=NULL, LPK=NULL) # # Tabla 2.1 Valores observados concesionarios # x1 x2 y1 y2 # A 8 8 14 20 # B 11 15 25 42 # C 14 12 8 30 # D 12 13 25 8 # E 11 18 40 22 # F 18 20 24 30 # link baixar o arquivo .csv # http://www.datafilehost.com/download-b7ad5202.html read.csv("Coll-Blasco_exemplo.csv", header=TRUE) coll.blasco <- read.csv("Coll-Blasco_exemplo.csv", header=TRUE) namesAF <- c("A", "B", "C", "D", "E", "F") namesX <-c("x1", "x2") namesY <- c("y1", "y2") namesXY <- c("x1", "x2", "y1", "y2") insumos <- matrix(c(coll.blasco$x1, coll.blasco$x2), nrow=6, ncol=2, byrow=FALSE, dimnames=list(namesAF, namesX)) is.matrix(insumos) insumos produtos <- matrix(c(coll.blasco$y1, coll.blasco$y2), nrow=6, ncol=2, byrow=FALSE, dimnames=list(namesAF, namesY)) is.matrix(produtos) produtos ######## CCR Input Orientado ####### # Em benchmarking modelo CCR (crs) input orientado, o problema primal é a minimização do insumos dada a # quantidade de produtos. Modelo envoltório (que usa lambda), ver notação página 72 Ferreira e Gomes (2009) ccr.in <- dea(insumos, produtos, RTS="crs", ORIENTATION="in", SLACK=TRUE) # eff CCR insumo orientado com folgas summary(ccr.in) # função que trás o resumo das medidas, como número e % de unidades eficientes em cada # faixa de valor, assim como folgas. Muito útil para grande quantidade de DMU's lambda(ccr.in) # equivalente a saída Benchmarking do SAID, indica quais DMU's estão sendo referência, # as que estão nas colunas, para as demais DMU's, que estão nas linhas. # em outras palavras, parceiros relevantes (peers) para as DMU's ineficientes, e o valor do # lambda indica o quanto a DMU eficiente, o benchmarking, é importante para a DMU ineficiente lambda(ccr.in, KEEPREF = TRUE) # quando a opção "KEEPREF = TRUE" é utilizada todas as DMU's são mostradas # nas colunas, não penas as eficiêntes. print(ccr.in$eff, digits=2) # para mostrar a eficiência das DMU's which( ccr.in$eff == 1 & !ccr.in$slack) # para mostrar apenas as DMU's eficientes e sem folga data.frame("eff"=ccr.in$eff, "ins"=insumos, "rad"=ccr.in$eff * insumos, "fol"=ccr.in$sx, "alv"=insumos - ccr.in$sx) # i) eff, ii) insumos observados ou atual, iii) movimento radial, iv) folga e v) alvo # # iii) movimento radial é o cálculo da redução dos insumos em direção a fronteira eficiente # e é obtido pela multiplicação dos insumos observados pela eficiência das respectivas # unidades, ou pela multiplicação da unidade do insumo da unidade ineficiente pelo # lambda do(s) seu(s) benchmarks. (ver página 101 de Ferreira e Gomes, 2009) # iv) mesmo projetando DMU em direção a fronteira eficiente devido a possibilidade # de existir alguns seguimentos da fronteira poliangular linear paralelos aos eixos # coordenadas é possível ocorrer folgas nesses pontos. Ou seja, mesmo que o movimento # radial tenha projetado a DMU para a fronteira eficiente é possivel existir alguma # ineficiencia, que caracterizamos como folga. (ver página 102 de Ferreira e Gomes, 2009) # v) o alvo é o movimento radial dimunuido das possíveis folgas data.frame("eff"=ccr.in$eff,"pro"=produtos,"rad"=ccr.in$eff*produtos,"fol"=ccr.in$sy,"alv"=produtos+ccr.in$sy) # i) eff, ii) produtos observados ou atual, iii) movimento radial, iv) folga e v) alvo # o mesmo do exemplo anterior mas agora aplicado para aos produtos ## Problema da Dualidade ## # no problema dual do modelo CCR (crs) insumos orientado [maximização] é expressa a forma multiplicada desse modelo, # onde os lambdas são substituidos pelos peso insumo (u) e peso produto (v). ver notação pág 72 Ferreira e Gomes (2009) # # os pesos, peso insumo (u) e peso produto (v), que permitem calcular os insumos e produtos virtuais conforme são # obtidos na saída do SAID e no modelo insumo orientado multiplicadores (primal) de Coll e Blasco (2006) só são # encontrados pelo problema dual no pacote Benchmarking ccr.in.dual <- dea.dual(insumos, produtos, RTS="crs", ORIENTATION="in") names(ccr.in.dual) print(cbind("eff"=ccr.in.dual$eff, ccr.in.dual$u, ccr.in.dual$v), digits=5) # o valor eficiência é exatamente igual aos do SAID e do Excel, mas os pesos apresentam valores diferentes # para detalhes ver tabela 4.8 na página 130 de Ferreira e Gomes (2009) ou nas páginas 117 e 118 tabelas 4.2 e 4.3 ######## CCR Output Orientado ####### # Em benchmarking modelo CCR (crs) output orientado, o problema primal é a maximização do produto dada a # quantidade de insumos. Modelo envoltório (que usa lambda), ver notação página 130 Ferreira e Gomes (2009) ccr.out <- dea(insumos, produtos, RTS="crs", ORIENTATION="out", SLACK=TRUE) summary(ccr.out) lambda(ccr.out) # apresentou uma pequena mudança no valor dos lambdas comparado com o resultado do SAID, mas como # o modelo é CCR acredito não deveria acontecer isso lambda(ccr.out, KEEPREF = TRUE) print(1/ccr.out$eff, digits=2) which(1/ccr.out$eff == 1 & !ccr.out$slack) data.frame("eff"=1/ccr.out$eff, "ins"=insumos, "rad"=(1/ccr.out$eff) * insumos, "fol"=ccr.out$sx, "alv"=insumos - ccr.in$sx) data.frame("eff"=1/ccr.out$eff, "pro"=produtos, "rad"=1/ccr.out$eff * produtos, "fol"=ccr.out$sy, "alv"=produtos + ccr.out$sy) # apresentou pequenas mudanças diante dos resultodos do SAID # # o problema dual do modelo CCR (crs) output orientado é a minimização, modelo dos multiplicadores, que considera # os pesos insumos (u) e produtos (v) ccr.out.dual <- dea.dual(insumos, produtos, RTS="crs", ORIENTATION="out") # problema dual output orientado names(ccr.out.dual) print(cbind("eff"=ccr.out.dual$eff, ccr.out.dual$u, ccr.out.dual$v), digits=5) ######## BCC Input Orientado ####### # Em benchmarking modelo BCC (vrs) input orientado, o problema primal é a minimização do insumos dada a # quantidade de produtos. Modelo envoltório (que usa lambda), ver notação página 130 Ferreira e Gomes (2009) bcc.in <- dea(insumos, produtos, RTS="vrs", ORIENTATION="in", SLACK=TRUE) # eff BCC insumo orientado com folgas summary(bcc.in) lambda(bcc.in) lambda(bcc.in, KEEPREF = TRUE) print(bcc.in$eff, digits=2) which(bcc.in$eff == 1 & !bcc.in$slack) data.frame("eff"=bcc.in$eff, "ins"=insumos, "rad"=bcc.in$eff * insumos, "fol"=bcc.in$sx, "alv"=insumos - bcc.in$sx) data.frame("eff"=bcc.in$eff, "ins"=produtos, "rad"=bcc.in$eff * produtos,"fol"=bcc.in$sy, "alv"=produtos + bcc.in$sy) ## RENDIMENTOS DE ESCALA ## # na saida do SAID, além dos pesos do CCR, tem v0 que corresponde ao k da abordagem de Coll e Blasco (2006) bcc.in.irs <- dea(insumos, produtos, RTS="irs", ORIENTATION="in", SLACK=TRUE) # eff BCC ins orientado rend crescentes bcc.in.drs <- dea(insumos, produtos, RTS="drs", ORIENTATION="in", SLACK=TRUE) # eff BCC ins orientado rend decrescentes # Ferreira e Gomes (2009) página 198 # CCR = BCC rendimentos constantes de escala # DRS = RVE rendimentos decrescentes, se DRS != RVE rendimentos crescentes # IRS = RVE rendimentos crescentes, se IRS != RVE rendimentos decrescentes data.frame("CRS"=ccr.in$eff,"VRS"=bcc.in$eff, "IRS"=bcc.in.irs$eff, "DRS"=bcc.in.drs$eff, "E_ESC"=ccr.in$eff/bcc.in$eff, "REND"=ifelse(ccr.in$eff == bcc.in$eff | bcc.in$eff == bcc.in.irs$eff & bcc.in$eff == bcc.in.drs$eff, "constante", ifelse(bcc.in$eff == bcc.in.drs$eff & bcc.in$eff != bcc.in.irs$eff, "decrescen", "crescente"))) # de modo equivalente pode-se fazer data.frame("CRS"=ccr.in$eff,"VRS"=bcc.in$eff, "IRS"=bcc.in.irs$eff, "DRS"=bcc.in.drs$eff, "E_ESC"=ccr.in$eff/bcc.in$eff, "REND"=ifelse(ccr.in$eff == bcc.in$eff | bcc.in$eff == bcc.in.irs$eff & bcc.in$eff == bcc.in.drs$eff, "constante", ifelse(bcc.in$eff == bcc.in.irs$eff & bcc.in$eff != bcc.in.drs$eff, "crescente", "decrescen"))) # no problema dual do modelo BCC (vrs) insumos orientado [maximização] é expressa na forma multiplicada, onde os # lambdas são substituidos pelos peso insumo (u) e peso produto (v). ver notação pág 117 Ferreira e Gomes (2009) bcc.in.dual <- dea.dual(insumos, produtos, RTS="vrs", ORIENTATION="in") names(bcc.in.dual) print(cbind("eff"=bcc.in.dual$eff, bcc.in.dual$u, bcc.in.dual$v), digits=3) ######## BCC Output Orientado ####### # Em benchmarking modelo BCC (vrs) output orientado, o problema primal é a maximização do produto dada a # quantidade de insumos. Modelo envoltório (que usa lambda), ver notação página 118 Ferreira e Gomes (2009) bcc.out <- dea(insumos, produtos, RTS="vrs", ORIENTATION="out", SLACK=TRUE) summary(bcc.out) lambda(bcc.out) lambda(bcc.out, KEEPREF = TRUE) print(1/bcc.out$eff, digits=2) which(1/bcc.out$eff == 1 & !bcc.out$slack) data.frame("eff"=1/bcc.out$eff, "ins"=insumos, "rad"=(1/bcc.out$eff) * insumos,"fol"=bcc.out$sx, "alv"=insumos - bcc.in$sx) data.frame("eff"=1/bcc.out$eff, "pro"=produtos, "rad"=1/bcc.out$eff * produtos,"fol"=bcc.out$sy, "alv"=produtos + bcc.out$sy) bcc.out.irs <- dea(insumos, produtos, RTS="irs", ORIENTATION="out", SLACK=TRUE) # eff BCC out orientado rend crescentes bcc.out.drs <- dea(insumos, produtos, RTS="drs", ORIENTATION="out", SLACK=TRUE) # eff BCC out orientado rend decrescentes data.frame("CRS"=ccr.out$eff,"VRS"=bcc.out$eff, "IRS"=bcc.out.irs$eff, "DRS"=bcc.out.drs$eff, "E_ESC"=ccr.out$eff/bcc.out$eff, "REND"=ifelse(ccr.out$eff == bcc.out$eff | bcc.out$eff == bcc.out.irs$eff & bcc.out$eff == bcc.out.drs$eff, "constante", ifelse(bcc.out$eff == bcc.out.irs$eff & bcc.out$eff != bcc.out.drs$eff,"crescente","decrescen"))) # de modo equivalente pode-se escrever data.frame("CRS"=ccr.out$eff,"VRS"=bcc.out$eff, "IRS"=bcc.out.irs$eff, "DRS"=bcc.out.drs$eff, "E_ESC"=ccr.out$eff/bcc.out$eff, "REND"=ifelse(ccr.out$eff == bcc.out$eff | bcc.out$eff == bcc.out.irs$eff & bcc.out$eff == bcc.out.drs$eff, "constante", ifelse(bcc.out$eff == bcc.out.drs$eff,"decrescen","crescente"))) ### CUIDADO!! quando exite diferença entre as variáves selecionadas, considerando muitas casa ### decimais, pode ocorrer equívocuo na análise de rendimentos de escala. para evitar esse erro é necessário ### limitar o número de casa decinais com a função round(x, digits=n) print(cbind("CRS"=ccr.out$eff, "VRS"=bcc.out$eff, "IRS"=bcc.out.irs$eff,"DRS"=bcc.out.drs$eff), digits=18) ee <- data.frame(round((cbind("CRS"=ccr.out$eff, "VRS"=bcc.out$eff, "IRS"=bcc.out.irs$eff,"DRS"=bcc.out.drs$eff)), digits=6)) data.frame(ee, "E_ESC"=ee$CRS/ee$VRS,"REND"=ifelse(ee$CRS == ee$VRS | ee$VRS == ee$IRS & ee$VRS == ee$DRS, "constante", ifelse(ee$VRS == ee$IRS & ee$VRS != ee$DRS,"crescente","decrescen"))) # o problema dual BCC (vrs) output orientado é de minimização, na forma multiplicada que considera os pesos dos # insumos (u) e produtos (v) bcc.out.dual <- dea.dual(insumos, produtos, RTS="vrs", ORIENTATION="out") # problema dual names(bcc.out.dual) print(cbind("eff"=bcc.out.dual$eff, bcc.out.dual$u, bcc.out.dual$v), digits=5) ######## Análise de Supereficiência ####### # ver Coll e Blasco (2006) capítulo 4 página 135 # ver Ferreira e Gomes (2009) capítulo 4 página 136 s.ccr.in <- sdea(insumos, produtos, RTS="crs", ORIENTATION="in") data.frame("CRS"=ccr.in$eff, "CCR_SUPER"=s.ccr.in$eff) s.ccr.out <- sdea(insumos, produtos, RTS="crs", ORIENTATION="out") data.frame("CCR"=1/ccr.out$eff, "CCR_SUPER"=1/s.ccr.out$eff) ######## Modelo FHD ####### # ver Ferreira e Gomes (2009) capítulo 4 página 143 bcc.in.fhd <- dea(insumos, produtos, RTS="fdh", ORIENTATION="in") data.frame("CRS"=ccr.in$eff,"VRS"=bcc.in$eff, "IRS"=bcc.in.irs$eff, "DRS"=bcc.in.drs$eff, "FHD"=bcc.in.fhd$eff) bcc.out.fhd <- dea(insumos, produtos, RTS="fdh", ORIENTATION="out") data.frame("CRS"=1/ccr.out$eff,"VRS"=1/bcc.out$eff, "IRS"=1/bcc.out.irs$eff, "DRS"=1/bcc.out.drs$eff, "FHD"=1/bcc.out.fhd$eff) ######## Seleção de variáveis ######## cor(coll.blasco, use = "all.obs", method = c("spearman")) # teste de correlação de Spearman cor(coll.blasco, use = "all.obs", method = c("kendall")) # teste de correlação de Kendall cor(coll.blasco, use = "all.obs", method = c("pearson")) # teste de correlação de Pearson ######## Eficiência custo (econômica), alocativa, receita e lucro ####### # ver Ferreira e Gomes (2009) capítulo 5 página 213 ins <- data.frame(insumos) p.ins <- data.frame("px1"=c(2,2,2,2,2,2), "px2"=c(6,4,3,4,3,2)) eff.cost <- cost.opt(ins, produtos, p.ins, RTS="vrs") print(cbind(ins,p.ins, "custo_min"= eff.cost$cost, "eff_econ"=eff.cost$cost/ (ins$x1 * p.ins$px1 + ins$x2 * p.ins$px2), "eff_tec"=bcc.in$eff, "eff_aloc"=(eff.cost$cost/ (ins$x1 * p.ins$px1 + ins$x2 * p.ins$px2))/bcc.in$eff),digits=3) pro <- data.frame(produtos) p.pro <- matrix(1,nrow=dim(pro)[1],ncol=2) eff.renevue <- revenue.opt(ins, pro, p.pro, RTS="vrs") eff.profit <- profit.opt(ins, pro, p.ins, p.ins, RTS="vrs")

2 1

Re: [R-br] stepAIC
by Daniela Recchia 16 Ago '11

16 Ago '11

Olá Walmes, eu usei no stepAIC(modelo,na.action=na.exclude) isso ja deveria bastar né? Mas mesmo assim o erro no stepAIC continua. -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "It is better to have an inexact answer to the right question than a precise one to the wrong question!" John Tuckey ________________________________ De: Walmes Zeviani <walmeszeviani(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br; Daniela Recchia <daniela_recchia(a)yahoo.com.br> Enviadas: Segunda-feira, 15 de Agosto de 2011 15:07 Assunto: Re: [R-br] stepAIC Daniela, A princípio é só remover os NA, deixando a base completa para todas as descritoras. O fato de ser glm binomial e descritoras categóricas não limita o uso do step(). À 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 ==========================================================================

2 1

Dados sem distribuição
by Adriano Carvalho Costa 15 Ago '11

15 Ago '11

Boa noite, Estou com a base de dados acima e estou fazendo transformação box-cox e os dados não atenderam a distribuição normal, voces tem alguma rotina de análise não paramétrica? Script dados<-read.table("final.txt", h=T) dados names(dados) attach(dados) head(dados) dim(dados) Transformação Box-Cox require(MASS) library(MASS) motil Trat boxcox(motil+0.01~Trat, data=dados) lambida=0.35 motili=(motil^(0.39) - 1)/0.39 motili plot(motil~Trat) plot(motili~Trat) anava <- aov(motili ~ Trat, data = dados) plot(anava) shapiro.test(anava$res)

3 7

[OFF-Topic] Curso on line do R
by Gilson Sanchez 15 Ago '11

15 Ago '11

Só repasando, este foi o mail da semana pasada, mas curso ainda esta em andamento. Walmes Haber se você se anima fazer um curso on line do R parecido. Espero que aproveitem, Abç’s ---------- Forwarded message ---------- From: Paul Geissler <GoToWebinar.Notifications(a)citrixonline.com> Date: Thu, Aug 11, 2011 at 1:19 PM Subject: Reminder : Your Webinar is on Thursday, August 11, 2011 12:00 PM - 1:00 PM MDT To: Gilson Sanchez Chia <gilsonsch(a)gmail.com> Learn the R Statistical Package <https://www1.gotomeeting.com/join/655473105/106560608> Join us on Thursday, August 11, 2011 12:00 PM - 1:00 PM MDT Dear Gilson, This message is to remind you that the following Webinar will take place Thursday, August 11, 2011 12:00 PM - 1:00 PM MDT. *Learn the R Statistical Package* 1. Click here to join: https://www1.gotomeeting.com/join/655473105/106560608 This link *should not be shared* with others; it is unique to you. 2. You will be connected to audio using your computer's microphone and speakers (VoIP). A headset is recommended. Or, you may select Use Telephone after joining the Webinar. Toll: +1 (805) 309-0021 Access Code: 234-132-963 Audio PIN: Shown after joining the Webinar Webinar ID: 655-473-105 Please send your questions, comments and feedback to: Paul_Geissler(a)usgs.gov. *System Requirements* PC-based attendees Required: Windows® 7, Vista, XP or 2003 Server Macintosh®-based attendees Required: Mac OS® X 10.5 or newer Read our Audio Checklist<https://www1.gotomeeting.com/default/help/g2w/g2w%5fuser%5fhelp.htm#Webinar…>for tips on using your computer's microphone and speakers with GoToWebinar. Add to your Outlook calendar<https://www1.gotomeeting.com/synchOutlook?wid=655473105&uid=106560608> You are receiving this email because you registered for this Webinar. Please cancel your registration<https://www1.gotomeeting.com/island/webinar/cancelRegistration.tmpl?Webinar…>if you cannot attend. You can also opt-out <https://www1.gotomeeting.com/unsubscribe/655473105/106560608> from receiving further emails from this Webinar's organizers. This event is powered by GoToWebinar, Webinars Made Easy™. View the GoToWebinar Privacy Policy <https://www1.gotomeeting.com/island/webinar/privacy.tmpl>. GoToWebinar is a service provided by : Citrix Online | 7414 Hollister Avenue | Goleta, CA 93117 -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

6 6

Re: [R-br] calculo de volume de função numérica (Paulo J Ribeiro Jr)
by Renato Assunção 15 Ago '11

15 Ago '11

Outra opcao e' usar tecnicas Monte Carlo. Joguem N pontos iid num quadrado (basta fazer uniforme em cada dimensao do quadrado. Selecione apenas aqueles pontos que cairem no seu dominio (o pacote splancs tem uma funcao para isto, verificar se um ponto esta num poligono) Tire a media de f(x_i) que e' uma boa aproximacao de integral de f(x) dividida pela AREA do poligono. Multiplique entao pela area do poligono (tambem tem uma funcao no splancs para isto). Renato

1 0

stepAIC
by Daniela Recchia 15 Ago '11

15 Ago '11

Boa tarde pessoal, estou tentando usar a funcao step ou stepAIC mas recebo um erro quando ela termina de rodar os modelos: o erro diz que o numero de linhas utilizadas foi alterado: Numeros faltantes devem ser apagados? Tenho dados categoricos, que vieram de um questionario que aplicamos com os alunos aqui da universidade. Queremos saber pq grande parte dos alunos interrompe ou desiste do curso de estatistica. A variavel resposta é Se a pessoa esta pensando em desistir do curso, se ja pensou ou se nunca pensou, e as descritoras sao, quantidade de disciplinas aprovadas, situacao financeira etc.. todas categoricas. Pensei entao em ajustar um modelo de regressao logistica com lm ou glm, family=binomial, nao sei o que seria mais correto. Mas para a selecao das variaveis, gostaria de usar o step, e nao estou conseguindo, por causa do erro. ALguem teria uma dice? Abracos, -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "It is better to have an inexact answer to the right question than a precise one to the wrong question!" John Tuckey

2 1

Re: [R-br] .Rprofile no Linux
by Pedro Rafael 15 Ago '11

15 Ago '11

Isso funciona sem alteração no Windows? Em 13 de agosto de 2011 20:34, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Isso funciona sem alteração no Windows? > > > Saudações, > Pedro Rafael > Em 13 de agosto de 2011 08:20, jcfaria [via R-br] < > ml-node+3741147-2070390175-223914(a)n4.nabble.com> escreveu: > > Pessoal, >> >> Encerrando esse tópico (que ficou bem grande), informo que, seguindo >> as sugestões do Paulo J., as intruções originais foram melhoradas e >> adicionadas nas R-idículas do >> LEG:http://www.leg.ufpr.br/doku.php/ridiculas?&#rprofile_no_linux >> >> >> Abs, >> -- >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> Jose Claudio Faria >> Estatistica - Prof. Pleno >> UESC/DCET/Brasil >> joseclaudio.faria at gmail.com >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> >> >> Em 13 de agosto de 2011 07:16, Jose Claudio Faria >> <[hidden email] <http://user/SendEmail.jtp?type=node&node=3741147&i=0>> >> escreveu: >> >> > Pedro, >> > >> > Não! >> > Seus objetos serão criados automaticamente a cada nova sessão, desde >> > que as instruções estejam localizadas no arquivo .Rprofile. >> > Dê uma olhada no post original. >> > >> > Abs, >> > -- >> > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> > Jose Claudio Faria >> > Estatistica - Prof. Pleno >> > UESC/DCET/Brasil >> > joseclaudio.faria at gmail.com >> > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> > >> > >> > >> > Em 12 de agosto de 2011 17:46, Pedro Rafael >> > <[hidden email] <http://user/SendEmail.jtp?type=node&node=3741147&i=1>> >> escreveu: >> >> Pessoal entendi perfeitamente a utilidade de criar novos workspaces, >> contudo >> >> quando se finaliza uma seção do R essa função não irá ser mais >> conhecida. O >> >> trabalho de rodar os comandos para criar um workspace oculto não é o >> mesmo >> >> para rodar as funções? >> >> Pedro >> >> Em 12 de agosto de 2011 17:45, Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=3741147&i=2>> >> >> >> escreveu: >> >>> >> >>> Pessoal entendi perfeitamente a utilidade de criar novos workspaces, >> >>> contudo quando se finaliza uma seção do R essa função não irá ser mais >> >> >>> conhecida. O trabalho de rodar os comandos para criar um workspace >> oculto >> >>> não é o mesmo para rodar as funções? >> >>> Pedro >> >>> >> >>> Em 12 de agosto de 2011 17:28, Leonard Assis [via R-br] >> >>> <[hidden email]<http://user/SendEmail.jtp?type=node&node=3741147&i=3>> >> escreveu: >> >>>> >> >>>> A grosso modo, o JavbaDoc pega seus comentários e transforma em >> >>>> Documentação, principalmente os cabeçalhos de metodos >> >>>> >> >>>> roxygen eu nunca usei, mas deve funcionar parecido >> >>>> lmassis <at> yahoo <dot> com <dot> br >> >>>> assis.leonard <at> gmail <dot> com >> >>>> >> >>>> >> >>>> 2011/8/12 Leandro Marino <[hidden email]> >> >>>>> >> >>>>> Hum... >> >>>>> >> >>>>> Tentei procurar na internet e não entendi muito bem. Para que o >> JavaDoc >> >>>>> e o roxygen servem?! >> >>>>> >> >>>>> Atenciosamente, >> >>>>> Leandro Marino >> >>>>> http://www.leandromarino.com.br (Fotógrafo) >> >>>>> http://est.leandromarino.com.br/Blog (Estatístico) >> >>>>> Cel.: <a href="tel:%2B%2055%2021%209845-7707" target="_blank" >> >>>>> value="+552198457707">+ 55 21 9845-7707 >> >>>>> Cel.: <a href="tel:%2B%2055%2021%208777-7907" target="_blank" >> >>>>> value="+552187777907">+ 55 21 8777-7907 >> >>>>> >> >>>>> >> >>>>> >> >>>>> Em 12 de agosto de 2011 16:08, Leonard Assis <[hidden email]> >> escreveu: >> >>>>>> >> >>>>>> Cara, não conhecia. vou olhar >> >>>>>> lmassis <at> yahoo <dot> com <dot> br >> >>>>>> assis.leonard <at> gmail <dot> com >> >>>>>> >> >>>>>> >> >>>>>> 2011/8/12 Gustavo Carvalho <[hidden email]> >> >>>>>>> >> >>>>>>> O roxygen é um pouco parecido, não? >> >>>>>>> http://roxygen.org/#documentation >> >>>>>>> >> >>>>>>> 2011/8/12 Leonard Assis <[hidden email]>: >> >>>>>>> > Ahhh, bem que o R poderia ter um treco parecido com o JavaDoc, >> seria >> >>>>>>> > perfeito. >> >>>>>>> > lmassis <at> yahoo <dot> com <dot> br >> >>>>>>> > assis.leonard <at> gmail <dot> com >> >>>>>>> > >> >>>>>>> > >> >>>>>>> > 2011/8/12 Jose Claudio Faria <[hidden email]> >> >>>>>>> >> >> >>>>>>> >> Leandro, >> >>>>>>> >> >> >>>>>>> >> Chegar ao nível de pacote (sem ou com documentação) para uso >> >>>>>>> >> (pessoal >> >>>>>>> >> ou coletivo), sem dúvidas é o ideal, embora mais trabalhoso! >> >>>>>>> >> >> >>>>>>> >> Todos os pacotes que tenho participação (bpca, TinnR, fdth, >> >>>>>>> >> ScottKnott) começaram como funções simples que foram sendo >> usadas e >> >>>>>>> >> testadas ao longo do tempo. Todas elas inicialmente foram >> usadas >> >>>>>>> >> por >> >>>>>>> >> longo tempo de uma das 4 formas que citei no post original. >> >>>>>>> >> >> >>>>>>> >> Para várias soluções, chega um dia em que você decide dar uma >> >>>>>>> >> polida >> >>>>>>> >> final, empacotar e disponibilizar na vitrine da CRAN. >> >>>>>>> >> >> >>>>>>> >> Eu particularmente não gosto muito dos pacotes Xmisc, prefiro >> as >> >>>>>>> >> soluções mais categorizadas, pois acredito que o futuro seja >> por >> >>>>>>> >> ai. >> >>>>>>> >> >> >>>>>>> >> Abs, >> >>>>>>> >> -- >> >>>>>>> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> >>>>>>> >> Jose Claudio Faria >> >>>>>>> >> Estatistica - Prof. Pleno >> >>>>>>> >> UESC/DCET/Brasil >> >>>>>>> >> joseclaudio.faria at gmail.com >> >>>>>>> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> >>>>>>> >> >> >>>>>>> >> Em 11 de agosto de 2011 21:26, Leandro Marino >> >>>>>>> >> <[hidden email]> escreveu: >> >>>>>>> >> > Gosto da possibilidade de pacotes, apesar de ainda não ter >> feito >> >>>>>>> >> > nenhum. >> >>>>>>> >> > nos >> >>>>>>> >> > meus primeiros testes todos estao totalmente NÃO >> documentados... >> >>>>>>> >> > >> >>>>>>> >> > Enfim... Como são funções quase auto-explicativas (só eu e o >> >>>>>>> >> > Prof. Ruben >> >>>>>>> >> > entendemos) =P >> >>>>>>> >> > >> >>>>>>> >> > O problema é que ainda temos que generalizar muitas delas... >> >>>>>>> >> > >> >>>>>>> >> > Atenciosamente, >> >>>>>>> >> > Leandro Marino >> >>>>>>> >> > http://www.leandromarino.com.br (Fotógrafo) >> >>>>>>> >> > http://est.leandromarino.com.br/Blog (Estatístico) >> >>>>>>> >> > Cel.: <a href="tel:%2B%2055%2021%209845-7707" target="_blank" >> >> >>>>>>> >> > value="+552198457707">+ 55 21 9845-7707 >> >>>>>>> >> > Cel.: <a href="tel:%2B%2055%2021%208777-7907" target="_blank" >> >> >>>>>>> >> > value="+552187777907">+ 55 21 8777-7907 >> >>>>>>> >> > >> >>>>>>> >> > >> >>>>>>> >> > >> >>>>>>> >> > Em 11 de agosto de 2011 21:05, Jose Claudio Faria >> >>>>>>> >> > <[hidden email]> escreveu: >> >>>>>>> >> >> >> >>>>>>> >> >> Eu venho cedendo a essa tentação a tantos anos... ;) >> >>>>>>> >> >> (Mas minhas funções de uso pessoal são bem documentadas) >> >>>>>>> >> >> >> >>>>>>> >> >> Abs, >> >>>>>>> >> >> -- >> >>>>>>> >> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> >>>>>>> >> >> Jose Claudio Faria >> >>>>>>> >> >> Estatistica - Prof. Pleno >> >>>>>>> >> >> UESC/DCET/Brasil >> >>>>>>> >> >> joseclaudio.faria at gmail.com >> >>>>>>> >> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> >>>>>>> >> >> >> >>>>>>> >> >> Em 11 de agosto de 2011 20:17, Benilton Carvalho >> >>>>>>> >> >> <[hidden email]> escreveu: >> >>>>>>> >> >> > Exatamente, Leonard, esse e' o modo que eu considero >> "correto" >> >>>>>>> >> >> > para >> >>>>>>> >> >> > gerenciar coisas assim. >> >>>>>>> >> >> > >> >>>>>>> >> >> > Mas ha' vezes nas quais cedemos `a tentacao... :-) >> >>>>>>> >> >> > _______________________________________________ >> >>>>>>> >> >> > R-br mailing list >> >>>>>>> >> >> > [hidden email] >> >>>>>>> >> >> > 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 >> >>>>>>> >> >> [hidden email] >> >>>>>>> >> >> 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 >> >>>>>>> >> > [hidden email] >> >>>>>>> >> > 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 >> >>>>>>> >> [hidden email] >> >>>>>>> >> 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 >> >>>>>>> > [hidden email] >> >>>>>>> > 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 >> >>>>>>> [hidden email] >> >>>>>>> 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 >> >>>>>> [hidden email] >> >>>>>> 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 >> >>>>> [hidden email] >> >>>>> 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 >> >>>> [hidden email] >> >>>> 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. >> >>>> >> >>>> ________________________________ >> >>>> If you reply to this email, your message will be added to the >> discussion >> >>>> below: >> >>>> >> >>>> >> http://r-br.2285057.n4.nabble.com/R-br-Rprofile-no-Linux-tp3736241p3740402.… >> >>>> To unsubscribe from R-br, click here. >> >>> >> >>> >> >>> -- >> >>> Saudações, >> >>> Pedro Rafael Diniz Marinho - Estatístico SES-PB >> >>> Currículo >> >>> Lattes: >> http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 >> >> >> >> >> >> >> >> -- >> >> Saudações, >> >> Pedro Rafael Diniz Marinho - Estatístico SES-PB >> >> Currículo >> >> Lattes: >> http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 >> >> >> >> _______________________________________________ >> >> R-br mailing list >> >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3741147&i=4> >> >> 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 >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3741147&i=5> >> 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. >> >> >> ------------------------------ >> If you reply to this email, your message will be added to the discussion >> below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Rprofile-no-Linux-tp3736241p3741147.… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

3 3

Sweave
by jorgeas＠ufcg.edu.br 15 Ago '11

15 Ago '11

Pessoal tenho tentado usar o Sweave, mas não dá certo onde estou errando, o arquivo que estou tentando rodar é este utilizando o TexnicCenter e R 2.13.1 \documentclass[a4paper]{article} \title{Sweave Example 1} \author{Friedrich Leisch} \begin{document} \maketitle In this example we embed parts of the examples from the \texttt{kruskal.test} help page into a \LaTeX{} document: <<>>= data(airquality) library(stats) kruskal.test(Ozone ~ Month, data = airquality) @ which shows that the location parameter of the Ozone distribution varies significantly from month to month. Finally we include a boxplot of the data: \begin{center} <<fig=TRUE,echo=FALSE>>= boxplot(Ozone ~ Month, data = airquality) @ \end{center} \end{document} e resulta neste anexo

2 2

calculo de volume de função numérica
by julio cesar araujo 14 Ago '11

14 Ago '11

Boa tarde, Tenho uma função, uma forma gráfica que não tem forma analítica, que é avaliada num grid (construi de 100 mas posso fazer mais fino caso queira). Estou precisando calcular o volume do gráfico dessa função para comparar com outro gráfico que possuo similar. Alguem tem alguma sugestão??? Existe alguma função que faça algum tipo de interpolação...etc? Desde já agradeço Julio

3 3

Re: [R-br] .Rprofile no Linux
by Pedro Rafael 13 Ago '11

13 Ago '11

Pessoal entendi perfeitamente a utilidade de criar novos workspaces, contudo quando se finaliza uma seção do R essa função não irá ser mais conhecida. O trabalho de rodar os comandos para criar um workspace oculto não é o mesmo para rodar as funções? Pedro Em 12 de agosto de 2011 17:45, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Pessoal entendi perfeitamente a utilidade de criar novos workspaces, > contudo quando se finaliza uma seção do R essa função não irá ser mais > conhecida. O trabalho de rodar os comandos para criar um workspace oculto > não é o mesmo para rodar as funções? > > Pedro > > Em 12 de agosto de 2011 17:28, Leonard Assis [via R-br] < > ml-node+3740402-1505413747-223914(a)n4.nabble.com> escreveu: > > A grosso modo, o JavbaDoc pega seus comentários e transforma em >> Documentação, principalmente os cabeçalhos de metodos >> >> roxygen eu nunca usei, mas deve funcionar parecido >> >> lmassis <at> yahoo <dot> com <dot> br >> assis.leonard <at> gmail <dot> com >> >> >> 2011/8/12 Leandro Marino <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=0> >> > >> >>> *Hum... >>> >>> Tentei procurar na internet e não entendi muito bem. Para que o JavaDoc e >>> o roxygen servem?! >>> >>> * >>> Atenciosamente, >>> Leandro Marino >>> http://www.leandromarino.com.br (Fotógrafo) >>> http://est.leandromarino.com.br/Blog (Estatístico) >>> Cel.: <a href="tel:%2B%2055%2021%209845-7707" target="_blank" value=" >>> +552198457707">+ 55 21 9845-7707 >>> Cel.: <a href="tel:%2B%2055%2021%208777-7907" target="_blank" value=" >>> +552187777907">+ 55 21 8777-7907 >>> >>> >>> >>> Em 12 de agosto de 2011 16:08, Leonard Assis <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=1> >>> > escreveu: >>> >>> Cara, não conhecia. vou olhar >>>> >>>> lmassis <at> yahoo <dot> com <dot> br >>>> assis.leonard <at> gmail <dot> com >>>> >>>> >>>> 2011/8/12 Gustavo Carvalho <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=2> >>>> > >>>> >>>>> O roxygen é um pouco parecido, não? >>>>> http://roxygen.org/#documentation >>>>> >>>>> 2011/8/12 Leonard Assis <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=3> >>>>> >: >>>>> > Ahhh, bem que o R poderia ter um treco parecido com o JavaDoc, seria >>>>> > perfeito. >>>>> > lmassis <at> yahoo <dot> com <dot> br >>>>> > assis.leonard <at> gmail <dot> com >>>>> > >>>>> > >>>>> > 2011/8/12 Jose Claudio Faria <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=4> >>>>> > >>>>> >> >>>>> >> Leandro, >>>>> >> >>>>> >> Chegar ao nível de pacote (sem ou com documentação) para uso >>>>> (pessoal >>>>> >> ou coletivo), sem dúvidas é o ideal, embora mais trabalhoso! >>>>> >> >>>>> >> Todos os pacotes que tenho participação (bpca, TinnR, fdth, >>>>> >> ScottKnott) começaram como funções simples que foram sendo usadas e >>>>> >> testadas ao longo do tempo. Todas elas inicialmente foram usadas por >>>>> >> longo tempo de uma das 4 formas que citei no post original. >>>>> >> >>>>> >> Para várias soluções, chega um dia em que você decide dar uma polida >>>>> >> final, empacotar e disponibilizar na vitrine da CRAN. >>>>> >> >>>>> >> Eu particularmente não gosto muito dos pacotes Xmisc, prefiro as >>>>> >> soluções mais categorizadas, pois acredito que o futuro seja por ai. >>>>> >> >>>>> >> Abs, >>>>> >> -- >>>>> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >>>>> >> Jose Claudio Faria >>>>> >> Estatistica - Prof. Pleno >>>>> >> UESC/DCET/Brasil >>>>> >> joseclaudio.faria at gmail.com >>>>> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >>>>> >> >>>>> >> Em 11 de agosto de 2011 21:26, Leandro Marino >>>>> >> <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=5>> >>>>> escreveu: >>>>> >> > Gosto da possibilidade de pacotes, apesar de ainda não ter feito >>>>> nenhum. >>>>> >> > nos >>>>> >> > meus primeiros testes todos estao totalmente NÃO documentados... >>>>> >> > >>>>> >> > Enfim... Como são funções quase auto-explicativas (só eu e o Prof. >>>>> Ruben >>>>> >> > entendemos) =P >>>>> >> > >>>>> >> > O problema é que ainda temos que generalizar muitas delas... >>>>> >> > >>>>> >> > Atenciosamente, >>>>> >> > Leandro Marino >>>>> >> > http://www.leandromarino.com.br (Fotógrafo) >>>>> >> > http://est.leandromarino.com.br/Blog (Estatístico) >>>>> >> > Cel.: <a href="tel:%2B%2055%2021%209845-7707" target="_blank" >>>>> value="+552198457707">+ 55 21 9845-7707 >>>>> >> > Cel.: <a href="tel:%2B%2055%2021%208777-7907" target="_blank" >>>>> value="+552187777907">+ 55 21 8777-7907 >>>>> >>>>> >> > >>>>> >> > >>>>> >> > >>>>> >> > Em 11 de agosto de 2011 21:05, Jose Claudio Faria >>>>> >> > <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=6>> >>>>> escreveu: >>>>> >> >> >>>>> >> >> Eu venho cedendo a essa tentação a tantos anos... ;) >>>>> >> >> (Mas minhas funções de uso pessoal são bem documentadas) >>>>> >> >> >>>>> >> >> Abs, >>>>> >> >> -- >>>>> >> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >>>>> >> >> Jose Claudio Faria >>>>> >> >> Estatistica - Prof. Pleno >>>>> >> >> UESC/DCET/Brasil >>>>> >> >> joseclaudio.faria at gmail.com >>>>> >> >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >>>>> >> >> >>>>> >> >> Em 11 de agosto de 2011 20:17, Benilton Carvalho >>>>> >> >> <[hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=7>> >>>>> escreveu: >>>>> >> >> > Exatamente, Leonard, esse e' o modo que eu considero "correto" >>>>> para >>>>> >> >> > gerenciar coisas assim. >>>>> >> >> > >>>>> >> >> > Mas ha' vezes nas quais cedemos `a tentacao... :-) >>>>> >> >> > _______________________________________________ >>>>> >> >> > R-br mailing list >>>>> >> >> > [hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=8> >>>>> >> >> > 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 >>>>> >> >> [hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=9> >>>>> >> >> 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 >>>>> >> > [hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=10> >>>>> >> > 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 >>>>> >> [hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=11> >>>>> >> 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 >>>>> > [hidden email]<http://user/SendEmail.jtp?type=node&node=3740402&i=12> >>>>> > 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 >>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3740402&i=13> >>>>> 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 >>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3740402&i=14> >>>> 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 >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3740402&i=15> >>> 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 >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3740402&i=16> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the discussion >> below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Rprofile-no-Linux-tp3736241p3740402.… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 2

Análise Envoltória de Dados (DEA)
by Daniel Dantas 13 Ago '11

13 Ago '11

Caros, Novamente busco auxílio. Coloquei nesta madrugada o R para calcular a eficiência DEA com uma base de dados de 5000 x 5 mas me surgiram os seguintes avisos: > spx_invert: the basis matrix is singular spx_simplex: numerical problems with basis matrix spx_simplex: sorry, basis recovery procedure not implemented yet Talvez possa ser problemas com os dados ou minha rotina está muito demorada e precisa ser revista ou algum outro problema desconhecido para os meus conhecimentos atuais. Abaixo segue a rotina: x = matriz de insumos (5592 x 4) y = matriz de produtos (5592 x 7) #------------------- # Calculando DEA #------------------- library("DEA") ### SBM BCC - Sem Orientacao write.csv(data.frame(base$base.COD_ID,as.data.frame(dea.sbm.bcc(x,y))),"/host/DEA/DEA_SBM_BCC.csv") Com esse comando deixei horas rodando e o que me retornou foram aqueles 3 avisos no terminal. PS: Usei o Rcmr no Ubuntu. O que pode ter acontecido? Obrigado, Daniel

1 0

Re: [R-br] Operações com matrizes
by Pedro Rafael 12 Ago '11

12 Ago '11

Quem é a? Pedro Em 12 de agosto de 2011 20:24, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > em é a? > > > Pedro > > Em 12 de agosto de 2011 18:33, Edmundo Caetano [via R-br] < > ml-node+3740510-1898737130-223914(a)n4.nabble.com> escreveu: > > Boa noite! Estou com um problema, pretendo determinar o valor de *a* >> de modo que a matriz A dada abaixo seja singular. Como tenho que proceder >> no R???? >> [,1] [,2] >> [1,] 1 8 >> [2,] 3 a >> >> >> >> *Edmundo do R. R. Caetano * >> >> >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3740510&i=0> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the discussion >> below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Operacoes-com-matrizes-tp3740510p374… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

teste wald glm
by Vinicius Brito Rocha 12 Ago '11

12 Ago '11

Senhores o padrão de teste de significância um ajuste glm com link binomial ou quasibinomial é um teste t Como posso realizar um teste de wald nos coeficientes? Att. -- 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.

2 2

Operações com matrizes
by Edmundo Caetano 12 Ago '11

12 Ago '11

Boa noite! Estou com um problema, pretendo determinar o valor de a de modo que a matriz A dada abaixo seja singular. Como tenho que proceder no R???? [,1] [,2] [1,] 1 8 [2,] 3 a Edmundo do R. R. Caetano

1 0

[Dúvidas] Vantagens de rodar o R no Linux.
by Pedro Rafael 12 Ago '11

12 Ago '11

Não estou discutindo qual é o melhor sistema, não quero cair em assuntos que Linux é melhor que Windows ou o contrário. Gostaria de saber quais as vantagens do Linux com base em um usuário que prioriza o uso do R em sua máquina. Quais vantagens computacionais eu teria no Linux que não tenho no Windows para executar um script em R? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

4 5

.Rprofile no Linux
by Jose Claudio Faria 12 Ago '11

12 Ago '11

Olá pessoal, Este "post" tem a finalidade de compartilhar algumas coisas que considero importantes na inicialização do R. É voltado para o Linux, mas os princípios, com os devidos cuidados, podem ser usados também no Windows. Muitas das opções importantes, do ponto de vista funcional (não relativos à aparência), podem ser feitas no arquivo .Rprofile. Esse arquivo deve ficar localizado no home do usuário: ~/.Rprofile e é um dos primeiros a ser lido quando uma sessão do R é iniciada. Tenho uma função (bem simples) que uso bastante em meu dia a dia: "cv" para calcular o coef. de variação de uma ANOVA: cv <- function(av) { if(is.null(av) || !inherits(av, 'aov')) stop('Please, check the parameter!') qmee <- with(av, sum(residuals^2) / df.residual) cv <- 100 * sqrt(qmee) / mean(av$fitted.values) return(round(cv, 2)) } Pois bem, a medida que vamos aumentando nossa intimidade com o R (inevitavelmente) iremos desenvolvendo nossas próprias funções (o R foi projetado para isso!). Ai vem o problema de ter que SEMPRE carregar a função quando for usar, o que pode se tornar uma chatisse. Pior ainda, ao limpar o workspace do usuário (.GlobalEnv) elas são removidas e precisam ser recarregadas. Tem como contornar? Sim, de várias formas: 1. Opção muito pouco prática: oldp <- getwd() setwd('/home/jcfaria/dados/r/funcoes/') source('cv.r') setwd(oldp) A função "cv" ficará disponível no workspace do usuário mas será removida com a instrução: > rm(list=ls()) muito usada por várias GUIs. 2. Opção "mais" prática: A mesma que a anterior, contudo, a função não deverá se chamar "cv", mas sim ".cv". Nesse caso ela permanecerá como um objeto oculto no workspace e não será removida com a intrução: > rm(list=ls()) Contudo poderá ser removida com a intrução: > rm(list=ls(all=TRUE)) 3. Colocando suas funções em algum ambiente (environment) do R (optei pelo base) oldp <- getwd() setwd('/home/jcfaria/dados/r/funcoes/') source('cv.r', local=baseenv()) setwd(oldp) 4. Criando seu próprio ambiente (acho a solução mais elegante): .env <- new.env() oldp <- getwd() setwd('/home/jcfaria/dados/r/funcoes/') source('cv.r', local=.env) setwd(oldp) Nesse último caso (4): - O objeto ".env" ficará oculto no seu workspace evitando ser deletado com rm(list=ls()) - Parar acessar a função "cv" será necessário > .env$cv Por exemplo: > av <- aov(Sepal.Length ~ Species, data=iris) > .env$cv(av) [1] 8.81 Esta forma de carregar as funções de forma permanente no R pode ser usado para qualquer objeto. No Windows bastava usar: source('cv', local=TRUE) que ela fica disponível no pacote base, contudo, não testei na versão mais nova. Proveito para disponibilizar meu .Rprofile que poderá ser útil para os que estão começando ## José Cláudio Faria/UESC/DCET ##¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬ ## General options library(utils) options(list(repos='http://cran-r.c3sl.ufpr.br/', width=50, editor='vim', browser='chromium')) ##¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬ ## Package debug: options options(list(debug.height=10, debug.width=125, debug.screen.pos='-1680+00')) # PARA DOIS MONITORES ##¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬ ## Functions oldp <- getwd() setwd('/home/jcfaria/dados/r/funcoes/') #source('fitreg.r', local=baseenv()) #source('fitregl.r', local=baseenv()) #source('fitrlcor.r', local=baseenv()) #source('plotreg.r', local=baseenv()) #source('plotrl.r', local=baseenv()) #source('plotlm3d.r', local=baseenv()) #source('cv.r', local=baseenv()) #source('rs.r', local=baseenv()) #source('sourcedir.r', local=baseenv()) .env <- new.env() source('fitreg.r', local=.env) source('fitregl.r', local=.env) source('fitrlcor.r', local=.env) source('plotreg.r', local=.env) source('plotrl.r', local=.env) source('plotlm3d.r', local=.env) source('cv.r', local=.env) source('rs.r', local=.env) source('sourcedir.r', local=.env) ##¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬ ## Removing variables setwd(oldp) rm(oldp) ##¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬¬ ## Packages: loading #library(fdth) #library(ScottKnott) library(bpca) library(reshape) library(debug) Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

6 22

Pacotes e C, Fortran ou outros
by Leandro Marino 12 Ago '11

12 Ago '11

Caros, Com a discussão anterior sobre o rprofile resolvi me aventurar em terreno completamente desconhecido. Alguém teria uma boa referência para iniciar a programação em qualquer outra linguagem que possa ser utilizada dentro do R? Tenho algumas funções que faço no R, no entanto, acredito que não estão eficientes como deveriam ser. Sendo assim, gostaria de levá-las para alguma outra linguagem e usá-las eventualmente dentro do R. Alguma sugestão de linguagem para começar meus estudos? E será que alguém tem material que possa ser enviado?! Grato, Leandro Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

5 6

Instalação do R no win, chave de registro e R2wd
by Pedro Emmanuel Alvarenga Americano do Brasil 12 Ago '11

12 Ago '11

Amigos de R, Acho que já mandei uma mensagem pareceida no passado, mas infelizmente ainda estou sofrendo do mesmo problema. Ja mandei uma mensagem para o desenvolvedor e ele recomendou mandar a mensagem pra uma outra lista que ninguém respondeu. Instalei o pacote R2wd, mas quando eu chamo o pacote com library('R2wd') ha um aviso que diz ser necessário instalar um programa chamado statconnDCOM. Depois da instalação que me parece bem sucedida, quando eu rodo o R de novo e chamo o R2wd novamente vem o mesmo aviso como se o programa não estivesse instalado e as funções do pacote não funcionam. Lendo a documentação do statconnDCOM, ha um interface de teste que parece estar ok. Mas duas recomendações na instalação desse programa que são ncessarias, a primeira é que as dependencias estejam instaldas dentro do diretorio do R e não em pastas pessoais e a segunda que durante a instalação do R, ha necessidade de marcar a opção de guardar a chave de registro da versão instalada do R. Eu rodei a instalação do R novamente e nao consegui perceber essa opção durante a instalação. A pergunta é, alguém saberia me dizer como eu faço pra marcar essa opção durante a instalação, ou como eu faço pra verificar que essa chave foi guardada no sistema? windows 7 home edition R 2.13.1 word 2003 , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Fwd: Eliminar NA de uma soma
by Pedro Rafael 12 Ago '11

12 Ago '11

Uma outra forma de fazer é: sum(na.omit(x)) Saudações, Pedro Rafael Diniz Marinho Em 12 de agosto de 2011 12:07, Carlos Mendonça [via R-br] < ml-node+3739556-550071832-223914(a)n4.nabble.com> escreveu: > > Em 12 de agosto de 2011 11:38, Carlos Mendonça <[hidden email]<http://user/SendEmail.jtp?type=node&node=3739556&i=0> > > escreveu: > > >> Gustavo e Benilton, muito obrigado pela ajuda. >> > > >> Carlos Mendonça. >> > > _______________________________________________ > R-br mailing list > [hidden email] <http://user/SendEmail.jtp?type=node&node=3739556&i=1> > 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. > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Eliminar-NA-de-uma-soma-tp3739464p37… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 1

apagar valores de um dataframe e colocar NA neles
by Luis Iván Ortiz Valencia 12 Ago '11

12 Ago '11

Galera Tenho um data frame com valores, e preciso apagar esses valores e colocar NA neles.....qual a melhor forma? Ivan

8 15

Eliminar NA de uma soma
by Carlos Mendonça 12 Ago '11

12 Ago '11

Caros, como faço para eliminar o valor NA do resultado de uma soma, pois se tenho algum variável com NA, quando somo essa variável o resultado obtido é NA, mesmo tendo vários registros com valores válidos. Carlos Mendonça

3 3

Re: [R-br] Digest R-br, volume 6, assunto 16
by Julio Bueno 12 Ago '11

12 Ago '11

Obrigado Benilton e demais que comentaram sobre o is.na() e o subsetting. A aplicação é importante pra programas atualizam registros com erro de atualização... matrizes do tipo mil linhas e dez mil colunas... Mas para além disso, a idéia do subsetting em outros programas sempre me surpreende, como na sua solução, mais elegante que a solução 1! Obrigado também pela dica de postagem... o usuário aqui postou duas vezes e não percebeu que o padrão de anexar texto do gmail era o problema :-). []s Júlio -- Júlio Sílvio de Sousa Bueno Filho Departamento de Ciências Exatas Universidade Federal de Lavras Caixa Postal 3037 Lavras, MG - Brasil - 37200 000 tel: (+55) 35 3829 1369 fax: (+55) 35 3829 1371 "Natureza da gente não cabe em nenhuma certeza." João G. Rosa - Riobaldo (Grande Sertão: Veredas) "Homo sum; humani nihil a me alienum puto." Terêncio - Chremes (Heautontimoroumenos)

1 0

Re: [R-br] Dúvida - Operadores %any%, :: e :::
by Pedro Rafael 11 Ago '11

11 Ago '11

Obrigado senhores pelas dicas, e desculpem a demora para responder. Estou aprendendo bastante sobre R nesta comunidade. Abraço à todos. Em 11 de agosto de 2011 17:05, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > > > Em 1 de agosto de 2011 18:52, jcfaria [via R-br] < > ml-node+3710920-1613572389-223914(a)n4.nabble.com> escreveu: > > Complementando as dicas do Benilton, >> >> Em programação orientada a objetos (R é uma delas e usa atualmente >> dois sistemas: S3 e S4) >> algumas vezes o desenvolvedor opta por encapsular um objeto. >> >> Encapsular um objeto (no caso uma função) significa ocultar do usuário >> detalhes usados no desenvolvimento dessa função. >> >> Por exemplo, as funções abaixo (plot.acf*, plot.data.frame*, >> plot.factor*, plot.formula*, etc) estão todas encapsuladas. >> Observe o marcador "*" em todas elas! >> >> > methods(plot) >> [1] plot.acf* plot.data.frame* plot.decomposed.ts* >> [4] plot.default plot.dendrogram* plot.density >> [7] plot.ecdf plot.factor* plot.formula* >> [10] plot.hclust* plot.histogram* plot.HoltWinters* >> [13] plot.isoreg* plot.lm plot.medpolish* >> [16] plot.mlm plot.ppr* plot.prcomp* >> [19] plot.princomp* plot.profile.nls* plot.spec >> [22] plot.spec.coherency plot.spec.phase plot.stepfun >> [25] plot.stl* plot.table* plot.ts >> [28] plot.tskernel* plot.TukeyHSD >> >> Non-visible functions are asterisked >> >> Se você digita o nome de uma delas para ver seu código fonte: >> > plot.data.frame >> >> Recebe a mensagem: Erro: objeto 'plot.data.frame' não encontrado! >> Ou seja, esta função está encapsulada (não visível ao usuário comum). >> >> O operador ":::" lhe permite ter acesso ao código fonte de qualquer >> função encapsulada (não visível): >> >> > graphics:::plot.data.frame >> function (x, ...) >> { >> plot2 <- function(x, xlab = names(x)[1L], ylab = names(x)[2L], >> ...) plot(x[[1L]], x[[2L]], xlab = xlab, ylab = ylab, >> ...) >> if (!is.data.frame(x)) >> stop("'plot.data.frame' applied to non data frame") >> if (ncol(x) == 1) { >> x1 <- x[[1L]] >> cl <- class(x1) >> if (cl %in% c("integer", "numeric")) >> stripchart(x1, ...) >> else plot(x1, ...) >> } >> else if (ncol(x) == 2) { >> plot2(x, ...) >> } >> else { >> pairs(data.matrix(x), ...) >> } >> } >> <environment: namespace:graphics> >> >> >> Posso estar equivocado e gostaria de escutar a opinião da comunidade, >> mas não vejo muito sentido em usar encapsulamento numa linguagem com >> as características do R. >> >> Abs, >> -- >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> Jose Claudio Faria >> Estatistica - Prof. Pleno >> UESC/DCET/Brasil >> joseclaudio.faria at gmail.com >> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3710920&i=0> >> 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. >> >> >> ------------------------------ >> If you reply to this email, your message will be added to the discussion >> below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Operadores-any-e-tp3710751p37… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 1

formato Eps, Wmf
by Gilson Sanchez 11 Ago '11

11 Ago '11

Pessoal do R, mais uma vez eu procurando ajuda de vc's, esteve procurando informação ni historial do R-Br mas nao encontrei muita informação, então como eu faço para mudar a extensão do formato dos graficos para eps o wmf. agradeço desde já! gilson -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

2 1

[Dúvida] Rodando Script R no terminal sem entrar no prompt.
by Carlos Mendonça 11 Ago '11

11 Ago '11

A minha dúvida é em relação ao Windows, como faço para executar um script do R sem ser necessário "entar no R"? -- Um abraço, Carlos Mendonça

4 7

Re: [R-br] Digest R-br, volume 6, assunto 12
by Julio Bueno 11 Ago '11

11 Ago '11

Pessoal, Escrevo porque uma combinação de mensagens me chamou a atenção No problema da substituição (mensagem 15) a solução é a função específica is.na() apresentada pelo Paulo na mensagem 17,,, e alguém falou sobre o uso do subsetting, que eu tenho usado com muita frequência. Um uso desta função com "subsetting" seria: x[is.na(x)==1] <- 0.1 Alguém sabe dizer a diferença em termo sde número de cálculos e tempo entre fazer o subsetting e o laço ifelse? []s Júlio 2011/8/10 <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 > > Par: > a 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. Re: Obter coordenadas a partir de endereço > (Pedro Emmanuel Alvarenga Americano do Brasil) > 2. Re: Obter coordenadas a partir de endereço (Benilton Carvalho) > 3. automatização de ajuste de modelos de regressao (Samuel Carvalho) > 4. Re: automatização de ajuste de modelos de regressao > (Ivan Bezerra Allaman) > 5. Re: automatização de ajuste de modelos de regressao > (Benilton Carvalho) > 6. Re: automatização de ajuste de modelos de regressao > (Henrique Dallazuanna) > 7. [Dúvida] Funções para alto filtros. (Pedro Rafael) > 8. apagar valores de um dataframe e colocar NA neles > (Luis Iván Ortiz Valencia) > 9. Accuracy 2012 - First Call for Papers - July 10-13, 2012, > Florianópolis, SC - Brazil (RODRIGO LILLA MANZIONE) > 10. Re: apagar valores de um dataframe e colocar NA neles > (Paulo Justiniano) > 11. Re: apagar valores de um dataframe e colocar NA neles > (Luis Iván Ortiz Valencia) > 12. Re: funcao para valores iguais (Luis Iván Ortiz Valencia) > 13. Usar peso para tabela cruzada (Clayton Santos Delfino) > 14. Substituir NA's (Daniel Dantas) > 15. Re: Substituir NA's (Gustavo Carvalho) > 16. Re: Substituir NA's (Daniel Dantas) > 17. Re: Substituir NA's (Paulo Justiniano) > 18. Re: Substituir NA's (Edson Lira) > 19. Re: Substituir NA's (Eder David Borges da Silva) > 20. Re: Substituir NA's (Daniel Dantas) > 21. Re: Substituir NA's (eder(a)leg.ufpr.br) > 22. Re: [Dúvida] Funções para alto filtros. (Marcos Silva) > 23. Re: Substituir NA's (Paulo Justiniano) > 24. Re: Usar peso para tabela cruzada (Jakson Alves de Aquino) > 25. Material Novo no CRAN (Marcos Silva) > 26. Re: [Dúvida] Rodando Script R no terminal sem entrar no > prompt. (Carlos Mendonça) > 27. Re: Accuracy 2012 - First Call for Papers - July 10-13, 2012, > Florianópolis, SC - Brazil (Leonard Assis) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 10 Aug 2011 12:11:17 -0300 > From: Pedro Emmanuel Alvarenga Americano do Brasil > <emmanuel.brasil(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Obter coordenadas a partir de endereço > Message-ID: > <CAFfGvy+rMiS0DrNEY2GLCKU9F9Bq5b7wojGSXrDQXJ6wNr5wPg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Amigos de R, > > Nao trabalho mais com espacial ha muitos e muitos anos, porem fiquei > curioso > com essa função, será que essa função exige que o os endereços seja > formatados sempre do mesmo jeito? Estou lembrando de um discussão que tive > quase dez anos atras por conta de muita falta de dados em > georeferenciamento > para fazer mapas de dengue no rio de janeiro. QUando havia dados de > endereços eles eram irregualres. Ou seja, nomes de ruas soletrados de forma > incorreta, as vezes o numero da casa era separado com virgula outras vezes > com ponto do nome da rua, e muitas vezes não havia separador entre o numero > do apartamento e o numero da predio. Fiquei pensando se talvez o R consiga > buscar o CEP do endereço no site dos correios, ou se é possivel baixar dos > correios um banco de dados com CEP, e tnetar corresponder esses endereços > com CEP. Nao seria mais facil? > > > Abraço forte e que a força esteja com você, > > Dr. Pedro Emmanuel A. A. do Brasil > Instituto de Pesquisa Clínica Evandro Chagas > Fundação Oswaldo Cruz > Rio de Janeiro - Brasil > Av. Brasil 4365 > Tel 55 21 3865-9648 > email: pedro.brasil(a)ipec.fiocruz.br > email: emmanuel.brasil(a)gmail.com > > ---Apoio aos softwares livres > www.zotero.org - gerenciamento de referências bibliográficas. > www.broffice.org ou www.libreoffice.org - textos, planilhas ou > apresentações. > www.epidata.dk - entrada de dados. > www.r-project.org - análise de dados. > www.ubuntu.com - sistema operacional > > > Em 9 de agosto de 2011 08:49, Thiago Veloso <thi_veloso(a)yahoo.com.br > >escreveu: > > > Luís Gustavo, > > > > Muito obrigado pela sua dica. Com ela consegui fazer a conversão que > > necessitava. > > > > Agradeço também aos demais participantes que responderam a minha > dúvida. > > > > Um abraço, > > > > Thiago Veloso. > > > > --- On *Fri, 5/8/11, Luís Gustavo <lgsilvaesilva(a)gmail.com>* wrote: > > > > > > From: Luís Gustavo <lgsilvaesilva(a)gmail.com> > > > > Subject: Re: [R-br] Obter coordenadas a partir de endereço > > To: r-br(a)listas.c3sl.ufpr.br > > Date: Friday, 5 August, 2011, 12:04 > > > > > > Prezado Thiago, > > > > Se não estou enganado você consegue as coordendas a partir do CEP, da > > seguinte forma: > > > > > > *require(XML) > > > > coordenadas<- function(cep) { > > url_lat_lon <- paste(sprintf(" > > http://maps.google.com/maps/api/geocode/xml?address=%s,", > > cep),"%20Brasil&sensor=false", sep="") > > lat_lon=xmlApply(xmlRoot(xmlTreeParse( > > readLines(url_lat_lon)))[['result']][['geometry']][['location']], "[[", > 1) > > return(lat_lon) > > } > > > > cep=36026300 > > coordenadas(cep) > > > > Obs.: Existe posts passados abordando estes assunto, lembro que salvei > esse > > script em alguns desses posts, porém não > > me lembro exatamente quem postou, portanto fica as minhas desculpas pela > a > > falta de referência. > > * > > Em 5 de agosto de 2011 00:19, Daniel Marcelino <dmsilva.br(a)gmail.com< > http://mc/compose?to=dmsilva.br@gmail.com>> > > escreveu: > > > > > > Olá Thiago, > > > Deve haver mais opções para fazer isso, eu conheço o PBSmapping. Há > > inclusive um material vendido pelo O' Reilly's com o passo-a-passo para > > conseguir isso. > > > Daniel > > > > > > 2011/8/4 Thiago Veloso <thi_veloso(a)yahoo.com.br< > http://mc/compose?to=thi_veloso@yahoo.com.br> > > > > > >> > > >> Olá pessoal, > > >> > > >> Tenho quase certeza que esse assunto já foi discutido aqui na lista, > > porém não estou conseguindo encontrar o tópico. > > >> > > >> Possuo uma lista com cerca de 40 endereços em Porto Alegre. Gostaria > de > > obter a localização geográfica destes endereços para mostrá-la em um > mapa. > > Apesar de ter somente o endereço, com um pouco de esforço é possível > > conseguir o CEP de cada um deles. > > >> > > >> Sendo assim, há algum pacote combinado R/google maps que permita > > conhecer as coordenadas de um endereço a partir do seu CEP? > > >> > > >> Grato desde já, > > >> > > >> Thiago Veloso. > > >> _______________________________________________ > > >> R-br mailing list > > >> R-br(a)listas.c3sl.ufpr.br< > http://mc/compose?to=R-br@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. > > > > > > > > > > > > -- > > > Daniel Marcelino > > > http://danielmarcelino.com > > > Skype: d_marcelino > > > > > > _______________________________________________ > > > R-br mailing list > > > R-br(a)listas.c3sl.ufpr.br < > http://mc/compose?to=R-br@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. > > > > > > > > -- > > Luís Gustavo Silva e Silva > > > > > > -----Inline Attachment Follows----- > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br <http://mc/compose?to=R-br@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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/04ae926b/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 10 Aug 2011 16:39:25 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Obter coordenadas a partir de endereço > Message-ID: > <CAO-arWNLQq1SATXKNcTq7yUpbdg5b-fCNJm0tBru=rQ=2FDKkw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > talvez, inclusive, alguma interacao com o pessoal da pc2: > > http://ceplivre.pc2consultoria.com/ > > b > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/2416ff6e/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 10 Aug 2011 09:23:51 -0700 (PDT) > From: Samuel Carvalho <samukajm(a)yahoo.com.br> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] automatização de ajuste de modelos de regressao > Message-ID: > <1312993431.87925.YahooMailNeo(a)web161201.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros(as) > Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e para > isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar > este modelo. > O que consigo, é fazer este ajuste inserindo o efeito da parcela na > composição do modelo linear. Abaixo segue um CMR > > dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15)) > mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da > parcela no modelo > ou > mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de > fazer de uma maneira mais automatizada > mod2 <- lm(y ~ x, subset(dados, parcela==2)) > mod3<- lm(y ~ x, subset(dados, parcela==3)) > > Pela atenção > Obrigado > > ==================================== > Samuel P. C. Carvalho > Mestre em Ciências Florestais [UFLA] > Doutorando em Recursos Florestais [ESALQ/USP] > ============================================= > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/172ab19a/atta… > > > > ------------------------------ > > Message: 4 > Date: Wed, 10 Aug 2011 10:01:42 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] automatização de ajuste de modelos de regressao > Message-ID: > <1312995702.55404.YahooMailNeo(a)web161820.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > De um modo bem menos elegante do que os nossos amigos da programação, > pode-se ter: > > > length(levels(factor(dados$parcela))) > results <- list() > for(i in 1:length(levels(factor(dados$parcela)))){ > results[[i]] <- lm(y ~ x, subset(dados, parcela==i)) > } > results > lapply(results,summary) > lapply(results,anova) > > Allaman > (S,f,P) > > > > > > Ivan Bezerra Allaman > Doutor em Zootecnia/UFLA > email e msn - ivanalaman(a)yahoo.com.br > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/1e7c4bcc/atta… > > > > ------------------------------ > > Message: 5 > Date: Wed, 10 Aug 2011 18:04:27 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] automatização de ajuste de modelos de regressao > Message-ID: > <CAO-arWOvRoGoHO0yFQhA55ejqsSHZt0spv=UDKnMjXpD+OxVMQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > by(dados, dados[['parcela']], function(z) lm(y~x, data=z)) > > > ------------------------------ > > Message: 6 > Date: Wed, 10 Aug 2011 14:11:19 -0300 > From: Henrique Dallazuanna <wwwhsd(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Samuel Carvalho <samukajm(a)yahoo.com.br> > Subject: Re: [R-br] automatização de ajuste de modelos de regressao > Message-ID: > <CAPvBnPEvmGhh4cOgwi-e0XrfnHVB3kObwqTxOqz-9iMSB8dsoQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Tente assim: > > mod <- lapply(split(dados, dados$parcela), lm, formula = x ~ y) > > 2011/8/10 Samuel Carvalho <samukajm(a)yahoo.com.br>: > > Caros(as) > > Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e > para > > isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar > > este modelo. > > O que consigo, é fazer este ajuste inserindo o efeito da parcela na > > composição do modelo linear. Abaixo segue um CMR > > > > dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15)) > > mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da > parcela > > no modelo > > ou > > mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de > > fazer de uma maneira mais automatizada > > mod2 <- lm(y ~ x, subset(dados, parcela==2)) > > mod3<- lm(y ~ x, subset(dados, parcela==3)) > > > > Pela atenção > > Obrigado > > > > ==================================== > > Samuel P. C. Carvalho > > Mestre em Ciências Florestais [UFLA] > > Doutorando em Recursos Florestais [ESALQ/USP] > > ============================================= > > _______________________________________________ > > 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. > > > > > > -- > Henrique Dallazuanna > Curitiba-Paraná-Brasil > 25° 25' 40" S 49° 16' 22" O > > > ------------------------------ > > Message: 7 > Date: Wed, 10 Aug 2011 15:38:05 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [Dúvida] Funções para alto filtros. > Message-ID: > <CAKwavrm+GG6vLrGpGcbN1Mvskp3ZLBp0TCv3_cUMDAeHNvQxHA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pessoa existe funções analogas ao subset? Achei muito interessante a função > subset para filtrar grandes massas de dados segundo algumas características > desejadas. O R fornece alguma outra função que tenha esse mesmo papel > (fazer > filtros em massa) ou algum outro recurso que seja interessante para fazer > filtros em data.frames enormes? Para mim só basta o me informarem o nome da > função, detalhamentos eu procuro no help. > > Pedro Rafael > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/5a917572/atta… > > > > ------------------------------ > > Message: 8 > Date: Wed, 10 Aug 2011 16:21:41 -0300 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] apagar valores de um dataframe e colocar NA neles > Message-ID: > <CAJ+S-MH=a86GH79wa=Sq4HnDndSVtE5tvC0Cav-kbfp4encjtw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Galera > > Tenho um data frame com valores, e preciso apagar esses valores e colocar > NA > neles.....qual a melhor forma? > > Ivan > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/05450618/atta… > > > > ------------------------------ > > Message: 9 > Date: Wed, 10 Aug 2011 16:42:56 -0300 > From: RODRIGO LILLA MANZIONE <rlmanzione(a)ig.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Accuracy 2012 - First Call for Papers - July 10-13, > 2012, Florianópolis, SC - Brazil > Message-ID: > <CAE2kf_5DdKc_KR+4SC1PumXxhRT6oPRjB+YaEvDQaAYxr9034w(a)mail.gmail.com > > > Content-Type: text/plain; charset="windows-1252" > > *ACCURACY 2012** > International Spatial Accuracy Research Association (ISARA) > The Tenth International Symposium on Spatial Accuracy Assessment in Natural > Resources and Environmental Sciences > July 10-13, 2012 > Florianópolis, SC - Brazil* > > http://2012.spatial-accuracy.org/ > > TOPICS: All topics concerned with spatial accuracy and uncertainty in a > natural resources and environmental sciences context are appropriate, > for example: > > * Semantic uncertainty and vagueness > * Modelling uncertainty using geostatistics > * Propagation of uncertainty in GIS > * Visualizing spatial uncertainty > * Uncertainty in Remote Sensing > * Spatio-temporal uncertainty > * Accuracy and uncertainty of DEMs > * Modelling scale in environmental systems > * Positional uncertainty > > PUBLICATION: We welcome the submission of both oral papers and posters. > For both modes of presentation authors are invited to submit abstracts > of no more than 500 words in English on original research.Those accepted > will be presented as oral papers/posters at the symposium and published as > short papers in the symposium proceedings (ISBN listed). All accepted > authors will also be > invited to submit a full paper for consideration towards an special > issue of an international journal. > > MODES OF PARTICIPATION: Oral papers and posters. Pre-conference > workshop proposals also welcome. > > DEADLINES: > > Abstract submission: *December** 9th, ****2011* > > Abstract acceptance notification: *March** 1st, ****2012* > > Short final paper submission: *May** 11th, ****2012* > > Workshop proposals: *October** 18th, **2011* > > Workshop acceptance notification: *October** 31st, **2011* > It is important to mention that this will be the first time the Spatial > Accuracy Group co-organizes a conference with The Commission on Modelling > Geographical System of the International Geographical Union > (IGU-CMGS)<http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research > area of common interest. > > Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) > Mariane Dal Santo (Secretary) > > Contact: accuracy2012(a)cfh.ufsc.br< > http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf… > > > **************************************** > Download a pdf of the full call for papers at: > http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg > > > > PS: Sorry for e-mail list crossed. > > -- > Prof. Carlos Antonio Oliveira Vieira > Departamento de Geociências ? CFH/UFSC > Trindade - Florianópolis - SC - Brasil > CEP 88040-900 Cx. Postal 476 > Fone: ++55(48)3721-8593 Cel: (48)9915-3653 > E-mails: carlos.vieira(a)ufv.br > carlos.vieira(a)ufsc.br > > Visite o site: www.accuracy2012.ufsc.br > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/03a37f1f/atta… > > > > ------------------------------ > > Message: 10 > Date: Wed, 10 Aug 2011 16:41:53 -0300 (BRT) > From: Paulo Justiniano <paulojus(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] apagar valores de um dataframe e colocar NA neles > Message-ID: <alpine.DEB.2.00.1108101641400.3064(a)pataxo.est.ufpr.br> > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > todos os valores? > > df[,] <- NA > > > On Wed, 10 Aug 2011, Luis Iván Ortiz Valencia wrote: > > > > > Galera > > > > Tenho um data frame com valores, e preciso apagar esses valores e colocar > NA neles.....qual a melhor forma? > > > > Ivan > > > > > > > > > > ------------------------------ > > Message: 11 > Date: Wed, 10 Aug 2011 16:51:23 -0300 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] apagar valores de um dataframe e colocar NA neles > Message-ID: > <CAJ+S-MHyQ3pixjTgc1oPou3X=Z+=GnG0FRBJXdzh0gjUhKQVpA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > sim e mantendo a estrutura do data frame , funcionou, obrigado Paulo. > > Em 10 de agosto de 2011 16:41, Paulo Justiniano <paulojus(a)leg.ufpr.br > >escreveu: > > > todos os valores? > > > > df[,] <- NA > > > > > > > > On Wed, 10 Aug 2011, Luis Iván Ortiz Valencia wrote: > > > > > >> Galera > >> > >> Tenho um data frame com valores, e preciso apagar esses valores e > colocar > >> NA neles.....qual a melhor forma? > >> > >> Ivan > >> > >> > >> > >> > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/1d05403a/atta… > > > > ------------------------------ > > Message: 12 > Date: Wed, 10 Aug 2011 16:52:48 -0300 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] funcao para valores iguais > Message-ID: > <CAJ+S-MEfpPfbzJjT5WTc-xAuZdipD32FYZBCNWVyt0-js_+m-A(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > obrigado....isso que precisava... > > Ivan > > Em 10 de agosto de 2011 11:21, Henrique Dallazuanna <wwwhsd(a)gmail.com > >escreveu: > > > Tente assim: > > > > !length(unique(v)) > 1 > > > > Onde v é o seu vetor > > > > 2011/8/10 Luis Iván Ortiz Valencia <liov2067(a)gmail.com>: > > > Bom dia > > > > > > existe uma funcao no R que de um valor logico (T/F) se um conjunto de > > > valores ( como datas) sao iguais? > > > > > > atte > > > > > > IVAN > > > > > > > > > > > > _______________________________________________ > > > 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. > > > > > > > > > > > -- > > Henrique Dallazuanna > > Curitiba-Paraná-Brasil > > 25° 25' 40" S 49° 16' 22" 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. > > > > > > > ............................................... > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/6780e97b/atta… > > > > ------------------------------ > > Message: 13 > Date: Wed, 10 Aug 2011 17:00:31 -0300 > From: "Clayton Santos Delfino" <csd(a)arandanet.com.br> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Usar peso para tabela cruzada > Message-ID: > <4DF023B51C70864891A31FDC6281AA6404F11170(a)arandaad.arandanet.com.br > > > Content-Type: text/plain; charset="iso-8859-1" > > Olá pessoal, boa tarde. > > Estou usando a função CrossTable do pacote "descr" para gerar uma tabulação > sobre dois fatores: marcas x regiao, tudo certo. > > O problema é que eu preciso que os percentuais sejam dados sobre o número > de questionários enviados e não sobre o número de marcas apontadas. > > Explico: na minha pesquisa, havia a opção do respondente incluir mais de > uma marca na resposta para um item. Neste caso, eu tenho que considerar > um "voto", por assim dizer, para cada marca que ele indicou mas preciso > considerar apenas um questionário quando fizer a tabela cruzada, ou seja, > O universo tem que ser o número de questionarios e não o número de votos. > > Se algum puder me indicar alguma função ou algum caminho para isso eu > agradeço. > > Segue um fragmento dos dados para ajudar. > > Valeu. > > Clayton > > > que marca > regiao > 1 OLIVO NE > 1 ABALUX NE > 2 LUMIBRAS CO > 2 RCG CO > 2 SKYLUX CO > 3 DAUTEC GSP > 3 ECP GSP > 3 LUMIBRAS GSP > 3 DECADA GSP > 5 RESMINI RS > 6 LUMICENTER ISP > 6 INTRAL ISP > 6 ITAIM ISP > 7 ABALUX CO > 7 MAEL CO > 8 TASCHIBRA SC > 9 ABALUX ISP > 9 BLAN ISP > 11 ABALUX SC > 12 LUMIMUNDI PR > 12 ABALUX PR > 12 ECP PR > 13 ECP PR > 13 ABALUX PR > 14 INDELPA ISP > 14 INTRAL ISP > 15 ART-LUZ ISP > 15 CAROLINO ISP > 15 ITAIM ISP > 15 LUMISOFT ISP > 16 RCG PR > 16 ABALUX PR > 17 SKYLUX MG > 17 LUMILUZ MG > 17 BLUMENAU MG > 17 LUMIBRAS MG > 18 ITAIM GSP > 19 ABALUX SC > 19 TASCHIBRA SC > 20 LUMA RJ > 20 ABALUX RJ > 20 BARSOTTI RJ > 20 ITAIM RJ > > > > ------------------------------ > > Message: 14 > Date: Wed, 10 Aug 2011 17:09:08 -0300 > From: Daniel Dantas <daniel.dantas(a)hotmail.com> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Substituir NA's > Message-ID: <COL120-W184A367D82C4BA229E3003E7230(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > Boa tarde pessoal, > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os > que não forem NA manter os valores de x? > > Lógica: Se x = NA, então 0.1, caso contrário x. > > Tentei da forma abaixo mas não obtive exito: > > > x > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > y<- ifelse(x==NA,0.1,x) > > y > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] NA NA NA NA NA NA NA > [2,] NA NA NA NA NA NA NA > [3,] NA NA NA NA NA NA NA > [4,] NA NA NA NA NA NA NA > [5,] NA NA NA NA NA NA NA > [6,] NA NA NA NA NA NA NA > > Tentei também de outra forma mas também sem sucess: > y<- ifelse(x=="NA",0.1,x) > > Obrigado, > Daniel > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/5531199f/atta… > > > > ------------------------------ > > Message: 15 > Date: Wed, 10 Aug 2011 17:13:15 -0300 > From: Gustavo Carvalho <gustavo.bio(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > Message-ID: > <CAKBHXWtfGf_zecqXWDF4RTHV2z=ubYY4iryQTVKKX8Q5WUJuRQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=UTF-8 > > ifelse(is.na(x), 0.1, x) funciona? > > 2011/8/10 Daniel Dantas <daniel.dantas(a)hotmail.com>: > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e > os > > que não forem NA manter os valores de x? > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > Tentei da forma abaixo mas não obtive exito: > > > >> x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > >> y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > > > _______________________________________________ > > 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. > > > > ------------------------------ > > Message: 16 > Date: Wed, 10 Aug 2011 17:15:29 -0300 > From: Daniel Dantas <daniel.dantas(a)hotmail.com> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Substituir NA's > Message-ID: <COL120-W4362FC8215D1B945B327EFE7230(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-2" > > > Melhor impossível! heheeh > > Muito obrigado Gustavo!!!!!! > > > From: gustavo.bio(a)gmail.com > > Date: Wed, 10 Aug 2011 17:13:15 -0300 > > To: r-br(a)listas.c3sl.ufpr.br > > Subject: Re: [R-br] Substituir NA's > > > > ifelse(is.na(x), 0.1, x) funciona? > > > > 2011/8/10 Daniel Dantas <daniel.dantas(a)hotmail.com>: > > > Boa tarde pessoal, > > > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 > e os > > > que n?o forem NA manter os valores de x? > > > > > > Lógica: Se x = NA, ent?o 0.1, caso contrário x. > > > > > > Tentei da forma abaixo mas n?o obtive exito: > > > > > >> x > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > > > y<- ifelse(x==NA,0.1,x) > > >> y > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] NA NA NA NA NA NA NA > > > [2,] NA NA NA NA NA NA NA > > > [3,] NA NA NA NA NA NA NA > > > [4,] NA NA NA NA NA NA NA > > > [5,] NA NA NA NA NA NA NA > > > [6,] NA NA NA NA NA NA NA > > > > > > Tentei também de outra forma mas também sem sucess: > > > y<- ifelse(x=="NA",0.1,x) > > > > > > Obrigado, > > > Daniel > > > > > > > > > > > > _______________________________________________ > > > 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. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/6884f07b/atta… > > > > ------------------------------ > > Message: 17 > Date: Wed, 10 Aug 2011 17:15:42 -0300 (BRT) > From: Paulo Justiniano <paulojus(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > Message-ID: <alpine.DEB.2.00.1108101714050.3064(a)pataxo.est.ufpr.br> > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > > veja a seguir, troque -1 por 0.1 > > > ap <- cbind(c(1, NA, 2), 1:3, c(NA, 1, 2)) > > ap > [,1] [,2] [,3] > [1,] 1 1 NA > [2,] NA 2 1 > [3,] 2 3 2 > > ifelse(is.na(ap), -1, ap) > [,1] [,2] [,3] > [1,] 1 1 -1 > [2,] -1 2 1 > [3,] 2 3 2 > > > > > On Wed, 10 Aug 2011, Daniel Dantas wrote: > > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e > os que não forem NA manter os valores de x? > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > Tentei da forma abaixo mas não obtive exito: > > > > > x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > > > y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > > > > > ------------------------------ > > Message: 18 > Date: Wed, 10 Aug 2011 16:20:00 -0400 > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Substituir NA's > Message-ID: <D2C27C38-9D36-41B4-8986-927DC071B541(a)yahoo.com.br> > Content-Type: text/plain; charset=utf-8 > > If(var="",0.1,var) > > Edson Lira > Estatístico > Ma-Am > > Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> > escreveu: > > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e > os que não forem NA manter os valores de x? > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > Tentei da forma abaixo mas não obtive exito: > > > > > x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > > > y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > _______________________________________________ > > 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. > > > ------------------------------ > > Message: 19 > Date: Wed, 10 Aug 2011 17:25:55 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > Message-ID: > <CALKkXXqLYrunKHOAJ0+dpVBGmueAVii-MDrUSAHtiDt6KbN1Fw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Daniel, > Pode ser que resolva > x[is.na(x)]<- 0.1 > Att > > Em 10 de agosto de 2011 17:09, Daniel Dantas > <daniel.dantas(a)hotmail.com>escreveu: > > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e > os > > que não forem NA manter os valores de x? > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > Tentei da forma abaixo mas não obtive exito: > > > > > x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > > > y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/39a8ecdb/atta… > > > > ------------------------------ > > Message: 20 > Date: Wed, 10 Aug 2011 17:30:52 -0300 > From: Daniel Dantas <daniel.dantas(a)hotmail.com> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Substituir NA's > Message-ID: <COL120-W236336DEC4913359EB1A33E7230(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o > Prof. Paulo Justiniano! > > E agora me surge uma outra dúvida: > > ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a linha) > ? > > Como escrevo isso acima em linguagem R? > > Obrigado, > Daniel > > > From: edinhoestat(a)yahoo.com.br > > Date: Wed, 10 Aug 2011 16:20:00 -0400 > > To: r-br(a)listas.c3sl.ufpr.br > > Subject: Re: [R-br] Substituir NA's > > > > If(var="",0.1,var) > > > > Edson Lira > > Estatístico > > Ma-Am > > > > Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> > escreveu: > > > > > Boa tarde pessoal, > > > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 > e os que não forem NA manter os valores de x? > > > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > > > Tentei da forma abaixo mas não obtive exito: > > > > > > > x > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > > > y<- ifelse(x==NA,0.1,x) > > > > y > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] NA NA NA NA NA NA NA > > > [2,] NA NA NA NA NA NA NA > > > [3,] NA NA NA NA NA NA NA > > > [4,] NA NA NA NA NA NA NA > > > [5,] NA NA NA NA NA NA NA > > > [6,] NA NA NA NA NA NA NA > > > > > > Tentei também de outra forma mas também sem sucess: > > > y<- ifelse(x=="NA",0.1,x) > > > > > > Obrigado, > > > Daniel > > > > > > > > > _______________________________________________ > > > 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. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/8dad73c6/atta… > > > > ------------------------------ > > Message: 21 > Date: Wed, 10 Aug 2011 17:32:28 -0300 > From: eder(a)leg.ufpr.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > Message-ID: > <d80ddd8615de5b2e7273f0b59fe32ed2.squirrel(a)www.leg.ufpr.br> > Content-Type: text/plain;charset=iso-8859-1 > > Daniel, > na.omit(x) > Att > > > > Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o > > Prof. Paulo Justiniano! > > > > E agora me surge uma outra dúvida: > > > > ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a > linha) > > ? > > > > Como escrevo isso acima em linguagem R? > > > > Obrigado, > > Daniel > > > >> From: edinhoestat(a)yahoo.com.br > >> Date: Wed, 10 Aug 2011 16:20:00 -0400 > >> To: r-br(a)listas.c3sl.ufpr.br > >> Subject: Re: [R-br] Substituir NA's > >> > >> If(var="",0.1,var) > >> > >> Edson Lira > >> Estatístico > >> Ma-Am > >> > >> Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> > >> escreveu: > >> > >> > Boa tarde pessoal, > >> > > >> > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 > >> e os que não forem NA manter os valores de x? > >> > > >> > Lógica: Se x = NA, então 0.1, caso contrário x. > >> > > >> > Tentei da forma abaixo mas não obtive exito: > >> > > >> > > x > >> > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > >> > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > >> > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > >> > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > >> > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > >> > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > >> > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > >> > > >> > y<- ifelse(x==NA,0.1,x) > >> > > y > >> > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > >> > [1,] NA NA NA NA NA NA NA > >> > [2,] NA NA NA NA NA NA NA > >> > [3,] NA NA NA NA NA NA NA > >> > [4,] NA NA NA NA NA NA NA > >> > [5,] NA NA NA NA NA NA NA > >> > [6,] NA NA NA NA NA NA NA > >> > > >> > Tentei também de outra forma mas também sem sucess: > >> > y<- ifelse(x=="NA",0.1,x) > >> > > >> > Obrigado, > >> > Daniel > >> > > >> > > >> > _______________________________________________ > >> > 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. > > > > > ------------------------------ > > Message: 22 > Date: Wed, 10 Aug 2011 17:37:00 -0300 > From: Marcos Silva <marcosfs2006(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Funções para alto filtros. > Message-ID: > <CAHKdobzMmKAgfnVGWXfqN5Gd0_Rmdv368oJoNODuxrAyBTHPWQ(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Acho que talvez o pacote data.table lhe seja útil. Eu ainda não o utilizei, > mas pelo que li parece ser bem rápido. > > > http://cran.r-project.org/web/packages/data.table/vignettes/datatable-intro… > > Abs. > > Em 10 de agosto de 2011 15:38, Pedro Rafael > <pedro.rafael.marinho(a)gmail.com>escreveu: > > > Pessoa existe funções analogas ao subset? Achei muito interessante a > função > > subset para filtrar grandes massas de dados segundo algumas > características > > desejadas. O R fornece alguma outra função que tenha esse mesmo papel > (fazer > > filtros em massa) ou algum outro recurso que seja interessante para fazer > > filtros em data.frames enormes? Para mim só basta o me informarem o nome > da > > função, detalhamentos eu procuro no help. > > > > Pedro Rafael > > > > _______________________________________________ > > 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. > > > > > > -- > Marcos F. Silva > http://sites.google.com/site/marcosfs2006 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/a650d702/atta… > > > > ------------------------------ > > Message: 23 > Date: Wed, 10 Aug 2011 17:37:08 -0300 (BRT) > From: Paulo Justiniano <paulojus(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > Message-ID: <alpine.DEB.2.00.1108101736310.3064(a)pataxo.est.ufpr.br> > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > complete.cases() > retem apenas as linahs completas do objeto > > > > > > On Wed, 10 Aug 2011, Daniel Dantas wrote: > > > Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o > Prof. Paulo Justiniano! > > > > E agora me surge uma outra dúvida: > > > > ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a > linha) ? > > > > Como escrevo isso acima em linguagem R? > > > > Obrigado, > > Daniel > > > > > From: edinhoestat(a)yahoo.com.br > > > Date: Wed, 10 Aug 2011 16:20:00 -0400 > > > To: r-br(a)listas.c3sl.ufpr.br > > > Subject: Re: [R-br] Substituir NA's > > > > > > If(var="",0.1,var) > > > > > > Edson Lira > > > Estatístico > > > Ma-Am > > > > > > Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> > escreveu: > > > > > > > Boa tarde pessoal, > > > > > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor > 0.1 e os que não forem NA manter os valores de x? > > > > > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > > > > > Tentei da forma abaixo mas não obtive exito: > > > > > > > > > x > > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > > > > > y<- ifelse(x==NA,0.1,x) > > > > > y > > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > > [1,] NA NA NA NA NA NA NA > > > > [2,] NA NA NA NA NA NA NA > > > > [3,] NA NA NA NA NA NA NA > > > > [4,] NA NA NA NA NA NA NA > > > > [5,] NA NA NA NA NA NA NA > > > > [6,] NA NA NA NA NA NA NA > > > > > > > > Tentei também de outra forma mas também sem sucess: > > > > y<- ifelse(x=="NA",0.1,x) > > > > > > > > Obrigado, > > > > Daniel > > > > > > > > > > > > _______________________________________________ > > > > 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. > > > > > > ------------------------------ > > Message: 24 > Date: Wed, 10 Aug 2011 17:39:42 -0300 > From: Jakson Alves de Aquino <jalvesaq(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usar peso para tabela cruzada > Message-ID: > <CAGBu4CMqxYior8fjuT2fHqS7V5VsWB-fAjbAzfsoXP8dTu=ugQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=UTF-8 > > 2011/8/10 Clayton Santos Delfino <csd(a)arandanet.com.br>: > > Estou usando a função CrossTable do pacote "descr" para gerar uma > > tabulação sobre dois fatores: marcas x regiao, tudo certo. > > > > O problema é que eu preciso que os percentuais sejam dados sobre o > > número de questionários enviados e não sobre o número de marcas > > apontadas. > > > > Explico: na minha pesquisa, havia a opção do respondente incluir > > mais de uma marca na resposta para um item. Neste caso, eu tenho que > > considerar um "voto", por assim dizer, para cada marca que ele > > indicou mas preciso considerar apenas um questionário quando fizer a > > tabela cruzada, ou seja, O universo tem que ser o número de > > questionarios e não o número de votos. > > > > Se algum puder me indicar alguma função ou algum caminho para isso > > eu agradeço. > > > > Segue um fragmento dos dados para ajudar. > > Supondo que o data.frame se chama "b", veja se o código abaixo faz o > que você quer: > > peso <- 1 / table(b$que) > peso <- data.frame(que = as.numeric(names(peso)), peso = as.numeric(peso)) > b2 <- merge(b, peso) > > library(descr) > crosstab(b2$marca, b2$regiao, b2$peso) > > -- > Jakson > > > ------------------------------ > > Message: 25 > Date: Wed, 10 Aug 2011 17:42:56 -0300 > From: Marcos Silva <marcosfs2006(a)gmail.com> > To: R-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Material Novo no CRAN > Message-ID: > <CAHKdobx9GgTJ9koG1bmNb6E2TaFFM8RqhdJtwHAVbKfjAW=kJQ(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Só para avisar que existe um novo material sobre o R em portugues no CRAN. > > http://cran.r-project.org/doc/contrib/Provete-Estatistica_aplicada.pdf > > Valeu. > > -- > Marcos F. Silva > http://sites.google.com/site/marcosfs2006 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/cd1ded06/atta… > > > > ------------------------------ > > Message: 26 > Date: Wed, 10 Aug 2011 18:49:46 -0300 > From: Carlos Mendonça <csaeslpv(a)centroin.com.br> > To: Forum_Linguagem R <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] [Dúvida] Rodando Script R no terminal sem entrar > no prompt. > Message-ID: > <CA+sjXFrF7oM9eKWbv2d5MxrLhowwCYZpW9c0Bk0rBY508ZJ6cg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Guilherme, mais uma vez, obrigado pela ajuda. Eu fiz o > que você postou, porém, mas quando eu entro com o > > segundo comando (no seu exemplo: "c:\R\tarefasC.r" ) é aberto o Tinn-R, sem > que seja executado o script automaticamente, > > Será que teria uma outra alternativa? > > Um abraço, > > Carlos Mendonça. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/ebdcf292/atta… > > > > ------------------------------ > > Message: 27 > Date: Wed, 10 Aug 2011 18:50:48 -0300 > From: Leonard Assis <assis.leonard(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Accuracy 2012 - First Call for Papers - July > 10-13, 2012, Florianópolis, SC - Brazil > Message-ID: > <CAEG0FK_juNtjXmZ9iHYYyc1VzsYFx8vXn7mofY7jcBpQA0kwnw(a)mail.gmail.com > > > Content-Type: text/plain; charset="utf-8" > > Mais uma da série "Coisas interessantes que acontecem em Floripa depois que > eu me mudei de lá" > > > lmassis <at> yahoo <dot> com <dot> br > assis.leonard <at> gmail <dot> com > > > 2011/8/10 RODRIGO LILLA MANZIONE <rlmanzione(a)ig.com.br> > > > *ACCURACY 2012** > > International Spatial Accuracy Research Association (ISARA) > > The Tenth International Symposium on Spatial Accuracy Assessment in > Natural > > Resources and Environmental Sciences > > July 10-13, 2012 > > Florianópolis, SC - Brazil* > > > > http://2012.spatial-accuracy.org/ > > > > TOPICS: All topics concerned with spatial accuracy and uncertainty in a > > natural resources and environmental sciences context are appropriate, > > for example: > > > > * Semantic uncertainty and vagueness > > * Modelling uncertainty using geostatistics > > * Propagation of uncertainty in GIS > > * Visualizing spatial uncertainty > > * Uncertainty in Remote Sensing > > * Spatio-temporal uncertainty > > * Accuracy and uncertainty of DEMs > > * Modelling scale in environmental systems > > * Positional uncertainty > > > > PUBLICATION: We welcome the submission of both oral papers and posters. > > For both modes of presentation authors are invited to submit abstracts > > of no more than 500 words in English on original research.Those accepted > > will be presented as oral papers/posters at the symposium and published > as > > short papers in the symposium proceedings (ISBN listed). All accepted > > authors will also be > > invited to submit a full paper for consideration towards an special > > issue of an international journal. > > > > MODES OF PARTICIPATION: Oral papers and posters. Pre-conference > > workshop proposals also welcome. > > > > DEADLINES: > > > > Abstract submission: *December** 9th, ****2011* > > > > Abstract acceptance notification: *March** 1st, ****2012* > > > > Short final paper submission: *May** 11th, ****2012* > > > > Workshop proposals: *October** 18th, **2011* > > > > Workshop acceptance notification: *October** 31st, **2011* > > It is important to mention that this will be the first time the Spatial > > Accuracy Group co-organizes a conference with The Commission on Modelling > > Geographical System of the International Geographical Union (IGU-CMGS)< > http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research area of common > interest. > > > > Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) > > Mariane Dal Santo (Secretary) > > > > Contact: accuracy2012(a)cfh.ufsc.br< > http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf… > > > > **************************************** > > Download a pdf of the full call for papers at: > > http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg > > > > > > > > PS: Sorry for e-mail list crossed. > > > > -- > > Prof. Carlos Antonio Oliveira Vieira > > Departamento de Geociências ? CFH/UFSC > > Trindade - Florianópolis - SC - Brasil > > CEP 88040-900 Cx. Postal 476 > > Fone: ++55(48)3721-8593 Cel: (48)9915-3653 > > E-mails: carlos.vieira(a)ufv.br > > carlos.vieira(a)ufsc.br > > > > Visite o site: www.accuracy2012.ufsc.br > > > > > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/9e2ccd14/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 6, assunto 12 > **************************************** > -- Júlio Sílvio de Sousa Bueno Filho Departamento de Ciências Exatas Universidade Federal de Lavras Caixa Postal 3037 Lavras, MG - Brasil - 37200 000 tel: (+55) 35 3829 1369 fax: (+55) 35 3829 1371 "Natureza da gente não cabe em nenhuma certeza." João G. Rosa - Riobaldo (Grande Sertão: Veredas) "Homo sum; humani nihil a me alienum puto." Terêncio - Chremes (Heautontimoroumenos)

3 5

Verificar Matrizes
by Edmundo Caetano 11 Ago '11

11 Ago '11

Bom dia! Gostaria que me ajudassem várias questões relacionadas com matrizes: 1. Como verificar se uma matriz é ortogonal. 2. Como verificar se uma matriz é singular 3. Calcular o posto de uma matriz. ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 10 de Agosto de 2011 18:51 Assunto: Digest R-br, volume 6, assunto 12 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. Re: Obter coordenadas a partir de endereço (Pedro Emmanuel Alvarenga Americano do Brasil) 2. Re: Obter coordenadas a partir de endereço (Benilton Carvalho) 3. automatização de ajuste de modelos de regressao (Samuel Carvalho) 4. Re: automatização de ajuste de modelos de regressao (Ivan Bezerra Allaman) 5. Re: automatização de ajuste de modelos de regressao (Benilton Carvalho) 6. Re: automatização de ajuste de modelos de regressao (Henrique Dallazuanna) 7. [Dúvida] Funções para alto filtros. (Pedro Rafael) 8. apagar valores de um dataframe e colocar NA neles (Luis Iván Ortiz Valencia) 9. Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil (RODRIGO LILLA MANZIONE) 10. Re: apagar valores de um dataframe e colocar NA neles (Paulo Justiniano) 11. Re: apagar valores de um dataframe e colocar NA neles (Luis Iván Ortiz Valencia) 12. Re: funcao para valores iguais (Luis Iván Ortiz Valencia) 13. Usar peso para tabela cruzada (Clayton Santos Delfino) 14. Substituir NA's (Daniel Dantas) 15. Re: Substituir NA's (Gustavo Carvalho) 16. Re: Substituir NA's (Daniel Dantas) 17. Re: Substituir NA's (Paulo Justiniano) 18. Re: Substituir NA's (Edson Lira) 19. Re: Substituir NA's (Eder David Borges da Silva) 20. Re: Substituir NA's (Daniel Dantas) 21. Re: Substituir NA's (eder(a)leg.ufpr.br) 22. Re: [Dúvida] Funções para alto filtros. (Marcos Silva) 23. Re: Substituir NA's (Paulo Justiniano) 24. Re: Usar peso para tabela cruzada (Jakson Alves de Aquino) 25. Material Novo no CRAN (Marcos Silva) 26. Re: [Dúvida] Rodando Script R no terminal sem entrar no prompt. (Carlos Mendonça) 27. Re: Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil (Leonard Assis) ---------------------------------------------------------------------- Message: 1 Date: Wed, 10 Aug 2011 12:11:17 -0300 From: Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Obter coordenadas a partir de endereço Message-ID: <CAFfGvy+rMiS0DrNEY2GLCKU9F9Bq5b7wojGSXrDQXJ6wNr5wPg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Amigos de R, Nao trabalho mais com espacial ha muitos e muitos anos, porem fiquei curioso com essa função, será que essa função exige que o os endereços seja formatados sempre do mesmo jeito? Estou lembrando de um discussão que tive quase dez anos atras por conta de muita falta de dados em georeferenciamento para fazer mapas de dengue no rio de janeiro. QUando havia dados de endereços eles eram irregualres. Ou seja, nomes de ruas soletrados de forma incorreta, as vezes o numero da casa era separado com virgula outras vezes com ponto do nome da rua, e muitas vezes não havia separador entre o numero do apartamento e o numero da predio. Fiquei pensando se talvez o R consiga buscar o CEP do endereço no site dos correios, ou se é possivel baixar dos correios um banco de dados com CEP, e tnetar corresponder esses endereços com CEP. Nao seria mais facil? Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional Em 9 de agosto de 2011 08:49, Thiago Veloso <thi_veloso(a)yahoo.com.br>escreveu: > Luís Gustavo, > > Muito obrigado pela sua dica. Com ela consegui fazer a conversão que > necessitava. > > Agradeço também aos demais participantes que responderam a minha dúvida. > > Um abraço, > > Thiago Veloso. > > --- On *Fri, 5/8/11, Luís Gustavo <lgsilvaesilva(a)gmail.com>* wrote: > > > From: Luís Gustavo <lgsilvaesilva(a)gmail.com> > > Subject: Re: [R-br] Obter coordenadas a partir de endereço > To: r-br(a)listas.c3sl.ufpr.br > Date: Friday, 5 August, 2011, 12:04 > > > Prezado Thiago, > > Se não estou enganado você consegue as coordendas a partir do CEP, da > seguinte forma: > > > *require(XML) > > coordenadas<- function(cep) { > url_lat_lon <- paste(sprintf(" > http://maps.google.com/maps/api/geocode/xml?address=%s,", > cep),"%20Brasil&sensor=false", sep="") > lat_lon=xmlApply(xmlRoot(xmlTreeParse( > readLines(url_lat_lon)))[['result']][['geometry']][['location']], "[[", 1) > return(lat_lon) > } > > cep=36026300 > coordenadas(cep) > > Obs.: Existe posts passados abordando estes assunto, lembro que salvei esse > script em alguns desses posts, porém não > me lembro exatamente quem postou, portanto fica as minhas desculpas pela a > falta de referência. > * > Em 5 de agosto de 2011 00:19, Daniel Marcelino <dmsilva.br(a)gmail.com<http://mc/compose?to=dmsilva.br@gmail.com>> > escreveu: > > > > Olá Thiago, > > Deve haver mais opções para fazer isso, eu conheço o PBSmapping. Há > inclusive um material vendido pelo O' Reilly's com o passo-a-passo para > conseguir isso. > > Daniel > > > > 2011/8/4 Thiago Veloso <thi_veloso(a)yahoo.com.br<http://mc/compose?to=thi_veloso@yahoo.com.br> > > > >> > >> Olá pessoal, > >> > >> Tenho quase certeza que esse assunto já foi discutido aqui na lista, > porém não estou conseguindo encontrar o tópico. > >> > >> Possuo uma lista com cerca de 40 endereços em Porto Alegre. Gostaria de > obter a localização geográfica destes endereços para mostrá-la em um mapa. > Apesar de ter somente o endereço, com um pouco de esforço é possível > conseguir o CEP de cada um deles. > >> > >> Sendo assim, há algum pacote combinado R/google maps que permita > conhecer as coordenadas de um endereço a partir do seu CEP? > >> > >> Grato desde já, > >> > >> Thiago Veloso. > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br<http://mc/compose?to=R-br@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. > > > > > > > > -- > > Daniel Marcelino > > http://danielmarcelino.com > > Skype: d_marcelino > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br <http://mc/compose?to=R-br@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. > > > > -- > Luís Gustavo Silva e Silva > > > -----Inline Attachment Follows----- > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br <http://mc/compose?to=R-br@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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/04ae926b/atta…> ------------------------------ Message: 2 Date: Wed, 10 Aug 2011 16:39:25 +0100 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Obter coordenadas a partir de endereço Message-ID: <CAO-arWNLQq1SATXKNcTq7yUpbdg5b-fCNJm0tBru=rQ=2FDKkw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" talvez, inclusive, alguma interacao com o pessoal da pc2: http://ceplivre.pc2consultoria.com/ b -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/2416ff6e/atta…> ------------------------------ Message: 3 Date: Wed, 10 Aug 2011 09:23:51 -0700 (PDT) From: Samuel Carvalho <samukajm(a)yahoo.com.br> To: r-br <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] automatização de ajuste de modelos de regressao Message-ID: <1312993431.87925.YahooMailNeo(a)web161201.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Caros(as) Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e para isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar este modelo. O que consigo, é fazer este ajuste inserindo o efeito da parcela na composição do modelo linear. Abaixo segue um CMR dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15)) mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da parcela no modelo ou mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de fazer de uma maneira mais automatizada mod2 <- lm(y ~ x, subset(dados, parcela==2)) mod3<- lm(y ~ x, subset(dados, parcela==3)) Pela atenção Obrigado ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] ============================================= -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/172ab19a/atta…> ------------------------------ Message: 4 Date: Wed, 10 Aug 2011 10:01:42 -0700 (PDT) From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> To: R Brasil <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] automatização de ajuste de modelos de regressao Message-ID: <1312995702.55404.YahooMailNeo(a)web161820.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" De um modo bem menos elegante do que os nossos amigos da programação, pode-se ter: length(levels(factor(dados$parcela))) results <- list() for(i in 1:length(levels(factor(dados$parcela)))){ results[[i]] <- lm(y ~ x, subset(dados, parcela==i)) } results lapply(results,summary) lapply(results,anova) Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/1e7c4bcc/atta…> ------------------------------ Message: 5 Date: Wed, 10 Aug 2011 18:04:27 +0100 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Subject: Re: [R-br] automatização de ajuste de modelos de regressao Message-ID: <CAO-arWOvRoGoHO0yFQhA55ejqsSHZt0spv=UDKnMjXpD+OxVMQ(a)mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1 by(dados, dados[['parcela']], function(z) lm(y~x, data=z)) ------------------------------ Message: 6 Date: Wed, 10 Aug 2011 14:11:19 -0300 From: Henrique Dallazuanna <wwwhsd(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Samuel Carvalho <samukajm(a)yahoo.com.br> Subject: Re: [R-br] automatização de ajuste de modelos de regressao Message-ID: <CAPvBnPEvmGhh4cOgwi-e0XrfnHVB3kObwqTxOqz-9iMSB8dsoQ(a)mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1 Tente assim: mod <- lapply(split(dados, dados$parcela), lm, formula = x ~ y) 2011/8/10 Samuel Carvalho <samukajm(a)yahoo.com.br>: > Caros(as) > Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e para > isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar > este modelo. > O que consigo, é fazer este ajuste inserindo o efeito da parcela na > composição do modelo linear. Abaixo segue um CMR > > dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15)) > mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da parcela > no modelo > ou > mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de > fazer de uma maneira mais automatizada > mod2 <- lm(y ~ x, subset(dados, parcela==2)) > mod3<- lm(y ~ x, subset(dados, parcela==3)) > > Pela atenção > Obrigado > > ==================================== > Samuel P. C. Carvalho > Mestre em Ciências Florestais [UFLA] > Doutorando em Recursos Florestais [ESALQ/USP] > ============================================= > _______________________________________________ > 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. > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O ------------------------------ Message: 7 Date: Wed, 10 Aug 2011 15:38:05 -0300 From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] [Dúvida] Funções para alto filtros. Message-ID: <CAKwavrm+GG6vLrGpGcbN1Mvskp3ZLBp0TCv3_cUMDAeHNvQxHA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Pessoa existe funções analogas ao subset? Achei muito interessante a função subset para filtrar grandes massas de dados segundo algumas características desejadas. O R fornece alguma outra função que tenha esse mesmo papel (fazer filtros em massa) ou algum outro recurso que seja interessante para fazer filtros em data.frames enormes? Para mim só basta o me informarem o nome da função, detalhamentos eu procuro no help. Pedro Rafael -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/5a917572/atta…> ------------------------------ Message: 8 Date: Wed, 10 Aug 2011 16:21:41 -0300 From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] apagar valores de um dataframe e colocar NA neles Message-ID: <CAJ+S-MH=a86GH79wa=Sq4HnDndSVtE5tvC0Cav-kbfp4encjtw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Galera Tenho um data frame com valores, e preciso apagar esses valores e colocar NA neles.....qual a melhor forma? Ivan -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/05450618/atta…> ------------------------------ Message: 9 Date: Wed, 10 Aug 2011 16:42:56 -0300 From: RODRIGO LILLA MANZIONE <rlmanzione(a)ig.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil Message-ID: <CAE2kf_5DdKc_KR+4SC1PumXxhRT6oPRjB+YaEvDQaAYxr9034w(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" *ACCURACY 2012** International Spatial Accuracy Research Association (ISARA) The Tenth International Symposium on Spatial Accuracy Assessment in Natural Resources and Environmental Sciences July 10-13, 2012 Florianópolis, SC - Brazil* http://2012.spatial-accuracy.org/ TOPICS: All topics concerned with spatial accuracy and uncertainty in a natural resources and environmental sciences context are appropriate, for example: * Semantic uncertainty and vagueness * Modelling uncertainty using geostatistics * Propagation of uncertainty in GIS * Visualizing spatial uncertainty * Uncertainty in Remote Sensing * Spatio-temporal uncertainty * Accuracy and uncertainty of DEMs * Modelling scale in environmental systems * Positional uncertainty PUBLICATION: We welcome the submission of both oral papers and posters. For both modes of presentation authors are invited to submit abstracts of no more than 500 words in English on original research.Those accepted will be presented as oral papers/posters at the symposium and published as short papers in the symposium proceedings (ISBN listed). All accepted authors will also be invited to submit a full paper for consideration towards an special issue of an international journal. MODES OF PARTICIPATION: Oral papers and posters. Pre-conference workshop proposals also welcome. DEADLINES: Abstract submission: *December** 9th, ****2011* Abstract acceptance notification: *March** 1st, ****2012* Short final paper submission: *May** 11th, ****2012* Workshop proposals: *October** 18th, **2011* Workshop acceptance notification: *October** 31st, **2011* It is important to mention that this will be the first time the Spatial Accuracy Group co-organizes a conference with The Commission on Modelling Geographical System of the International Geographical Union (IGU-CMGS)<http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research area of common interest. Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) Mariane Dal Santo (Secretary) Contact: accuracy2012(a)cfh.ufsc.br<http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf…> **************************************** Download a pdf of the full call for papers at: http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg PS: Sorry for e-mail list crossed. -- Prof. Carlos Antonio Oliveira Vieira Departamento de Geociências ? CFH/UFSC Trindade - Florianópolis - SC - Brasil CEP 88040-900 Cx. Postal 476 Fone: ++55(48)3721-8593 Cel: (48)9915-3653 E-mails: carlos.vieira(a)ufv.br carlos.vieira(a)ufsc.br Visite o site: www.accuracy2012.ufsc.br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/03a37f1f/atta…> ------------------------------ Message: 10 Date: Wed, 10 Aug 2011 16:41:53 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] apagar valores de um dataframe e colocar NA neles Message-ID: <alpine.DEB.2.00.1108101641400.3064(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" todos os valores? df[,] <- NA On Wed, 10 Aug 2011, Luis Iván Ortiz Valencia wrote: > > Galera > > Tenho um data frame com valores, e preciso apagar esses valores e colocar NA neles.....qual a melhor forma? > > Ivan > > > > ------------------------------ Message: 11 Date: Wed, 10 Aug 2011 16:51:23 -0300 From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] apagar valores de um dataframe e colocar NA neles Message-ID: <CAJ+S-MHyQ3pixjTgc1oPou3X=Z+=GnG0FRBJXdzh0gjUhKQVpA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" sim e mantendo a estrutura do data frame , funcionou, obrigado Paulo. Em 10 de agosto de 2011 16:41, Paulo Justiniano <paulojus(a)leg.ufpr.br>escreveu: > todos os valores? > > df[,] <- NA > > > > On Wed, 10 Aug 2011, Luis Iván Ortiz Valencia wrote: > > >> Galera >> >> Tenho um data frame com valores, e preciso apagar esses valores e colocar >> NA neles.....qual a melhor forma? >> >> Ivan >> >> >> >> > _______________________________________________ > 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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/1d05403a/atta…> ------------------------------ Message: 12 Date: Wed, 10 Aug 2011 16:52:48 -0300 From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] funcao para valores iguais Message-ID: <CAJ+S-MEfpPfbzJjT5WTc-xAuZdipD32FYZBCNWVyt0-js_+m-A(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" obrigado....isso que precisava... Ivan Em 10 de agosto de 2011 11:21, Henrique Dallazuanna <wwwhsd(a)gmail.com>escreveu: > Tente assim: > > !length(unique(v)) > 1 > > Onde v é o seu vetor > > 2011/8/10 Luis Iván Ortiz Valencia <liov2067(a)gmail.com>: > > Bom dia > > > > existe uma funcao no R que de um valor logico (T/F) se um conjunto de > > valores ( como datas) sao iguais? > > > > atte > > > > IVAN > > > > > > > > _______________________________________________ > > 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. > > > > > > -- > Henrique Dallazuanna > Curitiba-Paraná-Brasil > 25° 25' 40" S 49° 16' 22" 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. > ............................................... -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/6780e97b/atta…> ------------------------------ Message: 13 Date: Wed, 10 Aug 2011 17:00:31 -0300 From: "Clayton Santos Delfino" <csd(a)arandanet.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Usar peso para tabela cruzada Message-ID: <4DF023B51C70864891A31FDC6281AA6404F11170(a)arandaad.arandanet.com.br> Content-Type: text/plain; charset="iso-8859-1" Olá pessoal, boa tarde. Estou usando a função CrossTable do pacote "descr" para gerar uma tabulação sobre dois fatores: marcas x regiao, tudo certo. O problema é que eu preciso que os percentuais sejam dados sobre o número de questionários enviados e não sobre o número de marcas apontadas. Explico: na minha pesquisa, havia a opção do respondente incluir mais de uma marca na resposta para um item. Neste caso, eu tenho que considerar um "voto", por assim dizer, para cada marca que ele indicou mas preciso considerar apenas um questionário quando fizer a tabela cruzada, ou seja, O universo tem que ser o número de questionarios e não o número de votos. Se algum puder me indicar alguma função ou algum caminho para isso eu agradeço. Segue um fragmento dos dados para ajudar. Valeu. Clayton que marca regiao 1 OLIVO NE 1 ABALUX NE 2 LUMIBRAS CO 2 RCG CO 2 SKYLUX CO 3 DAUTEC GSP 3 ECP GSP 3 LUMIBRAS GSP 3 DECADA GSP 5 RESMINI RS 6 LUMICENTER ISP 6 INTRAL ISP 6 ITAIM ISP 7 ABALUX CO 7 MAEL CO 8 TASCHIBRA SC 9 ABALUX ISP 9 BLAN ISP 11 ABALUX SC 12 LUMIMUNDI PR 12 ABALUX PR 12 ECP PR 13 ECP PR 13 ABALUX PR 14 INDELPA ISP 14 INTRAL ISP 15 ART-LUZ ISP 15 CAROLINO ISP 15 ITAIM ISP 15 LUMISOFT ISP 16 RCG PR 16 ABALUX PR 17 SKYLUX MG 17 LUMILUZ MG 17 BLUMENAU MG 17 LUMIBRAS MG 18 ITAIM GSP 19 ABALUX SC 19 TASCHIBRA SC 20 LUMA RJ 20 ABALUX RJ 20 BARSOTTI RJ 20 ITAIM RJ ------------------------------ Message: 14 Date: Wed, 10 Aug 2011 17:09:08 -0300 From: Daniel Dantas <daniel.dantas(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Substituir NA's Message-ID: <COL120-W184A367D82C4BA229E3003E7230(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Boa tarde pessoal, A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? Lógica: Se x = NA, então 0.1, caso contrário x. Tentei da forma abaixo mas não obtive exito: > x [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 [3,] NA 5.333333 1.000000 2.000000 1 0 1 [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA y<- ifelse(x==NA,0.1,x) > y [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] NA NA NA NA NA NA NA [2,] NA NA NA NA NA NA NA [3,] NA NA NA NA NA NA NA [4,] NA NA NA NA NA NA NA [5,] NA NA NA NA NA NA NA [6,] NA NA NA NA NA NA NA Tentei também de outra forma mas também sem sucess: y<- ifelse(x=="NA",0.1,x) Obrigado, Daniel -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/5531199f/atta…> ------------------------------ Message: 15 Date: Wed, 10 Aug 2011 17:13:15 -0300 From: Gustavo Carvalho <gustavo.bio(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Substituir NA's Message-ID: <CAKBHXWtfGf_zecqXWDF4RTHV2z=ubYY4iryQTVKKX8Q5WUJuRQ(a)mail.gmail.com> Content-Type: text/plain; charset=UTF-8 ifelse(is.na(x), 0.1, x) funciona? 2011/8/10 Daniel Dantas <daniel.dantas(a)hotmail.com>: > Boa tarde pessoal, > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os > que não forem NA manter os valores de x? > > Lógica: Se x = NA, então 0.1, caso contrário x. > > Tentei da forma abaixo mas não obtive exito: > >> x > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > y<- ifelse(x==NA,0.1,x) >> y > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] NA NA NA NA NA NA NA > [2,] NA NA NA NA NA NA NA > [3,] NA NA NA NA NA NA NA > [4,] NA NA NA NA NA NA NA > [5,] NA NA NA NA NA NA NA > [6,] NA NA NA NA NA NA NA > > Tentei também de outra forma mas também sem sucess: > y<- ifelse(x=="NA",0.1,x) > > Obrigado, > Daniel > > > > _______________________________________________ > 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. > ------------------------------ Message: 16 Date: Wed, 10 Aug 2011 17:15:29 -0300 From: Daniel Dantas <daniel.dantas(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Substituir NA's Message-ID: <COL120-W4362FC8215D1B945B327EFE7230(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-2" Melhor impossível! heheeh Muito obrigado Gustavo!!!!!! > From: gustavo.bio(a)gmail.com > Date: Wed, 10 Aug 2011 17:13:15 -0300 > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > > ifelse(is.na(x), 0.1, x) funciona? > > 2011/8/10 Daniel Dantas <daniel.dantas(a)hotmail.com>: > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os > > que n?o forem NA manter os valores de x? > > > > Lógica: Se x = NA, ent?o 0.1, caso contrário x. > > > > Tentei da forma abaixo mas n?o obtive exito: > > > >> x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > >> y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > > > _______________________________________________ > > 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. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/6884f07b/atta…> ------------------------------ Message: 17 Date: Wed, 10 Aug 2011 17:15:42 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Substituir NA's Message-ID: <alpine.DEB.2.00.1108101714050.3064(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" veja a seguir, troque -1 por 0.1 > ap <- cbind(c(1, NA, 2), 1:3, c(NA, 1, 2)) > ap [,1] [,2] [,3] [1,] 1 1 NA [2,] NA 2 1 [3,] 2 3 2 > ifelse(is.na(ap), -1, ap) [,1] [,2] [,3] [1,] 1 1 -1 [2,] -1 2 1 [3,] 2 3 2 On Wed, 10 Aug 2011, Daniel Dantas wrote: > Boa tarde pessoal, > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? > > Lógica: Se x = NA, então 0.1, caso contrário x. > > Tentei da forma abaixo mas não obtive exito: > > > x > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > y<- ifelse(x==NA,0.1,x) > > y > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] NA NA NA NA NA NA NA > [2,] NA NA NA NA NA NA NA > [3,] NA NA NA NA NA NA NA > [4,] NA NA NA NA NA NA NA > [5,] NA NA NA NA NA NA NA > [6,] NA NA NA NA NA NA NA > > Tentei também de outra forma mas também sem sucess: > y<- ifelse(x=="NA",0.1,x) > > Obrigado, > Daniel > > > > ------------------------------ Message: 18 Date: Wed, 10 Aug 2011 16:20:00 -0400 From: Edson Lira <edinhoestat(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Substituir NA's Message-ID: <D2C27C38-9D36-41B4-8986-927DC071B541(a)yahoo.com.br> Content-Type: text/plain; charset=utf-8 If(var="",0.1,var) Edson Lira Estatístico Ma-Am Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> escreveu: > Boa tarde pessoal, > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? > > Lógica: Se x = NA, então 0.1, caso contrário x. > > Tentei da forma abaixo mas não obtive exito: > > > x > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > y<- ifelse(x==NA,0.1,x) > > y > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] NA NA NA NA NA NA NA > [2,] NA NA NA NA NA NA NA > [3,] NA NA NA NA NA NA NA > [4,] NA NA NA NA NA NA NA > [5,] NA NA NA NA NA NA NA > [6,] NA NA NA NA NA NA NA > > Tentei também de outra forma mas também sem sucess: > y<- ifelse(x=="NA",0.1,x) > > Obrigado, > Daniel > > > _______________________________________________ > 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. ------------------------------ Message: 19 Date: Wed, 10 Aug 2011 17:25:55 -0300 From: Eder David Borges da Silva <eder(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Substituir NA's Message-ID: <CALKkXXqLYrunKHOAJ0+dpVBGmueAVii-MDrUSAHtiDt6KbN1Fw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Daniel, Pode ser que resolva x[is.na(x)]<- 0.1 Att Em 10 de agosto de 2011 17:09, Daniel Dantas <daniel.dantas(a)hotmail.com>escreveu: > Boa tarde pessoal, > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os > que não forem NA manter os valores de x? > > Lógica: Se x = NA, então 0.1, caso contrário x. > > Tentei da forma abaixo mas não obtive exito: > > > x > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > y<- ifelse(x==NA,0.1,x) > > y > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > [1,] NA NA NA NA NA NA NA > [2,] NA NA NA NA NA NA NA > [3,] NA NA NA NA NA NA NA > [4,] NA NA NA NA NA NA NA > [5,] NA NA NA NA NA NA NA > [6,] NA NA NA NA NA NA NA > > Tentei também de outra forma mas também sem sucess: > y<- ifelse(x=="NA",0.1,x) > > Obrigado, > Daniel > > > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/39a8ecdb/atta…> ------------------------------ Message: 20 Date: Wed, 10 Aug 2011 17:30:52 -0300 From: Daniel Dantas <daniel.dantas(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Substituir NA's Message-ID: <COL120-W236336DEC4913359EB1A33E7230(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o Prof. Paulo Justiniano! E agora me surge uma outra dúvida: ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a linha) ? Como escrevo isso acima em linguagem R? Obrigado, Daniel > From: edinhoestat(a)yahoo.com.br > Date: Wed, 10 Aug 2011 16:20:00 -0400 > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Substituir NA's > > If(var="",0.1,var) > > Edson Lira > Estatístico > Ma-Am > > Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> escreveu: > > > Boa tarde pessoal, > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > Tentei da forma abaixo mas não obtive exito: > > > > > x > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > y<- ifelse(x==NA,0.1,x) > > > y > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > [1,] NA NA NA NA NA NA NA > > [2,] NA NA NA NA NA NA NA > > [3,] NA NA NA NA NA NA NA > > [4,] NA NA NA NA NA NA NA > > [5,] NA NA NA NA NA NA NA > > [6,] NA NA NA NA NA NA NA > > > > Tentei também de outra forma mas também sem sucess: > > y<- ifelse(x=="NA",0.1,x) > > > > Obrigado, > > Daniel > > > > > > _______________________________________________ > > 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. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/8dad73c6/atta…> ------------------------------ Message: 21 Date: Wed, 10 Aug 2011 17:32:28 -0300 From: eder(a)leg.ufpr.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Substituir NA's Message-ID: <d80ddd8615de5b2e7273f0b59fe32ed2.squirrel(a)www.leg.ufpr.br> Content-Type: text/plain;charset=iso-8859-1 Daniel, na.omit(x) Att > > Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o > Prof. Paulo Justiniano! > > E agora me surge uma outra dúvida: > > ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a linha) > ? > > Como escrevo isso acima em linguagem R? > > Obrigado, > Daniel > >> From: edinhoestat(a)yahoo.com.br >> Date: Wed, 10 Aug 2011 16:20:00 -0400 >> To: r-br(a)listas.c3sl.ufpr.br >> Subject: Re: [R-br] Substituir NA's >> >> If(var="",0.1,var) >> >> Edson Lira >> Estatístico >> Ma-Am >> >> Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> >> escreveu: >> >> > Boa tarde pessoal, >> > >> > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 >> e os que não forem NA manter os valores de x? >> > >> > Lógica: Se x = NA, então 0.1, caso contrário x. >> > >> > Tentei da forma abaixo mas não obtive exito: >> > >> > > x >> > [,1] [,2] [,3] [,4] [,5] [,6] [,7] >> > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA >> > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 >> > [3,] NA 5.333333 1.000000 2.000000 1 0 1 >> > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 >> > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 >> > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA >> > >> > y<- ifelse(x==NA,0.1,x) >> > > y >> > [,1] [,2] [,3] [,4] [,5] [,6] [,7] >> > [1,] NA NA NA NA NA NA NA >> > [2,] NA NA NA NA NA NA NA >> > [3,] NA NA NA NA NA NA NA >> > [4,] NA NA NA NA NA NA NA >> > [5,] NA NA NA NA NA NA NA >> > [6,] NA NA NA NA NA NA NA >> > >> > Tentei também de outra forma mas também sem sucess: >> > y<- ifelse(x=="NA",0.1,x) >> > >> > Obrigado, >> > Daniel >> > >> > >> > _______________________________________________ >> > 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. ------------------------------ Message: 22 Date: Wed, 10 Aug 2011 17:37:00 -0300 From: Marcos Silva <marcosfs2006(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] [Dúvida] Funções para alto filtros. Message-ID: <CAHKdobzMmKAgfnVGWXfqN5Gd0_Rmdv368oJoNODuxrAyBTHPWQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Acho que talvez o pacote data.table lhe seja útil. Eu ainda não o utilizei, mas pelo que li parece ser bem rápido. http://cran.r-project.org/web/packages/data.table/vignettes/datatable-intro… Abs. Em 10 de agosto de 2011 15:38, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Pessoa existe funções analogas ao subset? Achei muito interessante a função > subset para filtrar grandes massas de dados segundo algumas características > desejadas. O R fornece alguma outra função que tenha esse mesmo papel (fazer > filtros em massa) ou algum outro recurso que seja interessante para fazer > filtros em data.frames enormes? Para mim só basta o me informarem o nome da > função, detalhamentos eu procuro no help. > > Pedro Rafael > > _______________________________________________ > 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. > -- Marcos F. Silva http://sites.google.com/site/marcosfs2006 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/a650d702/atta…> ------------------------------ Message: 23 Date: Wed, 10 Aug 2011 17:37:08 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Substituir NA's Message-ID: <alpine.DEB.2.00.1108101736310.3064(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" complete.cases() retem apenas as linahs completas do objeto On Wed, 10 Aug 2011, Daniel Dantas wrote: > Maravilha! Deu certo aqui pessoal. Muito Obrigado ao Gustavo, Edinho e o Prof. Paulo Justiniano! > > E agora me surge uma outra dúvida: > > ifelse(existe pelo menos um NA na linha, Exclua a linha, Mantenha a linha) ? > > Como escrevo isso acima em linguagem R? > > Obrigado, > Daniel > > > From: edinhoestat(a)yahoo.com.br > > Date: Wed, 10 Aug 2011 16:20:00 -0400 > > To: r-br(a)listas.c3sl.ufpr.br > > Subject: Re: [R-br] Substituir NA's > > > > If(var="",0.1,var) > > > > Edson Lira > > Estatístico > > Ma-Am > > > > Em 10/08/2011, às 16:09, Daniel Dantas <daniel.dantas(a)hotmail.com> escreveu: > > > > > Boa tarde pessoal, > > > > > > A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? > > > > > > Lógica: Se x = NA, então 0.1, caso contrário x. > > > > > > Tentei da forma abaixo mas não obtive exito: > > > > > > > x > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA > > > [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 > > > [3,] NA 5.333333 1.000000 2.000000 1 0 1 > > > [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 > > > [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 > > > [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA > > > > > > y<- ifelse(x==NA,0.1,x) > > > > y > > > [,1] [,2] [,3] [,4] [,5] [,6] [,7] > > > [1,] NA NA NA NA NA NA NA > > > [2,] NA NA NA NA NA NA NA > > > [3,] NA NA NA NA NA NA NA > > > [4,] NA NA NA NA NA NA NA > > > [5,] NA NA NA NA NA NA NA > > > [6,] NA NA NA NA NA NA NA > > > > > > Tentei também de outra forma mas também sem sucess: > > > y<- ifelse(x=="NA",0.1,x) > > > > > > Obrigado, > > > Daniel > > > > > > > > > _______________________________________________ > > > 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. > > ------------------------------ Message: 24 Date: Wed, 10 Aug 2011 17:39:42 -0300 From: Jakson Alves de Aquino <jalvesaq(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Usar peso para tabela cruzada Message-ID: <CAGBu4CMqxYior8fjuT2fHqS7V5VsWB-fAjbAzfsoXP8dTu=ugQ(a)mail.gmail.com> Content-Type: text/plain; charset=UTF-8 2011/8/10 Clayton Santos Delfino <csd(a)arandanet.com.br>: > Estou usando a função CrossTable do pacote "descr" para gerar uma > tabulação sobre dois fatores: marcas x regiao, tudo certo. > > O problema é que eu preciso que os percentuais sejam dados sobre o > número de questionários enviados e não sobre o número de marcas > apontadas. > > Explico: na minha pesquisa, havia a opção do respondente incluir > mais de uma marca na resposta para um item. Neste caso, eu tenho que > considerar um "voto", por assim dizer, para cada marca que ele > indicou mas preciso considerar apenas um questionário quando fizer a > tabela cruzada, ou seja, O universo tem que ser o número de > questionarios e não o número de votos. > > Se algum puder me indicar alguma função ou algum caminho para isso > eu agradeço. > > Segue um fragmento dos dados para ajudar. Supondo que o data.frame se chama "b", veja se o código abaixo faz o que você quer: peso <- 1 / table(b$que) peso <- data.frame(que = as.numeric(names(peso)), peso = as.numeric(peso)) b2 <- merge(b, peso) library(descr) crosstab(b2$marca, b2$regiao, b2$peso) -- Jakson ------------------------------ Message: 25 Date: Wed, 10 Aug 2011 17:42:56 -0300 From: Marcos Silva <marcosfs2006(a)gmail.com> To: R-br(a)listas.c3sl.ufpr.br Subject: [R-br] Material Novo no CRAN Message-ID: <CAHKdobx9GgTJ9koG1bmNb6E2TaFFM8RqhdJtwHAVbKfjAW=kJQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Só para avisar que existe um novo material sobre o R em portugues no CRAN. http://cran.r-project.org/doc/contrib/Provete-Estatistica_aplicada.pdf Valeu. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/cd1ded06/atta…> ------------------------------ Message: 26 Date: Wed, 10 Aug 2011 18:49:46 -0300 From: Carlos Mendonça <csaeslpv(a)centroin.com.br> To: Forum_Linguagem R <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] [Dúvida] Rodando Script R no terminal sem entrar no prompt. Message-ID: <CA+sjXFrF7oM9eKWbv2d5MxrLhowwCYZpW9c0Bk0rBY508ZJ6cg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Guilherme, mais uma vez, obrigado pela ajuda. Eu fiz o que você postou, porém, mas quando eu entro com o segundo comando (no seu exemplo: "c:\R\tarefasC.r" ) é aberto o Tinn-R, sem que seja executado o script automaticamente, Será que teria uma outra alternativa? Um abraço, Carlos Mendonça. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/ebdcf292/atta…> ------------------------------ Message: 27 Date: Wed, 10 Aug 2011 18:50:48 -0300 From: Leonard Assis <assis.leonard(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil Message-ID: <CAEG0FK_juNtjXmZ9iHYYyc1VzsYFx8vXn7mofY7jcBpQA0kwnw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Mais uma da série "Coisas interessantes que acontecem em Floripa depois que eu me mudei de lá" lmassis <at> yahoo <dot> com <dot> br assis.leonard <at> gmail <dot> com 2011/8/10 RODRIGO LILLA MANZIONE <rlmanzione(a)ig.com.br> > *ACCURACY 2012** > International Spatial Accuracy Research Association (ISARA) > The Tenth International Symposium on Spatial Accuracy Assessment in Natural > Resources and Environmental Sciences > July 10-13, 2012 > Florianópolis, SC - Brazil* > > http://2012.spatial-accuracy.org/ > > TOPICS: All topics concerned with spatial accuracy and uncertainty in a > natural resources and environmental sciences context are appropriate, > for example: > > * Semantic uncertainty and vagueness > * Modelling uncertainty using geostatistics > * Propagation of uncertainty in GIS > * Visualizing spatial uncertainty > * Uncertainty in Remote Sensing > * Spatio-temporal uncertainty > * Accuracy and uncertainty of DEMs > * Modelling scale in environmental systems > * Positional uncertainty > > PUBLICATION: We welcome the submission of both oral papers and posters. > For both modes of presentation authors are invited to submit abstracts > of no more than 500 words in English on original research.Those accepted > will be presented as oral papers/posters at the symposium and published as > short papers in the symposium proceedings (ISBN listed). All accepted > authors will also be > invited to submit a full paper for consideration towards an special > issue of an international journal. > > MODES OF PARTICIPATION: Oral papers and posters. Pre-conference > workshop proposals also welcome. > > DEADLINES: > > Abstract submission: *December** 9th, ****2011* > > Abstract acceptance notification: *March** 1st, ****2012* > > Short final paper submission: *May** 11th, ****2012* > > Workshop proposals: *October** 18th, **2011* > > Workshop acceptance notification: *October** 31st, **2011* > It is important to mention that this will be the first time the Spatial > Accuracy Group co-organizes a conference with The Commission on Modelling > Geographical System of the International Geographical Union (IGU-CMGS)<http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research area of common interest. > > Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) > Mariane Dal Santo (Secretary) > > Contact: accuracy2012(a)cfh.ufsc.br<http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf…> > **************************************** > Download a pdf of the full call for papers at: > http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg > > > > PS: Sorry for e-mail list crossed. > > -- > Prof. Carlos Antonio Oliveira Vieira > Departamento de Geociências ? CFH/UFSC > Trindade - Florianópolis - SC - Brasil > CEP 88040-900 Cx. Postal 476 > Fone: ++55(48)3721-8593 Cel: (48)9915-3653 > E-mails: carlos.vieira(a)ufv.br > carlos.vieira(a)ufsc.br > > Visite o site: www.accuracy2012.ufsc.br > > > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110810/9e2ccd14/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 6, assunto 12 ****************************************

2 1

Re: [R-br] Modelos não lineares....
by Thiago De paula protásio 10 Ago '11

10 Ago '11

Obrigado Walmes pelos esclarecimentos e me desculpe por enviar os dados no anexo. Sua ajuda sempre nos salva.... Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Substituir NA's
by Daniel Dantas 10 Ago '11

10 Ago '11

Boa tarde pessoal, A minha dúvida é: Como substituir os NA's de uma matriz pelo valor 0.1 e os que não forem NA manter os valores de x? Lógica: Se x = NA, então 0.1, caso contrário x. Tentei da forma abaixo mas não obtive exito: > x [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.359375 2.593750 2.453125 0.531250 NA NA NA [2,] 5.095238 6.619048 3.857143 2.119048 11 0 11 [3,] NA 5.333333 1.000000 2.000000 1 0 1 [4,] 0.800000 7.850000 3.850000 2.250000 4 1 2 [5,] 5.833333 9.000000 4.833333 3.666667 2 0 0 [6,] 2.285714 8.857143 3.428571 2.714286 NA NA NA y<- ifelse(x==NA,0.1,x) > y [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] NA NA NA NA NA NA NA [2,] NA NA NA NA NA NA NA [3,] NA NA NA NA NA NA NA [4,] NA NA NA NA NA NA NA [5,] NA NA NA NA NA NA NA [6,] NA NA NA NA NA NA NA Tentei também de outra forma mas também sem sucess: y<- ifelse(x=="NA",0.1,x) Obrigado, Daniel

8 10

Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil
by RODRIGO LILLA MANZIONE 10 Ago '11

10 Ago '11

*ACCURACY 2012** International Spatial Accuracy Research Association (ISARA) The Tenth International Symposium on Spatial Accuracy Assessment in Natural Resources and Environmental Sciences July 10-13, 2012 Florianópolis, SC - Brazil* http://2012.spatial-accuracy.org/ TOPICS: All topics concerned with spatial accuracy and uncertainty in a natural resources and environmental sciences context are appropriate, for example: * Semantic uncertainty and vagueness * Modelling uncertainty using geostatistics * Propagation of uncertainty in GIS * Visualizing spatial uncertainty * Uncertainty in Remote Sensing * Spatio-temporal uncertainty * Accuracy and uncertainty of DEMs * Modelling scale in environmental systems * Positional uncertainty PUBLICATION: We welcome the submission of both oral papers and posters. For both modes of presentation authors are invited to submit abstracts of no more than 500 words in English on original research.Those accepted will be presented as oral papers/posters at the symposium and published as short papers in the symposium proceedings (ISBN listed). All accepted authors will also be invited to submit a full paper for consideration towards an special issue of an international journal. MODES OF PARTICIPATION: Oral papers and posters. Pre-conference workshop proposals also welcome. DEADLINES: Abstract submission: *December** 9th, ****2011* Abstract acceptance notification: *March** 1st, ****2012* Short final paper submission: *May** 11th, ****2012* Workshop proposals: *October** 18th, **2011* Workshop acceptance notification: *October** 31st, **2011* It is important to mention that this will be the first time the Spatial Accuracy Group co-organizes a conference with The Commission on Modelling Geographical System of the International Geographical Union (IGU-CMGS)<http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research area of common interest. Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) Mariane Dal Santo (Secretary) Contact: accuracy2012(a)cfh.ufsc.br<http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf…> **************************************** Download a pdf of the full call for papers at: http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg PS: Sorry for e-mail list crossed. -- Prof. Carlos Antonio Oliveira Vieira Departamento de Geociências – CFH/UFSC Trindade - Florianópolis - SC - Brasil CEP 88040-900 Cx. Postal 476 Fone: ++55(48)3721-8593 Cel: (48)9915-3653 E-mails: carlos.vieira(a)ufv.br carlos.vieira(a)ufsc.br Visite o site: www.accuracy2012.ufsc.br

2 1

[Dúvida] Funções para alto filtros.
by Pedro Rafael 10 Ago '11

10 Ago '11

Pessoa existe funções analogas ao subset? Achei muito interessante a função subset para filtrar grandes massas de dados segundo algumas características desejadas. O R fornece alguma outra função que tenha esse mesmo papel (fazer filtros em massa) ou algum outro recurso que seja interessante para fazer filtros em data.frames enormes? Para mim só basta o me informarem o nome da função, detalhamentos eu procuro no help. Pedro Rafael

2 1

funcao para valores iguais
by Luis Iván Ortiz Valencia 10 Ago '11

10 Ago '11

Bom dia existe uma funcao no R que de um valor logico (T/F) se um conjunto de valores ( como datas) sao iguais? atte IVAN

4 5

Accuracy 2012 - First Call for Papers - July 10-13, 2012, Florianópolis, SC - Brazil
by RODRIGO LILLA MANZIONE 10 Ago '11

10 Ago '11

*ACCURACY 2012** International Spatial Accuracy Research Association (ISARA) The Tenth International Symposium on Spatial Accuracy Assessment in Natural Resources and Environmental Sciences July 10-13, 2012 Florianópolis, SC - Brazil* http://2012.spatial-accuracy.org/ TOPICS: All topics concerned with spatial accuracy and uncertainty in a natural resources and environmental sciences context are appropriate, for example: * Semantic uncertainty and vagueness * Modelling uncertainty using geostatistics * Propagation of uncertainty in GIS * Visualizing spatial uncertainty * Uncertainty in Remote Sensing * Spatio-temporal uncertainty * Accuracy and uncertainty of DEMs * Modelling scale in environmental systems * Positional uncertainty PUBLICATION: We welcome the submission of both oral papers and posters. For both modes of presentation authors are invited to submit abstracts of no more than 500 words in English on original research.Those accepted will be presented as oral papers/posters at the symposium and published as short papers in the symposium proceedings (ISBN listed). All accepted authors will also be invited to submit a full paper for consideration towards an special issue of an international journal. MODES OF PARTICIPATION: Oral papers and posters. Pre-conference workshop proposals also welcome. DEADLINES: Abstract submission: *December** 9th, ****2011* Abstract acceptance notification: *March** 1st, ****2012* Short final paper submission: *May** 11th, ****2012* Workshop proposals: *October** 18th, **2011* Workshop acceptance notification: *October** 31st, **2011* It is important to mention that this will be the first time the Spatial Accuracy Group co-organizes a conference with The Commission on Modelling Geographical System of the International Geographical Union (IGU-CMGS)<http://www.science.mcmaster.ca/%7Eigu-cmgs/>on the research area of common interest. Carlos A. O. Vieira and Francisco Henrique de Oliveira (co-Chairs) Mariane Dal Santo (Secretary) Contact: accuracy2012(a)cfh.ufsc.br<http://quimera.cpd.ufv.br/webmail/src/compose.php?send_to=accuracy2012%40cf…> **************************************** Download a pdf of the full call for papers at: http://www.accuracy2012.ufsc.br/files/2011/07/Accuracy2012_Flyer.jpg PS: Sorry for e-mail list crossed. -- Prof. Carlos Antonio Oliveira Vieira Departamento de Geociências – CFH/UFSC Trindade - Florianópolis - SC - Brasil CEP 88040-900 Cx. Postal 476 Fone: ++55(48)3721-8593 Cel: (48)9915-3653 E-mails: carlos.vieira(a)ufv.br carlos.vieira(a)ufsc.br Visite o site: www.accuracy2012.ufsc.br

1 0

automatização de ajuste de modelos de regressao
by Samuel Carvalho 10 Ago '11

10 Ago '11

Caros(as) Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e para isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar este modelo. O que consigo, é fazer este ajuste inserindo o efeito da parcela na composição do modelo linear. Abaixo segue um CMR dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15)) mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da parcela no modelo ou mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de fazer de uma maneira mais automatizada mod2 <- lm(y ~ x, subset(dados, parcela==2)) mod3<- lm(y ~ x, subset(dados, parcela==3)) Pela atenção Obrigado ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

Re: [R-br] automatização de ajuste de modelos de regressao
by Ivan Bezerra Allaman 10 Ago '11

10 Ago '11

De um modo bem menos elegante do que os nossos amigos da programação, pode-se ter: length(levels(factor(dados$parcela))) results <- list() for(i in 1:length(levels(factor(dados$parcela)))){ results[[i]] <- lm(y ~ x, subset(dados, parcela==i)) } results lapply(results,summary) lapply(results,anova) Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Obter coordenadas a partir de endereço
by Thiago Veloso 10 Ago '11

10 Ago '11

Olá pessoal, Tenho quase certeza que esse assunto já foi discutido aqui na lista, porém não estou conseguindo encontrar o tópico. Possuo uma lista com cerca de 40 endereços em Porto Alegre. Gostaria de obter a localização geográfica destes endereços para mostrá-la em um mapa. Apesar de ter somente o endereço, com um pouco de esforço é possível conseguir o CEP de cada um deles. Sendo assim, há algum pacote combinado R/google maps que permita conhecer as coordenadas de um endereço a partir do seu CEP? Grato desde já, Thiago Veloso.

9 11

duvida
by Luis Iván Ortiz Valencia 10 Ago '11

10 Ago '11

boa tarde O R carrega arquivos CSV colocando as colunas (variaveis) todas como mode = numeric mesmo as do tipo carater como sexo, qual a melhor forma para mudar o mode das colunas? atte IVAN ...............................................

5 6

Linux - Compilando o R a partir dos fontes
by Jose Claudio Faria 09 Ago '11

09 Ago '11

Olá pessoal, Tenho percebido um ganho de velocidade muito grande do R compilando a partir dos fontes no Linux (tenho estado usando o LMDE/Xfce rolling - Linux Mint derivado do Debian com ambiente gráfico Xfce). Creio que todos que estão acompanhando as mudanças da versão 14 tem percebido isso. Quem não está acompanhando sugiro testarem e se informarem a respeito. Abaixo uma breve lista do que venho fazendo para baixar os fontes, atualizar, compilar e instalar. R-patched: --------------- 1. svn checkout https://svn.r-project.org/R/branches/R-2-13-branch/ R-pat 2. ./tools/rsync-recommended 3. ./configure --with-tcltk 4. make 5. sudo make install R-dev: -------- 1. svn checkout https://svn.r-project.org/R/trunk/ R-dev 2. ./tools/rsync-recommended 3. ./configure --with-tcltk 4. make 5. sudo make install Atualizações (R-pat ou R-dev): ----------------------------------------- 1. svn update 2. svn commit 3. make distclean 4. ./tools/rsync-recommended 5. ./configure --with-tcltk 6. make 7. sudo make install Este post tem duplo objetivo: 1- Dar umas dicas básicas (e críticas) para quem está começando; 2- Compartilhar minhas experiências e ver o que poderia fazer de uma forma melhor via sugestões. Para quem não conheçe o svn: http://pt.wikipedia.org/wiki/Subversion Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 3