Janeiro 2016 - R-br

extraindo dados de arquivo img
by Hélio Gallo Rocha 21 Jan '16

21 Jan '16

Boa tarde Listeiros, Alguém poderia me dar uma mão nesta tarefa: O site : http://spirits.jrc.ec.europa.eu/ dispõem arquivos climatológicos com as extenções : img e hdr Quem duipõem um CRM ou indicaria um pacote que extraia os dados e suas coordenadas? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Gerador automático de listagens
by Nei 21 Jan '16

21 Jan '16

Oi Pessoal, Uma ferramenta que pode ser útil: http://www.fakenamegenerator.com/order.php Um Abraço, Nei

2 3

Erro função lme
by Victor Fioreze 20 Jan '16

20 Jan '16

Olá! Estou rodando um script com a função lme do pacote nlme e esta dando um erro.Após eu rodar o modelo usando lme aparece: Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Revisei os dados e o script e esta tudo ok.O que pode estar acontecendo? Att, Victor Ionatan FiorezeMédico Veterinário e ZootecnistaMestrando em Nutrição de Ruminantes, PPGZ-UFPelBolsista Capes Integrante GERUMEN (53) 81032887lattes: http://lattes.cnpq.br/5026587958713372 Este e-mail foi enviado por um computador sem vírus e protegido pelo Avast. www.avast.com

2 1

Conversão de matriz em raster georeferenciado
by ASANTOS 19 Jan '16

19 Jan '16

Prezados Listeiros, Gostaria de realizar a conversão de uma matriz em raster georeferenciado com duas informações, a representação de plantas em uma matriz (B) e juntar as informações geográficas conhecendo-se o espaçamento entre plantas (d1), no entanto, gostaria de saber se conforme o meu CRM abaixo essa conversão esta correta, pois não estou bem certo de que realização de um cbind() entre os objetos B e d1 correspondam a real posição desses elementos. Alguém poderia me ajudar? ##Dados artificiais injuria<-rnorm(1550, mean=12) ##Crio a matriz com os dados artificiais B = matrix( data=injuria, nrow=25, ncol=62) # Inserindo coordenadas espaciais p0<-c(680234.482625,7803361.55128)# Ponto inicial c0 <- seq(p0[1],p0[1]+182,by=2) d0 <- seq(p0[2],p0[2]+50,by=3) d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2 d1<-d1[-(1:14),] ## Removo algumas corrdenadas para ficar com o mesmo número de dados plot(d1[,1],d1[,2]) length(d1[,1]) require(raster) FORM<-cbind(d1,B) colnames(FORM)<-c('x','y','area') r <- rasterFromXYZ(as.data.frame(FORM)[, c("x", "y", "area")]) extent(r) <- c(680234.5,680416.5,7803362,7803410) projection(r) <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") plot(r) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 4

Múltiplas funções em uma única função
by ASANTOS 15 Jan '16

15 Jan '16

Prezados Listeiros, Gostaria de criar uma função com múltiplas funções em seu interior, no meu caso gostaria de criar uma função denominada MYsummary, com três funções no interior dela, sendo que à partir de banco de dados artificial: ##Banco de dados artificial who<-rep(c("JOSE","CARLOS"),40) dates<-rep(as.Date(c("2007-06-22", "2004-02-13")),40) var1<-rnorm(80) tableFE<-data.frame(who, dates,var1) head(tableFE) Tenho uma primeira função que identifica o objeto who, ou seja, quem são as pessoas: ## Função Identificadora da pessoa people<-function(x, db = tableFE) { x<-NULL results <- x x<-unique(db[,1]) results <- x results } PER<-people(db = tableFE) #--# print(PER) A segunda função me fornece o período da coleta de dados: ## Função do período timeFE<-function(x, db = tableFE) { x<-NULL minT<-min(db[,2]) maxT<-max(db[,2]) results <- x x<-cbind(minT,maxT) colnames(x)<-c("start","end") results <- x results } TIM<-timeFE(db = tableFE) #--# print(TIM) A terceira função calcula a média: MD<-mean(tableFE[,3]) print(MD) Porém quando tento juntar as três funções na minha função MYsummary(), para que as três funções sejam executadas concomitantemente, exibindo os três prints ao mesmo tempo não funciona, sendo: MYsummary<-function(x, db = tableFE) { c( ## Função Identificadora da pessoa people<-function(x, db = tableFE) { x<-NULL results <- x x<-unique(db[,1]) results <- x results print(results) } , ## Função do período timeFE<-function(x, db = tableFE) { x<-NULL minT<-min(db[,2]) maxT<-max(db[,2]) results <- x x<-cbind(minT,maxT) colnames(x)<-c("start","end") results <- x results print(results) } , MD=mean(tableFE[,3]) ) } MYsummary(tableFE)## Não funciona Alguém poderia me dar uma luz? Muito obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

1 1

Classificação de pontos em cada polígono em um shapefile
by ASANTOS 14 Jan '16

14 Jan '16

Caros Listeiros Eu criei um shapefile POLY e gostaria de lançar alguns pontos aleatórios (objeto pts) e depois classifica-los com a função over() para saber em qual polígono esta contido cada ponto, no entanto, todas as modificações que eu faço geram em incompatibilidade na avaliação de argumentos entre ppp e cordinates e gostaria que alguém me desse uma dica de como otimizar essa operação sem ficar fazendo inúmeras transformações de objetos, segue CRM: #Pacotes require(shapefiles) require(rgdal) require(maptools) require(spatstat) #------------------------------------------------------------------------------- #Criação de 4 polígonos sr1=Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409, 180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676, 332618, 332413, 332349)))),'1') sr2=Polygons(list(Polygon(cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437, 179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683, 331133, 331623, 332152, 332357, 332373)))),'2') sr3=Polygons(list(Polygon(cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875, 179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110), c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004, 329783, 329665, 329720, 329933, 330478, 331062, 331086)))),'3') sr4=Polygons(list(Polygon(cbind(c(180304, 180403,179632,179420,180304), c(332791, 333204, 333635, 333058, 332791)))),'4') #Converte em polígonos espaciais sr=SpatialPolygons(list(sr1,sr2,sr3,sr4)) srdf=SpatialPolygonsDataFrame(sr, data.frame(row.names=c('1','2','3','4'), PIDS=1:4)) #Cria o shapefile writeOGR(srdf, getwd(), 'POLY', 'ESRI Shapefile') #Leitura do shapefile contorno_line_X <- readShapeSpatial("POLY.shp") plot(contorno_line_X) #Seleção de pontos contorno_line_spat <- as(contorno_line_X,"owin") pts <- rpoispp(lambda=0.000012, win=contorno_line_spat) points(pts, col="blue") #Classificação dos pontos em cada área pts$PIDS <- over(pts, contorno_line_X)$PIDS pontos<-as.data.frame(pts) pontos #------------------------------------------------------------------------------- Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 6

Large Merge
by Wagner Tassinari 14 Jan '16

14 Jan '16

Olá pessoal, tudo bem ? Estou com um problema, estou querendo fazer o merge entre dois bancos de dados: > dim(bancolink) [1] 1630328 8 > dim(fichaAlink) [1] 2954930 6 Minha variável - chave primária está com o formato character. Daí qunado eu faço: > CadUnico_fichaA_merge = merge(bancolink, fichaAlink, by.x=c('CPF'), by.y=c('NUMERO_CPF')) Error in vecseq(f__, len__, if (allow.cartesian || notjoin || !anyDuplicated(f__, : Join results in 150596694 rows; more than 4585258 = nrow(x)+nrow(i). Check for duplicate key values in i each of which join to the same group in x over and over again. If that's ok, try by=.EACHI to run j for each group to avoid the large allocation. If you are sure you wish to proceed, rerun with allow.cartesian=TRUE. Otherwise, please search for this error message in the FAQ, Wiki, Stack Overflow and datatable-help for advice. Tentei fazer usando o formato data.table, daí ele dá um erro também: > CadfAmerge = bancolink[fichaAlink] Error in `[.data.table`(bancolink, fichaAlink) : When i is a data.table (or character vector), x must be keyed (i.e. sorted, and, marked as sorted) so data.table knows which columns to join to and take advantage of x being sorted. Call setkey(x,...) first, see ?setkey. Desde já agradeço a todos, -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

4 3

Formatando datas
by Wagner Tassinari 14 Jan '16

14 Jan '16

Olá pessoal, Estou com um probleminha. Em meu banco de dados, minha variável data se encontra desta forma: > data [1] 4042009 4042009 27072009 27072009 27072009 6082010 6082010 6082010 [9] 6082010 6082010> library(chron)> data2 = chron(dates=data, format=c('dd/mm/yy'))> data2 [1] 24/Aug/13036 24/Aug/13036 23/Aug/76090 23/Aug/76090 23/Aug/76090 24/Dec/18621 [7] 24/Dec/18621 24/Dec/18621 24/Dec/18621 24/Dec/18621 O certo seria: [1] 04/04/2009 04/04/2009 27/07/2009 27/07/2009 06/08/2010 ... Alguém poderia me dar uma dica ? Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis) <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> Este e-mail foi enviado por um computador sem vírus e protegido pelo Avast. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> <#DDB4FAA8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

6 5

Seleção de dígitos dos valores de um vetor
by ASANTOS 13 Jan '16

13 Jan '16

Prezados Listeiros, Tenho um vetor de dados db<-c(201002,486003,486016,486014,486004,486020,486002,486023,486027) e gostaria de extrair somente os três últimos dígitos de cada valor. Isso de uma maneira que eu tenha um novo vetor de dados sendo newdb: 002 003 016 014 004 020 002 023 027 Para isto pergunto, se alguém conhece alguma função que faça isso? Muito obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 3

Critério de Seleção de Modelos
by Andressa Siroky 10 Jan '16

10 Jan '16

Boa tarde! Quais critérios de seleção são adequados quando os modelos não são encaixados? Vi que AIC e BIC não são. Obg Andressa Siroky

5 4