automaticamente apos atingir um limite
de bounces.
Por hora vou aumentar este limite mas não creio que vá resolver.
Meu unico conselho até momento é assinar a lista com outro email.
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa
André
tenho recebido notificações de muitos emails yahoo e hotmail sendo
removidos automaticamente da lista
parece que há algum porblema específico com estes domínios.
Vou tentar investigar com os administradores do mailman novamente
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de
---
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR - Brasil
Tel: (+55) 41 3361 3573
VOIP: (+55) (41) (3361 3600) 1053 1066
Fax: (+55) 41 3361 3141
e
strsplit() permite partir a string e produz uma lista de vetores
Para pegar o ultimo pode inverter o vetor:
foo <- "Manuel dos Santos de Jesus"
rev(strsplit(foo, split=" ")[[1]])[1]
[1] "Jesus"
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoin
Olá Fernanda
voce quer dizer de média 0 ?
ou sem intercepto na função de covariância?
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR - Brasil
Tel: (+55) 41 3361 3573
VOIP: (+55) (41
(pode ser aproximado)
Com isto veria quais os centroides de menuicipios que estariam dentro
deste poligono
coordinates(OBJETO)
te dá as goordenadas e basta um poligono aproximado
qua contenha os centroides dos municipios
e com iuso seleciona subset com over()
Paulo Justiniano Ribeiro Jr
LEG
options(OutDec=",")
formatC(8, digits=2, format="f")
[1] "8,00"
On Tue, 6 Oct 2015, Andre Oliveira wrote:
Pessoas,
pesquisei na lista e no google, mas não resolvi o que desejo. Como modificar de
forma global ou especifica de uma tabela e exibir números inteiros com duas casa
decimais (00),
uma possível solução
rootSolve::polyroot()
se puder escrever o probçlema como encotnrar as multiplas raizes de um
polinomio, o que em geral é possível
On Fri, 2 Oct 2015, Wagner Wolff wrote:
Olá pessoal!
Estou precisando de uma ajuda, preciso encontrar uma solução para uma variável
que
Wagner
1 unico ponto talvez não seja adequado para testar uma distribuição e sim
para avaliar o quanto e se ele é atipico de algums forma
há aogumas propostar para medir pontos influentes que foram investigadas
pelo grupo da Eng. Agricola de Cascavel (orientados do Prof. Miguel Opazo)
que
pac - packageStatus()
names(pac)
terá uma dista com istalados e dispobíveis
exporte em um arquivo e depois leir e passe o vetor
dos instalados para o install.packages()
On Fri, 28 Aug 2015, Fernando Antonio de souza wrote:
Caros colegas.
Preciso formatar meu computador. Eu utilizo Ubuntu
A explicação é devida ao fato de que os pontos iniciais (menores
distâncias) da correlação determinam predominantemente o resultados da
krigagem (ou os maiores pesos, se quiser)
Nesta proção os variogramas não são tão diferentes assim.
Entretanto note que seu comentário diz respeito às médias
Wagner
talvez queira dar uma olhada em osso material do sinape 2012
www.leg.ufpr.br/mcie
On Tue, 14 Jul 2015, Wagner Tassinari wrote:
Gostaria de saber se alguém poderia me enviar ou indicar algum material a
respeito de estimação de parametros via máxima verossimilhança no R.
Obrigado
O INLA nóa est[a dispon[ivel no CRAN
www.r-inla.org
On Thu, 11 Jun 2015, Mauro Sznelwar wrote:
Tentei instalar esta biblioteca INLA e não consegui, por acaso é exclusiva
do LINUX ou não está disponível no CRAN?
Como os dados são contínuos seria necessário considerar uma distribuição
log(0) = -Inf
Veja que voce tem muitos zeros, mais que 25% dos teus dados,
visto que o 1o quartil ainda é zero
Não tenhosugestão específica mas acredito que voce poderia:
- ver se seus zeros são zeros mesmo ou valores censurados abaixo de um
cento limite de detecção
- sendo zeros mesmo eu
de uma olhada da função ncf::correlog
que parece fezer algo neste linha (centered Mantel statisdics for
muultivariate data)
On Wed, 10 Jun 2015, Daniel Marcelino wrote:
Caros, estou procurando uma função que me retorne a estatística de
Moran-I bivariada. Alguém da lista que faça/estude
considerar os -inf como zero de algum modo? Ou `e o caso de
encontrar um modelo/distribuição mais adequado mesmo?
Agradecido,
Samuel
Em 10 de junho de 2015 10:07, Paulo Justiniano paulo...@leg.ufpr.br escreveu:
log(0) = -Inf
Veja que voce tem muitos zeros, mais que 25% dos teus
Felinto,
sugiro que ambos esteja no formato do pacote sp
O shape pode ser lido diretamente neste formato
e voce precisaria converter a saida da krige.conv
em um SpatialPixelDataFrame
On Wed, 3 Jun 2015, Felinto COSTA wrote:
Prezados.
Sabem alguma maneira de plotar um arquivo em formato
Helio
O que o Luiz Ivan mencionou é outra possibilidade
me passe inbox um CMR e dados ou dump da sessão que exaqmino com mais
cuidado
On Sat, 16 May 2015, Hélio Gallo Rocha wrote:
Caros da lista,
Com certeza tem algo nos dados, pois além do erro na krigagem apareceu este ao
estimar por
e patamar, ai rola, mas de qualquer forma,
fiquei intrigado com o erro.
Em 18 de maio de 2015 10:27, Paulo Justiniano [via R-br]
ml-node+s2285057n4664521...@n4.nabble.com escreveu:
Helio
O que o Luiz Ivan mencionou é outra possibilidade
me passe inbox um CMR e dados ou dump
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que está
usnado modelo com covariáveis (trend) e que estes possuem valores muito
grandes ou pequenos
Se suas coordenadas forem em UTM considere dividi-las por 100
quanto ao warning é inofensivo mas isto foi resolvido na ultima
meteorológicas no Paraná que serão
comparadas a um modelo ECMWF.
Quanto a sessionInfo() :
R - 3.1.3
geoR - 1.7-4.1
Em 15 de maio de 2015 17:02, Paulo Justiniano [via R-br]
ml-node+s2285057n4664494...@n4.nabble.com escreveu:
Quanto ao erro pode ser vevido a diversas coisas mas suspeito que
mais detalhadamente solve() resolve um sistema de equações lineares
Description:
This generic function solves the equation ‘a %*% x = b’ for ‘x’,
where ‘b’ can be either a vector or a matrix.
O uso para inverter matrizes é possivel mas NAO RECOMENDADO na maioria dos
casos
Por
vinheta
pense nas vinhetas da globo...
uma animação rápida (ou nem tanto) chamando atenção para algo
(no caso a funcionaldade do pacote
On Fri, 24 Apr 2015, Anselmo Alves de Sousa wrote:
Pessoal,
pergunta nível iniciante: existe algum termo em português para Vignettes.
Alguém poderia me
isto é pq certos municipios e estados podem ser definidos por mais de um
poligono
Por isto p polygon que defini a unidade é formada por uma lista de
Polygons, que é o desesho de cada (sub) área
Pense por exemplo em um estado com ilhas que pertençam a ele.
On Fri, 6 Mar 2015, Augusto Ribas
Eu faria:
1. um subset somente com o 14 e out4ro com o resto ordenado da forma
desejada
2. rbind para agrupar novamente com o 14 primeiro
On Wed, 4 Feb 2015, Alexandre Santos wrote:
Caros Listeiros,
Tenho um data frame constituído da variável y (valores aleatório) e da
variável
Em Terça-feira, 23 de Dezembro de 2014 7:47, Paulo Justiniano
paulo...@leg.ufpr.br escreveu:
Andre
Voce dver ter experimentado outros valores mas
4 no top me parece muito
eu frequentemente uso algo como:
mar=c(3.2, 3.2, 0.3, 0.3), mgp=c(1.8, 0.8, 0)
quando quero reduzir espaços em figuras
!
André Oliveira Souza.
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de
Educação, Ciência e Tecnologia do Espirito Santo.
IFES
Em Segunda-feira, 22 de Dezembro de 2014 14:39, Paulo Justiniano
paulo...@leg.ufpr.br escreveu:
voce pode diminuir no latex ou talvez
voce pode diminuir no latex ou talvez no proprio R
voce tentou tambem diminuir o espaco ao redor das figuras
com par()?
par(mar=...)
ou similar?
On Sun, 21 Dec 2014, Andre Oliveira wrote:
Pessoas,
boa noite. Depois de muita pesquisas eu recorro a vocês. O que posso colocar no
preâmbulo de
Vamos por partes:
Vou considerar apenas anisotropia geométrica (que é o caso)
1. Para a versão empírica (amostral) baseada nos dados voce deve
calcular as correlações entre pares de pontos
agrupar isto em um grid e mapear por exemplo com countour() ou
contourplot()
2. Para a teórica/ajustada
Helio
se voce salvas sua sessao ao final qdo reinicias vai estar tudo ali
todo conteúdo vai estar no arquivo .RData do seu workspace
este arquivo pode ser inclusive renomeado para sua organização
por exemplo analise01.RData
e posteriormente carregado com
source(analise01.RData)
para
Sim e sim
Existem expressoes analiticas para transformar
as predicoes e erros de predição para escala original
(procure log-normal kriging e volte a escrever se nao encontrar
rapidamente)
Na geoR se voce usa
krige.conv()
com os dados na escala original e define lambda = 0
(ou passa um objeto
__
On 26/11/2014 13:21, Paulo Justiniano wrote:
Sim e sim
Existem expressoes analiticas para transformar
as predicoes e erros de predição para escala original
(procure log-normal kriging e volte a escrever se nao
On Tue, 11 Nov 2014, Waldelene Moura wrote:
Olá
Alguem poderia me dar algumas dicas de como proceder para submeter ao CRAN um
pacote que implementei?
Já contrui o arquivo .tar.gz através do R CMD BUILD.
seria bom rodar o
R CMD check _seu_pacote_
para ver se ele passa nos testes
veja
help(predict.BGCCM)
## S3 method for class 'BGCCM'
predict(object, locations, borders,
variable.to.predict = 1, ...)
variable.to.predict: scalar with options for values or 2 indicating
which variable is to be predicted.
On Mon, 3 Nov 2014,
as rotinas da geoR fazem a predicao retornando na escala original
para transforma;óes da familia box-cox atraves do argumento lambda
para alguns valores exitem texprtessoes analiticas e de forma geral
pode-se transofrmar de volta utilizando as simulacoes condicionais
On Wed, 29 Oct 2014,
faça o gráfico com a opções legend=F e depois use o comando legend()
para posicionar onde desejado
On Tue, 28 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal,
Sempre que utilizo a função variog4 do pacote geoR as legendas estão ficando
sobrepostas às curvas dos semivariogramas
Existe uma confusão de terminologia aqui
likfit() não ajusta um modelo ao variograma, e sim os parametros uma
distribuição normal multivariada com determinada estrutura de covariancias
aos dados.
a ideia é exatamente estimar tudo de uma so ver sem procedimentos por
etapar
corrigir
Taynana
eu realmente nao entendi direito o que voce quer
mas
dimnames()
vai retornar uma listas com nomes de linhas e colunas
que voce pode transformar em vetor com unique()
e eleminar repeticoes com unique().
On Mon, 27 Oct 2014, Taynãna César Simões wrote:
Oi, Davi
Queria o nome
table(cut(dados, breaks=c(...))
On Tue, 21 Oct 2014, Flávio Fagundes wrote:
Prezados(as),Por um acaso existe alguma função semi-pronta no R que produz uma
tabela de frequência para intervalo de classes com amplitudes de
tamanhos desproporcionais?
Imagino ser simples a programação (vou
defina estas coordenadas como pontos a serem interpolados e roda a função
novamente:
myLOCS - cbind(c(5, 25, 15), c(8, 10, 15))
krige.conv(..., loc=myLOCS, ...)
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal,
Para código reproduzível desta dúvida pode-se utilizar o exemplo
qualquer dado do R nestes formatos
cbind(predLOCS, kr$predict)
write.table(...)
alternativamente voce pode exportar como raster, shapefile ou outros
formatos
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Olá Prof. Paulo Justiniano,
Funcionou por aqui, muito obrigado pela ajuda!
Outra dúvida: É
NOOOSSAA quanta coluna duplicada
será mesmo que voce quer/precisa fazer isto?
On Mon, 20 Oct 2014, Giselle Davi wrote:
Prezados,
Preciso juntar colunas no R !!!
Tenho um arquivo assim:
D01 D02 D03 D04 D05 D06 D07 D08 D09 D10 F01 F02 F03 F04 F05 F06 F07 F08 F09
1 G G G G G
se M é sua matrix neste caso
M[,sort(colnames(M)]
vai ordenar como desejado
On Fri, 17 Oct 2014, Giselle Davi wrote:
Prezados,
Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1,
A4, A2, B2, B1, B3 ?
O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2,
ops...
order() no lugar de sort()
On Fri, 17 Oct 2014, Giselle Davi wrote:
Prezados,
Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1,
A4, A2, B2, B1, B3 ?
O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3
Utilizei o comando order, mas o
Nao, isto nao está implementando
teria que escrever a funcao objetivo
realisations é usada na incomum situação de que voce tem varias
replicacoes (objervação) independentes entre si do mesmo processo.
Ou seja a variavelrealizations é indicadora de qual reALIZAÇÃO DO MEMSO
PROCESSO É cada
geoR:::BGCCM
ou os modelos do spBayes e RandomFields
sao algumas opcoes
On Wed, 27 Aug 2014, Jacqueline Cantú wrote:
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
Prezados,
Agradeço imensamente a todos que colaboraram. As propostas de Benilton e Ivan
funcionaram com muita velocidade, mesmo
quando solicitei uma simulação de 1534 dados. No entanto, quando apliquei o
teste de aderência (Shapiro-Wilk), o valor
de p indicou que os dados NÃO obedeciam a uma
se o seu data-frame só tiver as colunas y e iso*
baste fazer:
lm(y ~. , ...)
ond o . significa todas as demais colunas do data-frame exceto o y
On Mon, 28 Jul 2014, Alexandre Loures wrote:
Boa tarde!
Preciso rodar uma regressão em que há muitas variáveis cujo nome inicia da
mesma forma.
pq ovce nao cria um data-frame apenas com as variáveis que quer usar?
On Mon, 28 Jul 2014, Alexandre Loures wrote:
Obrigado Paulo!
Mas meu data-frame tem outras colunas além de y e iso1.
Em 28/07/2014 15:38, Paulo Justiniano escreveu:
se o seu data-frame só tiver as colunas y e
Uma certa forma é possivel sim
as funcoes likfitBGCCM() e o respectivo metodo para predict()
permitem o chamado Bivariate Gaussian Common Component Model
para estimacao (MV) e predicao em prblemas bivariados
Mas sen duvida a gstat possui metodos mais gerais (embora nao estimacao
por MV)
Nao sabemos o contexto mas eu tambem iria refletir se desconsiderar
tais amostras causa algum problema.
Voce teria como testar tal singularidade antes de chamar a optim?
neste caso um if() poderia reolver.
On Wed, 23 Jul 2014, walmes . wrote:
Use a função try() nas porções de código sujeitas
eu comecaria por representar o objeto como um spatialpixeldataframe do sp
e ver se o funcao que monta a matrix de vizinhanca tem as ocoes desejadas
se tiver fica moleza pq é so multiplicar a patriz pelo vetor de dados
On Wed, 9 Jul 2014, Wagner Bonat wrote:
Ola pessoal !
Estou com um pequeno
lembrando que o pacote ubuntu vai funcionar se o R foi instalado
também pelos binários
se o R fpoi compilado o plugin r-gedit deve ser instalado a partir do
disponivel na pagina do porojeto
Basicamente é só baixar o arquivvo (.tar.gz ou análogo)
descomprimir e copiar para os diretórios
A pergunta foi sobre o estimador da **variancia** da população e nao o
estimador de lambda.
No mais os colegas já responderam, e para conferir
é de fato notar que se trata de uma distribuicao uniforme. Note que a
densidade nao inclui o *x* e o espacao para metrico depende do parametro.
anova()
com apenas um objeto de modelo nada acrescenta ao summary neste caso ,as a
funcao permite voce pssar varios modelos e comparar o ajuste deles
?anova
On Sun, 22 Jun 2014, Jefferson Ferreira Ferreira wrote:
OI Fernando, obrigado pela resposta.
Na verdade eu não tenho uma dúvida
Ir de 240 para 1200 pontos de preodicao nao deve causar toda esta
diferenca no tempo de predicao se voce estiver simplesmente obtendo media
e variancia da preditiva.
Mesmo simulando de preditiva eu nao esperaria este tempo todo
nao tenho familiaridade com alguns objetos que voce está usando
residuos tem media zero
para calcular o CV vc está entao dividindo por zero
alias, para que mesmo calcular o CV? pense nisto pois em geral isto é de
pouca ou nhemuma valia
On Fri, 13 Jun 2014, Hélio Gallo Rocha wrote:
Boa noite listeiros.
Uma dúvida de estatística básica.
A questão é
é a mesca coisa a funcao da geoR é só um metodo para a do MASS
On Sat, 14 Jun 2014, Helder Gramacho wrote:
Boa tarde pessoal,
Observando a ajuda da Transformação Box-Cox por meio do help(boxcox), vi que o
objeto precisa ser da classe fórmula
como no exemplo abaixo:
boxcox(Volume ~
como no R de forma geral os numeros como estes são armazenados de double
mas exibidos com 7 caracters por default
veja opcoes de argumento digits em options() e print()
On Mon, 9 Jun 2014, Alex wrote:
Bom dia aos colegas da lista!
Após usar a função read.geodata, percebi que no resultado
para comecar veja o resultado de
summary(Modelo)
On Fri, 6 Jun 2014, Adriano Borges Costa wrote:
Boa Tarde,
Sou iniciante no uso do R e estou com uma dívida bastante simples. Estou
fazendo a análise de uma política pública e quero medir a relação entre o
grau de execução da política e
Sugerimos isto a eles mas acharam melhor ter liste propria...
On Fri, 6 Jun 2014, Éder Comunello wrote:
Walmes, bom dia!
Existiria alguma possibilidade de unificar as listas, criando uma grande
comunidade de usuários do R 'falantes' do português?
Éder Comunello comunello.e...@gmail.com
se M é o seu obheto desta forma
as.vector(t(as.matrix(M)))
On Wed, 4 Jun 2014, luc.so...@usp.br wrote:
Para reforçar minha dúvida:
tenho a planilha com várias linhas
linha 1- 1 2 3 4 5 6 7
linha 2- 0 5 6 4 3 3 2
linha 3- 8 6 4 0 3 2 1
gostaria de transformar para:
coluna 1 coluna
sim especialmente nas variancias de predicao (krigagem) que vao ser
menores
On Tue, 3 Jun 2014, Hélio Gallo Rocha wrote:
Caros da lista,
Uma dúvida conceitual
Comparando os estimadores de semivariancia de Materon e de Cressie, observei
que de há um acréscimo dos valores da
semivariancia
coloque o objeto que vai receber as importações como uma **lista*
defina
x - list()
e dentro do loop
x[[i]]
On Mon, 2 Jun 2014, Paulo Pimenta wrote:
Prezados,
Estou com dificuldades para abrir uma série de arquivos (88 arquivos .dbf) por
meio de um loop e armazenar suas informações numa
Para isto voce teria que esccrver o pacote o que pode ser mais do que voce
deseja para carregar uma função
porque nao dar simplesmente um sourve no arquivo dela?
On Sat, 31 May 2014, Luiz Roberto Martins Pinto wrote:
Caros listeiros,
Preparei uma função e desejo colocá-la 'residente' (isto
Estatística Computacional
Universidade Estadual de Santa Cruz
Ilhéus-Bahia
luizroberto.u...@gmail.com
skype: lrmpinto
http://lattes.cnpq.br/2732314327604831
Em 31 de maio de 2014 19:24, Paulo Justiniano paulo...@leg.ufpr.br escreveu:
Para isto voce teria que esccrver o pacote o que pode
O Elias Krainski vai passar uns dias aqui em Curitiba e no período
vai apresentar um curso sobre o INLA.
A enfase será mais na motivaçãoe entranhas do INLA (e nao apenas no uso
do pacote e exemplos de dados)
A seguir um resumo, tb em www.leg.ufpr.br/seminarios
todos interessados são bem
os dados (dataset JURA?) são EXTREMAMENTE assimétricos e o variograma
usual nao vai ter bom comportamento e não é surpresa que
ajuste o modelo com efeito pepita puro
seguem algums comandos com algumas ideias explorando uma possivel
trensformação de variáveis
2 comentarios adicionais:
1. se
Em 5 de maio de 2014 10:08, Paulo Justiniano paulo...@leg.ufpr.br escreveu:
os dados (dataset JURA?) são EXTREMAMENTE assimétricos e o variograma
usual nao vai ter bom comportamento e não é surpresa que
ajuste o modelo com efeito pepita puro
seguem algums comandos com algumas
Complementando o comentário do Vinicius
voce tem 2 opcoes
- 1. faça sua função retornar o negativo da (log)verossimilhnaça
- 2. chame a optim() com a opções control=list(fnscale=-1)
On Sat, 26 Apr 2014, Vinicius Brito Rocha wrote:
Como o Benilton e o Wagner pontuaram, problemas de
O AIC da verossimilhança NAO pode ser comparado com qq AIC de ajuste de
variogramas.
O 2o seria o AIC do ajuste de um modelo nao linear ao variograma,
supondo qe os dados sao os pontos do variograma
Portanto, apeesar de mesmo nome, sao coisas completamente distintas e
totalmente nao
De fato, faltou fechar o parentesis na ultima linha.
Bem percebido
Em Segunda-feira, 24 de Março de 2014 23:04, Mauro Sznelwar
sznel...@uol.com.br escreveu:
Este código não funciona.
O seu loop for tá sobre escrevendo resultados em
cor.mult- cor.test(x, y[,i])
e portanto somente o
O seu loop for tá sobre escrevendo resultados em
cor.mult- cor.test(x, y[,i])
e portanto somente o ultimo está sendo guardado
(note tb que sua matrz de dados esá reciclando valores)
Uma sugestão:
x - 1:10
y - matrix(rnorm(200), nrow=10)
apply(y, 2, function(y, x=x) cor.test(x,y)$p.value
On
R is available as Free Software under the terms of the Free Software
Foundation's GNU General Public License
http://www.r-project.org/COPYING
On Wed, 15 Jan 2014, Leonard Assis wrote:
É um projeto que estamos pegando
Uma das coisas que pediram foi informar a licença do software
O método de inferencia.
geoRglm implementa inferencia por
(i) maxima verossimilhança
(ii) bayesiana
glm + ajuste de residuos
seria um método ad hoc de estimação, sem bases em pricípios gerais com
caracteristicas e propriedades conhecidas.
Alem disto, na 2a estratégia, a predição tb não segue
Nao!
o ponto . siginifica todas as variaveis no data-frame
alem da resposta (Y)
Se o data frame contem apenas Y, T, e C entao as especificaqcoes sao
iguais
On Fri, 13 Dec 2013, Fátima Lima Paula wrote:
Ana, o modelo com o ponto é o modelo nulo.
O outro está ajustando pelo T e C.
Não
mais uma repositorio
http://archive.ics.uci.edu/ml/
On Wed, 11 Dec 2013, Jakson Alves de Aquino wrote:
2013/12/11 Andre Oliveira andreolso...@yahoo.com.br:
Boa noite,
estou para fazer uma apresentação sobre análise descritiva de dados no R e
não tenho um banco de dados com descrição das
http://cran-r.c3sl.ufpr.br/web/views/Multivariate.html
On Fri, 6 Dec 2013, Jodavid Ferreira wrote:
Alguém pode me ajudar a como realizar uma análise Multivariada no R?
--
Jodavid Ferreira
Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com
Membro da equipe
tudo parece correto a menos que o pacote nao estaja na sua maquina
para comecar ... certifique-se que ele aparece listado em:
search()
[1] .GlobalEnvpackage:pscl package:vcd
[4] package:grid package:gam package:splines
[7] package:coda package:lattice
depende de como as informações sobre os bairros estão disponiveis
para voce
por exemplo, suponha que o objeto
rj é um SpatialPolygonsataframe (possivelmente importado de um
shapefile)
coordinates(rj)
vai te dar as coordenadas pontuais dos controides do bairro
On Wed, 4 Dec 2013, João
se isto são os noems de colunas em uma planilha (data-frame) com os dados
entao parece que voce nao tem coordenadas ai e vai precisar obter o mapa
de bairros do rj de alguma fonte
On Wed, 4 Dec 2013, João Sousa wrote:
Os dados são do tipo
ID
ANO
TIPO
BAIRRO
LOGRADOURO
DIASEM
FXHORA
MES
se M é uma matriz quadrada
diag(M) - 1.02 * diag(M)
aumenta em 2% cada valor da diagonal
igualmente pode-se somar ctes etc
On Tue, 3 Dec 2013, Marcos Toebe wrote:
Prezados
Estou desenvolvendo um estudo durante o doutorado, no qual busco definir por
reamostragem no r, o dimensionamento
Tb nao entendi.
se o objeto já está ali entao ele já existe no R?
Voce quer sabver este objeto e ler depois?
melhor explicar com mais detalhes
On Mon, 2 Dec 2013, Fernanda De Bastiani wrote:
Caros,
alguém conhece alguma função no R que possa substituir a função scan()?
No sentido de que
Invertendo a pergunta... por que deveriam ser iguais??
likfit estimou um modelo gaussiano para os dados
com um termo espacial (efeito aleatório espacial)
No caso um modelo com 3 parametros, um de mdia e 2 de covariancia.
glm dos dados contra os preditos estimou uma regressao linear destes,
e
Fernanda
Imagino quie os warnings sejam quando voce utiliza isto dentro de uma
funcao sua que está escrevendo (em um pacote talvez) seria isto?
quais sao os warnings?
A funcao acessada é de outro pacote?
se for isto de uma olhada no manual Writing R extensions para ver se há
alguma
own objects:
Muito obrigada,
Fernanda
Em Quarta-feira, 6 de Novembro de 2013 12:10, Paulo Justiniano
paulo...@leg.ufpr.br escreveu:
Fernanda
Imagino quie os warnings sejam quando voce utiliza isto dentro de uma
funcao sua que está escrevendo (em um pacote talvez) seria isto?
quais sao os
to use ::: for its own objects:
Muito obrigada,
Fernanda
Em Quarta-feira, 6 de Novembro de 2013 12:10, Paulo Justiniano
paulo...@leg.ufpr.br escreveu:
Fernanda
Imagino quie os warnings sejam quando voce utiliza isto dentro de uma
funcao sua que está escrevendo (em um pacote talvez) seria isto
que cada nível tem uma
recomendação de aplicação do fertilizante.
Se souber qual a porcentagem de cada fatia sobre a área total do componente
fósforo poderei recomendar a quantidade do fertilizante na medida
certa.
Hélio
Em 6 de novembro de 2013 08:54, Paulo Justiniano [via R-br]
ml-node
Parece que o posto no blog estṕa incompleto, sugiro contatar o autor do
mesmo
TALVEZ o desejado necessite incluir
jura - prediction.dat
coordinates(jura) = ~Xloc+Yloc
x - krige(log(Pb)~1, jura, jura.grid, model = m)
On Sat, 2 Nov 2013, Helder Gramacho wrote:
Boa noite,
sou iniciante na
Eu prefiro usar o tempo
modelar a média
ao inves de retirar a tendencia
pois aco que o 1o expressa mais precisamente o que de fato é feito
On Thu, 31 Oct 2013, Hélio Gallo Rocha wrote:
Caros da lista,
Estudando outros dados, não houve problema na validação cruzada para retirada de
m1-lm(trh~mus+tipo, data=subset(b1, tipo==0))
On Tue, 29 Oct 2013, Edson Lira wrote:
Caros amigos, estou tentando fazer um estudo de regressão com duas variáveis
preditoras. Dados e comandos abaixo:
mus-c(2,6,8,3,2,7,9,8,4,6)
tr-c(e,m,e,m,e,e,m,m,e,e)
x - c(10, 12, 11, 5, 4, -5, -10, -11, -15, 4, 5, 6, 8, 10, 8, 5, 4, -6, -2, 4,
8)
x1 - ifelse(x 0, 1, 0)
x[which(diff(x1) == 1)+1]
[1] -5 -6
On Mon, 28 Oct 2013, Alisson Lucrecio wrote:
Caro Colegas da lista r-br,
Boa tarde.
Eu estou precisando encontrar o primeiro numero negativo
Ainda preciso dedicar um tempo para estudar adequadamente o que foi colocado,
mas por ora gostaria apenas de esclarecer se a 'retirada' do efeito
da tendência (por ajuste polinomial) é uma estratégia válida ou há alguma
ressalva a ser feita. Agradeço se os senhores, ou outro colega da lista,
...@gmail.com
escreveu:
Caros da lista,
Obrigado pelas dicas, se puderem me elucidar em dois pontos:
1. Quais outros pacotes trabalham com geoestatistica.
2. Com analisar o residuo se ele é gerado na validação?
grato
Hélio
Em 17 de outubro de 2013 09:20, Paulo Justiniano [via R-br]
ml-node
Comentário abaixo
Prezado Helio, boa tarde!
Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior
não faz muito sentido, já que a krigagem é um estimador exato e a
operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia
ser zero/nulo). O que
ver post anterior sobre patch a ser incorporado em breve
On Sun, 20 Oct 2013, Hélio Gallo Rocha wrote:
Eder,
valeu pela força.
Amanhã cedo vou analisar sua sugestão, mas de antemão adianto que sem a
validação cruzada não tem como ter o dado do erro, nem tampouco fazer a
krigagem, a não ser
Os comentários estao ininteligiveis.
A verossimilhanca nao gera mapa, e sim estima parâmetros.
nem faz a krigagem -- nao dá apra entender o que quis dizer
os codigos dos pacotes geoestatísticos (qq um deles) para obeter as
quantidades desejadas
sugiro ser mais rpeciso e eventualmente enviar
Caro Wenceslau
Aproveitando a discussão para aprender sobre a validação de modelos em
geoestatisticos. Dado que estest modelos e tecnicas são novas para mim.
A krigagem é um estimador exato se o numero de pontos estimados for o mesmo dos
avaliados, não ?
CVes posta enviado agora pouco sobre
Há um problema na geoR para fazer a validacao cruzada com trend (model de
media nao constante.
Algumas opcoes:
1. fazer um pequeno programa que remova uma observacao por vez e chame a
krige.conv()
2. Usar outros pacotes do R
3. Fazer a valizacao nos residuos
4. Pensar se a validacao
se os zeros sao reais entao a funcao gama nao é adequada ...
entretanto talvez seus zeros sejam dados abaixo de um limite de detecao,
portanto sao dados intervalares
entre zero e o limite de detecao
voce pode escrever a verossimilhanca para ajustar os dados assim
ou em uma solucao aproximada,
1 - 100 de 376 matches
Mail list logo