- R-br - DINF UFPR

modelo BCRM
by Jacqueline Cantú 27 Ago '14

27 Ago '14

Olá pessoal,Alguém sabe dizer se já existe algum pacote ou função implementada no R sobre o modelo espacial gaussiano bivariado de co-regionalização (BCRM)? Att, Jacqueline

2 1

Fit não linear com função que aprensenta singularidades.
by Michelle Baú Graczyk 27 Ago '14

27 Ago '14

Boa Tarde, Eu estou fazendo um fit nâo linear com um polinômio. Porém, a fiunção que eu estou usando : y= - {[a*(x^(2*n)) + b]/[(x-186)(x+186)]} apresenta duas assíntotas em +- 186 e por isso dá um erro no fit: Erro em numericDeriv(form[[3l]], names(ind), env): obtido valor faltante ou infinito quando avaliando o modelo. Pelo que eu entendi, o motivo é que a derivada também possui assíntotas nestes mesmo postos e daí na função vai para o infinito. Eu fiz D( (-(a*(x^(2*n)) + b)/((x-186)(x+186))), "x") e a mensagem de erro foi: Erro em x^(2*n)argumento não numérico para operador binário. Aí eu pensei que era por causa do n e deixei então só o x^2 e o erro foi o mesmo. Eu procurei na internet e não achei nada que explicasse o que está querendo dizer. O que seria um operador binário? Voltando às assintotas, eu achei uma função chamada ‘NLSstAsymptotic(xy)’ que é para assintotas verticais em gráficos mas eu tentei emplementá-la e nâo deu certo. Tem outra função que eu possa usar para que o fit seja feito mesmo que a função usada para o fit apresente singularidades? Abaixo está o código: file <- read.table(file="/Users/bau/Desktop/teste/DISparam.txt", col.names=c('x','y')) plot(x=file$x, y=file$y, xlab='Valores de x', ylab='Valores de y', lwd=0.5,) dados <- data.frame(x=file$x,y=file$y) #chute inicial guess <- list(a=1, b=3, n=3) # funcao usada para o fit # func <- function(x,n, a, b){ y <- (-(a*(x^(2*n)) +b)/((x-186)*(x+186))) return(y) } #executa o fit e armazena resultado na variavel 'fit' fit <- nls(y ~ func(x,n,a,b), data=dados, start=guess, lower=-200, upper=200) fit summary(fit) #gera valores de y a partir do resultado do fit, para plotar #o resultado no mesmo grafico que os dados yfit <- predict(fit) #plota a linha do fit lines(dados$x, yfit, col='blue',lwd=2) #carrega os parametros param <- coef(fit) Obrigada, Michelle

2 2

Re: [R-br] CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

1 0

Re: [R-br] CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

1 0

Dados em uma tabela
by jose gomes de souza filho Gomes 27 Ago '14

27 Ago '14

Boa Noite, Tenho o seguinte banco de dados em txt( em anexo): JOSÉ A................................. 2 B.....................................0 C...................................0 D.................................... 6 E.......................0 F............................. 1 G................................ 13 MARIA A................................. 3 B.....................................0 C...................................0 D.................................... 1 E.......................0 F............................. 1 G................................ 6 PEDRO A.................................0 B..................................... C...................................0 D...................................0 E.......................0 F.............................0 G................................ 3 Gostaria que ficasse assim: Nome A B C D E G F JOSÉ 2 0 0 6 0 1 13 MARIA 3 0 0 1 0 1 6 PEDRO 0 0 0 0 0 0 3 Como fazer isso? Att, Gomes

4 4

Sugestão de site para download de livros grátis.
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

9 13

Números grandes estourando!
by Jodavid Ferreira 25 Ago '14

25 Ago '14

Alguém sabe algum tipo númerico de armazenamento maior que o "as.double" no R. Estou fazendo uma simulação, porém os valores são tão altos que estão estourando -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Bacharelando em Estatística - Universidade Federal da Paraíba - www.de.ufpb.br Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://blog.jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

R lendo arquivos .spt
by Alisson Lucrécio 24 Ago '14

24 Ago '14

Caros colegas da r-br, Boa noite. Algum de vocês saberia se R consegue ler arquivos .spt? Obrigado. -- Alisson Lucrecio da Costa

1 0

kml - importação para análise de dados
by Leonard de Assis 24 Ago '14

24 Ago '14

1. Qual o melhor jeito de ler um arquivo kml contendo polígonos no R? 2. Como identificar em quais destes polígonos está contida uma determinada coordenada? É que eu tenho dados geoespaciais para vários polígonos parte armazenados em arquivos kml e parte eu tenho a referência (x,y) com os dados da ocorrência, preciso de uma maneira saudável para unir essas informações

3 5

Re: [R-br] Dúvida Simples - Histograma
by andrebvs＠bol.com.br 23 Ago '14

23 Ago '14

4 4