- R-br - DINF UFPR

Problema para mudar projeção em um raster
by ASANTOS 22 Jul '13

22 Jul '13

Boa tarde Pessoal, Criei um raster e referencie em UTM sendo: require(raster) require(rgdal) require(sp) #Criando raster r <- raster(nc=1000, nr=1000) set.seed(20130622) stackIm1 <- stack(lapply(1, function(x) setValues(r, round(runif(ncell(r))* 255))))## Simulation red band # # Definindo sistema de projeção r.geo <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") # geographical datum WGS84 proj4string(stackIm1) <- r.geo # mas quando tento converter para latlong com spTransform: stackIm1.ll<-spTransform(stackIm1, CRS("+proj=longlat+datum=WGS84")) Erro em spTransform(stackIm1, CRS("+proj=longlat+datum=WGS84")) : erro na avaliação do argumento 'CRSobj' na seleção do método para a função 'spTransform': 'Erro em CRS("+proj=longlat+datum=WGS84") : unknown projection id Obtenho este erro e não consigo achar a solução, alguém poderia dar um help, 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 ======================================================================

4 14

Duvida Concatenação
by Sérgio Henrique almeida da silva ju 22 Jul '13

22 Jul '13

Prezados Estou tentando concatenar uma variável com a seguinte condição: Se ela tiver tamanho 1, gostaria de concatenar 000, se ela tiver tamanho 2 concateno 00 e se ela tiver tamanho 3 concateno 0. Eu pensei na seguinte função: dados$len <- nchar(dados$var) func <- function(){ if (dados$len == 1){ dados$var1 = paste("000", dados$var, colapse="") } else if (dados$len == 2){ dados$var1 = paste("00", dados$var, colapse="")} else if (dados$len == 3){ dados$var1 = paste("0", dados$var, colapse="")} } func() Mas está dando a seguinte mensagem de erro: Mensagens de aviso perdidas: In if (dados$len == 1) { : a condição tem comprimento > 1 e somente o primeiro elemento será usado Vi em outro tópico que o Benilton sugere usar with, mas não sei como fazer. Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

3 9

Funcao L(K Ripley's transf.) com mais de 3000 pontos
by Felipe Sodre Barros 22 Jul '13

22 Jul '13

Boa noite. Estou fazendo uma análise de processos pontuais, utilizando o pacote spatstat. Ao tentar realizar a análise L, me deparei com um problema: Meus dados superam o 'nlarge' (defindo como 3.000). E por mais que eu tente aumentar o 'nlarge', tenho como resultado que os dados são muito grandes, sendo por isto impossível de alocar XXmb. Contudo, se eu rodar tal análise sem a geração de simulações (envelope), ele apenas me informa que o n dos dados que possuo são superiores a 3000 e por isto apenas calculará o 'bordar correction'. A minha duvida está, em que de fato representa o 'border correction' calculado. Posso utiliza-lo como se este representasse o padrão de distribuição do processo analisado? Desde já vos agradeço. Att

2 1

subset
by Fátima Lima Paula 21 Jul '13

21 Jul '13

Olá companheiros do R. Quero criar um dataframe a partir de algumas categorias de uma variável. Não sei o que estou fazendo de errado. Exemplo: nome=c("a","b","c","d","e","f","g","h","i","j","k") ide=c("1","2","3","1","2","4","6","6","9","2","1") df=cbind(nome,ide) df=as.data.frame(df) df #quero um dataframe apenas com os ide = 1, 3 e 6 teste=subset(df,ide==c("1","3","6")) teste Mensagens de aviso perdidas: 1: In is.na(e1) | is.na(e2) : comprimento do objeto maior não é múltiplo do comprimento do objeto menor 2: In `==.default`(ide, c("1", "3", "6")) : comprimento do objeto maior não é múltiplo do comprimento do objeto menor teste #Cria algo que não tem nada a ver: nome ide 1 a 1 4 d 1 Alguém pode me dizer onde estou errando? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 2

Problema no read.fwf.
by Rodrigo Coster 21 Jul '13

21 Jul '13

Caros, estou usando o read.fwf para ler um banco, e obtive a seguinte mensagem de erro: > len <- c(13, 11, 11, 2, 2, 1, 70, 8, 4, 8, 8, 2, 35, 7, 15, 11, 16, 5, 8, 2, 8, 7, 5, 8, 2, 13, 4, 4) > dados <- dadostemp <- read.fwf('05.txt', width=len, header=FALSE, stringsAsFactors=FALSE) Erro em scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : linha 28659 não tinha 28 elementos Pensei que na linha 28659 estava faltando algum caractere, mas verifiquei que todas as linhas tem o mesmo tamanho: > temp <- readLines('05.txt') > table(nchar(temp)) 290 386103 > length(temp) [1] 386103 Alguem tem alguma idéia do que pode ser? Ou de como nao usar o read.fwf e usar o carregado pelo readLines?

3 4

Eliminar o nome das colunas de um data.frame()
by Tiago Souza Marçal 20 Jul '13

20 Jul '13

Boa noite pessoal, gostaria de saber se é possível dar nomes em colunas alternadas de um data.frame() ou apagar o nome das colunas. CMR: x<-c("a","a","a","a") y<-c(3,3,3,3) data.frame(x,y) x y 1 a 3 2 a 3 3 a 3 4 a 3 # Este é o output normal do R Gostaria de saber se posso substituir os nomes x e y por somente um como por exemplo col1. col1 1 a 3 2 a 3 3 a 3 4 a 3 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 #################################################################

3 2

problemas com o símbolo %
by alanarocha＠sapo.pt 20 Jul '13

20 Jul '13

boa tarde, o meu ficheiro de dados tem o simbolo da % estou a utilizar stringAsFactor=False para o simbolo ser eliminado mas quando faço a leitura aparece-me na mesma como posso eliminar o simbolo % sem alterar a formatação dos numeros exemplo ter 70% e ao tirar o programa mudar para 0,7 o meu codigo é: setwd("f:/Ana_trabalho/indicadores_semestrais/1 semestre 2013") dados<-read.table("Hgb.csv",header=TRUE, sep=";",stringsAsFactors=FALSE,dec=",") #for (i in 1:7) { # Executa o comando a seguir para as colunas 2 a 6 (1º semestre de 2013) dados[,i] <- as.numeric(gsub('%','',dados[,i])) # Tira o símbolo de percentual e transforma texto em número #} colnames(dados) <- c("Clinica","Jan","Fev","Mar","Abr","Mai","Jun") dados os dados estão neste link: https://www.dropbox.com/sh/b9t2w2ljbi07ow4/8gONfv_4Sj/Hgb.csv Cumprimentos Ana

2 1

Matriz
by Tiago Souza Marçal 20 Jul '13

20 Jul '13

Boa tarde pessoal estou tentando construir uma matriz com letras e números, entretanto, gostaria de retirar as aspas de ambos que aparecem na matriz. Alguém tem alguma sugestão? Segue abaixo um CMR: x<-c("a","a","a","a") y<-c(2,2,2,2) matrix(c(y,x),nc=2) 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 #################################################################

2 1

Fatorial duplo com tratamento adicional
by ana paula coelho madeira 19 Jul '13

19 Jul '13

Boa noite! Tenho uns dados de um experimento em DIC em fatorial duplo com um tratamento adicional. Como devo entrar com os dados no R para usar a função fat2.ad.dic() do pacote ExpDes. Desde já agradeço Ana Paula

4 7

Analise de variância fator com níveis quantitativos
by Alisson Lucrecio 19 Jul '13

19 Jul '13

Caro Colegas r-br, Boa tarde. Nesse exemplo abaixo é certo transformar um dado quantitativo em fator para depois fazer a analise de variância? Obrigado Alisson Lucrécio da Costa adi <- expand.grid(cult=gl(1,5,la=LETTERS[1]), fert=101) fat <- expand.grid(cult=gl(5,5,la=LETTERS[2:6]), fert=seq(50,150,25)) da <- rbind(adi, fat) theta <- c(c(-194.29, -197.26, -197.85, -203.03, -190.20, -190.45), c(9.1797, 8.2686, 8.6437, 9.3438, 8.8773, 8.1872), c(-0.03382, -0.03479, -0.03632, -0.03341, -0.03597, -0.03675)) X <- model.matrix(~-1+cult/(fert+I(fert^2)), data=da) da$eta <- X%*%matrix(theta) da$y <- da$eta+rnorm(nrow(da),0,30) require(lattice) xyplot(y~fert|cult, data=da, type=c("p","a")) da$Fert <- factor(da$fert) #?????? levels(da$Fert) levels(da$cult) m0 <- aov(y~cult*Fert, data=da) anova(m0)

2 1