Re: [R-br] Fóruns e grupos
R Brasil no telegram Tito Conte Em 13 de fevereiro de 2017 22:31, Edson Lira via R-br < r-br@listas.c3sl.ufpr.br> escreveu: > Procure no Google > "Foruns de Estatística usando o software R" > > Boa noite. > > Edson Lira > Estatístico > Manaus-Amazonas > > > Em Sábado, 11 de Fevereiro de 2017 12:10, Roberto Guevara via R-br < > r-br@listas.c3sl.ufpr.br> escreveu: > > > Colegas usuários de R, estou em busca de fóruns na internet ou grupos de > whatsapp sobre Matemática e/ou Estatística (incluindo sobre R). Poderiam me > indicar algum? > > Obrigado > > -- > Roberto Guevara Ferreira Lima > > Biólogo - CRBio 6ª Região: 73146/06-D > Mestre em Zoologia > Laboratório de Ecologia e Zoologia de Vertebrados > Instituto de Ciências Biológicas > Universidade Federal do Pará > Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 > Currículo Lattes: http://lattes.cnpq.br/7363382159007600 > > ___ > R-br mailing list > 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@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@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.
Re: [R-br] Selecionar a extração de valores dos pixels em linhas a cada número de pixels arbitrário
Cara, não entendi direito o que você quer. mas no meu entendimento vc precisa: 1- obter o nx,ny do seu raster 2- obter o index do seu ponto um i,j e ai vc cria um vetor de 2 em 2 com seq (seq(i,ny,2)) E extrai pelo index Caso não seja possíve extrair pelo index vc pode criar um vetor de coordenadas e colocar o loop nelas. Eu nunca mexi muito com raster no R (geralmente uso Py pra isso por que a empresa usa mais py que R) mas acredito que funcione. Se tiver dificuldades com os comandos da um alô que eu te dou uma força Tito Conte Em 26 de janeiro de 2017 17:43, ASANTOS via R-br <r-br@listas.c3sl.ufpr.br> escreveu: > Prezados Membros, > > Tenho um raster hipotético r: > > require(raster); require(sp) > > ## Criando um raster > r <- raster(nc=10, nr=10) > r <- setValues(r, round(runif(ncell(r))* 255)) > > Consigo extrair os valores dos pixels de interesse dadas as > coordenadas centrais dos mesmos: > > x <- c(54,18,54,54) > y <- c(27,81,9,63) > > pontos <- SpatialPoints(cbind(x,y)) > > cells <- cellFromXY(r, pontos); cells > p.cen <- xyFromCell(r, cells); p.cen > > values <- extract(r,p.cen) > dados<-cbind(p.cen,cells,values) > > No entanto, eu gostaria de encontrar um modo para extrair o valor > dos pixels, no qual, eu selecionaria uma coordenada e a função realizaria a > extração dos pixels de linhas verticais inteiras do raster a cada dois > pixels à partir da coordenada dada, por exemplo? > > Alguém poderia me dar uma ideia de como chegar a esse resultado? > > Obrigado, > > -- > == > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > Campus Cáceres > Caixa Postal 244 > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) > e-mails:alexandresanto...@yahoo.com.br > alexandre.san...@cas.ifmt.edu.br > Lattes: http://lattes.cnpq.br/1360403201088680 > OrcID: orcid.org/-0001-8232-6722 - ResearcherID: A-5790-2016 > Researchgate: www.researchgate.net/profile/Alexandre_Santos10 > LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 > Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ > == > > ___ > R-br mailing list > 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 fornea cdigo > mnimo reproduzvel. ___ R-br mailing list 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.
Re: [R-br] Modificar um data frame em uma matriz de comparação múltipla específica
Alexandre, Eu te sugiro converter seu resultado para vetor e inserir os nomes out3 = list() for(i in seq(1,length(out2))) { out3[[i]]=as.vector(out2[[i]]) names= paste( rep(rownames(out2[[i]]),each=ncol(out2[[i]])), rep(colnames(out2[[i]]),nrow(out2[[i]])) ,sep='_') names(out3[[i]])=names } Tito Conte Em 11 de janeiro de 2017 21:28, ASANTOS via R-br <r-br@listas.c3sl.ufpr.br> escreveu: > Caros Membros, > > Eu gostaria de modificar um data frame em uma matriz de > comparação múltipla específica para tanto tenho um data frame DF com a nota > dada por dois indivíduos (Indv) para cinco diferentes produtos químicos, > sendo: > > Indv<-c(1,2) > deltametrina<-c(1,1) > fipronil<-c(5,3) > imidaclopride<-c(7,5) > sulfluramida<-c(3,7) > tiametoxam<-c(9,9) > DF<-cbind(Indv,deltametrina,fipronil,imidaclopride,sulfluram > ida,tiametoxam) > > DF > Indv deltametrina fipronil imidaclopride sulfluramida tiametoxam > [1,]115 7 3 9 > [2,]213 5 7 9 > > Depois eu montei uma matriz de comparação múltipla com a > seguinte regra, a variável (cada produto) com maior valor menos o de menor > valor e armazenei a nota de cada Indv em um list, sendo: > > df <- as.data.frame(t(DF[, -1])) > out <- lapply(df, function(x) outer(x, x, function(x, y) abs(x-y))) > out2 <- lapply(out, function(m) { > dimnames(m) <- list(rownames(df), rownames(df)) > m > }) > out2 > > $V1 > deltametrina fipronil imidaclopride sulfluramida tiametoxam > deltametrina 04 6 2 8 > fipronil 40 2 2 4 > imidaclopride62 0 4 2 > sulfluramida 22 4 0 6 > tiametoxam 84 2 6 0 > > $V2 > deltametrina fipronil imidaclopride sulfluramida tiametoxam > deltametrina 02 4 6 8 > fipronil 20 2 4 6 > imidaclopride42 0 2 4 > sulfluramida 64 2 0 2 > tiametoxam 86 4 2 0 > > No entanto, gostaria que meu output fosse: > > $V1 > - [deltametrina, fipronil, 4] > - [deltametrina, imidaclopride, 6] > - [deltametrina, sulfluramida, 2] > - [deltametrina, tiametoxam, 8] > - [fipronil, imidaclopride, 2] > - [fipronil, sulfluramida, 2] > - [fipronil, tiametoxam, 4] > - [imidaclopride, sulfluramida, 4] > - [imidaclopride, tiametoxam, 2] > - [sulfluramida, tiametoxam, 6] > > $V2 > - [deltametrina, fipronil, 2] > - [deltametrina, imidaclopride, 4] > - [deltametrina, sulfluramida, 6] > - [deltametrina, tiametoxam, 8] > - [fipronil, imidaclopride, 2] > - [fipronil, sulfluramida, 4] > - [fipronil, tiametoxam, 6] > - [imidaclopride, sulfluramida, 2] > - [imidaclopride, tiametoxam, 4] > - [sulfluramida, tiametoxam, 2] > > Alguém teria alguma sugestão para ajudar? > > Obrigado, > > Alexandre > > -- > == > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > Campus Cáceres > Caixa Postal 244 > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) > e-mails:alexandresanto...@yahoo.com.br > alexandre.san...@cas.ifmt.edu.br > Lattes: http://lattes.cnpq.br/1360403201088680 > OrcID: orcid.org/-0001-8232-6722 - ResearcherID: A-5790-2016 > Researchgate: www.researchgate.net/profile/Alexandre_Santos10 > LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 > Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ > == > > ___ > R-br mailing list > 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 fornea cdigo > mnimo reproduzvel. ___ R-br mailing list 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.
Re: [R-br] Pacotes do R
Edmar, o que você precisa fazer exatamente? É mais fácil buscar o que você precisa fazer na internet do que ficar consultando pacote por pacote. Tito Conte Em 14 de outubro de 2016 09:57, Edmar Caldas via R-br < r-br@listas.c3sl.ufpr.br> escreveu: > > Bom Dia! > > Preciso baixar vários pacotes do R. São 120 pacotes onde posso encontrar > porque ficar pesquisando um por um não dá. > > Edmar > > ___ > R-br mailing list > 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@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.
Re: [R-br] Problemas para instalar RNetCDF no Rstudio (linux)
Paulo, se o erro persistir, tente pelo conda! https://www.continuum.io/downloads Tito Conte 2016-10-11 18:19 GMT-03:00 Jônatan via R-br <r-br@listas.c3sl.ufpr.br>: > Paulo, o link abaixo deve ajudar. > > http://stackoverflow.com/questions/11319698/how-to- > install-r-packages-rnetcdf-and-ncdf-on-ubuntu > > On Tue, Oct 11, 2016 at 6:00 PM, Paulo Henrique Pimenta via R-br < > r-br@listas.c3sl.ufpr.br> wrote: > >> Olá a todos, >> >> Estou tentando instalar o RNetCDF no R do qual uso, num ubunutu 12.04 no >> RStudio, porém não consigo instalar o pacote aparecendo o seguinte erro: >> >> ### >> >> > install.packages("RNetCDF")Installing package into >> > ‘/home/notemaster/R/i686-pc-linux-gnu-library/3.0’ >> (as ‘lib’ is unspecified)--2016-10-11 17:49:52-- >> https://cran.rstudio.com/src/contrib/RNetCDF_1.8-2.tar.gzResolvendo >> cran.rstudio.com (cran.rstudio.com)... 54.192.56.107Conectando-se a >> cran.rstudio.com (cran.rstudio.com)|54.192.56.107|:443... conectado.A >> requisição HTTP foi enviada, aguardando resposta... 200 OKTamanho: 95693 >> (93K) [application/x-gzip]Salvando em: >> “/tmp/Rtmpbh1Jqr/downloaded_packages/RNetCDF_1.8-2.tar.gz” >> 0K .. .. .. .. .. 53% 68,7K 1s >> 50K .. .. .. .. ... 100% 50,2K=1,6s >> 2016-10-11 17:49:54 (58,7 KB/s) - >> “/tmp/Rtmpbh1Jqr/downloaded_packages/RNetCDF_1.8-2.tar.gz” salvo >> [95693/95693] >> * installing *source* package ‘RNetCDF’ ...** package ‘RNetCDF’ successfully >> unpacked and MD5 sums checkedchecking for nc-config... yes >> checking for gcc... gcc -std=gnu99 >> checking whether the C compiler works... yes >> checking for C compiler default output file name... a.out >> checking for suffix of executables... >> checking whether we are cross compiling... no >> checking for suffix of object files... o >> checking whether we are using the GNU C compiler... yes >> checking whether gcc -std=gnu99 accepts -g... yes >> checking for gcc -std=gnu99 option to accept ISO C89... none needed >> checking how to run the C preprocessor... gcc -std=gnu99 -E >> checking for grep that handles long lines and -e... /bin/grep >> checking for egrep... /bin/grep -E >> checking for ANSI C header files... yes >> checking for sys/types.h... yes >> checking for sys/stat.h... yes >> checking for stdlib.h... yes >> checking for string.h... yes >> checking for memory.h... yes >> checking for strings.h... yes >> checking for inttypes.h... yes >> checking for stdint.h... yes >> checking for unistd.h... yes >> checking netcdf linker flags... -L/usr/lib -lnetcdf >> checking netcdf pre-processor and compiler flags... -I/usr/include -DgFortran >> checking for library containing XML_StopParser... no >> checking for library containing utFree... no >> checking for library containing utScan... noconfigure: error: "unable to use >> udunits2 or udunits"ERROR: configuration failed for package ‘RNetCDF’ >> * removing ‘/home/notemaster/R/i686-pc-linux-gnu-library/3.0/RNetCDF’Warning >> in install.packages : >> installation of package ‘RNetCDF’ had non-zero exit status >> >> The downloaded source packages are in >> ‘/tmp/Rtmpbh1Jqr/downloaded_packages’> library(RNetCDF)Erro em >> library(RNetCDF) : there is no package called ‘RNetCDF’ >> >> >> >> >> >> Tentei usar o pacote netcdf, mas não consegui devido ao pacote ter sido >> descontinuado. Espero que alguém possa me ajudar, pois sou novo no uso de >> netcdf no R. >> >> Abraços a todos. >> >> >> Paulo Henrique de A. S. Pimenta. >> >> Graduando em Meteorologia (Bacharelado) - IAG/USP. >> Fone: +5511981318435. >> >> ___ >> R-br mailing list >> 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. >> > > > > -- > ### > ## Jônatan Dupont Tatsch > ## Professor do Departamento de Física > ## Centro de Ciências Exatas e Naturais (CCNE) > ## Universidade Federal de Santa Maria - UFSM > ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - > 97105-900 > ## Telefone: +55(55)33012083 > ## www.ufsm.br/meteorologia > ### > > ___ > R-br mailing list > 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@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.
Re: [R-br] Extração de dados climáticos
Ana, Que tal você usar o pacote netcdf? Ele é muito mais adequado para isso que você deseja fazer (além de ser mais simples). Segue o manual (é bem simples de fazer a extraçao de um ponto) https://cran.r-project.org/web/packages/RNetCDF/RNetCDF.pdf Tito Conte Em 9 de outubro de 2016 20:09, Mauro Sznelwar via R-br < r-br@listas.c3sl.ufpr.br> escreveu: > Poderia enviar o arquivo correspondente para poder rodar? > > > Prezados colegas, bom dia! > > > > Então, é o seguinte, no site http://berkeleyearth.org/ eu tenho dados de > temperatura média de toda a terra. Eu preciso extrair de acordo com os meus > pontos geográficos, a temperatura média anual para cada ponto. Só que eu > tenho os específicos que eu preciso. Comecei o script de forma que eu > pudesse extrair para um ponto específico e para todos os anos, para depois > arrumar o script para fazer para todos os meus pontos e anos de interesse. > > Segue o script utilizado! > > library(raster) > > b <- brick("D:/Documentos/Desktop/analises_cap2/dados_climaticos/ > precipt.nc", varname <- "precip") > print(b) > b > > idx <- getZ(b) > coords <- matrix(c(-55, -10), ncol = 2) > vals <- extract(b, coords, df = T) > df <- data.frame(idx, t(vals)[-1,]) > > rownames(df) <- NULL > names(df) <- c('date','value') > head(df) > > Seguindo o script no local onde deveria apresentar as variáveis > extraídas , aparece NA. Alguém já fez este tipo de extração? Peço ajuda. > > Desde já, muito obrigada. > > Ana Paula > > ___ > R-br mailing list > 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@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.
Re: [R-br] "alternativa não-paramétrica" para o teste de Dunnett
Achei isso: "I'm pretty sure Dunnett's Test is for parametric data only. You can use Mann Whitney follow up tests and just divide alpha by how many follow ups you need to run to take a conservative approach." Tito Conte 2016-10-03 13:24 GMT-03:00 Luiz Leal via R-br <r-br@listas.c3sl.ufpr.br>: > Boa tarde a todos. > Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de > Dunnett? > Desde já agradeço > Luiz > > ... > *Dunnett's test* is a multiple comparison > <https://en.wikipedia.org/wiki/Multiple_comparisons> procedure to compare > each of a number of treatments with a single control. > ... > > ___ > R-br mailing list > 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@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.
Re: [R-br] Processamento paralelo para um modelo de treinamento
Ele está rodando paralelo? Você verificou? Se unix use o comando top e aperte 1 Se windows abra o gerenciador de tarefas e verifique os núcleos de trabalho. Aproveite para ver se é a CPU do seu PC o problema. Pode ser a memória ou a escrita do disco. Uma solução de baixissimo nível, mas que funciona. é você quebrar o código. Em partes e rodar vários Rs, cada um com uma etapa, gerar uns arquivos intermediários e depois concatena tudo. Tito Conte Em 2 de outubro de 2016 08:41, Fernando Gama via R-br < r-br@listas.c3sl.ufpr.br> escreveu: > Pessoal, > > Estou tentando treinar um modelo com 192 atributos e meu objetivo é, no > pós-treinamento, identificar os atributos mais importantes. (seleção de > features). > > A questão é que estou tendo problemas para treinar o modelo porque o > processamento é extremamente lento. Pesquisei sobre o parallel e o > doParallel e coloquei no meu código mas aparentemente não obtive resultados > segue um trecho do código: > > library(caret) > library(doParallel) > > myControl <- trainControl(method = "repeatedcv", number = 10, repeats = 3, > allowParallel = TRUE) > > t<-proc.time() > > cl <- makeCluster(detectCores()) > > registerDoParallel(cl) > > model <- train(GENRE~., data=dtset_genres, method="lvq", preProcess = > "scale", trControl = myControl) > > stopCluster(cl) > > proc.time()-t > > > Alguma sugestão? > > -- > Att, > > | Fernando Gama da Mata | > | Database Specialist | Master's Degree UFPA | > > | Contacts: +55 91 99150 0365 | f.fabiogam...@gmail.com | Social > Networks: [ <https://www.facebook.com/fernando.gama.13>][ > <https://plus.google.com/+FernandoGama13>][ > <https://www.linkedin.com/in/fernandogama>] | > > > ___ > R-br mailing list > 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@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.
Re: [R-br] Título de gráficos em mais de uma linha
use strwrap ou conte caracteres e insira \n no espaço anterior (isso depende do tamanho do seu gráfico) Isso depende do que você precisa Tito Conte Em 30 de setembro de 2016 12:29, Guilherme Biz via R-br < r-br@listas.c3sl.ufpr.br> escreveu: > Boa tarde. > Estou tentando fazer um gráfico com título extenso e gostaria que fosse > dividido em duas ou mais linhas de forma automática, sem ser pelo comando > \n. > Ex: > > a=rnorm(100) > b=rnorm(100) > plot(a,b,main="Este título é muito grande para ficar só em uma linha e > gostaria de dividi-lo em duas ou mais linhas de forma automática") > > Sei que usando "\n" eu passo para a linha de baixo e no Rstudio se eu > escrever o título em duas linhas também ocorre isso, mas preciso de forma > automática, pois são muitos gráficos e estou fazendo de forma iterativa. > > Alguém tem uma solução? > > Agradeço desde já. > Att, > > Guilherme Biz > > > ___ > R-br mailing list > 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@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.
Re: [R-br] Fwd: Erro
Pedro, você esta rodando em linux ou windows? Tito Conte Em 29 de julho de 2015 09:27, Pedro Mestre pmes...@pesca.sp.gov.br escreveu: Prezad@s, boa tarde. Estou tentando atualizar o R pelo installr seguindo os passos disponíveis no link ( http://www.r-bloggers.com/a-step-by-step-screenshots-tutorial-for-upgrading-r-on-windows/), porém, está dando o seguinte erro: Erro em file(con, r) : não é possível abrir a conexão As linhas de comando que executei (inclusive com as explicações) foram essas: # installing/loading the latest installr package: install.packages(installr); library(installr) # install+load installr updateR() # updating R. Alguém sabe o que pode ser? PS. Desculpem se a dúvida é básica, mas não domino linguagem de programação. Abraços e obrigado, PS. Poderiam me confirmar o recebimento da mensagem? Obrigado de novo. -- Pedro Mestre Ferreira Alves Assistente Técnico de Pesquisa Instituto de Pesca Av. Bartolomeu de Gusmão, 192 Ponta da Praia - Santos - SP CEP: 11030-906 Fone/Fax: (13) 3261 2342 www.pesca.sp.gov.br/pmes...@pesca.sp.gov.br Sent with MailTrack https://mailtrack.io/install?source=signaturelang=enreferral=pedromestr...@gmail.comidSignature=22 -- Pedro Mestre Ferreira Alves Assistente Técnico de Pesquisa Instituto de Pesca Av. Bartolomeu de Gusmão, 192 Ponta da Praia - Santos - SP CEP: 11030-906 Fone/Fax: (13) 3261 2342 www.pesca.sp.gov.br/pmes...@pesca.sp.gov.br Sent with MailTrack https://mailtrack.io/install?source=signaturelang=enreferral=pedromestr...@gmail.comidSignature=22 ___ R-br mailing list 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@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.
Re: [R-br] Fwd: Erro
E em qual step você esta tendo este erro? Tito Conte Em 6 de agosto de 2015 10:55, Tito Conte tito.co...@gmail.com escreveu: Pedro, você esta rodando em linux ou windows? Tito Conte Em 29 de julho de 2015 09:27, Pedro Mestre pmes...@pesca.sp.gov.br escreveu: Prezad@s, boa tarde. Estou tentando atualizar o R pelo installr seguindo os passos disponíveis no link ( http://www.r-bloggers.com/a-step-by-step-screenshots-tutorial-for-upgrading-r-on-windows/), porém, está dando o seguinte erro: Erro em file(con, r) : não é possível abrir a conexão As linhas de comando que executei (inclusive com as explicações) foram essas: # installing/loading the latest installr package: install.packages(installr); library(installr) # install+load installr updateR() # updating R. Alguém sabe o que pode ser? PS. Desculpem se a dúvida é básica, mas não domino linguagem de programação. Abraços e obrigado, PS. Poderiam me confirmar o recebimento da mensagem? Obrigado de novo. -- Pedro Mestre Ferreira Alves Assistente Técnico de Pesquisa Instituto de Pesca Av. Bartolomeu de Gusmão, 192 Ponta da Praia - Santos - SP CEP: 11030-906 Fone/Fax: (13) 3261 2342 www.pesca.sp.gov.br/pmes...@pesca.sp.gov.br Sent with MailTrack https://mailtrack.io/install?source=signaturelang=enreferral=pedromestr...@gmail.comidSignature=22 -- Pedro Mestre Ferreira Alves Assistente Técnico de Pesquisa Instituto de Pesca Av. Bartolomeu de Gusmão, 192 Ponta da Praia - Santos - SP CEP: 11030-906 Fone/Fax: (13) 3261 2342 www.pesca.sp.gov.br/pmes...@pesca.sp.gov.br Sent with MailTrack https://mailtrack.io/install?source=signaturelang=enreferral=pedromestr...@gmail.comidSignature=22 ___ R-br mailing list 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@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.
Re: [R-br] MCA gráfico
Edson, da um help(par) lá vai ter boa parte das coisas que você precisa E sim, aumente os limites do eixo Tito Conte Em 26 de fevereiro de 2015 11:41, Edson Lira edinhoes...@yahoo.com.br escreveu: Bom dia caros amigos. Novamente retomo um problema que ainda não consegui resolver. Download MCA COMUNICACAO.docx http://www.datafilehost.com/d/8d884324 # nesse link tem o script o gráfico e as dúvidas [image: image] http://www.datafilehost.com/d/8d884324 Download MCA COMUNICACAO.docx http://www.datafilehost.com/d/8d884324 Download MCA COMUNICACAO.docx free. File hosted free by DataFileHost. Visualizar em www.datafileho... http://www.datafilehost.com/d/8d884324 Visualizado por Yahoo Download mj1.csv http://www.datafilehost.com/d/71e81372 # os dados [image: image] http://www.datafilehost.com/d/71e81372 Download mj1.csv http://www.datafilehost.com/d/71e81372 Download mj1.csv free. File hosted free by DataFileHost. Visualizar em www.datafileho... http://www.datafilehost.com/d/71e81372 Visualizado por Yahoo No geral, necessito melhorar a aparência e a escala do gráfico e ainda não consegui. [ ]'s. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list 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@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.
Re: [R-br] Extrair coordenadas de um raster dada uma condição
Talvez ajude valores=t(matrix(r$layer[,],10,10)) Tito Conte Em 1 de fevereiro de 2015 23:48, ASANTOS alexandresanto...@yahoo.com.br escreveu: Caros Listeiros, Fiz a seleção da coordenada dos pixels (p.cen) em um raster que continham pontos de interesse (pontos) no CRM abaixo, gostaria de saber como impor uma condição para selecionar as coordenadas dos pixels que não contém nenhum ponto. Tentei criar um novo objeto com todas as coordenadas no raster e remover as coordenadas p.cen sem sucesso, alguém poderia me ajudar? Obrigado, ### code r require(raster); require(sp) ## Criando um raster r - raster(nc=10, nr=10) r - setValues(r, round(runif(ncell(r))* 255)) ##Pontos selecionados x - c(44.25,33.06,63.22,71.35) y - c(22.00,81.90,15.18,71.03) pontos - SpatialPoints(cbind(x,y)) cells - cellFromXY(r, pontos); cells coordinates(r)[cells,] p.cen - xyFromCell(r, cells); p.cen plot(r) points(pontos, col=red) points(p.cen) values - extract(r,p.cen) ##Extrai as coordendas dados-cbind(p.cen,cells,values) dados ### /code -- == Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresanto...@yahoo.com.br alexandre.san...@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 == --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com ___ R-br mailing list 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@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.
Re: [R-br] Extração de coordenadas e tempo em múltiplas imagens
Alexandre, Fiz isso para imagens de satélite O procedimento que usei para criar o laço foi com um array de 3 dimensões, onde as duas primeiras eram o x e o y e coloquei um matriz a traz da outra como se fosse construindo um cubo então tinha um objeto com x,y e t . O registro do nome de cada figura, eu armazenei em um vetor Pra ser sincero, não entendi muito bem sua rotina pois não domino algumas das funções que vc usou mesmo assim mando um exemplo do que fiz, que espero que te ajude. # carregar os arquivos de interesse no diretorio my.dir = dir(/home/tito/tcc/tau,pattern=.nc,full.names=T) # criar a variável nula do array para não dar problema TAU =NULL # loop para ler cada um dos arquivos for(arquivo in my.dir){ # carregar o arquivo pelo nome nc = open.ncdf(arquivo) # extrair variavel de tau em x taux = get.var.ncdf(arquivo,x) # armazenar a variavel em um array segundo as dimenções e o numero da figura e de acordo com o nome do arquivo TAU = array(c(TAU,taux),dim=c(nrow(taux),ncol(tauy),which(my.dir==arquivo } Espero ter ajudado Tito Conte 2013/12/4 ASANTOS alexandresanto...@yahoo.com.br Caros Membros, Estou precisando de ajuda para criar um laço para extrair as coordenadas e o tempo (no nome do arquivo) em múltiplas figuras. No meu exemplo abaixo e usando apenas a figura 1, tenho uma regiões de pixels de interesse em branco e o nome da figura Fig1-0_00_00.png, preciso das coordenadas x e y das regiões em branco (Valor digital = 0), representado pelo objeto fig.data.pix e juntei ao tempo que é 0:00:00, pois o nome da Figura é 0_00_00. Conforme CRM abaixo, fazendo figura por figura esta OK, mas gostaria de fazer um laço para extrair essas informações das três figuras abaixo e armazená-las em um único objeto, alguém teria alguma sugestão? #--- # Rotina para extração de coordenadas e tempo em imagens #--- # Abertura das imagens require(RCurl) #Figura 1 -- get.dropbox - function(file, token) { require(RCurl, quietly=TRUE) bin - getBinaryURL(paste0(https://www.dropbox.com/s/ 8b99gn9ve1yo7is/Fig1-0_00_00.png, token, /, file), ssl.verifypeer = FALSE) con - file(file, open = wb) writeBin(bin, con) close(con) message(noquote(paste(file, read into, getwd( } # get.dropbox(Fig1-0_00_00.png, Fig1-0_00_00.png) #Figura 2 -- get.dropbox - function(file, token) { require(RCurl, quietly=TRUE) bin - getBinaryURL(paste0(https://www.dropbox.com/s/ 8b99gn9ve1yo7is/Fig2-0_20_00.png, token, /, file), ssl.verifypeer = FALSE) con - file(file, open = wb) writeBin(bin, con) close(con) message(noquote(paste(file, read into, getwd( } # get.dropbox(Fig2-0_20_00.png, Fig2-0_20_00.png) #Figura 3 - get.dropbox - function(file, token) { require(RCurl, quietly=TRUE) bin - getBinaryURL(paste0(https://www.dropbox.com/s/ bo6gbt3ca6xwxm0/Fig3-0_40_00.png, token, /, file), ssl.verifypeer = FALSE) con - file(file, open = wb) writeBin(bin, con) close(con) message(noquote(paste(file, read into, getwd( } # get.dropbox(Fig3-0_40_00.png, Fig3-0_40_00.png) # Extraido as coordenadas dos pixels para figura 1 - require(raster) r - raster(Fig1-0_00_00.png) image(r) r2-as.data.frame(r) # Cria as coordenadas -- grid0 - expand.grid(x=seq(0,length.out=r@nrows), y=seq(0,length.out=r@ncols)) # Coordenadas + valor dos pixels fig.data-cbind(grid0,r2) # Região de interesse - Número digital = 0 fig.data.pix-fig.data[fig.data[,3]==0,] # Extrai o tempo do nome do arquivo tempo2-r@data@names tempo1-gsub(Fig1., , tempo2) tempo-gsub(_, :, tempo1) # Cria o arquivo final dados-cbind(fig.data.pix[,1:2],rep(tempo)) head(dados) #--- Obrigado, -- == Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresanto...@yahoo.com.br alexandre.san...@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680
Re: [R-br] Gráfico de pontos
O que exatamente você quer na legenda? uma barra de cores degrade? Tito Conte 2013/12/4 Tiago Souza Marçal tiagosouzamar...@hotmail.com Nos envie um CMR para que possamos tentar ajudá-lo. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas # -- Date: Wed, 4 Dec 2013 08:12:29 -0200 From: alxc...@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Gráfico de pontos Ao grupo r-br, por favor, eu criei um gráfico de pontos onde o tamanho e cor de cada um é baseado no quartil ao qual ele pertence, o trecho do código está abaixo. Entretanto estou com dificuldades de criar a legenda. Alguém poderia me orientar? *attach(dadosT1)* *par(mfrow=c(1,1))* *ramp-colorRamp(c(yellow, red))* *points(dadosT1,xlab=Coordenada X (coluna),* *ylab=Coordenada Y (linha),* *pt.divide=quartiles, col=rgb(ramp(seq(0,1, length=5)), max=255))* Desde já agradeço, Alex Santos ___ R-br mailing list 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@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@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.
Re: [R-br] Arquivos hdf
Thiago, sua solução foi MUITO BOA! OBRIGADO! Tito Conte 2013/8/13 Eder Comunello ecomu...@gmail.com Caro Tito, bom dia! Apenas cogitando... Os arquivos podem ser HDF4 e não HDF5. Não sei se foi corrigido, mas no passado o pacote hdf5() não abria arquivos HDF4 diretamente e era necessário convertê-los primeiro. No endereço que segue abaixo tem inclusive o link pra um conversor: [R] hdf files https://stat.ethz.ch/pipermail/r-help/2010-January/225901.html -- Éder Comunello e comunello.e...@gmail.comcomu...@gmail.com Dourados, MS - [22 16.5'S, 54 49'W] ___ R-br mailing list 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@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.
Re: [R-br] Arquivos hdf
Benilton, o hdf5 não abre o arquivo Thiago os arquivos são do seawifs, ao que parece não tem que processar a imagem, o problema é abri-los mesmo Tito Conte ___ R-br mailing list 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] Arquivos hdf
Caros, estou processadno arquivos de imagens de satélite e muitos deles vem em formto .hdf, alguns de vocês faz ideia de como abri-los (o pacote hdf5, foi testado exaustivamente e não tive bons resultados) Tito Conte ___ R-br mailing list 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.
Re: [R-br] Coordenadas e distâncias
utilize um laço, o comando for ou while resolvem, dependendo do que você deseje fazer Tito Conte 2013/7/8 Roberto Guevara robertoguevara@gmail.com Olá pessoal do R-br. Meus dados consistem basicamente em coordenadas e o horário em que cada ponto foi obtido, como exemplificado a seguir: obs dia minxy 11400 425 475 21410 425 425 31 420 475 325 41 430 525 325 51435 525 275 61443 525 275 obs se refere a cada unidade de amostra, dia a cada dia de coleta (um total de 73), min é o horário em unidades de minutos e x e y cada elemento do ponto cartesiano. O que preciso obter desses dados é a distribuição de distâncias euclidianas entre um ponto escolhido aleatoriamente e todos os pontos seguintes de um mesmo dia. Ou seja, preciso escolher aleatoriamente um ponto por dia, calcular, para todos os dias, a euclidiana entre este ponto e todos os pontos posteriores do mesmo dia e relacionar a distância espacial com a diferença de tempo entre os pontos e o dia em que cada distância está associada. Eu já elaborei um código pra isso. Caso tenham sugestões para melhorar o código, ficarei agradecido. Porém, não consigo passar para uma próxima etapa. Como poderia fazer para repetir esse procedimento o número de vezes que fosse necessário? Por exemplo, existem milhares de combinações de pontos escolhidos aleatoriamente ao longo dos 73 dias. Eu gostaria de repetir o procedimento 1000 vezes para ter uma ideia da variação de distribuição de distâncias que posso obter em decorrência da escolha aleatória de um ponto de origem por dia. Código já elaborado: arquivo-read.table(arquivo.txt,head=T) fixo-numeric() #vetor com origens de cada dia ult-numeric() #vetor com últimos pontos de cada dia xy-list() #lista para armazenar a parte dos pontos de cada dia delimitadas por fixo e ult hor-list() #lista para armazenar a parte dos horarios de cada dia delimitadas por fixo e ult distancias-list()#distancias euclidianas tempo-list()#diferença de tempo entre os pontos dia-list()#dia de coleta for (i in 1:max(unique(arquivo$dia))){ #repetindo até o total de dias fixo[i]-sample(arquivo$obs[arquivo$dia==i],1)#sorteio ponto de origem ult[i]-max(arquivo$obs[arquivo$dia==i])#armazenando último ponto de cada dia xy[[i]]-arquivo[fixo[i]:ult[i],4:5]#armazenando parte dos pontos de cada dia hor[[i]]-arquivo[fixo[i]:ult[i],6]#armazenando parte dos horários de cada dia distancias[[i]]-dist(xy[[i]])[1:nrow(xy[[i]])-1]#somente as distâncias em relação ao primeiro ponto interessam tempo[[i]]-dist(horarios[[i]])[1:length(horarios[[i]])-1]#somente diferenças ao primeiro ponto dia[[i]]-rep(i,length(arquivo$dia[fixo[i]:ult[i]])-1)$cada dia de coleta } distancias-unlist(distancias) tempo-unlist(tempo) dia-unlist(dia) dtd-data.frame(distancias,tempo,dia) Este data.frame seria como uma unidade de amostra das milhares de distribuições que são possíveis de obter. Como otimizar para que qualquer número desejado não seja dispendioso de obter? Obrigado pela atenção -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600 ___ R-br mailing list 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@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.
Re: [R-br] #plotar nome mês/ano no eixo x
De uma olhada no pacote zoo Tito Conte On Mon, Jun 24, 2013 at 3:05 PM, walmes . walmeszevi...@gmail.com wrote: Uma forma de formatar é usar a função axis.POSIX() veja, #-- # gráficos de séries no tempo met - read.table(http://www.leg.ufpr.br/~walmes/data/metereologia.txt;, header=TRUE, sep=\t, stringsAsFactors=FALSE) str(met) names(met) - tolower(names(met)) # passa os nomes para minúsculas #-- # conversão com a função as.POSIX* met$data - as.POSIXct(met$data, format=%Y-%m-%d) plot(temp~data, data=met, type=l, xaxt=n) axis.POSIXct(1, at=seq(as.POSIXct(2009-10-01), as.POSIXct(2011-10-01), by=month), format=%B, col=2) plot(temp~data, data=met, type=l, xaxt=n) axis.POSIXct(1, at=seq(as.POSIXct(2009-10-01), as.POSIXct(2011-10-01), by=month), format=%m/%y, col=2) #-- # duas séries em um mesmo gráfico par(mar=c(5.1,4.1,4.1,4.1)) plot(temp~data, data=met, type=l, xaxt=n) par(new=TRUE) plot(urel~data, data=met, type=l, xaxt=n, yaxt=n, ylab=NA, col=2) axis(4, col=2) mtext(urel, side=4, line=3) axis.POSIXct(1, at=seq(as.POSIXct(2009-10-01), as.POSIXct(2011-10-01), by=month), format=%m/%y) À 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: wal...@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 == ___ R-br mailing list 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@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.
Re: [R-br] arquivo netcdf
Heloíse verifique os seus pacotes do windows para netcdf Tito Conte 2013/6/13 Thiago V. dos Santos thi_vel...@yahoo.com.br Heloíse, O pacote ncdf é antigo e deve ser descontinuado em breve. Instale o ncdf4 no lugar. Os binários para windows não existem no CRAN, mas estão disponíveis no site do autor: http://cirrus.ucsd.edu/~pierce/ncdf/. Esse pacote deve suportar o protocolo opendap. Caso ele não suporte, uma solução temporária é baixar o arquivo para testar localmente (download.file('url_to_file'http://dl.dropboxusercontent.com/u/27700634/transition.zip',%C2%A0'transition.zip' )), e em seguida escrever para os autores dos pacotes (incluindo o Robert Hijmans do raster) perguntando se eles têm algum plano de incluir o opendap em um futuro próximo. Alternativamente, existe uma opção mais complicada de compilar o pacote gdal com suporte a opendap. Porém, não sei se isso pode ser feito em Windows. Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm Phone: (612) 323 9898 -- *From:* Heloíse Pavanato helopavan...@gmail.com *To:* r-br@listas.c3sl.ufpr.br *Sent:* Wednesday, June 12, 2013 9:45 PM *Subject:* [R-br] arquivo netcdf Olá, Tenho tentado abrir um arquivo netcdf mas sem êxito. Tentei com o pacote ncdf para abrir um arquivo através de um link OPeNDAP dataset - open.ncdf(' http://dap.marinexplore.org/dap/u/911215431c07/02841d66-d38e-11e2-8d83-22000a1c5967 ') e através do RNetCDF dataset - open.nc(' http://dap.marinexplore.org/dap/u/911215431c07/02841d66-d38e-11e2-8d83-22000a1c5967 ') e sempre dá o Error in open.ncdf trying to open file. Estou trabalhando em ambiente Windows. Isso seria um problema? Alguém tem alguma dica para me dar? Grata, Heloise ___ R-br mailing list 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@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@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] média entre pontos de uma mesma matrix
Bom dia, Preciso calcular a média de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,]369 [2,]47 10 [3,]58 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,]22 NA [2,]222 [3,]22 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte ___ R-br mailing list 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.
Re: [R-br] média entre pontos de uma mesma matrix
Benilton, como criar um array 3d (nunca fiz isso)? Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com Crie um array 3d e use rowMeans com o argumento na.rm=TRUE On 12 Jun 2013 11:14, Tito Conte tito.co...@gmail.com wrote: Bom dia, Preciso calcular a média de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,]369 [2,]47 10 [3,]58 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,]22 NA [2,]222 [3,]22 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte ___ R-br mailing list 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@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@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.
Re: [R-br] Como ordenar os elementos das linhas de uma matriz de dados?
apply(suamatrix,1,sort) Tito Conte 2013/6/12 Gilenio Borges Fernandes gile...@ufba.br Prezados(as) Como ordenar os elementos das linhas de uma matriz de dados? Exemplo, 10 linhas e 7 colunas: C1 C2 C3 C4 C5 C6 C7 L1 28 33 36 69 35 45 62 L2 21 28 33 36 16 22 59 L3 23 37 40 61 22 26 30 L4 20 61 74 78 15 22 62 L5 13 51 61 637 41 45 L6 13 20 21 424 41 62 L7 13 21 28 4368 12 L8 32 42 43 44 45 48 50 L95 20 69 741 50 52 L10 20 36 51 78 12 49 79 Grato Gilenio -- Gilenio Borges Fernandes Professor Associado III Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860 ___ R-br mailing list 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@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.
Re: [R-br] média entre pontos de uma mesma matrix
ok Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com o comando eh array(), use o argumento dim= apropriadamente... qdo houver um exemplo reproduzivel, nao me incomodarei em dar as duas linhas de comando q vc procura. Em 12 de junho de 2013 11:36, Tito Conte tito.co...@gmail.com escreveu: Benilton, como criar um array 3d (nunca fiz isso)? Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com Crie um array 3d e use rowMeans com o argumento na.rm=TRUE On 12 Jun 2013 11:14, Tito Conte tito.co...@gmail.com wrote: Bom dia, Preciso calcular a média de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,]369 [2,]47 10 [3,]58 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,]22 NA [2,]222 [3,]22 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte ___ R-br mailing list 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@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@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@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@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.
Re: [R-br] média entre pontos de uma mesma matrix
usei um a=cbind(matrix1,matrix2,matrix3) dim(a)=c(3,3,3) rowMeans(a,dims=2) Muito obrigado pessoal Tito Conte 2013/6/12 andre...@bol.com.br Mais simples do que isso, impossível, rs Benilton resumiu o que tinha feito em uma única linha de comando. Então Benilton, o atributo dims=2 na função rowMeans, foi pegar cada elemento de cada matrix em igual posição e fazer a média? Pelo resultado eu vi que sim, agora, é dificil de enxergar isso num plano tridimensional. No caso para dim=1 ou 3, o que se teria como resultado? *Att.* *André* -- Em 12/06/2013 14:55, *Benilton Carvalho beniltoncarva...@gmail.com *escreveu: mais simples que isso, dado seu array 3d 'M', M rowMeans(M, dims=2, na.rm=T) (bem mais simples qdo ha' um exemplo q pode ser reproduzido) Em 12 de junho de 2013 14:51, andre...@bol.com.br escreveu: Olá Tito, acredito que seja essa solução a qual procuras: # Matriz Tridimensional: M A nM nE S1 S2 for (i in 1:(nM*nE)) { A[i] } Result Result Espero ter ajudado colega. Att. André BVS Em 12/06/2013 14:31, Tito Conte tito.co...@gmail.com escreveu: ok Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com o comando eh array(), use o argumento dim= apropriadamente... qdo houver um exemplo reproduzivel, nao me incomodarei em dar as duas linhas de comando q vc procura. Em 12 de junho de 2013 11:36, Tito Conte tito.co...@gmail.com escreveu: Benilton, como criar um array 3d (nunca fiz isso)? Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com Crie um array 3d e use rowMeans com o argumento na.rm=TRUE On 12 Jun 2013 11:14, Tito Conte tito.co...@gmail.com wrote: Bom dia, Preciso calcular a m édia de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,] 3 6 9 [2,] 4 7 10 gt; [3,] 5 8 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,] 2 2 NA [2,] 2 2 2 [3,] 2 2 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de gt ; cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte ___ R-br mailing list R-br@listas.c3sl.u fpr.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@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. g t; ___ R-br mailing list 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@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@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ív el. ___ R-br mailing list 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@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@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.
Re: [R-br] média entre pontos de uma mesma matrix
Isso pq estava dentro de um loop que cria cada matriz Tito Conte 2013/6/12 Tito Conte tito.co...@gmail.com usei um a=cbind(matrix1,matrix2,matrix3) dim(a)=c(3,3,3) rowMeans(a,dims=2) Muito obrigado pessoal Tito Conte 2013/6/12 andre...@bol.com.br Mais simples do que isso, impossível, rs Benilton resumiu o que tinha feito em uma única linha de comando. Então Benilton, o atributo dims=2 na função rowMeans, foi pegar cada elemento de cada matrix em igual posição e fazer a média? Pelo resultado eu vi que sim, agora, é dificil de enxergar isso num plano tridimensional. No caso para dim=1 ou 3, o que se teria como resultado? *Att.* *André* -- Em 12/06/2013 14:55, *Benilton Carvalho beniltoncarva...@gmail.com *escreveu: mais simples que isso, dado seu array 3d 'M', M rowMeans(M, dims=2, na.rm=T) (bem mais simples qdo ha' um exemplo q pode ser reproduzido) Em 12 de junho de 2013 14:51, andre...@bol.com.br escreveu: Olá Tito, acredito que seja essa solução a qual procuras: # Matriz Tridimensional: M A nM nE S1 S2 for (i in 1:(nM*nE)) { A[i] } Result Result Espero ter ajudado colega. Att. André BVS Em 12/06/2013 14:31, Tito Conte tito.co...@gmail.com escreveu: ok Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com o comando eh array(), use o argumento dim= apropriadamente... qdo houver um exemplo reproduzivel, nao me incomodarei em dar as duas linhas de comando q vc procura. Em 12 de junho de 2013 11:36, Tito Conte tito.co...@gmail.com escreveu: Benilton, como criar um array 3d (nunca fiz isso)? Tito Conte 2013/6/12 Benilton Carvalho beniltoncarva...@gmail.com Crie um array 3d e use rowMeans com o argumento na.rm=TRUE On 12 Jun 2013 11:14, Tito Conte tito.co...@gmail.com wrote: Bom dia, Preciso calcular a m édia de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,] 3 6 9 [2,] 4 7 10 gt; [3,] 5 8 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,] 2 2 NA [2,] 2 2 2 [3,] 2 2 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de gt ; cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte ___ R-br mailing list R-br@listas.c3sl.u fpr.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@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. g t; ___ R-br mailing list 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@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@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ív el. ___ R-br mailing list 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@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@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.
Re: [R-br] Coordenadas UTM para Lat Lon
Obrigado a todos! Tito Conte 2013/6/4 Luciano F. Sgarbi luciano.f.sga...@gmail.com Veja se isso funciona: UTM.to.Lat_Long-function(coord){ b=6356752.3142 a=6378137 e=0.0818191909 eisq=0.0067394968 k0=0.9996 ei=(1-(1-e*e)^(1/2))/(1+(1-e*e)^(1/2)) C1=3*ei/2-27*ei^3/32 C2=21*ei^2/16-55*ei^4/32 C3=151*ei^3/96 C4=1097*ei^4/512 polo=coord[,1] zone=coord[,2] lon=coord[,3] lat=coord[,4] Corrected.Northing=ifelse(polo==N,lat,1000-lat) East.Prime=50-lon Arc.Length=lat/k0 mu=Arc.Length/(a*(1-e^2/4-3*e^4/64-5*e^6/256)) phl=mu+C1*sin(2*mu)+C2*sin(4*mu)+C3*sin(6*mu)+C4*sin(8*mu) C=eisq*cos(phl)^2 T=tan(phl)^2 N=a/(1-(e*sin(phl))^2)^(1/2) R=a*(1-e*e)/(1-(e*sin(phl))^2)^(3/2) D=East.Prime/(N*k0) F1=N*tan(phl)/R F2=D*D/2 F3=(5+3*T+10*C-4*C*C-9*eisq)*D^4/24 F4=(61+90*T+298*C+45*T*T-252*eisq-3*C*C)*D^6/720 LF1=D LF2=(1+2*T+C)*D^3/6 LF3=(5-2*C+28*T-3*C^2+8*eisq+24*T^2)*D^5/120 delta.lon=(LF1-LF2+LF3)/cos(phl) zone.cm=6*zone-183 raw.lat=180*(phl-F1*(F2+F3+F4))/pi LAT=ifelse(polo==N,raw.lat,-1*raw.lat) LON=zone.cm-delta.lon*180/pi lon.deg=as.integer(LON) lon.min=abs(as.integer(60*(LON-lon.deg))) lon.sec=3600*(abs(LON)-abs(lon.deg)-lon.min/60) lat.deg=as.integer(LAT) lat.min=abs(as.integer(60*(LAT-lat.deg))) lat.sec=3600*(abs(LAT)-abs(lat.deg)-lat.min/60) df=data.frame(coord,LON,LAT,lon.deg,lon.min,lon.sec,lat.deg,lat.min,lat.sec) return(df) } coord=data.frame(polo=sample(c(N,S),100,r=T),zone=sample(0:23,100,r=T),lon=sample(0:9,100,r=T),lat=sample (0:9,100,r=T)) UTM.to.Lat_Long(coord) Em 4 de junho de 2013 09:21, Elias Teixeira Krainski eliaskrain...@yahoo.com.br escreveu: Vc precisa informar a projeção no formato Proj4 na função project() On 03/06/13 19:08, Tito Conte wrote: Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte __**_ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -- *Luciano F. Sgarbi* * Mestre em Ecologia e Evolução Dep. Ecologia, Instituto de Ciências Biológicas I* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil* ___ R-br mailing list 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@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.
Re: [R-br] teste de Jonckheere
google it teste Jonckheere r project Tito Conte 2013/6/4 Cristina d...@ufscar.br ** Olá pessoal! Gostaria de saber se alguém já usou, no R, o teste de Jonckheere. tentei usar o teste mas consigo encontrar qual pacote está instalado. Obrigada Teresa Cristina ___ R-br mailing list 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@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.
Re: [R-br] MObilidade urbana
Edson, não entendo de mobilidade urbana, então vou perguntar se o que você esta dizendo é o que estou entendendo voce tem placas de carro no ponto X, estes carros podem ir pra Y1,Y2,Y3,Y4. Neste caso você quer saber quais carros de X foram para cada um dos Ys certo? se eu estiver certo o comando which pode ajudar Tito Conte 2013/5/31 Edson Lira edinhoes...@yahoo.com.br Caros amigos, estou tentando estabalecer uma amostra para o seguinte problema: Tenho uma origem e quatro destinos. Na origem consigo coletar 1326 placas de carros em 60 minutos, São 4 destinos e consigo identificar 620 placas desta origem, como estabelecer uma amostra nos 4 destinos para identificar o origem? [ ]'s. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list 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@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.
Re: [R-br] Plotagem 3D
da uma olhada na função cloud do lattice Tito Conte 2013/6/3 Tiago Souza Marçal tiagosouzamar...@hotmail.com Bom dia pessoal, estou querendo plotar um gráfico tridimensional com os escores canônicos que estão abaixo: dados-structure(list(V1 = c(-18.7175, -20.1915, -21.1557, -17.5965, -18.142, -13.857, -18.0086, -14.0436, -12.571), V2 = c(38.3004, 36.4955, 34.4964, 39.9303, 40.4819, 38.3763, 34.9869, 35.0723, 35.9022), V3 = c(-20.9237, -20.12, -17.8586, -17.0183, -17.9961, -20.8457, -19.5551, -17.5703, -18.7393)), .Names = c(V1, V2, V3), class = data.frame, row.names = c(NA, -9L)) require(lattice) zlab - VC2 xlab - VC1 ylab - VC3 wireframe(V2~V1*V3, scales=list(arrows=FALSE), zlab=list(zlab, rot=90), xlab=list(xlab, rot=24), ylab=list(ylab, rot=-37),xlim=c(-12,-20),ylim=c(-18,-22),type=c(on,bottom)) No entanto, não estou conseguindo fazer com que os pontos apareçam no gráfico e gostaria de eliminar as 3 arestas frontais. Agradeço desde já pela ajuda. att. Tiago # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas # ___ R-br mailing list 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@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] Coordenadas UTM para Lat Lon
Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte ___ R-br mailing list 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.
Re: [R-br] DE - PARA
Não entendi o que você quer edmar Tito Conte 2013/5/28 Edmar Caldas edmar_cal...@yahoo.com.br Bom Dia! Alguém poderia me dizer se dá pra fazer um de - para. Ex: Tenho uma tabela: Cod_item Produto JanFev mAR ABRIL maiojunhojlhoagostoset out nov dez 001 AAA 1015 108 900000 000 002 BBB 00 0 0 0 54 5 4 444 agora quero fazer que o cod_item 002 que tem valor diferente de ZERO paraa para a serie do produto AAA fica assim então. Cod_item Produto Jan Fev mAR ABRIL maio junho jlho aosto set out nov dez 001AAA 10 15 10 8 9 5 4 5 4 4 4 4 Edmar ___ R-br mailing list 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@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.
Re: [R-br] Detecçção de pontos externos
Então na verdade o problema com estes dados foram levantados por terceiros e repassados. O problema é que estou questionando os resultados, e quera eu mesmo interpola-los, mas ao que parece esses já são os interpolados. Pra esse caso solicitei os brutos. Estou aguardando. Agradeço a sua disposição e discussão, me ajudaram a criticar um pouco o que tenho recebido desse pessoal. Tito Conte 2013/5/27 Hélio Gallo Rocha heliogalloro...@gmail.com Tito, Andei olhando, parece que o comando over extrai dados como uma seleção, Talvez tenha me expresso mal, o que tava pensando era detectar os pontos mais externos da área, para fazer a borda na krigagem. Como estão os estudos da krigagem? Consegui ver seus dados usando o pacote RgoogleMaps, é em são luis do maranhão? mas pra te dizer a verdade, não tem o que krigar ali, já tem informação demais ali. Hélio Em 25 de maio de 2013 11:54, TitoConte [via R-br] ml-node+s2285057n4659465...@n4.nabble.com escreveu: Funcao over do pacote sp ( desculpe sofri correção automática) On Saturday, May 25, 2013, Tito Conte wrote: Overdose do pacote sp On Saturday, May 25, 2013, Hélio Gallo Rocha wrote: Bom dia a todos da lista, Pegando o exemplo do Tito onde seus dados possuem 25000 coordenadas, teria uma função onde se possa separar apenas os pontos externos para fazer a borda? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho -- Tito Conte -- Tito Conte ___ R-br mailing list [hidden email] http://user/SendEmail.jtp?type=nodenode=4659465i=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-Detecccao-de-pontos-externos-tp4659463p4659465.html To unsubscribe from R-br, click herehttp://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_codenode=3357982code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4 . NAMLhttp://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewerid=instant_html%21nabble%3Aemail.namlbase=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespacebreadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho ___ R-br mailing list 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@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.
Re: [R-br] Detecçção de pontos externos
Overdose do pacote sp On Saturday, May 25, 2013, Hélio Gallo Rocha wrote: Bom dia a todos da lista, Pegando o exemplo do Tito onde seus dados possuem 25000 coordenadas, teria uma função onde se possa separar apenas os pontos externos para fazer a borda? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho -- Tito Conte ___ R-br mailing list 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.
Re: [R-br] Detecçção de pontos externos
Funcao over do pacote sp ( desculpe sofri correção automática) On Saturday, May 25, 2013, Tito Conte wrote: Overdose do pacote sp On Saturday, May 25, 2013, Hélio Gallo Rocha wrote: Bom dia a todos da lista, Pegando o exemplo do Tito onde seus dados possuem 25000 coordenadas, teria uma função onde se possa separar apenas os pontos externos para fazer a borda? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho -- Tito Conte -- Tito Conte ___ R-br mailing list 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.
Re: [R-br] Cluster com var categorica e continua
Princípios de Estatística Em Ecologia - Gotelli, Nicholas J./ Ellison é um bom começo Tito Conte 2013/5/23 nani nani.des...@bol.com.br Conhece alguma apostila/ livro um pouco mais didático sobre o assunto? Obrigada! N. -- Em 22/05/2013 17:43, *Tito Conte tito.co...@gmail.com * escreveu: pacote vegan de uma olhada nos métodos Tito Conte 2013/5/22 nani nani.des...@bol.com.brhttp://../../../undefined/compose?to=nani.des...@bol.com.br Olá, Alguém tem algum exemplo de clusterização utilizando variáveis categóricas conjuntamente com variáveis contínuas? Ou alguma sugestão de material? Obr, N. ___ R-br mailing list R-br@listas.c3sl.ufpr.brhttp://../../../undefined/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@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@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.
Re: [R-br] Dúvida sobre Krigagem
Os dados pertencem a levantamentos batimétricos SObre o SO eu uso o Unix então esse não é um problema Tito Conte 2013/5/21 Benilton Carvalho beniltoncarva...@gmail.com Note que nao adianta ter 3 ou 8 ou XYZ 2 gigas de RAM, se o sistema operacional for 32 bits Em 21 de maio de 2013 12:19, Hélio Gallo Rocha heliogalloro...@gmail.com escreveu: Ok Elias, testei numa máquina com 3 gb, vou testar no not que tem 8 Abraço Hélio Em 21 de maio de 2013 10:48, Elias Krainski [via R-br] ml-node+s2285057n4659371...@n4.nabble.com escreveu: On 05/20/2013 09:31 PM, Hélio Gallo Rocha wrote: Tito, krigagem sempre me interessa, pretendo trabalhar com soft livre e agricultura de precisão, Olhei seus dados, tem tanta informação que parecem ser provenientes de krigagem, o R não conseguiu processar os dados, parece que o R tem algum limite de memória, não sei ao certo, talvez o Elias saiba... O R não possui limite de memória. Num bom sistema operacional você pode usar toda a memória RAM disponível no seu computador. O problema é usar um algoritmo em que há necessidade em se trabalhar em memória e se seus dados tiverem tamanho em memória maior que a memória do seu computador. Sobre krigagem para grande número de pontos, o problema é a necessidade de inverter uma matriz n por n. Você pode evitar isso se usar low rank kriging or the SPDE approach. Att. Elias ___ 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-Duvida-sobre-Krigagem-tp4659319p4659371.html To unsubscribe from R-br, click here. NAML -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho ___ R-br mailing list 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@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@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.
Re: [R-br] Dúvida sobre Krigagem
Tito Conte 2013/5/22 Benilton Carvalho beniltoncarva...@gmail.com Ainda assim, se seu *nix for 32 bits, vc podera encontrar problemas no enderecamento de memoria. Em 22 de maio de 2013 14:01, Tito Conte tito.co...@gmail.com escreveu: Os dados pertencem a levantamentos batimétricos SObre o SO eu uso o Unix então esse não é um problema Tito Conte 2013/5/21 Benilton Carvalho beniltoncarva...@gmail.com Note que nao adianta ter 3 ou 8 ou XYZ 2 gigas de RAM, se o sistema operacional for 32 bits Em 21 de maio de 2013 12:19, Hélio Gallo Rocha heliogalloro...@gmail.com escreveu: Ok Elias, testei numa máquina com 3 gb, vou testar no not que tem 8 Abraço Hélio Em 21 de maio de 2013 10:48, Elias Krainski [via R-br] ml-node+s2285057n4659371...@n4.nabble.com escreveu: On 05/20/2013 09:31 PM, Hélio Gallo Rocha wrote: Tito, krigagem sempre me interessa, pretendo trabalhar com soft livre e agricultura de precisão, Olhei seus dados, tem tanta informação que parecem ser provenientes de krigagem, o R não conseguiu processar os dados, parece que o R tem algum limite de memória, não sei ao certo, talvez o Elias saiba... O R não possui limite de memória. Num bom sistema operacional você pode usar toda a memória RAM disponível no seu computador. O problema é usar um algoritmo em que há necessidade em se trabalhar em memória e se seus dados tiverem tamanho em memória maior que a memória do seu computador. Sobre krigagem para grande número de pontos, o problema é a necessidade de inverter uma matriz n por n. Você pode evitar isso se usar low rank kriging or the SPDE approach. Att. Elias ___ 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-Duvida-sobre-Krigagem-tp4659319p4659371.html To unsubscribe from R-br, click here. NAML -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho ___ R-br mailing list 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@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@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@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@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.
Re: [R-br] Dúvida sobre Krigagem
Hélio, Estes dados eu recebi do pessoal que fez o levantamento pra mim, Será que eles mesmo já submeteram os dados brutos a krigagem pela cara dos dados. Se for este o caso isso justificaria por que o plot com o image para teste está tão estranho Tito Conte 2013/5/22 Tito Conte tito.co...@gmail.com Tito Conte 2013/5/22 Benilton Carvalho beniltoncarva...@gmail.com Ainda assim, se seu *nix for 32 bits, vc podera encontrar problemas no enderecamento de memoria. Em 22 de maio de 2013 14:01, Tito Conte tito.co...@gmail.com escreveu: Os dados pertencem a levantamentos batimétricos SObre o SO eu uso o Unix então esse não é um problema Tito Conte 2013/5/21 Benilton Carvalho beniltoncarva...@gmail.com Note que nao adianta ter 3 ou 8 ou XYZ 2 gigas de RAM, se o sistema operacional for 32 bits Em 21 de maio de 2013 12:19, Hélio Gallo Rocha heliogalloro...@gmail.com escreveu: Ok Elias, testei numa máquina com 3 gb, vou testar no not que tem 8 Abraço Hélio Em 21 de maio de 2013 10:48, Elias Krainski [via R-br] ml-node+s2285057n4659371...@n4.nabble.com escreveu: On 05/20/2013 09:31 PM, Hélio Gallo Rocha wrote: Tito, krigagem sempre me interessa, pretendo trabalhar com soft livre e agricultura de precisão, Olhei seus dados, tem tanta informação que parecem ser provenientes de krigagem, o R não conseguiu processar os dados, parece que o R tem algum limite de memória, não sei ao certo, talvez o Elias saiba... O R não possui limite de memória. Num bom sistema operacional você pode usar toda a memória RAM disponível no seu computador. O problema é usar um algoritmo em que há necessidade em se trabalhar em memória e se seus dados tiverem tamanho em memória maior que a memória do seu computador. Sobre krigagem para grande número de pontos, o problema é a necessidade de inverter uma matriz n por n. Você pode evitar isso se usar low rank kriging or the SPDE approach. Att. Elias ___ 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-Duvida-sobre-Krigagem-tp4659319p4659371.html To unsubscribe from R-br, click here. NAML -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho ___ R-br mailing list 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@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@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@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@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.
Re: [R-br] Tratamento de imagens no R
Gilbert, tanto é possível que estudantes de ecologia se utilizam de imagens de ossos de ouvido de peixe, para estimar a idade através do R Tito Conte 2013/5/21 Alejandro C. Frery acfr...@gmail.com Bom dia Gilbert e demais redistas, Recentemente Springer publicou o seguinte livro @BOOK{IntroImageProcessingR, title = {Introduction to Image Processing with R: Learning by Examples}, publisher = {Springer}, year = {2013}, author = {A. C. Frery and T. Perciano}, quality = {1}, url = { http://www.springer.com/computer/image+processing/book/978-1-4471-4949-1}} que, embora não trate especificamente de imagens satelitais, trata do uso de R em processamento de imagens. Alejandro 2013/4/6 Gilbert Queiroz gilbert_quei...@yahoo.com.br: Pessoal, Semana passada, um amigo me disse que vium livro sobre tratamento de imagens de satelite no R. É possível isto??? Outra pergunta: sobre aqruivos raster - é possível utilizá-los no R e fazer um tratamento vetorial (vetorização)?. Abs. Gilbert. ___ R-br mailing list 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. -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008 ___ R-br mailing list 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@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.
Re: [R-br] Tratamento de imagens no R
QUanto a imagens de satélite eu estou fazendo isso no meu TCC, dependendo do que você precisar posso te ajudar com isso Tito Conte 2013/5/22 Tito Conte tito.co...@gmail.com Gilbert, tanto é possível que estudantes de ecologia se utilizam de imagens de ossos de ouvido de peixe, para estimar a idade através do R Tito Conte 2013/5/21 Alejandro C. Frery acfr...@gmail.com Bom dia Gilbert e demais redistas, Recentemente Springer publicou o seguinte livro @BOOK{IntroImageProcessingR, title = {Introduction to Image Processing with R: Learning by Examples}, publisher = {Springer}, year = {2013}, author = {A. C. Frery and T. Perciano}, quality = {1}, url = { http://www.springer.com/computer/image+processing/book/978-1-4471-4949-1 }} que, embora não trate especificamente de imagens satelitais, trata do uso de R em processamento de imagens. Alejandro 2013/4/6 Gilbert Queiroz gilbert_quei...@yahoo.com.br: Pessoal, Semana passada, um amigo me disse que vium livro sobre tratamento de imagens de satelite no R. É possível isto??? Outra pergunta: sobre aqruivos raster - é possível utilizá-los no R e fazer um tratamento vetorial (vetorização)?. Abs. Gilbert. ___ R-br mailing list 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. -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008 ___ R-br mailing list 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@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.
Re: [R-br] concatenar dados de muitas estações
? Tito Conte 2013/5/17 Helena Turon helena.tu...@gmail.com prec é a quantidade de precipitação diária, que deveria ser um vetor numérico, você acha que é lista? Como posso resolver isto? Obrigada, Helena ___ R-br mailing list 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@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.
Re: [R-br] Dúvida sobre Krigagem
Alias tem como avaliar a qualidade destes dados? Tito Conte 2013/5/22 Tito Conte tito.co...@gmail.com Hélio, Estes dados eu recebi do pessoal que fez o levantamento pra mim, Será que eles mesmo já submeteram os dados brutos a krigagem pela cara dos dados. Se for este o caso isso justificaria por que o plot com o image para teste está tão estranho Tito Conte 2013/5/22 Tito Conte tito.co...@gmail.com Tito Conte 2013/5/22 Benilton Carvalho beniltoncarva...@gmail.com Ainda assim, se seu *nix for 32 bits, vc podera encontrar problemas no enderecamento de memoria. Em 22 de maio de 2013 14:01, Tito Conte tito.co...@gmail.com escreveu: Os dados pertencem a levantamentos batimétricos SObre o SO eu uso o Unix então esse não é um problema Tito Conte 2013/5/21 Benilton Carvalho beniltoncarva...@gmail.com Note que nao adianta ter 3 ou 8 ou XYZ 2 gigas de RAM, se o sistema operacional for 32 bits Em 21 de maio de 2013 12:19, Hélio Gallo Rocha heliogalloro...@gmail.com escreveu: Ok Elias, testei numa máquina com 3 gb, vou testar no not que tem 8 Abraço Hélio Em 21 de maio de 2013 10:48, Elias Krainski [via R-br] ml-node+s2285057n4659371...@n4.nabble.com escreveu: On 05/20/2013 09:31 PM, Hélio Gallo Rocha wrote: Tito, krigagem sempre me interessa, pretendo trabalhar com soft livre e agricultura de precisão, Olhei seus dados, tem tanta informação que parecem ser provenientes de krigagem, o R não conseguiu processar os dados, parece que o R tem algum limite de memória, não sei ao certo, talvez o Elias saiba... O R não possui limite de memória. Num bom sistema operacional você pode usar toda a memória RAM disponível no seu computador. O problema é usar um algoritmo em que há necessidade em se trabalhar em memória e se seus dados tiverem tamanho em memória maior que a memória do seu computador. Sobre krigagem para grande número de pontos, o problema é a necessidade de inverter uma matriz n por n. Você pode evitar isso se usar low rank kriging or the SPDE approach. Att. Elias ___ 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-Duvida-sobre-Krigagem-tp4659319p4659371.html To unsubscribe from R-br, click here. NAML -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho ___ R-br mailing list 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@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@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@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@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.
Re: [R-br] Dúvida sobre Krigagem
Não, obrigado pelo documento Tito Conte On Fri, May 17, 2013 at 5:01 PM, Elias Krainski eliaskrain...@yahoo.com.brwrote: Ja viu os exemplos em http://www.leg.ufpr.br/geoR/geoRdoc/tutorialsBR.html ? On 05/17/2013 04:32 PM, Tito Conte wrote: Sim essa parte eu sei, mas não sei como utilizar a função de krigagem, em ambos os casos a minha dificuldade é q o R faça a interpolação, coisa que não estou conseguindo Tito Conte 2013/5/17 Elias Krainski eliaskrain...@yahoo.com.br Basicamente, vc precisa de um modelo geoestatístico para isso. Mas há métodos de interpolação que não usam geoestatística/krigagem. ___ R-br mailing list 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 listr...@listas.c3sl.ufpr.brhttps://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@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@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.
Re: [R-br] teste de Jonckheere
http://lmgtfy.com/?q=r+project+Jonckheere Tito Conte 2013/5/20 Cristina d...@ufscar.br ** Olá pessoal! Gostaria de saber se alguém já usou no R o teste de Jonckheere. tentei usar o teste mas não sei qual pacote está instalado. Obrigada Teresa Cristina ___ R-br mailing list 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@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.
Re: [R-br] Dúvida sobre Krigagem
Interpolar por krigagem os dados do arquivo em anexo txt Tito Conte 2013/5/16 Elias Krainski eliaskrain...@yahoo.com.br krigar o arquivo em anexo Como seria krigar um arquivo? :) Sobre krigagem, veja exemplos em http://www.leg.ufpr.br/geoR/**geoRdoc/tutorialsBR.htmlhttp://www.leg.ufpr.br/geoR/geoRdoc/tutorialsBR.html Basicamente, vc precisa de um modelo geoestatístico para isso. __**_ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. ___ R-br mailing list 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.
Re: [R-br] concatenar dados de muitas estações.
o que seria prec? e você criou uma matrix e esta adicionando como lista? Tito Conte 2013/5/17 Helena Turon helena.tu...@gmail.com Desculpe-me, o problema é: prect - matrix(nrow= 15000,ncol=i) prect[[,t]] - prec Erro em prect[[, t]] - prec : mais elementos fornecidos do que podem ser substituídos. Eu penso que é por conta do número de observações ser diferente nas estações (por conta dos dados faltantes), será? Obrigada, Helena ___ R-br mailing list 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@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.
Re: [R-br] construção de um objecto para calcular os outliers
Ana, você em primeiro lugar você esta colocando uma lista dentro de uma matrix se você quer contabilizar o número de outliers veja o tamanho do vetor que você obteve dado as condições para tal faça as alterações em vermelho porano-apply(dados[-c(19,23,3 6,37),-1],2,function(x)c(summary(x),0)[1:7]) desviopadrao-apply(dados[-c(19,23,36,37),-1],2,sd,na.rm=TRUE) variancia-apply(dados[-c(19,23,36,37),-1],2,var,na.rm=TRUE) r-apply(dados[-c(19,23,36,37),-1],2,range,na.rm=TRUE) diffr-diff(r) iqr-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x,na.rm=TRUE)[4]-quantile(x,na.rm=TRUE)[2]}) besq-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x, na.rm=TRUE)[2]})-1.5*iqr bdt-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x, na.rm=TRUE)[4]})+1.5*iqr outliers-rep(0,ncol(dados)-1) for(i in names(dados[-c(19,23,36,37),-1])){ tmp-dados[-c(19,23,36,37),i] outliers[i]-length(tmp[tmp bdt[i] | dados[-c(19,23,36,37),i] besq[i]]) } row.names(diffr) - Diferen?a do range # Como o resultado do comando diff ? uma matriz, voc? precisa colocar o nome da linha antes de juntar todos os resultados rownames(porano)-c(Mínimo,Primeiro Quartil,Mediana,Média,Terceiro Quartil,Máximo,NAs) todos_resultados_dados-rbind(porano,Desvio padr?o=desviopadrao, Vari?ncia=variancia, diffr, iqr,besq,bdt,outliers) todos_resultados_dados Tito Conte 2013/5/17 alanaro...@sapo.pt Boa tarde estou a criar objectos para poder fazer o calculo dos outliers mas está a dar tudo NULL, eu queria saber porque.. dados lidos: 2008 2009 2010 2011 2012 149 54 57 54 56 253 54 55 53 59 352 51 52 53 52 457 59 59 60 63 563 60 63 64 62 658 57 59 58 63 estatísticas: #estatísticas descritivas porano-apply(dados[-c(19,23,**36,37),-1],2,function(x)c(** summary(x),0)[1:7]) desviopadrao-apply(dados[-c(**19,23,36,37),-1],2,sd,na.rm=**TRUE) variancia-apply(dados[-c(19,**23,36,37),-1],2,var,na.rm=**TRUE) r-apply(dados[-c(19,23,36,37)**,-1],2,range,na.rm=TRUE) diffr-diff(r) iqr-apply(dados[-c(19,23,36,**37),-1],2,function(x){** quantile(x,na.rm=TRUE)[4]-**quantile(x,na.rm=TRUE)[2]}) besq-apply(dados[-c(19,23,36,**37),-1],2,function(x){**quantile(x, na.rm=TRUE)[2]})-1.5*iqr bdt-apply(dados[-c(19,23,36,**37),-1],2,function(x){**quantile(x, na.rm=TRUE)[4]})+1.5*iqr outliers-list() for(i in names(dados[-c(19,23,36,37),-**1])){ tmp-dados[-c(19,23,36,37),i] outliers[[i]]-tmp[tmp bdt[i] | dados[-c(19,23,36,37),i] besq[i]] } row.names(diffr) - Diferen?a do range # Como o resultado do comando diff ? uma matriz, voc? precisa colocar o nome da linha antes de juntar todos os resultados rownames(porano)-c(Mínimo,**Primeiro Quartil,Mediana,Média,**Terceiro Quartil,Máximo,NAs) todos_resultados_dados-rbind(**porano,Desvio padr?o=desviopadrao, Vari?ncia=variancia, diffr, iqr,besq,bdt,outliers) todos_resultados_dados porano resultado: 2008 2009 2010 2011 2012 Mínimo 4852 5655.3358 Primeiro Quartil 4.708149 4.198439 4.554801 5.255589 4.6888 Mediana22.16667 17.62689 20.74621 27.62121 21.98485 Média 1919 1922 18 Terceiro Quartil 6 57 67 Máximo 4344.5 41.5 44 41.5 NAs6764.5 69.5 68 69.5 Desvio padr?o Numeric,0 44 Numeric,0 43 Numeric,0 Vari?ncia NULL NULL NULL NULL NULL Diferen?a do range NULL NULL NULL NULL NULL iqrNULL NULL NULL NULL NULL besq NULL NULL NULL NULL NULL bdtNULL NULL NULL NULL NULL outliers NULL NULL NULL NULL NULL porquê NULL? Cumprimentos Ana Rocha __**_ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. ___ R-br mailing list 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] Krigagem
Bom dia, Alguém poderia me ajudar com krigagem? Estou usando o fields, e tenho uma matriz com 57703, com x, y, z pontos, quando uso o comando Krig ele da vários erros. O comando é K = Krig(xy,z) xy são as posições e z é a altura Tito Conte ___ R-br mailing list 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.
Re: [R-br] Legenda no eixo X entre dois gráficos
De um ?mtext Lá vai ter tudo o que você precisa! Tito Conte On Thu, Jan 17, 2013 at 12:37 PM, Humberto Hazin hgha...@hotmail.comwrote: Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo ( https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Specieslci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5Sailfish 0.1257 0.25 0.498 6White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci oruci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5Sailfish 0.16000 0.20 0.4000 6White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a-read.csv(tableMetaGW.csv, header = TRUE, sep = ,, quote=\, dec=.,fill =TRUE) tableMat1b-read.csv(tableMetaGB.csv, header = TRUE, sep = ,, quote=\, dec=.,fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis - c(length(tableMat1a$or):1) plot(tableMat1a$or,y.axis, type = p, axes = F, xlab = , ylab = , pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = r, main = Green vs White,yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = n) par(mar=c(5,5, 3, 5)) y.axis - c(length(tableMat1b$or):1) plot(tableMat1b$or,y.axis, type = p, axes = F, xlab = , ylab = , pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = r, main = Green vs Blue,yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = n) dev.off() Enviado do Email do Windows ___ R-br mailing list 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@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.
Re: [R-br] recorte
você pode colocar NA nos pontos que estão interligados suponha o seguinte: duas áreas limitadas por quadrados a=c(1,1,3,3,1,4,4,5,5,4) b=c(2,4,4,2,2,6,7,7,6,6) f=matrix(c(a,b),ncol=2) plot(f,type='l') # caso vc queira remover a linha basta adicionar um NA a=c(a[1:5],NA,a[6:10]) b=c(b[1:5],NA,b[6:10]) f=data.frame(a,b) plot(f,type='l') Seria isso? ___ R-br mailing list 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.
Re: [R-br] recorte
Caso você quira fazer um recorte no mapa, faça uma máscara Criar uma matriz de zero e um (onde zero é o que você quer omitir e 1 os seus dados), através de uma equação da forma geométrica que você deseja multiplique essa máscara pelos seus dados, e plote novamente segue um exemplo abaixo de função que remove as bordas de um circulo mask - function (dados){ # tamanho da coluna Y=c(1:dim(dados)[1]) # calcular o raio r=floor(length(Y)/2) # centro do círculo cp = round(length(Y)/2,0) # matriz com zeros TT=matrix(0,nrow=dim(dados)[1],ncol=dim(dados)[2]) # fazer uma máscara de um e zero for(X in c(1:dim(dados)[2])){ # calcula o raio de todos os Y e X com dados R=sqrt((X-cp)^2+(Y-cp)^2) #o maior valor do raio do circulo zeros=which(R=r) # remove os pontos de fora do circulo TT[X,-zeros]=1 } return(TT) } aplicando aos meus dados fica dados.recortados-dados*mask(dados) espero ter ajudado Tito Conte 2013/1/18 Tito Conte tito.co...@gmail.com você pode colocar NA nos pontos que estão interligados suponha o seguinte: duas áreas limitadas por quadrados a=c(1,1,3,3,1,4,4,5,5,4) b=c(2,4,4,2,2,6,7,7,6,6) f=matrix(c(a,b),ncol=2) plot(f,type='l') # caso vc queira remover a linha basta adicionar um NA a=c(a[1:5],NA,a[6:10]) b=c(b[1:5],NA,b[6:10]) f=data.frame(a,b) plot(f,type='l') Seria isso? ___ R-br mailing list 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.
Re: [R-br] ajuda loop
Heloíse, experimente fazer assim: lat - read.table(file=lat.txt) lon - read.table(file=lon.txt) library(R.matlab) # para ler .mat sx - readMat('SX.mat',maxLength=NULL, fixNames=TRUE,verbose=F,sparseMatrixClass='SparseM') # limites de longitude lonSX - sx$SX[3,1,1] lonmin - min(as.numeric(lonSX$Lon[1:130])) lonmax - max(as.numeric(lonSX$Lon[1:130])) # limites de latitude latSX - sx$SX[4,1,1] latmin - min(as.numeric(latSX$Lat[1:130])) latmax - max(as.numeric(latSX$Lat[1:130])) indice=which(latlatmax) #indices das latitudes menores do que a máxima indice=which(lat[indice]latmin) #indices das latitudes maiores do que a minima indice=which(lon[indice]lonmax) #indices das longitudes menores do que a máxima indice=which(lon[indice]lonmin) #indices das longitudes maiores do que a minima no final vc terá os indices das coordenadas que encerram a área de interesse Atenciosamente, Tito Conte Instituto Oceanográfico Universidade de São Paulo LabPesq - Laboratório Oceanografia Pesqueira ___ R-br mailing list 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.
Re: [R-br] regressão polinomial
Paulo mas o polinomio corresponde a uma borda por isso dentro fora o trabalho que estou realizando corresponde a saber se as partículas jogadas ao acaso em uma determinada área ultrapassam ou não as bordas. Os dados que passei para vocês corresponde a uma dessas bordas, Porém os dados não são responsáveis pela construção de nenhuma dessas bordas ___ R-br mailing list 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.
Re: [R-br] regressão polinomial
Paulo e Walmes, encontrei a seguinte forma para resolver o problema, gostaria que dessem uma olhada #dados do limite da borda direitra y=c(-26.,-25.9862,-25.9343,-25.8822,-25.8433,-25.8054,-25.7948,-25.7872,-25.7668,-25.7284,-25.7015,-25.6282,-25.4612,-25.3016,-25.2564,-25.2412,-25.2412,-25.2232,-25.0869,-25.,-25.,-24.9856,-24.9397,-24.8976,-24.8533,-24.6587,-24.6373,-24.5740,-24.5406,-24.4379,-24.3934,-24.3628,-24.3277,-24.3042,-24.2972,-24.2973,-24.3020,-24.3068,-24.3068,-24.2956,-24.2920) x=c(-45.8529,-45.8302,-45.7575,-45.6680,-45.6107,-45.5050,-45.4574,-45.4058,-45.3538,-45.2913,-45.2634,-45.2264,-45.1667,-45.0719,-45.0213,-45.,-45.,-44.9748,-44.8726,-44.8405,-44.8405,-44.8351,-44.8259,-44.7972,-44.7764,-44.6867,-44.6674,-44.5912,-44.5587,-44.4927,-44.4515,-44.4115,-44.3314,-44.2664,-44.2171,-44.1729,-44.1288,-44.0820,-44.0509,-44.0118,-44.) # partícula 1 amostra=c(x=-45.3,y=-25.111) # inserir dados das amostras no limite x[length(x)+1]=amostra['x'] y[length(y)+1]=amostra['y'] #criar polinomio m=lm(y~poly(x,3)) xpred - data.frame(x = seq(min(x), max(x), len=250)) #plotar dados plot(y,x) lines(xpred$x, predict(m, newdata=xpred)) if(abs(resid(m[length(m)]))0.1 or resid(m[length(m)]))0 ) i='dentro' else i='fora' ___ R-br mailing list 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.
Re: [R-br] regressão polinomial
Paulo agora outra dúvida, tenho dez mil pontos, como comparar os que estão a direita e a esquerda dessa curva (como se fosse um dentro ou fora do sistema), já que não consigo extrair o coeficiente de da função poly? Necessito dos dados numericamente pois estes pontos serão submetidos a novos tratamentos posteriormente ___ R-br mailing list 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.
Re: [R-br] regressão polinomial
Walmes vou criar um exemplo para facilitar o entendimento o que eu preciso é comparar pontos amostrados com os resultados de um modelo de interpolação polinomial por exemplo #minhas amostras y=seq(1:9) x=seq(3:11) #pontos do modelo xm=seq(2,11) ym=seq(0:8) #tirando os valores do modelo da reta r=lm(ym~poly(xm,1)) #extraindo coeficientes a=r$coefficients[1] b=r$coefficients[0] validacao=vector() #comparando cada ponto obtido pelo modelo ajustado com a amostra for(i in c(1:length(y)){ # calculando valor do modelo m=a*x[i]+b # comparando modelo com amostra ifelse(my[i],a='fora',a='dentro') validacao[i]=a } é algo deste estilo porém com um polinomio de grau elevado, Quando ploto os resultados eles estão bonitos, mas quando calculo ele como o método abaixo os valores extrapolam muito da realidade e não entendo o porque, veja o exemplo aqui y=c(-26.,-25.9862,-25.**9343,-25.8822,-25.8433,-25.** 8054,-25.7948,-25.7872,-25.**7668,-25.7284,-25.7015,-25.** 6282,-25.4612,-25.3016,-25. 2564,-25.2412,-25.2412,-25.**2232,-25.0869,-25.,-25.** ,-24.9856,-24.9397,-24.**8976,-24.8533,-24.6587,-24.** 6373,-24.5740,-24.5406,-24 .4379,-24.3934,-24.3628,-24.**3277,-24.3042,-24.2972,-24.** 2973,-24.3020,-24.3068,-24.**3068,-24.2956,-24.2920) x=c(-45.8529,-45.8302,-45.**7575,-45.6680,-45.6107,-45.** 5050,-45.4574,-45.4058,-45.**3538,-45.2913,-45.2634,-45.** 2264,-45.1667,-45.0719,-45. 0213,-45.,-45.,-44.**9748,-44.8726,-44.8405,-44.** 8405,-44.8351,-44.8259,-44.**7972,-44.7764,-44.6867,-44.** 6674,-44.5912,-44.5587,-44 .4927,-44.4515,-44.4115,-44.**3314,-44.2664,-44.2171,-44.** 1729,-44.1288,-44.0820,-44.**0509,-44.0118,-44.) r=lm(y~poly(x,5)) b0=r$coefficient[1] b1=r$coefficient[2] b2=r$coefficient[3] b3=r$coefficient[4] b4=r$coefficient[5] b5=r$coefficient[6] modelo=b0+b1*x+b2*x^2+b3*x^3+b4*x^4+b5*x^5 modelo veja que os resultados estão muito fora. mas deveriam ser os mesmos, como proceder? ___ R-br mailing list 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.
Re: [R-br] Geoprocessamento
Geologia ou geografia? ___ R-br mailing list 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] regressão polinomial
Pessoal tenho o seginte conjunto de dados: f=c(-26.,-25.9862,-25.9343,-25.8822,-25.8433,-25.8054,-25.7948,-25.7872,-25.7668,-25.7284,-25.7015,-25.6282,-25.4612,-25.3016,-25.2564,-25.2412,-25.2412,-25.2232,-25.0869,-25.,-25.,-24.9856,-24.9397,-24.8976,-24.8533,-24.6587,-24.6373,-24.5740,-24.5406,-24.4379,-24.3934,-24.3628,-24.3277,-24.3042,-24.2972,-24.2973,-24.3020,-24.3068,-24.3068,-24.2956,-24.2920) w=c(-45.8529,-45.8302,-45.7575,-45.6680,-45.6107,-45.5050,-45.4574,-45.4058,-45.3538,-45.2913,-45.2634,-45.2264,-45.1667,-45.0719,-45.0213,-45.,-45.,-44.9748,-44.8726,-44.8405,-44.8405,-44.8351,-44.8259,-44.7972,-44.7764,-44.6867,-44.6674,-44.5912,-44.5587,-44.4927,-44.4515,-44.4115,-44.3314,-44.2664,-44.2171,-44.1729,-44.1288,-44.0820,-44.0509,-44.0118,-44.) se vocês plotarem eles, obteram uma curva, gostaria de obter a função desta curva ou uma regressão polinomial de grau elevado dela já tentei usar a função lm e não deu um bom resultado, quando chego em graus elevados e calculo utilizando os coeficientes obtidos os valores dão 10 vezes maiores (não sei bem porque) mas quando ploto ela aparece bonitinha no lugar que deveria aprecer em outro caso tentei na raça mesmo por multiplicação matricial mas também deu errado pois os valores da multiplicação da inversa são muito pequenos alguém pode me dar uma sugestão ou uma luz do que eu estou fazendo de errado? Tito Conte ___ R-br mailing list 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.
Re: [R-br] regressão polinomial
Segue o código # pontos y=c(-26.,-25.9862,-25.**9343,-25.8822,-25.8433,-25.** 8054,-25.7948,-25.7872,-25.**7668,-25.7284,-25.7015,-25.** 6282,-25.4612,-25.3016,-25. 2564,-25.2412,-25.2412,-25.**2232,-25.0869,-25.,-25.** ,-24.9856,-24.9397,-24.**8976,-24.8533,-24.6587,-24.** 6373,-24.5740,-24.5406,-24 .4379,-24.3934,-24.3628,-24.**3277,-24.3042,-24.2972,-24.** 2973,-24.3020,-24.3068,-24.**3068,-24.2956,-24.2920) x=c(-45.8529,-45.8302,-45.**7575,-45.6680,-45.6107,-45.** 5050,-45.4574,-45.4058,-45.**3538,-45.2913,-45.2634,-45.** 2264,-45.1667,-45.0719,-45. 0213,-45.,-45.,-44.**9748,-44.8726,-44.8405,-44.** 8405,-44.8351,-44.8259,-44.**7972,-44.7764,-44.6867,-44.** 6674,-44.5912,-44.5587,-44 .4927,-44.4515,-44.4115,-44.**3314,-44.2664,-44.2171,-44.** 1729,-44.1288,-44.0820,-44.**0509,-44.0118,-44.) # regressão polinomial m=lm(y~x+I(x^2)+I(x^3)+I(x^4)+I(x^5)) # criar polinomio de ordem 5 # extraindo os coeficientes b0=coefficient(m)[1] b1=coefficient(m)[2] b2=coefficient(m)[3] b3=coefficient(m)[4] b4=coefficient(m)[5] b5=coefficient(m)[6] # recalculando m_teste=b0+b1*x+b2*x^2+b3*x^3+b4*x^4+b5*x^5 A idéia seria que os pontos de m_teste dessem semelhantes a y mas isso não acontece ___ R-br mailing list 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.
Re: [R-br] Como precificar um treinamento de R básico?
Curso de 4 horas por uma semana 100 reais por aluno ___ R-br mailing list 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.
Re: [R-br] Como precificar um treinamento de R básico?
Passei um treinamento pra um grupo de 10 pessoas, a idéia foi saber o nível de conhecimento de R e programação primeiro, depois sondei o que esperavam do curso e as apostilas que tinha. Elaborei um roteiro direcionado atendendo as espectativas do grupo e termieni dando o conselho do 'use o help ___ R-br mailing list 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.
Re: [R-br] Correlação entre vetores de comprimentos diferentes
De forma genérica a lista gerada, corresponde a x e y de uma reta quero saber quais pontos do arquivo estão acima ou abaixo dessa reta porém o maior vetor eh o da reta e não o meu conjunto de dados ___ R-br mailing list 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.
Re: [R-br] Problemas com comando witch
para esse caso preciso saber quais posições de um vetor de 1 elementos são maiores do que os valores de um vetor de 203900 elementos (há NANs e elementos repetidos) e comparar um outro vetor de 1 com outro vetor de 203900 e ver se há posições iguais entre os resultados obtidos. como uso o match ou o which nesse casou ou outra função que resolva meu problema ___ R-br mailing list 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] Problemas com comando witch
tenho dois vetores: p=c(2,3,4,5,6,7,8,9) r=c(3,4,5) quando uso which(p==r) ou which(r==p) ele me retorna as posições onde r ocorre em p como fazer para saber onde p ocorre em r? obrigado (já dei um google nisso e não encontrei ) Tito Conte ___ R-br mailing list 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.
Re: [R-br] Escalas
se quiser ir mais a fundo vá em axis ___ R-br mailing list 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.
Re: [R-br] Problemas com comando witch
e para maior e menor? ___ R-br mailing list 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.
Re: [R-br] Problemas com comando witch
e não da certo para vetores de tamanhos diferentes ___ R-br mailing list 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.
Re: [R-br] Problemas com comando witch
e para o caso de maior e menor? ___ R-br mailing list 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.
Re: [R-br] Problemas com R 2.15
Instalei alguns pacotes NetCDF e deu certo, valeu pessoal ___ R-br mailing list 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] Problemas com R 2.15
Acabei de instalar o R para Ubuntu e não estou conseguindo instalar nenhum pacote novo, fiz isso para o pacote ncdf, e tentei pra mais alguns e nada todos dão o mesmo erro (achava que era problema com o NetCDF) a mensagem que aparece segue abaixo installing *source* package ‘ncdf’ ... ** package ‘ncdf’ successfully unpacked and MD5 sums checked checking for nc-config... /usr/bin/nc-config configure: creating ./config.status config.status: creating src/Makevars ** libs gcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -I/usr/include -fpic -O3 -pipe -g -c ncdf.c -o ncdf.o ncdf.c:3:20: fatal error: netcdf.h: No such file or directory compilation terminated. make: *** [ncdf.o] Error 1 ERROR: compilation failed for package ‘ncdf’ * removing ‘/home/tito/R/i686-pc-linux-gnu-library/2.15/ncdf’ The downloaded source packages are in ‘/tmp/Rtmpn8hiQ4/downloaded_packages’ Warning message: In install.packages(ncdf) : installation of package ‘ncdf’ had non-zero exit status Obrigado Tito ___ R-br mailing list 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.
Re: [R-br] Loop varrendo quantidade enorme de arquivos
Cara quando eu executo vários arquivos faço o seguinte a-dir(diretório dos seus arquivos) # aqui cria um vetor com as strings dos diretórios (nome dos arquivos no seu caso) tabela-matrix(NA,2,2) resultados-tabela s=0 e para dar nome para eles do jeito que vc está querendo e criar a variável vc pode pegar o nome while(s length(a)){ s=1+s # conta +1 na leitura do arquivo nome_do_arquivo-substr(a[s],0,nchar(a[s])-4) # retira o nome da variável do nome do arquivo (subtrai o .tif) arquivo-raster('lai{ano}{diajuliano}.tif')#carrega seu tif multiplicado-arquivo*0.001 # multiplica por 0.001 tabela[s,]-assing(nome_do_arquivo,multiplicado) # atribui o resultado da multiplicação a vaiavel gerada e armazena em uma tabela na linha escolhida lá } s=0 f=0 while(f length(a)){ s=s+1 f=4*s g=f-3 resultados[s,]-tabela[c(g:f),] colnames(resultados)-cut(seq(ISOdate(2001,1,1),to=ISOdate(2011,12,31),by=month),month) } deve funcionar qualquer coisa de uma olhada no que as funções fazem ___ R-br mailing list 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.
Re: [R-br] precisão numérica e which.max com empates
Tito Conte 2012/5/2 Pedro Emmanuel Alvarenga Americano do Brasil emmanuel.bra...@gmail.com Amigos de R, Estou tentando comparar alguns métodos de estimação e um dos passos necessários é encontrar valores maximos ou minimos de vetores numéricos. O problema numero 1 é que quando comparo visualmente os vetores ocorre que ocasionalmente acontece 2 ou 3 empates nos valores máximos ou mínimos como no vetor abaixo. Reparem que o valor .7 aparece na posição 21 e 23. Mas quando utilizo a função which.max ou max eu so encontro a posição 23. POssivlemente ha alguma coisa relacionado com precisão nume´rica apesar de não fazer muito sentido no momento. test.values$Youden [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.55 0.60 0.65 0.60 0.65 0.60 0.65 [21] 0.70 0.65 0.70 0.65 0.60 0.65 0.60 0.65 0.60 0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 max(test.values$Youden)[1] 0.7 which.max(test.values$Youden)[1] 23 which(test.values$Youden==max(test.values$Youden)) [1] 23 test.values$Youden[21]==test.values$Youden[23] [1] FALSE O problema 2 é que quando de fato ocorrem empates a função which.max retorna somente a primeira posição. Então, a questão seria como retornar todas as posições com valores máximos ou encontrar uma mediana dos valores empatados? Qualquer luz é bem vinda. Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.bra...@ipec.fiocruz.br email: emmanuel.bra...@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 ___ R-br mailing list 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@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.
Re: [R-br] precisão numérica e which.max com empates
Tente isso: which(test.values$Youden==max(test.values$Youden),arr.ind=TRUE) ___ R-br mailing list 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.