Flávia, uso os pacotes FactomineR e BiplotGUI, dynGraph, para análise
multivariada e obtive bons resultados e gráficos muito bem definidos, você
encontra video aulas no youtube e na própria biblioteca do R tem arquivos pdf
com exemplos. Espero que ajude.
Att,
Ari Clecius Alves de Lima
o pacote spdep tem uma implementacao para o car normal
o R pode interfacear com turos programas que estimam este modelo (*BUGS
por exemplo)
o pacote INLA tb ajusta o CAR via Inf Bayesiana
Veja mais na sessão Spatial Task View da página do R que revisa os
recursos disponíveis nesta área
$ locate Rscript
No meu caso com R compilado das fontes:
paulojus@pj:~$ locate Rscript
/usr/local/bin/Rscript
/usr/local/lib/R/bin/Rscript
/usr/local/share/man/man1/Rscript.1
On Fri, 13 Apr 2012, Leandro Marino wrote:
Caros,
alguem saberia me dizer aonde fica localizada a pasta Rscript
http://cran.r-project.org/web/views/Multivariate.html
[]s
Leonard de Assis
assis dot leonard at gmail dot com
Em 12/04/2012 16:17, Flávia Souza escreveu:
Boa tarde a todos. Alguém pode me ajudar com material/link sobre
análise multivariada utilizando o R? Gostaria de saber
principalmente
Paulo,
No meu caso coloquei o comando locate Rscript e não retornou nada.
Ubuntu Server 10.04 LTS
R 2.15.0 instalado via apt-get
Alguma sugestão?
-Mensagem original-
De: r-br-boun...@listas.c3sl.ufpr.br
[mailto:r-br-boun...@listas.c3sl.ufpr.br] Em nome de Paulo Justiniano
Enviada em:
rbind não serve?
#exemplo:
tab3-rbind(tab1,tab2)
tab3
Em 13 de abril de 2012 09:36, ASANTOS alexandresanto...@yahoo.com.brescreveu:
Bom dia Pessoal,
Criei um objeto com duas variáveis:
v1-1:11
v2-rnorm(11)
tab1-cbind(v1,v2)
e outro:
v1-1:11
v2-rnorm(11)
tab2-cbind(v1,v2)
e
On 04/13/2012 10:59 AM, Augusto Ribas wrote:
Wikipedia não resolve?
http://en.wikipedia.org/wiki/Tukey%27s_range_test
O Wikipedia não resolveria o meu problema.
No livro do Crawley, o R book, na pagina 482 ele da uma explicação bem
legal sobre Comparações multiplas.
O R book pode ser que
Tem algum motivo específico pra vc usar header=F nas chamadas de read.table()?
Se vc colocar header=T em ambos os casos, o rbind() vai funcionar
---
Fernando Mayer
Universidade Federal de Santa Catarina - UFSC
Departamento de Ecologia e Zoologia - ECZ/CCB
URL:
o problema eh q vc esta' escrevendo os arquivos *com* cabecalho
(header) e lendo *sem*...
o q vc quer nao eh:
result1=read.table(tab1.txt,header=T)
result2=read.table(tab2.txt,header=T)
tab12=rbind(result1, result2)
?
2012/4/13 ASANTOS alexandresanto...@yahoo.com.br:
Augusto,
Fiquei
Humm, eu não entendi, mas qd vc da read.table vc não pode usar header=T? O
cabeçalho repete varias vezes no dataframe? Assim o r nem ta lendo numero
str(result1)
Em 13 de abril de 2012 10:19, ASANTOS alexandresanto...@yahoo.com.brescreveu:
Augusto,
Fiquei ate com vergonha da minha
Augusto,
Fiquei ate com vergonha da minha pergunta, mais na verdade eu
elaborei ela de maneira errada para o propósito da minha dúvida, na
realidade eu tenho:
v1-1:11
v2-rnorm(11)
tab1-cbind(v1,v2)
write.table(tab1, file =tab1.txt, col.names=T, row.names=F)
v1-1:11
v2-rnorm(11)
Não seria dessa forma?
v1-1:11
v2-rnorm(11)
tab1-data.frame(cbind(v1,v2))
v1-1:11
v2-rnorm(11)
tab2-data.frame(cbind(v1,v2))
tab12-rbind(tab1,tab2)
tab12
Edson Lira
Estatístico
Manaus-Amazonas
De: ASANTOS alexandresanto...@yahoo.com.br
Para:
No livro Designs and Analysis of experiments do Montgomery na página 96 tem um
tópico que mostra a matemática do teste de Tuckey. Tem também alguns exemplos
usando a LSD.
Espero ter ajudado.
Edson Lira
Estatístico
Manaus-Amazonas
De: Augusto Ribas
Olá,
Alguém sabe qual a funçao para rodar uma ancova hierárquica?
O modelo que quero testar é o crescimento de plantulas (RGR) em funçao do
tipo de capoeira (Standtype), espécies de plantulas (spp), tratamento de
abertura do dossel (treatment) e tendo como covariáveis a idade da capoeira
(age) e
Pessoal,
estou usando o comando cut, para criar a legenda de uma mapa do brasil que
estou plotando.
O problema é que minha variável tem uma escala muito alta e por isso ao
executar o comando cut estou ficando com algo parecido com isso:
cut(mapa$dados),10)
[1] (397,1.32e+05]
Eu faria algo como:
set.seed(1)
x = runif(100)
y = runif(100, 1e5, 1e6)
z = c(x, y)
qs = quantile(z)
lbls = prettyNum(qs)
lbls = paste((, lbls[-length(lbls)], ,, lbls[-1], ], sep='')
w = cut(z, qs, labels=lbls)
b
___
R-br mailing list
Caros,
estou postando abaixo um loop com o qual estou trabalhando. É longo, mas se
vocês olharem para as linhas finais, verão que estou querendo fazer um plot com
os resultados. Mas fazendo dessa maneira, enquanto o loop corre, o R vai
produzindo um gráfico com só um ponto (referente ao
Como seu exemplo nao e' reproduzivel, darei apenas indicacoes gerais
que vc deve adaptar a seu codigo:
1) antes de iniciar o while(), inicie a janela grafica ja' usando os
limites para os eixos X e Y, por exemplo:
plot(1, type='n', xlim=c(0, 1000), ylim=c(0, 2))
## vc tera um grafico vazio
On 04/13/2012 01:34 PM, Edson Lira wrote:
No livro Designs and Analysis of experiments do Montgomery na página
96 tem um tópico que mostra a matemática do teste de Tuckey. Tem
também alguns exemplos usando a LSD.
Era exatamente isso o que eu estava procurando!
Muito obrigado Edson!
Abraços!
No meu conjunto de dados têm muitos dados em branco e o R não reconhece
como NA, como posso fazer para que ele transforme o em branco para NA?
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 -
Sem a mensagem de erro fica difícil ajudar! Poste a mensagem de erro, ou
melhor, nos forneça um CRM para que a ajuda seja mais efetiva.
(S,f,P)
Allaman
\begin{signature}
=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Tá importando com que função? De que extensão? Cadê CMR? Peça o help da
função, se for das usuais deve ter um argumento na.strig=. Comece por ele.
À disposição
Walmes.
==
Walmes Marques Zeviani
LEG (Laboratório de
Oi Walmes
Estou importando do SPSS e uso a foreign.
Em 13 de abril de 2012 16:19, Walmes Zeviani walmeszevi...@gmail.comescreveu:
Tá importando com que função? De que extensão? Cadê CMR? Peça o help da
função, se for das usuais deve ter um argumento na.strig=. Comece por ele.
À disposição
Antes de tudo deve-se consultar o help da função. help(cut) te mostra que
argumento dig.lab= controla isso.
À disposição.
Walmes.
==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S,
Estou usando o seguinte comando:
library(foreign)
dados - read.spss(tmp.sav, use.value.labels=TRUE, max.value.labels=Inf,
to.data.frame=TRUE)
Em 13 de abril de 2012 16:21, Sérgio Henrique almeida da silva ju
sergio.edfis...@gmail.com escreveu:
Oi Walmes
Estou importando do SPSS e uso a
No arquivo de origem qual o caracter representador de NA, é ponto, é
espaço, é missing, é asterisco? Passe para o argumento na.strig= (ou
similar a este na função de leitura que você usa) o caratere
correspondente. Se não resolver providencie um CMR.
À disposição.
Walmes.
Cara.
Depois que vc fez:
library(foreign)
dados - read.spss(tmp.sav, use.value.labels=TRUE, max.value.labels=Inf,
to.data.frame=TRUE)
#faça o seguinte
exemplo-dput(dados)
#se o arquivo for muito gigante faz o dput de um pedaçinho so, so garanta
que no pedacinho vc mostre esses brancos que ta
Oi
Reduzi o banco e rodei o comando:
exemplo-dput(dados[1:10,])
structure(list(nu_notific = structure(c(9L, 1L, 13L, 16L, 4L,
5L, 2L, 14L, 11L, 3L), .Label = c(0010777, 0014623, 0036055,
0180019, 0180280, 0180567, 0180573, 0181106, 0181303,
0841583, 1634645, 1778689, 2082159, 2082161, 2087955,
#abrindo seus dados como exemplo
exemplo-structure(list(nu_notific = structure(c(9L, 1L, 13L, 16L, 4L,
5L, 2L, 14L, 11L, 3L), .Label = c(0010777, 0014623, 0036055,
0180019, 0180280, 0180567, 0180573, 0181106, 0181303,
0841583, 1634645, 1778689, 2082159, 2082161, 2087955,
9138695), class =
Quais as implicações de fazer isso para futuro uso dos dados?
Em 13 de abril de 2012 17:05, Augusto Ribas ribas@gmail.com escreveu:
#abrindo seus dados como exemplo
exemplo-structure(list(nu_notific = structure(c(9L, 1L, 13L, 16L, 4L,
5L, 2L, 14L, 11L, 3L), .Label = c(0010777, 0014623,
Boa tarde pessoal,
Tenho os centróides em UTM de 11 parcelas quadradas de 40x40m e
estou tentando calcular os vértices, inclusive com a repetição do
primeiro vértice para fechar o polígono. Ao final, quero armazenar os
resultados. Para tanto fiz:
Olá a todos,
Estou com o seguinte problema, quero calcular a integral de 0 até HT( valor da
altura total dde cada arvore presente na base de dados em anexo).
Rotininha:
rm(list=ls(all=T))
setwd(C:/Documents and Settings/Daniela/Desktop/Ajustes no R)
mdados-read.table(INTEGRAL.txt, h=T,
Sandro,
usando format com nsmall = 20, eu tenho um resultado de 1.00.
Mas eu tenho certeza que esse número é na realidade alguma coisa menor que 1.
Eu não sei se isso ocorre pelo modo como o R mostra o número (com 22 casas) ou
se é um problema em como os valores estão guardados
exemplo reproduzivel.
bastaria vc dizer quais sao os valores de sum(X2) e os graus de liberdade...
b
2012/4/13 Vitor Aguiar vitor.agu...@me.com:
Sandro,
usando format com nsmall = 20, eu tenho um resultado de 1.00.
Mas eu tenho certeza
esses valores são altos. É um banco de dados com muitas variáveis.
Veja:
[1] D3S1358
Chi square
481252.8
degrees of freedom
276
p-value
1
[1] D5S818
Chi square
7583.649
degrees of freedom
105
p-value
1
[1] D7S820
Chi square
72249.05
só pra explicar, o Chi square é o sum(X2).
On Apr 13, 2012, at 2:46 PM, Vitor Aguiar wrote:
esses valores são altos. É um banco de dados com muitas variáveis.
Veja:
[1] D3S1358
Chi square
481252.8
degrees of freedom
276
p-value
1
[1] D5S818
Chi
Evite dados em anexo. Onde não tá dando certo? Dê mais
informações/contextualize melhor o seu problema. Não precisamos dos dados
para integrar, só de uma particular função. Envie um CMR usando apenas o
corpo da mensagem, para digamos apenas 3 árvores, que acredito ser
suficiente. Leia o guia de
Eu vou apenas usar o primeiro marcador como exemplo (isso
representa o Minimo de Codigo Minimo Reproduzivel). Vou tbm assumir
o usual de teste qui-quadrado, no qual usa-se pvalor = P(Chi2
valorObservado)
sx2 = 481252.8
dfreedom = 276
pvalor1 = pchisq(sx2, df=dfreedom, lower.tail=FALSE) ## vc
atualize a sua base antes do locate:
sudo updatedb
e repita o locate
On Fri, 13 Apr 2012, Leandro Marino wrote:
Paulo,
No meu caso coloquei o comando locate Rscript e não retornou nada.
Ubuntu Server 10.04 LTS
R 2.15.0 instalado via apt-get
Alguma sugestão?
-Mensagem original-
Ótimo,
mas uma vez sem palavras Benilton.
Obrigado,
Vitor
On Apr 13, 2012, at 3:12 PM, Benilton Carvalho wrote:
Eu vou apenas usar o primeiro marcador como exemplo (isso
representa o Minimo de Codigo Minimo Reproduzivel). Vou tbm assumir
o usual de teste qui-quadrado, no qual usa-se
Se eu quiser listar arquivos de outras extensões, como doc, ppt, ele também
faz? Eu vi que o list.files() lista tudo, mas se eu quiser só algumas extensões?
Marcos,
Se vc tem certeza de que cada RData contém somente um data.frame, acho que
isso vai funcionar:
arquivos - list.files()
veja a ajuda do comando, em particular detalhes sobre o argumento 'pattern'...
adicionalmente, ao fazer uma pergunta associada a outro problema, crie
um novo topico.
b
2012/4/11 Mauro Sznelwar sznel...@uol.com.br:
Se eu quiser listar arquivos de outras extensões, como doc, ppt, ele também
Caros colegas do r-br,
Alguém saberia me dizer se existe algum problema do local onde
esta hospedado o r-br@listas e o yahoo mail, não sei se alguém esta
passando por isso, mas estou recebendo as 22h mensagens postadas as 11h
e algumas mensagens minhas estão chegando a lista duplicadas.
43 matches
Mail list logo