- R-br - DINF UFPR

Rotação de matrizes com coordenadas geográficas em UTM
by ASANTOS 08 Abr '13

08 Abr '13

Bom dia pessoal, Estou cá de volta com o problema de rotacionar matrizes para rotacionar um grid criado com expand.grid(), o exemplo numérico que havia postado na lista esta OK, porém agora quando emprego coordenadas geográficas em UTM não estou conseguindo plotar, pois os valores no eixo x estão estranhos e o eixo y esta negativo, alguém saberia me dar alguma sugestão, segue CRM abaixo, obrigado # p0<-c(524237,7978856)### Primeiro ponto (Coordenada geográfica conhecida em UTM) points(p0[1],p0[2], col="red") # # c0 <- seq(p0[1],p0[1]+400,by=100)### Criando o grid d0 <- seq(p0[2],p0[2]+500,by=100) d1 <- expand.grid(x=c0, y=d0) plot(d1[,1],d1[,2]) # d2<-as.matrix(d1) ang<- 60 #### Angulo a ser rotacionado Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang), cos(ang))) ## Matriz de rotação drt<-d2 %*% Q # Multiplicação de matrizes points(drt[,1],drt[,2], col="red") # -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

3 6

Modelo não linear- correção
by Fernando A. Souza 08 Abr '13

08 Abr '13

Caros colegas, Estou rodando uma série de equações não lineares, utilizando o modelo de gompertz. A minha dúvida é a seguinte: Quais estimativas eu devo fornecer para avaliação do modelo ? Eu penso em fornecer o "erro padrão residual" e o "log-likelihood". Quais a melhores informações a serem fornecidas? Abraço a todos -- Fernando Antônio de Souza Zootecnista, Dsc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

2 3

Criando uma mascara com shapefile sobre raster
by ASANTOS 06 Abr '13

06 Abr '13

Boa tarde pessoal, Estou com um problema e não consigo criar uma mascara com um arquivo no formato shapefile sobre raster usando a função crop() do pacote raster, por algum motivo meu objeto contorno que é um SpatialLinesDataFrame não esta servido para oferecer o polígono limítrofe da área, alguém poderia me dar um help, segue CRM: require("raster") A <- stack(c("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif"))## Imagem plotRGB(A) ## Plota a imagem contorno<- readShapeLines("Catas_Altas.shp") ## Contorno da área lines(contorno, col="red")Plota o contorno C <- crop(A,contorno) ## Recorta o que esta contido no contorno na imagem Erro em .local(x, y, ...) : nenhum slot de nome "legend" para esse objeto de classe "RasterStack" -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

2 13

Expandir valor!
by Sérgio Henrique almeida da silva ju 06 Abr '13

06 Abr '13

Prezados Como faço pra expandir um valor em uma coluna? Por exemplo, no banco abaixo, ele só me dá o valor da especie para a primeira linha que é o chefe do domicilio, mas eu gostaria de expandir esse valor para todos do cod 5001102003001. cod especie 5001102003001 1 5001102003001 5001102003001 5001102003001 5001102003001 5006309003002 2 5006309003002 5006309003002 5006309003002 alguém tem alguma ideia? 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 / 94429486

2 3

Modelo não linear
by Fernando A. Souza 06 Abr '13

06 Abr '13

Caros colegas, Estou rodando uma série de equações não lineares, utilizando o modelo de gompertz. A minha dúvida é a seguinte: Quais estimativas eu devo fornecer para avaliação do modelo, uma vez que o Eu penso em fornecer o "erro padrão residual" e o "log-likelihood". Quais a melhores informações a serem fornecidas? Abraço a todos -- Fernando Antônio de Souza Zootecnista, Dsc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

1 0

importar dados SRTM usando getData() do pacote raster
by Adriano S. Melo 06 Abr '13

06 Abr '13

Prezados colegas, Estou tendo dificuldade com isto: getData("SRTM", download=T, lat=-16.68, lon=-49.256) tentando a URL ' http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_27_16.ZIP' Content type 'application/x-zip-compressed' length 35447524 bytes (33.8 Mb) URL aberta ================================================== downloaded 33.8 Mb Erro em .SRTM(..., download = download, path = path) : file not found Já tentei especificar o argumento 'path' mas não funcionou. O problema acontece no Ubuntu12.04 (testei em 2 computadores). No windows 7 não tive problemas. Por outro lado, este funciona: getData('alt', country="BRA") Alguma sugestão? Abraços, Adriano S. Melo Ecologia - UFG

1 1

Utilizar o Hershey do R
by Jodavid Ferreira 05 Abr '13

05 Abr '13

Alguém saberia informar como posso utilizar o Hershey do R. Rodando o comando: demo(Hershey) estou precisando utilizar o \\#H0907 que se encontra na ultima tabela que se chama: "More Raw Hershey Escape Sequences" mas não sei como fazer! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

atualização de pacotes
by Rui Svensson 05 Abr '13

05 Abr '13

Não é muito inteligente.Quando vc aciona a atualização de pacotes, deveria estar subentendido que os pacotes já atualizados não precisam ser baixados de novo. Tal como as dependências.Isso já existe no update do antivirus: se está uptodate, não precisa atualizar.Deu pra entender ? RUI SVENSSON FONSECA analista de sistemas & analista de custos & especialista em planilhas financeiras analista de O&M & crítico estrutural compulsivoZEN

6 8

exportar um vetor do tipo list
by Aline Lipsky 04 Abr '13

04 Abr '13

Boa tarde, Eu tenho um vetor do tipo list com vetores de tamanhos diferentes. Preciso exportar do r no formato "csv" ou "txt", mas não estou conseguindo.Alguém sabe se tem como fazer essa exportação? Grata, Aline F. Lipsky

5 5

Boxplots
by Edson Lira 04 Abr '13

04 Abr '13

Caros amigos, tenho uma variiável chamada número de casos e uma outra (fator) chamada tipo de doação. Usando a rotina: boxplot(casos ~ doação,data=dc,xlab="Tipo de Doação", ylab="No. de casos",cex.axis=0.8,cex.lab=1) Tenho os dois boxplots (lado a lado) por tipo de doação (voluntário e depósito). A rotina points(mean(dc$casos),col="blue") me mostra o ponto da média no box somente em um box, como fazer para mostrar nos dois plots Edson Lira Estatístico Manaus-Amazonas

3 5