Re: [R-br] Amostragem em um vetor

[FHRB Toledo via R-br]

Isso tem bastante cara de um exercício de classe. De qualquer forma, segue
uma versão funcionando que se aproxima muito da sua explicação:

> library(plyr)
> library(ggplot2)

> R <- 1E4 # tamanho da populacao
> S <- 1E3 # quantidade de amostras
> T <- 2:50 # tamanho das amostras
> names(T) <- T # nomeando vetor

> pop <- rnorm(R) # populacao

> ## cavalor de forca... : gera os dados
> sampler <- function(t, s, p) ldply(t, .id = 'size',
 function(ts) t(sapply(lapply(1:S,

function(ss) sample(x = p,


size = t,


replace = TRUE)),
   function(x) c(mu
= mean(x), sigma = sd(x)

> ## resume por tamanho de amostras
> dados <- mutate(ddply(sampler(T, S, pop), .(size),
summarize, muBar = mean(mu), sigmaBar =
mean(sigma)),
  size = as.numeric(as.character(size)))

> ## media das medias
> ggplot(dados, mapping = aes(x = size, y = muBar)) + geom_point()

> ## media dos desvios
> ggplot(dados, mapping = aes(x = size, y = sigmaBar)) + geom_point()

Vale a pena revisar e testar tudo novamente, sobre as saídas como matriz,
deixo isso para suas adaptações.

att,
FH

2017-09-27 19:24 GMT-05:00 Andre Oliveira via R-br :

> Vou tentar explicar melhor!
>
> " para cada  tam repetir  nam (vezes)"
>
>
> Seria repetir as amostragens:
>
> amostras de 2 em 2  um número k vezes. (Teria aqui k desvios-padrão e k
> medias). salvar media dos desvios e media das media!
>
> amostras de 3 em 3  um número k vezes. (Teria aqui k desvios-padrão e k
> medias). salvar media dos desvios e media das media!
>
>
> amostras de 4 em 4  um número k vezes. (Teria aqui k desvios-padrão e k
> medias). salvar media dos desvios e media das media!
>
> 
> 
> 
> 
>
>
> amostras de j em j  um número k vezes. (Teria aqui k desvios-padrão e k
> medias). salvar media dos desvios e media das media!
>
>
> Plotar os vetores salvos!
>
> 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 Quarta-feira, 27 de Setembro de 2017 15:02, Fernando Souza <
> nandodeso...@gmail.com> escreveu:
>
>
> Acredito que isso deva te ajudar. Não entendi muito bem o que quis dizer
> com  " para cada  tam repetir  nam (vezes)" , por isso o código abaixo
> não faz nenhuma repetição dos valores amostrados , porém acredito que é
> fácil de ajustar.
>
> populacao <- rchisq(10, df = 10)
> tamanho <- c(6,14,16,18,20) # diferentes tamanhos para amostrar
>
> amostra <- function(x,n){
> a <- list()
>
> for(i in 1:length(n)){
>
> a[[i]] <- sample(x,size = n[i])
>
>
> }
>return(a)
>
> }
> amostrados <- amostra(populacao,tamanho)
>
> #install.packages("plyr")
> library(plyr)
>
> resultado <- ldply (amostrados,function(x){
>media <- mean(x)
>desvio <- sd(x)
>
>resposta <- data.frame(media,desvio)
>return(resposta)
>
> })
>
> ##você pode converter o dataframe para matrix
> as.matrix(resultado)
>
>
> Em 26 de setembro de 2017 22:42, Andre Oliveira via R-br <
> r-br@listas.c3sl.ufpr.br> escreveu:
>
> Boa noite,
> preciso de amostrar em um vetor de nome populacao amostras de tamanho n,
> com n variando de tamanhos  [2; tam] e para cada  tam repetir  nam
> (vezes) e assim pegar a estimativas de media e dos desvios -padrão das
> amostras.
>
>
> Para cada amostra obtida extrair desvio padrão, média e guardar em uma
> matriz com vazia com NAs!  Ao final plotar a média dos desvios e a media
> das medias em função de cada n.
>
> Tentei fazer,  mas sem sucesso! Caso alguém pude ajudar ficarei grato!
>
>
> ## ##
> #
>
> n=10
> populacao <- rchisq(10, df = n)
> mu=n
> s=sqrt(2*n)
> mu
> s
>
> par(mfrow = c(1, 3))
> hist(populacao, prob=T, main="Y ~ Quiquadrado(GL=10)",xlab=" ",
> yla="Densidade",col=" limegreen")
>
> amostragem<-function( populacao, tam, nam)
>   {
>SD<-NULL
>media<-NULL
> X <-matrix(NA, tam, 2)
> {
> for(i in 2:nam)
> amostra<-sample(populacao,tam)
> SD[i]=sd(amostra)
> media[i]<-mean(amostra)
>X[i,] <- c(mean(SD), mean(media))
> }
> print(cbind(mean(SD),mean( media))
> print(X)
> plot.ts(SD, pch=19, type = "p", xlab="Número de folhas amostrada",
> ylab="Estimativa da Variabiliade", lty=2,xaxp=c(0,150,10))
> plot.ts(media, pch=19, type = "p", xlab="Número de folhas amostrada",
> ylab="Estimativa da média", lty=2,xaxp=c(0,150,10))
> }
>
>
> amostragem(populacao, 150, 5)
>  # ##
> ## ##
>
> obg
>
> __ _
> R-br mailing list
> 

Re: [R-br] reshape

[FHRB Toledo via R-br]

Antonio,

Seu subset não está definido corretamente... deve ser: variable %in%
c("ninfa", "pulpa").

att,
FH

2017-08-16 12:11 GMT-05:00 Antonio Moita via R-br 
:

> # Suponha o sequinte experimento em blocos ao acaso
> # com tamanho da parcela experimental de 4 plantas
> # no qual se conta o número pulpas ninfas e ovos.
>
> # GOSTARIA DE OBTER O CONJUNTO DE DADOS QUE FOSSE A SOMA DE PLANTAS POR
> CULTIVAR BLOCOS,
> # QUANDO USO TODO CONJUNTO DE DADOS FUNCIONA BEM
>
>
>
> require(reshape)
>
>
> x<-factor(LETTERS[1:3])
> cultivar<-rep(x, each=6)
> bloco<-rep(1:2,each=3)
> planta<-rep(1:3)
>
> ninfa<-c(0,0,0,1,1,1,0,0,0,1,2,3,0,1,2,0,0,1)
> pulpa<-c(1,0,1,1,0,1,0,0,0,0,0,0,0,2,1,0,0,4)
> ovos<-c(2,2,0,3,0,0,0,0,4,0,0,0,0,2,0,0,0,0)
>
> df<-data.frame(cultivar, bloco, planta, ninfa, pulpa, ovos)
>
>
> df.m<-melt(df, id=c("cultivar","bloco","planta"))
> df.m
>
> df.c<-cast(df.m, cultivar + bloco ~ variable,   sum,   na.rm=TRUE )
> df.c
>
> # MAS QUANDO SELECIONO SÓ UMA PARTE
> #   OS VALORES NÃO BATEM COM O QUADRO ANTERIOR, ALGUEM SABERIA EXPLICAR
> 
> df.c1<-cast(df.m, cultivar + bloco ~ variable, 
> subset=variable==c("ninfa","pulpa"),
> sum)
> df.c1
>
> # a SOMA DE ninfa e pulpa tambem esta errada
> df.c0<-cast(df.m, cultivar + bloco ~ ., subset=variable==c("ninfa","pulpa"),
> sum)
> df.c0
>
>
> Att
> Antonio W Moita
> Embrapa Hortaliças
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Salvar saídas de um loop

[FHRB Toledo via R-br]

Yuri,

Veja se isso resolve seu problema... :

## meses
meses <- c('setembro', 'outubro', 'novembro', 'dezembro', 'janeiro',
'fevereiro', 'marco', 'abril')
names(meses) <- meses # precisam ser nomeados

## funcao que le um certo aquivo
readFile <- function(mes) read.table(paste0(mes, '/final.txt'), header =
TRUE, sep = '\t')

## carrega todos arquivos em 1 so
data_raw <- ldply(meses, readFile, .id = 'mes')

## formato longo
data_long <- melt(data_raw, id.vars = .(mes, cidades), variable.name =
'ano', value.name = 'prod')

## altera niveis de ano
levels(data_long$ano) <- 1980:2012

## uma pilha com as matrizes por cidade
stack <- dlply(data_long, .(cidades), function(set) dcast(set, mes ~ ano,
value.var = 'prod'))

Quanto a salvar cada cidade fica como desafio para você :)

2017-08-09 15:10 GMT-05:00 Yury Duarte via R-br :

> Obrigado Fernando!
>
> Vou dar uma olhada nessas sugestões, obrigado!
> Mas meu interesse é em executar a tarefa através do comando for() apenas.
> Gostaria de entender porque meu código não é reproduzível, já que forneci
> ele na íntegra, juntamente com os arquivos necessários para alimenta-lo.
>
> Att
>
> Yury Duarte
> Engenheiro Agrônomo - ESALQ/USP
>
> Em 9 de agosto de 2017 16:57, Fernando Antonio de souza <
> nandodeso...@gmail.com> escreveu:
>
>> Ah! Seu código não é reproduzível. Envie um código reproduzível para que
>> possa receber ajuda.
>>
>> att
>>
>> =
>> Fernando Souza
>> Celular: (31)99796-8781 (Vivo)
>> E-mail:nandodeso...@gmail.com
>> 
>> ==
>>
>> On Ago 9 2017, at 4:48 pm, Yury Duarte via R-br 
>> wrote:
>>
>>> Boa tarde colegas listeiros!
>>>
>>> Estou tendo um pequeno problema para salvar as saídas de um loop que fiz
>>> para preencher uma matriz com dados de diferentes arquivos. Meus arquivos
>>> base possuem informações anuais de produtividade (colunas) de 15 locais
>>> (linhas) para um determinado mês. Construí o loop com intenção de gerar um
>>> arquivo por local, contemplando as produtividades de cada mês (nas linhas)
>>> e de cada ano (nas colunas). Conferi o output do único arquivo salvo pelo
>>> script e os valores correspondem ao arranjo desejado, então acredito que a
>>> lógica para o preenchimento da matriz esteja correta. Os arquivos para
>>> executar a simulação estão em anexo e o script está no corpo do email.
>>>
>>> Desde já, agradeço pela colaboração de todos!
>>>
>>> rm(list = ls())
>>>
>>> raiz = 'C:\\Users\\Yury\\Desktop\\Mestrado\\1_TESE\\Model_Data\\Mai
>>> ze\\PREVISAO\\PA\\'
>>>
>>> estrategia = c('E1', 'E2', 'E3', 'E4', 'E5')
>>> meses = c('setembro', 'outubro', 'novembro', 'dezembro', 'janeiro',
>>> 'fevereiro', 'marco', 'abril')
>>> cidade = c("abelardoluz", "altamira", "brasilia", "caarapo", "catalao",
>>> "lagarto", "machado", "mateiros",
>>>"muitoscapoes", "piracicaba", "pontagrossa",
>>> "primaveradoleste", "srmangab", "sdesiderio", "urucui")
>>> anos = c(1980:2012)
>>> previsao = matrix(NA, nrow = length(meses), ncol = length(1980:2012))
>>> colnames(previsao) = 1980:2012
>>>
>>> #for(e in 1:length(estrategias)){
>>>
>>>   for(m in 1:length(meses)){
>>>
>>> a = read.table(paste0(raiz, 'E1\\', meses[m],'\\', 'final.txt'),
>>> header = T, sep = '\t')
>>> for(city in 1:length(cidade)){
>>>
>>>   for(i in 1:length(anos)){
>>>
>>> previsao[m,i] = a[city,i]
>>>   }
>>> }
>>> write.table(previsao, paste0(raiz, 'E1\\', cidade[city],'.txt'),
>>> row.names = F, col.names = T, sep = '\t')
>>>   }
>>> #}
>>>
>>> Yury Duarte
>>> Engenheiro Agrônomo - ESALQ/USP
>>>
>>> ___
>>> R-br mailing 
>>> listr...@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
>>> m�nimo reproduz�vel.
>>>
>>>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Bug R in Debian 8,

[FHRB Toledo via R-br]

Conforme o post relacionado "Instalação da lme4", e dependendo do pacote
que queira instalar, você pode optar por instalar o pacte *.deb ao invés do
CRAN, para isso (considerando ser o nlme) faça no terminal:

$ sudo apt-get install r-cran-nlme

HTH

FH

2017-07-05 18:23 GMT-05:00 Fernando Souza via R-br :

> Caros colegas
>
> Possuo R instalado em uma máquina Debian 8. Recentemente ao utilizar o
> pacote nome recebi algumas mesgs estranhas de erro. Fui pesquisar e parece
> estar relacionado a falhas de compilação devido ao período de congelamento
> de debian8 para passaram para debian9. Já testei e o mesmo acontece com
> outros pacote s. O erro pare estar relacionado aos computadores c e fortran
>
> https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=861333
>
> Minha dúvida é. Como resolver isso? Consegui corrigir o nlme(não sei como)
> após reinstalar o r e após o nlme. Mas como saber quais são os pacotes
> problemáticos e como corrigir? As dicas do link acima não funcionaram para
> mim.
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Instalação da lme4

[FHRB Toledo via R-br]

Não sei sua distribuição qual é, mas se for alguma baseada em Debian (com
aptitude), instale pelo terminal o pacote *.deb ao invés do pacote CRAN, ou
seja:

$ sudo apt-get install r-cran-lme4

Pela mensagem de erro, o que te falta é a dependência do RcppEigen.

HTH

FH

2017-07-05 17:45 GMT-05:00 Andre Oliveira via R-br :

> Boa noite, estou usando o terminal Karina.Vinícius os avisos são estes!
>
> > install.packages("lme4")
> Installing package into ‘/usr/local/lib/R/site-library’
> (as ‘lib’ is unspecified)
> tentando a URL 'https://cran.cnr.berkeley.edu/src/contrib/lme4_1.1-13.
> tar.gz'
> Content type 'application/x-gzip' length 3624472 bytes (3.5 MB)
> ==
> downloaded 3.5 MB
>
> * installing *source* package ‘lme4’ ...
> ** package ‘lme4’ successfully unpacked and MD5 sums checked
> ** libs
> g++  -I/usr/share/R/include -DNDEBUG  
> -I"/usr/local/lib/R/site-library/Rcpp/include"
> -I"/usr/local/lib/R/site-library/RcppEigen/include"   -DNDEBUG
> -DEIGEN_DONT_VECTORIZE -fpic  -g -O2 -fstack-protector
> --param=ssp-buffer-size=4 -Wformat -Werror=format-security
> -D_FORTIFY_SOURCE=2 -g  -c external.cpp -o external.o
> In file included from external.cpp:8:0:
> predModule.h:12:23: fatal error: RcppEigen.h: Arquivo ou diretório não
> encontrado
>  #include 
>^
> compilation terminated.
> make: ** [external.o] Erro 1
> ERROR: compilation failed for package ‘lme4’
> * removing ‘/usr/local/lib/R/site-library/lme4’
> Warning in install.packages :
>   installation of package ‘lme4’ had non-zero exit status
>
> The downloaded source packages are in
> ‘/tmp/RtmpNtJzRo/downloaded_packages’
>
>
> 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 Quarta-feira, 5 de Julho de 2017 17:48, Vinícius Lionel Mateus via R-br
>  escreveu:
>
>
> André,
>
> Já tive problemas similares com o lme4 usando o openSUSE.
> No meu caso, havia dependências não atendidas.
>
> Existe mais informação antes do output que você postou?
>
> Abs,
> Vinícius
>
> On 07/05/2017 01:51 AM, Andre Oliveira via R-br wrote:
>
>  Bom dia,
> alguém está com mesmo problema na instalação da lme4?
>
> compilation terminated.
> make: ** [external.o] Erro 1
> ERROR: compilation failed for package ‘lme4’
> * removing ‘/usr/local/lib/R/site-library/lme4’
>
> The downloaded source packages are in
> ‘/tmp/RtmpRFLqvs/downloaded_packages’
> Warning message:
> In install.packages("lme4") :
>   installation of package ‘lme4’ had non-zero exit status
> > library(lme4)
> Error in library(lme4) : there is no package called ‘lme4’
>
>
> Grato
>
> 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
>
>
> ___
> R-br mailing 
> listr...@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
> m�nimo reproduz�vel.
>
>
> --
> Best regards,
>
> Vinícius Lionel Mateus, PhD (http://lattes.cnpq.br/6501001637020665)
> Atmospheric Chemistry Laboratory - Dep. Chemistry
> PUC-Rio - Pontifical Catholic University of Rio de Janeiro
> Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brazil CEP: 
> 22451-900
> Phone:   (+55) (21) 982 - 889 - 221
> Skype: vinicius.lionelhttp://www.qui.puc-rio.br/index.html
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a
> c�digo m�nimo reproduz�vel.
>
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Interface para usuário realizar chamada de função em R

[FHRB Toledo via R-br]

Infinitas possibilidades podem ser feitas... minha sugestão seria você
tentar algo como:

http://dirk.eddelbuettel.com/code/rinside.html

Att,
FH

2017-04-20 12:42 GMT-05:00 Paulo Nogueira Starzynski via R-br <
r-br@listas.c3sl.ufpr.br>:

> Olá colegas de lista.
> Gostaria da opinião de vocês com relação a algo que envolve R, mas não
> apenas ele.
>
> Desenvolvi em R um código que automatiza alguns processos. Todo o trabalho
> foi encapsulado dentro de uma função principal (aqui chamada de
> *realiza.consulta*), que recebe um parâmetro com as configurações
> informando o que deve ser feito.
>
> Criei algo parecido para exemplificar:
>
> #
> # "Painel de Controle" do Usuário
> #
> configuracoes <- list(
>
>   idUser <- 'usuario',
>   idPassword <- 'senha123',
>
>   consultas <- list(c("info1", "todos"),
> c("info2", "apenas alguns"),
> c("info3", NA),
> c("info4", "opcao1", "opcao2")
> )
> )
>
>
> #
> # Excucao do PROGRAMA PRINCIPAL
> source("/codigos.R")
> realiza.consulta(configuracoes)
>
>
> Como vocês podem ver, as configurações ficam armazenadas em um objeto do
> tipo lista, que é passado para a função principal. É de responsabilidade do
> usuário criar esse objeto e chamar a função realiza.consulta.
>
> A questão é que eu gostaria de desenvolver uma "máscara" para esse
> processo, onde os usuários informariam tais configurações e clicariam em um
> botão para iniciar o procedimento, sem ter de abrir o R ou RStudio.
>
> Vocês já fizeram algo semelhante para poder dividir o conhecimento?
>
> Abraços,
> Paulo
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] código para construir tabela com dados pareados

[FHRB Toledo via R-br]

help(xtabs)

2017-04-07 12:59 GMT-05:00 Mauricio Cardeal via R-br <
r-br@listas.c3sl.ufpr.br>:

> Boa tarde!
>
> Por favor, como faço para construir uma tabela como a do modelo do final
> do script,
>
> onde o conteúdo das 4 caselas são contagens do número de pares que combinam
>
> melhora e não melhora de 2 tratamentos aplicados simultaneamente no mesmo
>
> individuo (perna direita e perna esquerda), por exemplo ?
>
> # CMR:
>
> par <- c(1,1,2,2,3,3,4,4,5,5,6,6,7,7)  # refere-se ao número do individuo
> trata <-  c(1,0,1,0,1,0,1,0,1,0,1,0,1,0) # tratamento 1 (A) ou 0 (B)
> melhora <-  c(1,0,1,0,1,0,0,1,0,1,0,0,1,1) # resultado do tratamento 1
> (melhora) 0 (não melhora)
> tratac <- factor(trata)
> levels(tratac)<-list("A"=1,"B"=0 )
> melhorac <- factor(melhora)
> levels(melhorac)<-list("+"=1,"-"=0)
>
> data.frame(par,tratac,melhorac)
>
>par tratac melhorac
> 11  A+
> 21  B-
> 32  A+
> 42  B-
> 53  A+
> 63  B-
> 74  A-
> 84  B+
> 95  A-
> 10   5  B+
> 11   6  A-
> 12   6  B-
> 13   7  A+
> 14   7  B+
>
> # COMO FAÇO PARA GERAR ESSA TABELA A PARTIR DOS DADOS ACIMA?
>
> #  tratamento B
> #  melhora melhora
> # +   -
> #tratamento A melhora +   1   3
> # melhora -   2   1
>
> Obrigado,
>
> Maurício Cardeal
>
> UFBA
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo
> mnimo reproduzvel.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Ajuda: Append

[FHRB Toledo via R-br]

Talvez o que você pretenda fazer seja um merge(), veja o help.

2017-04-04 9:52 GMT-05:00 Leonardo Aguirre via R-br <
r-br@listas.c3sl.ufpr.br>:

> rbind.fill()
>
> Em 4 de abril de 2017 11:51, Edmar Caldas via R-br <
> r-br@listas.c3sl.ufpr.br> escreveu:
>
>> Pessoal estou com outra duvida. Estou tentando unir dois banco com a
>> função append. tenho arquivo A que contém a variável id, nome, lazer e
>> arquivo B que contém a variável id, nome, trabalho. quero unir os 02
>> arquivos e ficar com id, nome, lazer, trabalho mesmo que alguns casos
>> fiquem NULL.
>>
>> append <- rbind(LAZER_1,LAZER_2) dá erro de numero de colunas ou
>> argumentos não coincidem. estou tentando fazer isso:
>>
>>
>> LAZER_1$lazer <-NULL
>> LAZER_1
>>
>> LAZER_2$trabalho <-NULL
>> LAZER_2
>>
>> depois
>>
>> append <- rbind(LAZER_1,LAZER_2)
>>
>> continua com o mesmo erro?
>>
>> Edmar
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
>
> --
> *Leonardo Aguirre*
> +55 (61) 98133-8373
> leoaguirr...@gmail.com
> [image: Facebook]  [image: LinkedIn]
> 
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Função aggregate.

[FHRB Toledo via R-br]

Edmar,

Alternativa não testada:

> plyr::ddply(hospital, plyr::.(paciente), summarize,
   first  = min(data_int, na.rm = TRUE),
   last  = max(data_int, na.rm = TRUE),
   valor = mean(valor_int, na.rm = TRUE))

Espero ter ajudado.

att,
FH

2017-03-24 10:36 GMT-06:00 Edmar Caldas via R-br :

>
>
> Estou com um arquivo neste formato.
>
> [image: Alinhar imagem]
>
>
> Estou querendo agregar por paciente. Informando a média da coluna
> valor_int e quero também pegar a menor e a maior data da coluna
> data_int.Quero que o arquivo fique assim;
>
> [image: Alinhar imagem]
>
>
>
>
> Desenvolvi esse codigo.
>
> [image: Alinhar imagem]
>
>
> agora estou tentando colocar as datas.
>
> [image: Alinhar imagem]
>
>
>
>
>
>
>
> Edmar
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Salvar e ler de bando de dados

[FHRB Toledo via R-br]

Use tanto no write.table() quanto no read.table() a função paste0() para
gerar o nome variando conforme i:

> write.table(Y, paste0(i, '"txt"))

E para ler:

BD_i <- read.table(paste0(i, ".txt") header = TRUE)

Espero ter ajudado.

att,
FH

2017-03-23 5:56 GMT-06:00 Adriele Giaretta Biase via R-br <
r-br@listas.c3sl.ufpr.br>:

>
> Olá pessoal,
>
>
>
> eu estou trabalhando no R com conjunto de dados muito grande (numa fase de
> teste) e a tendência é que essa banco de dados aumente ainda mais, numa
> próxima etapa que não será teste. E em algumas de minhas programações, o
> R retornou mensagem de aviso dizendo que não podia alocar o vetor de
> tamanho de 3.7 GB.
>
> Dessa forma, ao invés de concatenar uma matriz em baixo da outra num
> objeto criado (onde ele me retorna que não foi capaz de armazenar tudo), eu
> gostaria de ir salvando cada matriz gerada de resultado no diretório local.
>
> Minha dúvida, é como criar esse algoritmo de forma que eu consiga salvar
> dentro de um loop, data.frames com nomes diferentes, sem que ocorra
> substituição de um data.frame pelo outro.
>
> Eu criei um código mínimo (de forma que vocês possam me dar alguma ajuda),
> e eu possa adaptar para a situação que estou trabalhando. Esse código está
> substituindo o data.frame a cada rodada dentro do loop. Eu queria que o
> nome do data.frame alterasse para que guardasse todos os resultados com
> nomes diferentes, acompanhando a variação do loop (1,2, ..., N):
>
>
>
> library(MASS)
>
> salva_BD_i <- function(N){
>
>   for (i in 1:N){
>
> p  <- 30# numero de variáveis a serem geradas (p)
>
> ME <- rep(1, p)
>
> rho<- 0.5 # correlacao
>
> sigma2 <- 1
>
> sigma  <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p))
>
> n <- 1
>
> Y  <- mvrnorm(n, ME, sigma)
>
> write.table(Y, "i.txt")
>
> }
>
> }
>
>
>
> N=3
>
> salva = salva_BD_i(N)
>
>
>
>
>
> Numa segunda etapa, eu precisaria de ler (puxar todos os data.frame para
> dentro de um objeto do R, num loop. Também criei um código mínimo
> executável, mas acontece o mesmo, o índice do loop não é processado dentro
> do read.table:
>
> > setwd('C:\\Users\\Adriele\\ teste')
>
>
>
> > dir()
>
>  [1] "1.txt"   "2.txt" "3.txt" "4.txt" "5.txt" "6.txt"
> "7.txt" "8.txt" "9.txt"   "10.txt" "Adriele.R"
>
>
>
> Leia_BD_i <- function(N){
>
>   aux= c()
>
>   for (i in 1:N){
>
> BD_i <- read.table("i.txt",  header = TRUE)
>
> media = apply(BD_i, 2, mean)
>
> aux = rbind(aux, media)
>
>   }
>
>   return(medias = aux)
>
> }
>
>
>
> N=10
>
> Leitura = Leia_BD_i(N)
>
>
> Error in file(file, "rt") : cannot open the connection
>
> In addition: Warning message:
>
> In file(file, "rt") : cannot open file 'i.txt': No such file or directory
>
> Agradeço se alguém puder contribuir,
>
> --
> Adriele Giaretta Biase.
> Mestre em  Estatística e Experimentação Agropecuária - UFLA.
> Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP
> Contato: (19) 98861-0619.
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Emacs

[FHRB Toledo via R-br]

https://ess.r-project.org/

2017-02-24 9:42 GMT-06:00 Alexandre Loures via R-br <
r-br@listas.c3sl.ufpr.br>:

> Boa tarde!
>
> Alguém saberia informar algum material para configurar o Emacs para
> trabalhar com os scripts, ou seja, "destacar" os códigos do R?
>
> Desde já muito obrigado!
>
> --
> *Alexandre Rodrigues Loures*
> Doutorando em Economia Aplicada
> Universidade Federal da Paraíba - UFPB
> Centro de Ciências Sociais Aplicadas - CCSA
> Programa de Pós-Graduação em Economia - PPGE
> Site: www.ccsa.ufpb.br/ppge
> [image: orcid] www.orcid.org/-0002-1288-0135
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Ajuda para incluir valores

[FHRB Toledo via R-br]

Rafael,

> help(merge)

Espero ter ajudado.

att,
FH

2016-12-06 13:06 GMT-06:00 Rafael Garcia Cunha via R-br <
r-br@listas.c3sl.ufpr.br>:

> Boa tarde,
>
> eu tenho o seguinte problema. Eu preciso juntar dois conjuntos de dados
> que possuem 2 variáveis em comum, "uor" e "mes". Só que 1 variável não
> possui valores para determinado mês. E o que eu queria era que para os
> meses que não possuem informação, essa determinada variável, assumisse o
> valor 0 e fosse incluídos esses meses faltantes. Abaixo eu envio um código
> para exemplificar.
> A ideia é que seja incluída o valor 0 na variável x para os meses 4, 5 e 6
> da uor "B".
>
> ## exemplo de codigo ##
>
> library(data.table)
> uor <- c(rep("A", 8), rep("B", 5), rep("C", 8), rep("D", 8))
> uor2 <- c(rep("A", 8), rep("B", 8), rep("C", 8), rep("D", 8))
> mes1 <- c(1:8, 1,2,3,7,8, 1:8, 1:8)
> mes2 <- rep(1:8, 4)
>
> DT <- data.table(uor, mes1, x = rnorm(length(uor)))
> DT2 <- data.table(uor2, mes2, x1 = rnorm(length(uor2)), x2 =
> rpois(length(uor2), 4))
>
> ##
>
> Atenciosamente,
>
> --
> Rafael Garcia Cunha
>
> YNWA
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Loop script

[FHRB Toledo via R-br]

Fernando,

> long_shape <- reshape2::melt(dados, id.vars = 'Trat')
> model <- function(df) lm(value ~ Trat, data = df)
> fits <- plyr::ldply(long_shape, .(variable), model)

Deste modo, "fits" é uma lista que armazena os ajustes para todas as
variáveis. Minha sugestão e usar outro **ply para tratar tais ajustes e
obter os testes de média e demais analíses necessárias.

Att,
FH



2016-12-06 13:29 GMT-06:00 Fernando Rodrigo Bortolozo via R-br <
r-br@listas.c3sl.ufpr.br>:

> Eu estou fazendo anova e com posterior teste de médias de um experimento
> com muitas variáveis. Então eu tenho que repetir o script muitas vezes o
> mesmo script para cada variável.
>
> Os colegas sabem me dizer se existe alguma maneira de fazer um Loop para
> um script, ou mesmo criar um elemento com todas variáveis para analisar de
> uma só vez?
>
> Aqui vai uma pequena parte das variáveis que vou analisar, tenho 7
> tratamentos, 4 repetições.
>
> [image: pasted1]
>
> Aqui esta o script (Agricolae Package)
>
> dados <- read.csv("~/Documents/R/Data306090DAPv2.csv", header = TRUE)
> head(dados)
> model<-aov(Cas30 ~Trat, data=dados)
> out <- HSD.test(model,"Trat", group=TRUE,console=TRUE, main="Calcio")
> par(mfrow=c(1,2))
> out$means
> par (mar = c (3,3,2,0),cex=0.9)
> bar.group(out$groups,ylim=c(0, 360),density=50,border="blue", col="blue",
> las=1)
> bar.err(out$means,variation="SE",horiz=FALSE, ylim = c(0, 360),
> col=colors()[15],space=0.3,bar=TRUE, las=1)
> out<-HSD.test(model,"Trat", group=FALSE)
> means<-out$means
> df<-df.residual(model)
> MSerror<-deviance(model)/df
> with(dados,HSD.test(Cas30,Trat,df,MSerror, group=TRUE,console=TRUE, main=
> "Cálcio"))
>
>
> Muito obrigado desde já.
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Inputando registros

[FHRB Toledo via R-br]

banco$var3 <- with(banco, ifelse(is.na(var2), var1, var2))

2016-06-30 13:08 GMT-05:00 Mauricio Cardeal via R-br <
r-br@listas.c3sl.ufpr.br>:

> banco$var3 [is.na(banco$var2)] <-banco$var1
> Em 30 de jun de 2016 13:55, "Wagner Tassinari via R-br" <
> r-br@listas.c3sl.ufpr.br> escreveu:
>
>> Olá pessoal, tenho o seguinte banco de dados:
>>
>> var1 var2
>> 20 25
>> 40
>> 40 45
>> 40
>> 50 55
>> 50
>> 60 60
>>
>>
>> Quero criar uma var3 que seja igual a var2 mas onde está em branco impute
>> os valores de var1, ex:
>>
>> var3
>> 25
>> 40
>> 45
>> 40
>> 55
>> 50
>> 60
>>
>> Estou fazendo da seguinte forma:
>>
>> banco$var1=as.character(banco$var1)
>> banco$var2=as.character(banco$var2)
>>
>> banco$var3= banco$var2
>> banco$var3[banco$var2 == " "] = banco$var1
>>
>> Mas não esta dando certo !!
>>
>> Obrigado pessoal
>>
>> -
>> Wagner S. Tassinari
>> Departamento de Matemática
>> Universidade Federal Rural do Rio de Janeiro.
>> BR-465, Km 7 - Seropedica, RJ - Brasil
>> CEP: 23890-000
>> Cel: (21) 96488-5982 (WhatsApp)
>> Skype: wagner.tassinari
>> wtassin...@gmail.com
>> --
>>
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Teste de diferça entre proporções via loop

[FHRB Toledo]

Meu centavo:

> dados <- read.table(file = 'data_test.txt', header = TRUE, sep= '\t')
> dados <- plyr::ddply(.data = dados, .variables = .(AREA), .fun = mutate,
   pvalor = prop.test(x = c(X1, X2), n = c(N, N),
correct = FALSE)$p.value)

att,
FH

2016-03-23 14:33 GMT-06:00 Cesar Rabak :

> Como sugestão para um código mais enxuto, o objetivo poderia ser atingido
> com a aplicação da função apply:
>
> > dados <- read.table("data_test.txt", header=T, sep= "\t")
> > dados$p_valor <- apply(dados[,2:4], 1, function(linha)
> prop.test(x=c(linha[1],linha[2]),n=c(linha[3],linha[3]), correct=F)$p.value)
> > dados
>   AREA   X1   X2 N   p_valor
> 11 1982 9416 11398  0,00e+00
> 22 6546 4084 10630 5,362948e-250
> 33 3418  667 40850  0,00e+00
> 44 1503 4408 59110  0,00e+00
> 55  966  339 30500  6,472911e-69
>
> Ademais, se o OP tiver que resolver problemas dessa natureza com
> frequência, pode ser que valha a pena examinar alguns pacotes que têm
> funções para facilitar esse tipo de resumos, como, por exemplo, os pacotes
> *descr* e *reporttools*.
>
> HTH
> --
> Cesar Rabak
>
>
> On Wed, Mar 23, 2016 at 1:46 PM, Daniel Tiezzi  wrote:
>
>> Boa tarde Wagner
>>
>> Supondo que seus dados esteja uma tabela data_test (em anexo)
>>
>> acho que você pode resolver assim
>>
>>
>> data<- read.table(“data_test.txt”, header=T, sep= “\t”)
>>
>> data_df <- as.data.frame(data)
>>
>> res<- list()
>>
>> for (i in 1:nrow(data)) {
>>   y<- prop.test(x = c(data$X1[i], data$X2[i]), n = c(data$N[i],
>> data$N[i]), correct = FALSE)
>>   x= y$p.value
>>   res[[paste0(i)]] <- x
>> }
>>
>> result <- t(as.data.frame(res))
>>
>> data_df1 <- cbind(data_df, result)
>>
>> data_df1
>>
>>
>> Espero que funcione
>>
>>
>> daniel
>>
>>
>>
>> Daniel Tiezzi, MD, PhD
>> Professor Associado
>> Departamento de Ginecologia e Obstetrícia
>> Setor de Mastologia e Oncologia Ginecológica
>> Faculdade de Medicina de Ribeirão Preto - USP
>> Tel.: 16 3602-2488
>> e-mail: dtie...@fmrp.usp.br
>>
>>
>>
>>
>>
>> On Mar 23, 2016, at 3:50 PM, Wagner Tassinari 
>> wrote:
>>
>> AREAX1  X2   N   p_valor
>> 11982941611398  ?
>> 26546408410630  ?
>> 33418667   4085   ?
>> 415034408  5911   ?
>> 5966339   305   ?
>>
>>
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Problemas com criação de listas

[FHRB Toledo]

Veja se isso te ajuda:

> lista <- doBy::splitBy(formula = ~ V2, data = arq_1)

att,
FH

On 23 October 2015 at 13:31, Paulo Henrique Pimenta <
paulopimen...@hotmail.com> wrote:

> Olá Listeiros,
>
> Estou com um problema na criação de novas listas no R. Abaixo segue parte
> do meu código:
>
> arq_1 <-
> read.table("C:/Users/Paulo/Documents/Alturas_Ilha_Fiscal-2015-10-22/Alturas_Ilha_Fiscal/50140002560102196123121961ALT_MOD.xls")
>
> 
> Fazendo a mao mes a mes 
>
> jan <- arq_1[which(arq_1$V2==1),]
> fev <- arq_1[which(arq_1$V2==2),]
> mar <- arq_1[which(arq_1$V2==3),]
> abril <- arq_1[which(arq_1$V2==4),]
> maio <- arq_1[which(arq_1$V2==5),]
> junho <- arq_1[which(arq_1$V2==6),]
> julho <- arq_1[which(arq_1$V2==7),]
> agosto <- arq_1[which(arq_1$V2==8),]
> setembro <- arq_1[which(arq_1$V2==9),]
> outubro <- arq_1[which(arq_1$V2==10),]
> novembro <- arq_1[which(arq_1$V2==11),]
> dezembro <- arq_1[which(arq_1$V2==12),]
>
> automatizando#
>
> for (i in 1:12){
>
> mes[i] <- arq_1[which(arq_1$V2==i),]
>
>}
>
> Eu gostaria de poder automatizar o primeiro processo, a captação de um
> bloco da lista arq_1 para a variável mes[i], porém eu recebo a mensagem:
>
> There were 12 warnings (use warnings() to see them)
>
> e ao abrir os warnings:
>
> Warning messages:
> 1: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 2: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 3: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 4: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 5: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 6: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 7: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 8: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 9: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 10: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 11: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 12: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
>
> Fazendo pelo processo manual dá certo, mas pq fazendo pelo
> processo “automatizado" não é possível. Peço ajuda de vocês para encontrar
> uma saída!
>
> Um grande abraço a todos
>
> Paulo.
>
>
>
> Sent from Windows Mail
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Problemas com criação de listas

[FHRB Toledo]

Nesse caso,

Uma "automatização" bem naif seria:

lista <- lapply(1:12, function(mes) subset(arq_1, V2 == mes))

P.S.: Considere solicitar ao administrador da sua máquina atualizar sua
versão, bem como instalar pacotes úteis!

att,
FH

On 23 October 2015 at 16:19, Paulo Henrique Pimenta <
paulopimen...@hotmail.com> wrote:

> Olá FHRB,
>
> A solução  lista <- doBy::splitBy(formula = ~ V2, data = arq_1) parece
> interessante, no entanto o computador que eu uso não possui para mim acesso
> para realizar atualizações, assim o pacote doBy não pode ser usado, ele não
> esta disponível para a versão 2.14.1 do R
>
> Abraços.
>
> Paulo Henrique de A. S. Pimenta.
>
> Graduando em Meteorologia (Bacharelado) - IAG/USP.
> Fone: +5511981318435.
>
>
> --
> From: paulopimen...@hotmail.com
> To: r-br@listas.c3sl.ufpr.br
> Subject: Problemas com criação de listas
> Date: Fri, 23 Oct 2015 18:31:56 +
>
>
> Olá Listeiros,
>
> Estou com um problema na criação de novas listas no R. Abaixo segue parte
> do meu código:
>
> arq_1 <-
> read.table("C:/Users/Paulo/Documents/Alturas_Ilha_Fiscal-2015-10-22/Alturas_Ilha_Fiscal/50140002560102196123121961ALT_MOD.xls")
>
> 
> Fazendo a mao mes a mes 
>
> jan <- arq_1[which(arq_1$V2==1),]
> fev <- arq_1[which(arq_1$V2==2),]
> mar <- arq_1[which(arq_1$V2==3),]
> abril <- arq_1[which(arq_1$V2==4),]
> maio <- arq_1[which(arq_1$V2==5),]
> junho <- arq_1[which(arq_1$V2==6),]
> julho <- arq_1[which(arq_1$V2==7),]
> agosto <- arq_1[which(arq_1$V2==8),]
> setembro <- arq_1[which(arq_1$V2==9),]
> outubro <- arq_1[which(arq_1$V2==10),]
> novembro <- arq_1[which(arq_1$V2==11),]
> dezembro <- arq_1[which(arq_1$V2==12),]
>
> automatizando#
>
> for (i in 1:12){
>
> mes[i] <- arq_1[which(arq_1$V2==i),]
>
>}
>
> Eu gostaria de poder automatizar o primeiro processo, a captação de um
> bloco da lista arq_1 para a variável mes[i], porém eu recebo a mensagem:
>
> There were 12 warnings (use warnings() to see them)
>
> e ao abrir os warnings:
>
> Warning messages:
> 1: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 2: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 3: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 4: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 5: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 6: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 7: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 8: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 9: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 10: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 11: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
> 12: In mes[i] <- arq_1[which(arq_1$V2 == i), ] :
>   number of items to replace is not a multiple of replacement length
>
> Fazendo pelo processo manual dá certo, mas pq fazendo pelo
> processo “automatizado" não é possível. Peço ajuda de vocês para encontrar
> uma saída!
>
> Um grande abraço a todos
>
> Paulo.
>
>
>
> Sent from Windows Mail
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] RES: Uso da função sort()

[FHRB Toledo]

Michelle,

Minha sugestão seria:

plyr::arrange(df, data)

Em que df é um data frame e data é uma coluna do mesmo... ao colocar -data
a ordem é inversa... !

att,
FH

2015-09-02 7:02 GMT-05:00 Rodrigo Coster :

> Para não ter que mudar o formato da data no arquivo original, o melhor é
> declarar ela como data (através do as.Date()):
>
> datas <- c("10/11/2011", "17/03/2010", "15/05/2011")
> valor <- c(34, 27, 54)
> df <- data.frame(data = as.Date(datas, '%d/%m/%Y'), valor)
>
> df[order(df$data), ]
>
> 2015-09-02 1:14 GMT-03:00 Diego Miro :
>
>> Apenas complementando a solução do Jobenil, observe que se a coluna data
>> tem a forma dd/mm/, cria-se um problema na hora de ordenar.
>> O ideal seria ter a data em um formato mmdd.
>>
>> Em 1 de setembro de 2015 22:18, Jobenil - Gmail 
>> escreveu:
>>
>>> datas <- c("10/11/2011", "17/03/2010", "15/05/2011")
>>>
>>> valor <- c(34, 27, 54)
>>>
>>> df <- data.frame(data, valor)
>>>
>>>
>>>
>>> df[order(df$datas),]
>>>
>>> datasvalor
>>>
>>> 1 10/11/201134
>>>
>>> 3 15/05/201154
>>>
>>> 2 17/03/201027
>>>
>>> > df[order(df$valor),]
>>>
>>> Datasvalor
>>>
>>> 2 17/03/201027
>>>
>>> 1 10/11/201134
>>>
>>> 3 15/05/201154
>>>
>>>
>>>
>>> Dê um ?order
>>>
>>>
>>>
>>>
>>>
>>> *De:* R-br [mailto:r-br-boun...@listas.c3sl.ufpr.br] *Em nome de *Michelle
>>> Bau Graczyk
>>> *Enviada em:* terça-feira, 1 de setembro de 2015 21:28
>>> *Para:* R-br@listas.c3sl.ufpr.br
>>> *Assunto:* [R-br] Uso da função sort()
>>>
>>>
>>>
>>> Boa noite,
>>>
>>>
>>>
>>> Eu preciso ordenar uma tabela de dados na forma crescente de datas. Na
>>> coluna 1 eu tenho as datas e nas outras 400 colunas eu tenho números
>>> correspondente a valores de venda de um certo produto.
>>>
>>> O problema é que quando eu aplico o sort() na coluna de datas, ele só
>>> ordena a mesma. Para ser mais clara, ela não leva junto a linha de dados
>>> com as outros 400 valores para a nova posição.
>>>
>>>
>>>
>>> Por favor, Vocês sabem me dizer como eu posso arrumar isto?
>>>
>>>
>>>
>>> Muito obrigada,
>>>
>>>
>>>
>>> Michelle
>>>
>>> ___
>>> R-br mailing list
>>> R-br@listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Loop para mudar nomes de colunas em data.frames

[FHRB Toledo]

André,

Minha sugestão não exatamente sua resposta, mas pode ajudar:

data(iris)

list.data - lapply(X = list('1' = iris, '2' = iris, '3' = iris),
FUN = function(data) with(data, tapply(Sepal.Length,
Species, length)))

Att,
FH

2015-08-20 15:09 GMT-05:00 André Vieira andre...@gmail.com:

 Pessoal,

 Quero rodar um loop que, em um dos comandos, adicione números em sequência
 (e.g. 1, 2, 3) apenas à segunda coluna de várias bases de dados. No meu
 caso, imaginei que precisaria usar um loop dentro de outro, mas não
 funcionou. (Ainda estou começando no R, como vocês podem perceber).

 Segue um exemplo mais simples com a base iris:

 data1 - iris
 data2 - iris
 data3 - iris

 list_data - list(data1, data2, data3)
 list_emp - vector(list, length(list_data))

 for(i in 1:length(list_data)){
 file - list_data[[i]]
 Sepal.Length - tapply(file$Sepal.Length, file$Species, length)
 newData - as.data.frame(Sepal.Length)
 newData - cbind(Species = rownames(newData), newData)
 rownames(newData) - NULL
 list_emp[[i]] - newData
 for(j in 1:3){
 colnames(list_emp[[i]][2]) - paste(Sepal.Length, j, sep = _)
 # Nao funciona
 }
 rm(Sepal.Length, file, i, newData)
 }

 O loop deveria retornar uma lista com data frames na forma:

 # Species| Sepal.Length_1
 # setosa |  50
 # versicolor |  50
 # virginia   |  50

 # Species| Sepal.Length_2
 # setosa |  50
 # versicolor |  50
 # virginia   |  50

 # Species| Sepal.Length_3
 # setosa |  50
 # versicolor |  50
 # virginia   |  50

 Alguém pode me ajudar a corrigir a linha que não funciona? Além disso,
 aceito sugestões que não sejam com for loops.

 Abs.,

 André

 --
 *André Vieira*
 Mestre, Sociologia, UFMG (M.Sc., Sociology, UFMG)
 Bacharel, Sociologia, UnB (B.Sc., Sociology, UnB)
 Licenciado, Sociologia, UnB (B.Ed., Sociology, UnB)
 GitHub http://github.com/ahpvieira
 Currículo Lattes http://lattes.cnpq.br/6721713322923268
 Skype: andre.hpv

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Geração de números aleatórios com amplitude determinada

[FHRB Toledo]

Não seria o caso de você obter esses valores analiticamente, ao invés de
usar simulação? Veja a documentação da função pnorm().

att,
FH


2014-07-28 21:14 GMT-05:00 Marcos Bissoli mbiss...@gmail.com:

 Boa noite a todos,

 Preciso gerar uma sequência de números aleatórios, de uma distribuição
 normal. No entanto, a função rnorm não permite definir limites mínimo e
 máximo para o vetor de dados. Assim, criei uma função para tal fim:

 Gerador.Aleatorio - function(n,menor,maior,x,s)
 +{
 +  dados - rnorm(n,x,s)
 +  while ((min(dados)menor)|(max(dados)maior))
 +  {
 +dados - rnorm(n,x,s)
 +  }
 +  dados - dados
 +  print(dados)
 +}

 Onde:
 - n é o número de valores que quero gerar
 - menor é o limite mínimo do vetor
 - maior é o limite máximo do vetor
 - x é a média do vetor
 - s é o desvio padrão do vetor

 O problema é que, como usei um laço (loop), o custo de tempo
 computacional dessa função tornou-a simplesmente impraticável. Até 100
 valores eu consegui gerar. Mas agora tentei com 284 e já jantei, fumei,
 tratei de assuntos familiares e... Nada! E minha demanda, na verdade, é
 para gerar mais de 1500 valores, simulando estudo epidemiológico
 observacional, a partir de valores já descritos em um estudo piloto que
 realizei (ou seja, eu tenho já estimativas de n, mínimo, máximo, média e
 desvio padrão).

 Há alguma função, pacote ou sintaxe que torne isso possível? E qual seriam
 esses limites do custo computacional?

 Grato desde já por qualquer ajuda.

 Abraços fraternos,

 --
 MARCOS BISSOLI
 Faculdade de Nutrição
 Universidade Federal de Alfenas

 E-mail: mbiss...@gmail.com
 Twitter: #mbissoli
 Facebook: https://www.facebook.com/MarcosBissoli

 Alfenas, Minas Gerais, Brasil


 *Pense na Natureza antes de Imprimir*
 Divulgue ON-LINE

 Eu apoio a ENEN na luta por um Brasil sem fome

 por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan
 (para cada povo sua própria língua, para todos os povos o Esperanto)

 E nunca votarei no PSDB/DEM!

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Selecionar vários níveis dentro de um fator com for()

[FHRB Toledo]

x - paste('operario', 1:20, sep = '')

dados_caSUB - subset(dados_ca, dados_ca[,3] %in% x)


2014-07-16 16:29 GMT-05:00 ASANTOS alexandresanto...@yahoo.com.br:

 Boa noite pessoal,

 Tenho um fator (indiv) no objeto dados_cal[,3], composto de vários
 níveis sendo soldado1,  ate soldado 20 e operario1, ... ate operario20,
 gostaria de selecionar só os operário, para isto eu fiz:

 #Operários
 x- c(operario1,operario2,operario3,operario4,
 operario5,operario6,
 operario7 ,operario8,operario9,operario10,operario11,
 operario12,
 operario13,operario 14,operario15,operario16,
 operario17,operario18,
 operario19,operario20)

 RES=NULL
   for(l in 1:length(x)){# Rayon des cercles

 teste-dados_cal[dados_cal[,3]==x[l]]
 RES=rbind(RES,teste)
 }
 #

   Mas o meu for() aplicado ao objeto x não funcionou, 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:alexandresanto...@yahoo.com.br
 alexandre.san...@cas.ifmt.edu.br
 Lattes: http://lattes.cnpq.br/1360403201088680
 ==

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] função BiCopKPlot

[FHRB Toledo]

Consulte o arquivo fonte:

CDVine/R/BiCopChiPlot.R

Na linha 154 é definida a função H() e na linha 170 a função W()!

att,
FH


On 24 March 2014 16:36, Marcus Nunes marcus.nu...@gmail.com wrote:

 Natália,

 tenta rodar

 library(CDVine)
 CDVine:::H
 CDVine:::W

 Att,


 On Mon, Mar 24, 2014 at 4:30 PM, Natalia Martins nsmbarr...@gmail.comwrote:

 Prezados  boa tarde,
 estou estudando a função BiCopKPlot utilizada para a modelagem de
 cópulas.
 Nesta função tem-se:

 function (u1, u2, PLOT = TRUE, ...)
 {
 if (is.null(u1) == TRUE || is.null(u2) == TRUE)
 stop(u1 and/or u2 are not set or have length zero.)
 if (length(u1) != length(u2))
 stop(Lengths of 'u1' and 'u2' do not match.)
 if (length(u1)  2)
 stop(Number of observations has to be at least 2.)
 if (PLOT != TRUE  PLOT != FALSE)
 stop(The parameter 'PLOT' has to be set to 'TRUE' or 'FALSE'.)
 Wi - W(u1, u2)
 Hi - H(u1, u2)
 Hi.sort - sort(Hi)
 n - length(u1)
 W.in - rep(NA, n)
 for (i in 1:n) {
 f = function(w) {
 w * (-log(w)) * (w - w * log(w))^(i - 1) * (1 - w +
 w * log(w))^(n - i)
 }
 W.in[i] - n * choose(n - 1, i - 1) * (integrate(f, lower = 0,
 upper = 1)$value)
 }
 g - function(w) {
 w - w * log(w)
 }
 if (PLOT) {
 plot(g, xlim = c(0, 1), ylim = c(0, 1), pch = x, xlab =
 expression(W[1:n]),
 ylab = H, ...)
 points(W.in, Hi.sort, pch = x, cex = 0.4, ...)
 abline(a = 0, b = 1)
 }
 else {
 kendall.plot.output - list(W.in, Hi.sort)
 names(kendall.plot.output) - c(W.in, Hi.sort)
 return(kendall.plot.output)
 }
 }

 Porem, não encontrei no pdf da library(CDVine) como é definida a função H
 e W (em destaque).
 Por favor se alguem souber e puder me ajudar.

 Muito obrigada!

 --

 Natália Martins
 Contato: (35) 91812482

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Marcus Nunes
 http://marcusnunes.me/

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Pacote synbreed

[FHRB Toledo]

Giselle,

Note que tem dois pacotes...

synbreed e synbreedData!

As funções estão no synbreed enquanto que synbreedData contém apenas dados!

Leia:

http://cran.r-project.org/web/packages/synbreed/ e
http://cran.r-project.org/web/packages/synbreedData/index.html


2014-01-30 Giselle Davi giselle_d...@yahoo.com.br:

 Prezados,

 alguém saberia me dizer o que aconteceu com as funções do pacote synbreed
 ???
 Grande parte delas não podem ser encontradas quando se roda a rotina
 conforme o manual.

 Abaixo encontra-se o início da rotina descrita no manual do pacote
 synbreed, onde já existe o erro de função não encontrada.

 Desde já agradeço.


 library(synbreedData)
 data (maize)
 newDHpheno - data.frame(Trait=1000,row.names=newDH)
 newDHpheno

 #simulating genotypic data
 newDHgeno - matrix(sample(c(0,1), ncol(maize$geno), replace=TRUE), nrow=1)
 rownames(newDHgeno) - newDH

 #new pedigree
 newDHpedigree - data.frame(ID=newDH, Par1=0, Par2=0, gener=0)

 #new covar information
 newDHcovar - data.frame(family=NA, DH=1, tbv=1000, row.names=newDH)

 # add individual
 maize2 -
 add.individuals(maize,newDHpheno,newDHgeno,newDHpedigree,newDHcovar)
 summary(maize2)
 ##erro em add.individuals: função não pode ser encontrada


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Como colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R

[FHRB Toledo]

expression (Lesão (*mm^2*))


2014-01-30 Ederson Civardi civardi@gmail.com:

 Bom dia!
 Gostaria de uma ajuda básica para a formatação de gráficos no R.
 Preciso colocar caracteres, letras, números em subscrito ou sobrescrito em
 gráfico no R, porém não estou conseguindo fazer isso.
 Tentei vários comandos, porém sempre aparece como por ex. (M^2).
  O comandos para o gráfico em 3D é:
 persp(xr,yr,z, main = Banco, theta = 30, phi = 15, expand = 0.7, col =
 darkorchid4  , xlab=Temperatura (ºC), ylab=PMF (h), zlab= expression
 (Lesão (mm^2)),scale=T,ticktype = detailed)
 Outra dificuldade é inserir algumas linhas nos eixos x, y, z no gráfico
 que correspondam as escalas de cada eixo, para facilitar a visualização
 para que for ver.
 Alguém pode me ajudar?
 Desde já agradeço.
 --
 *_*
 *Ederson A. Civardi*
 *Doutorando Produção Vegetal - UFG - PPGA - Goiânia*
 Cel. (Claro): (62) 9237-7817
 

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Como colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R

[FHRB Toledo]

Edson, veja exatamente o comando que você postou!

Lembre-se que textos entre aspas são strings e não são avaliados pela
função expression()...

expression(Lesão ( * mm^2 * ) )

Veja: Lesão ( é uma string... * mm^2 * não é uma string... e a aspas
final deve estar entre aspas tb )!

Espero que ajude!

att,
FH


2014-01-30 Ederson Civardi civardi@gmail.com:


 Obrigado Fernando Toledo por responder. Porém o comando que me passou
 ainda não funciona.
 Da erro.
 Erro: símbolo inesperado in:
 persp(xr,yr,z, main = Banco, theta = 30, phi = 15, expand = 0.7, col =
 darkorchid4  , xlab=Temperatura (ºC), ylab=PMF (h), zlab= expression
 (Lesão (*mm^2*),scale=T,ticktype = detai
 persp

 --
 *_*
 *Ederson A. Civardi*
 *Doutorando Produção Vegetal - UFG - PPGA - Goiânia*
 Cel. (Claro): (62) 9237-7817
 

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Parcela subdividida: problema na declaração do erro associado a parcela

[FHRB Toledo]

Cláudio... isso não seria um warnings!

É apenas para avisar que o a solução do sistema está vindo por uma inversa
generalizada pois sua matriz de incidência não tem posto completo!


2013/12/4 Jose Claudio Faria joseclaudio.fa...@gmail.com

 Verdade!

 Como deu erro eu não tentei rodar o summary.

 Muito obrigado Tiago!

 2013/12/4 Tiago Souza Marçal tiagosouzamar...@hotmail.com:
  Mesmo com o aviso use o comando summary(m0) e o quadro será formado.
 
  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
 
  #
 
 
  Date: Wed, 4 Dec 2013 18:20:19 -0200
  From: joseclaudio.fa...@gmail.com
  To: r-br@listas.c3sl.ufpr.br
  Subject: [R-br] Parcela subdividida: problema na declaração do erro
  associado a parcela
 
 
  Pessoal,
 
  O CMR abaixo é da apostila do Walmes.
 
 
 
 #--
  # dados
  ps - expand.grid(BL=c(I, II, III, IV),
  ES=c(e1, e2),
  AD=c(a1, a2, a3))
  ps$alt -
 
 c(58,77,38,52,44,59,30,34,85,90,73,77,59,68,45,55,66,93,67,64,54,75,53,48)
  str(ps)
  #
 
 
 #--
  # análise de variância (erro A = BL x AD)
  m0 - aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data=ps) #
  termo Error para declarar erro A
  m0 - aov(alt ~ BL + AD * ES + Error(BL:AD), data=ps) # o mesmo com
  fórmula comprimida
 
  summary(m0)
 
 
 #--
 
  Encontrei muitos outros exemplos (na web e em livros) com essa
  notação: Error(BL:AD)
 
  Contudo o R (3.0.2patched) não está aceitando essa notação!
  Fornece sempre a mensagem de erro:
  In aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data = ps) :
  modelo Error() é singular
 
  Contudo o R está aceitando esse modelo:
  m0 - aov(alt ~ BL + AD + Error(BL/AD) + ES + AD:ES, data=ps)
 
  Mas nesse caso o resíduo A (associado a bloco e parcela - AD) não fica
  correto.
 
  Alguém que anda trabalhando com esses modelos poderia ajudar?
 
  Ab,
  --
  ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
  Jose Claudio Faria
  Estatistica
  UESC/DCET/Brasil
  joseclaudio.faria at gmail.com
  Telefones:
  55(73)3680.5545 - UESC
  55(73)9100.7351 - TIM
  55(73)8817.6159 - OI
  ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
  código mínimo reproduzível.
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.



 --
 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
 Jose Claudio Faria
 Estatistica
 UESC/DCET/Brasil
 joseclaudio.faria at gmail.com
 Telefones:
 55(73)3680.5545 - UESC
 55(73)9100.7351 - TIM
 55(73)8817.6159 - OI
 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Lattice segplot + panel.xblocks

[FHRB Toledo]

Valeu Walmes... na médida!

Obrigado!


On 2 December 2013 07:51, walmes . walmeszevi...@gmail.com wrote:

 Tenta assim Fernando,

 segplot(reorder(tratamento, media)~inferior+superior|fator, data=dados,
 centers=media, draw.bands=FALSE, segments.fun=panel.arrows,
 ends='both', angle=90, length=1, unit='mm',
 panel=function(x, y, z, centers, subscripts, ...) {
 wp - which.packet()
 panel.rect(dadosAUX[wp,xleft], min(as.numeric(z))-1,
dadosAUX[wp,xright], max(as.numeric(z))+1,
fill=rgb(0.5,0.5,0.5,0.7), border=FALSE)
 panel.segplot(x, y, z, centers=centers, subscripts=subscripts,
 ...)
 })

 À disposição.
 Walmes.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Lattice segplot + panel.xblocks

[FHRB Toledo]

Senhores,

Estou tentando fazer um gráfico usando a segplot()... até aí tudo perfeito!

Entretanto gostaria de hachurar uma área do background do gráfico e
tentei de várias formas, mas não consegui! O problema reside é que cada
painel terá uma área diferente para hachurar. Procurei e achei as funções
panel.xblocks() e pane.rect(). Acho que o caminho é por aí!

Gostaria da ajuda de vocês para tentar implementar isso... segue o CMR:


dados - dput(dados)
structure(list(fator = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L), .Label = c(a, b), class = factor), tratamento = c(1L,
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L,
16L, 17L, 18L, 19L, 20L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L,
10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L), media =
c(56.7273408697765,
21.6789539145706, 48.9147297325716, 16.128713312085, 72.2770042662147,
76.4078177803679, 50.9077623429051, 73.4185754594888, 6.10335920092642,
55.0889792680129, 73.2888965805648, 42.4138621293343, 61.1580652843309,
60.6674652959473, 61.985833396912, 59.1674199053476, 64.5810153586145,
49.9847839597647, 43.5403475266414, 73.2055926430283, 78.0443736554009,
29.0590364242541, 67.0013026876617, 21.7957799074036, 94.6993295797361,
97.3586875521169, 70.5384358034064, 84.0236784048615, 6.13251698933115,
72.2706008690223, 95.925565877041, 56.3692012642583, 79.8817867665571,
78.7296774114649, 78.3993769394509, 76.0775293205963, 80.0321741107855,
66.4321920049995, 54.1984087410355, 98.1778312021352), inferior =
c(43.1836043224278,
17.1528746401404, 37.1996905806528, 12.9475504916982, 46.7714226580457,
45.3692479773409, 39.7789554679844, 29.8329460853507, 2.56413580548224,
35.8089238584928, 47.4171926197379, 29.7410223171652, 39.6810438781391,
38.6961342195907, 37.1835630240529, 37.067425724059, 35.4724912266013,
35.2427573484592, 25.0802453164008, 51.397401326238, 67.626248555891,
18.7695537973605, 56.7158673051006, 14.4447744494442, 85.018794282225,
86.9126759092797, 60.2303195514873, 71.7768506217856, 0.319225404377661,
64.8335865483444, 85.7060969685811, 47.3588023848516, 70.3807330972353,
68.9839724287896, 65.992091623, 65.2925804505429, 69.4612476374327,
55.1810503876307, 44.5545294532709, 86.576660755252), superior =
c(62.3108540885626,
23.8615642963601, 53.7260587383325, 17.774181934749, 79.7146304606838,
84.9807494296331, 55.9924261554983, 86.5989631071634, 7.15646923884448,
60.7368137544209, 80.8316127544867, 46.572404145686, 67.4355492078165,
66.989163863718, 68.859173168856, 65.4348577926561, 72.5430140364942,
54.8771182917571, 48.6009761134763, 80.3851368007263, 87.1231257331536,
38.7305442725622, 76.1037592450441, 28.7278361811028, 102.689316368094,
106.025905852459, 79.6868529319033, 94.2065097903592, 11.6427132483976,
78.4162162468943, 104.39108823494, 64.3280982909318, 87.9131833084866,
86.9760720156632, 89.0620270643251, 85.2808640058557, 88.7701994563184,
76.4111720930247, 62.5587781906989, 107.948178577293)), .Names = c(fator,
tratamento, media, inferior, superior), row.names = c(NA,
-40L), class = data.frame)

segplot(reorder(tratamento, media) ~ inferior + superior | fator,
data = dados,
centers = media,
draw.bands = FALSE,
segments.fun = panel.arrows,
ends = 'both',
angle = 90,
length = 1,
unit = 'mm',
panel = function(x, y, z, centers, subscripts, ...) {
  panel.segplot(x, y, z, centers = centers, subscripts =
subscripts, ...)
})

Minha intenção é que o backgrund de cada painel tenha uma área preenchida
de vermelho, por exemplo, mas cada intervalo é diferente por painel... acho
que a única forma é informar essas áreas externamente, com outr data.frame,
seja por exemplo:

dadosAUX - dput(dadosAUX)
structure(list(fator = structure(1:2, .Label = c(a, b), class =
factor),
xleft = c(30, 50), xright = c(60, 80)), .Names = c(fator,
xleft, xright), row.names = c(NA, -2L), class = data.frame)

Assim... em cada painel, preciso que essa área seja preenchida por toda
extensão do eixo Y.

Qualquer ajuda é bem vinda.

Obrigado a todos!

Boa semana e que venha Dezembro!

att,
FH
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Plotar intervalo de confiança

[FHRB Toledo]

Fátima,

O post que compartilho não tem quase nada a ver com o que você pergunta não
fosse o fato de que nele você pode ver uma bela demonstração do uso do
latticeExtra::segplot()

Acho que é um ótimo ponto de partida para você fazer os gráficos com os
intervalos de confiança!

https://ridiculas.wordpress.com/tag/segplot/

Todo o crédito ao Walmes!

att,
FH


2013/11/26 Fátima Lima Paula fatima.lima.pa...@gmail.com

 Prezados, tenho a taxa de mortalidade por idade e queria plotar o
 intervalo de confiança. Alguém pode me ajudar?

 pop=c(37374,  54783,  86049, 126048, 287384)

 obfaF=c(6,15,27,55,239)

 txobfaF=(obfaF*365/popfaF)*100

 txobfaF

 idade=c(1,2,3,4,5)
 plot(idade, txobfaF,type=b,main=Sexo Feminino,xlab=Idade,ylab=Taxa
 de Mortalidade,lty=2,ylim=c(1,40),xaxt='n')

 Obrigada
 Fátima

 --
 Minha felicidade depende da qualidade dos meus pensamentos

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Parametrização com lme4

[FHRB Toledo]

Adriano e Walmes,

Obrigado pelas sugestões, e por compartilhar o livro!

O modelo ajustou corretamente, perfeito!

Obrigado mesmo.

abraço,
FH


2013/11/20 walmes . walmeszevi...@gmail.com

 Fernando,

 Acho que é assim:

 ## fatores cruzados
 xtabs(~local+trat, dados)

 ## fatores aninhados
 xtabs(~trat+eras, dados)

 ##  local
 ## +bloco dentro de local
 ## +tratamento
 ## +interação local tratamento
 mf1 - lm(terms(resp~local/rep+trat+local:trat, keep=TRUE),
   data=dados)
 anova(mf1)

 ##  local
 ## +bloco dentro de local
 ## +eras
 ## +tratamento dentro de eras
 ## +interação local eras
 ## +interação local tratamento dentro de eras
 mf2 - lm(terms(resp~local/rep+eras/trat+local:eras/trat, keep=TRUE),
   data=dados)
 anova(mf2)

 require(lme4)

 mm1 - lmer(resp~(1|local/rep)+(1|trat)+(1|local:trat), data=dados)
 summary(mm1)

 mm2 - lmer(resp~(1|local/rep)+(1|eras/trat)+
 (1|local:eras)+(1|local:eras:trat), data=dados)
 summary(mm2)

 À disposição.
 Walmes.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Extrair p-valor de um objeto lmer

[FHRB Toledo]

Sugiro também ler o help da última versão do pacote de 25 de outubro, há um
tópico especialmente nesse assunto

http://cran.r-project.org/web/packages/lme4/lme4.pdf

att,
FH


2013/11/18 walmes . walmeszevi...@gmail.com

 A lmer() não retorna p-valores e isso é proposital. Douglas Bates, o
 desenvolvedor, esclarece isso nessa mensagem

 https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html

 Além do mais, uma reunião das discussões pode ser acompanhada aqui com
 algumas soluções parciais

 http://rwiki.sciviews.org/doku.php?id=guides:lmer-tests


 À disposição.
 Walmes.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Parametrização com lme4

[FHRB Toledo]

Senhores,

Seguindo no assunto proposto pelo Alisson, no thread sobre pvalor com a
lme4, estou tentando parametrizar um modelo misto com o intuito de fazer a
inferência pelo perfil da verossimilhança.

Entretanto, estou tentando traduzir um modelo da lm() pela parametrização
da lme4 e enfrento grandes dificuldades, assim gostaria da ajuda de vocês
com esse intento...

Seguem dois modelos da lm() que são equivalentes, visto que o fator trat
tem 9 níveis e eras tem um certo confundimento com trat, assim procedo:

mf1 - lm(terms(resp ~ local/rep + trat + local:trat, keep = TRUE), data =
dados)
mf2 - lm(terms(resp ~ local/rep + eras/trat + local:eras/trat, keep =
TRUE), data = dados)

Esses dois modelo tem o mesmo erro e gostaria de reproduzi-los usando a
lmer(). Quando faço o primeiro modelo, sem problemas:

mm1 - lmer(resp ~ ( 1 | local/rep ) + ( 1 | trat ) + ( 1 | local:trat ),
data = dados)

Recupero a mesma variância do resíduo da lm(), todavia não consgigo fazer o
segundo modelo recuperar essa mesma variância... Alguém tem alguma sugestão?

P.S.: Se, por acaso alguém tiver outra sugestão, que use a lme(), ou lme2()
gostaria de saber... como disse, meu objetivo a obter os perfis da
verossimilhança!

Seja os dados de exemplo:

dados - dput(dados)
structure(list(local = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c(1, 2, 3,
4), class = factor), rep = structure(c(1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c(1,
2, 3, 4), class = factor), trat = structure(c(1L, 6L,
2L, 3L, 5L, 8L, 7L, 4L, 9L, 1L, 3L, 7L, 6L, 5L, 4L, 2L, 8L, 9L,
1L, 8L, 4L, 6L, 3L, 9L, 2L, 5L, 7L, 1L, 5L, 9L, 6L, 8L, 7L, 2L,
3L, 4L, 9L, 8L, 3L, 2L, 7L, 4L, 1L, 6L, 5L, 9L, 2L, 1L, 8L, 7L,
6L, 3L, 4L, 5L, 9L, 4L, 6L, 8L, 2L, 5L, 3L, 7L, 1L, 9L, 7L, 5L,
8L, 4L, 1L, 3L, 2L, 6L, 4L, 9L, 2L, 3L, 8L, 6L, 5L, 1L, 7L, 4L,
3L, 5L, 9L, 8L, 1L, 2L, 6L, 7L, 4L, 6L, 1L, 9L, 3L, 7L, 2L, 8L,
5L, 4L, 8L, 7L, 9L, 6L, 5L, 2L, 3L, 1L, 9L, 2L, 7L, 8L, 4L, 5L,
6L, 3L, 1L, 9L, 8L, 6L, 2L, 4L, 3L, 7L, 5L, 1L, 9L, 5L, 3L, 2L,
8L, 1L, 7L, 4L, 6L, 9L, 4L, 1L, 2L, 5L, 6L, 7L, 8L, 3L), .Label = c(a,
b, c, d, e, f, g, h, i), class = factor), eras =
structure(c(2L,
4L, 2L, 3L, 3L, 1L, 4L, 3L, 1L, 2L, 3L, 4L, 4L, 3L, 3L, 2L, 1L,
1L, 2L, 1L, 3L, 4L, 3L, 1L, 2L, 3L, 4L, 2L, 3L, 1L, 4L, 1L, 4L,
2L, 3L, 3L, 1L, 1L, 3L, 2L, 4L, 3L, 2L, 4L, 3L, 1L, 2L, 2L, 1L,
4L, 4L, 3L, 3L, 3L, 1L, 3L, 4L, 1L, 2L, 3L, 3L, 4L, 2L, 1L, 4L,
3L, 1L, 3L, 2L, 3L, 2L, 4L, 3L, 1L, 2L, 3L, 1L, 4L, 3L, 2L, 4L,
3L, 3L, 3L, 1L, 1L, 2L, 2L, 4L, 4L, 3L, 4L, 2L, 1L, 3L, 4L, 2L,
1L, 3L, 3L, 1L, 4L, 1L, 4L, 3L, 2L, 3L, 2L, 1L, 2L, 4L, 1L, 3L,
3L, 4L, 3L, 2L, 1L, 1L, 4L, 2L, 3L, 3L, 4L, 3L, 2L, 1L, 3L, 3L,
2L, 1L, 2L, 4L, 3L, 4L, 1L, 3L, 2L, 2L, 3L, 4L, 4L, 1L, 3L), .Label =
c(0,
1, 2, 3), class = factor), resp = c(1878L, 3417L, 3546L,
2284L, 3220L, 2692L, 2243L, 2088L, 1871L, 1793L, 1983L, 2087L,
2628L, 2909L, 1928L, 3065L, 2650L, 2390L, 2428L, 3125L, 3010L,
2181L, 2478L, 2575L, 2348L, 2492L, 2295L, 1892L, 2257L, 2003L,
3123L, 2673L, 2449L, L, 3040L, 3212L, 1767L, 2514L, 2503L,
2547L, 2055L, 1591L, 1898L, 2748L, 2376L, 1458L, 1709L, 1649L,
1784L, 1724L, 1956L, 1706L, 1043L, 1056L, 1346L, 1657L, 1970L,
1483L, 1103L, 985L, 779L, 915L, 800L, 1155L, 1074L, 1435L, 988L,
1052L, 932L, 1112L, 1728L, 2746L, 5150L, 3101L, 5020L, 5828L,
3927L, 5123L, 5765L, 3620L, 4901L, 5859L, 5050L, 5383L, 2920L,
4002L, 3850L, 4727L, 4998L, 5702L, 5154L, 5098L, 4039L, 3493L,
5417L, 5097L, 4954L, 3964L, 4371L, 5681L, 4090L, 4644L, 3370L,
5529L, 5491L, 5145L, 5158L, 3864L, 2777L, 2697L, 2739L, 1911L,
2103L, 2959L, 3241L, 2457L, 2321L, 2410L, 2085L, 3511L, 2657L,
3349L, 3157L, 2948L, 3491L, 2951L, 2626L, 2396L, 1997L, 2895L,
2244L, 2807L, 3408L, 3350L, 4104L, 2613L, 2589L, 2617L, 2791L,
3063L, 3665L, 2918L, 2092L, 2583L)), .Names = c(local, rep,
trat, eras, resp), class = data.frame, row.names = c(NA,
-144L))

att,
FH

[R-br] Círculo

[FHRB Toledo]

Senhores,

Segue um CMR extremamente simples... a pergunta o resultado está certo?

require(plotrix)
plot(c(-1.2, 1.2), c(-1.2, 1.2), type = 'n', axes = FALSE, xlab = '', ylab
= '')
rect(-1, -1, 1, 1, lwd = 1.5)
draw.circle(0, 0, 1)

Por que será que obtenho um gráfico em que o círculo ultrapassa o limite
superior do quadrado?

att,
FH
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Círculo

[FHRB Toledo]

Obtive uma resposta correta...

plot(c(-1.2, 1.2), c(-1.2, 1.2), type = 'n', axes = FALSE, xlab = '', ylab
= '')
rect(-1, -1, 1, 1, lwd = 1.5)
x - cos(seq(0, 2 * pi, length = 1000))
y - sin(seq(0, 2 * pi, length = 1000))
polygon(x, y, lwd = 1.5)

att,
FH


2013/11/14 FHRB Toledo fernandohtol...@gmail.com

 Senhores,

 Segue um CMR extremamente simples... a pergunta o resultado está certo?

 require(plotrix)
 plot(c(-1.2, 1.2), c(-1.2, 1.2), type = 'n', axes = FALSE, xlab = '', ylab
 = '')
 rect(-1, -1, 1, 1, lwd = 1.5)
 draw.circle(0, 0, 1)

 Por que será que obtenho um gráfico em que o círculo ultrapassa o limite
 superior do quadrado?

 att,
 FH

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Materiais de programacao gratuitos

[FHRB Toledo]

Muita leitura pela frente... não se esqueçam de comer e escovar os dentes!


On 18 October 2013 12:10, Benilton Carvalho beniltoncarva...@gmail.comwrote:


 https://github.com/vhf/free-programming-books/blob/master/free-programming-books.md

 divirtam-se


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Remover NULL de um output da função lapply

[FHRB Toledo]

THiago,

substitua o 0 ao seu bel prazer:

ifelse(x == 'a', 1, 0)

att,
FH


2013/10/14 Tiago Souza Marçal tiagosouzamar...@hotmail.com

 Bom dia pessoal,
 gostaria de saber como remover os NULL do CMR abaixo.

 x=c(b,a,c)
 u=lapply(1:3,function(i){
   if(x[i]==a){
   1}
   }
 )

 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

 #

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Aggregate

[FHRB Toledo]

Fatima,

Se entendi seu problema

Segue uma solução...

Note, que o vetor tcont deve ser numérico e foi preciso recodificar os s
e n, como 0 e 1!

sexo - c(1,1,2,2,1,2,1,1,1,1)
tcont - c(100,320,24,256,134,290,18,34,15,121)
read - c(s,n,n,n,s,s,s,n,n,s)

dados - data.frame(sexo,tcont,read)
dados

dados$read - as.numeric(ifelse(read == 's', 1, 0))

dados

dados.resumo - aggregate(dados[,2:3],
  by = list(sexo = dados$sexo),
  FUN = sum)

att,
FH


2013/10/2 Fátima Lima Paula fatima.lima.pa...@gmail.com

 Pessoal, estou querendo agregar umas variáveis, mas não sei.
 Imaginemos o seguinte data frame:
 sexo=c(1,1,2,2,1,2,1,1,1,1)
 tcont=c(100,320,24,256,134,290,18,34,15,121)
 read=c(s,n,n,n,s,s,s,n,n,s)
 df=data.frame(sexo,idade,read)
 df
 Quero agregar essas variáveis de forma que fique a primeira coluna sexo, a
 segunda coluna a soma do tcont e a terceira coluna a contagem dos read =
 sim.
 Alguém pode me ajudar, por favor?
 Obrigada
 Fátima

 --
 Minha felicidade depende da qualidade dos meus pensamentos

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Teste de Hartley

[FHRB Toledo]

Uai,

Mas obter a probabilidade da distribuição dada a sua razão entre as
variâncias não é o mesmo que fazer o teste?

Com certeza, essa função não é o teste, mas é toda a família funções p, q,
d e r da distribuição de máximos de F. Assim, um pouco de conhecimento
sobre o que realmente é o Teste de Hartley e está feito.

FH


On 20 September 2013 13:33, Luiz Henrique da Conceicao Leal 
richfi...@ig.com.br wrote:

 Já tinha visto, mas ela dá a distribuição e não o teste em sim

 Luiz

 Em 19 de setembro de 2013 15:29, FHRB Toledo 
 fernandohtol...@gmail.comescreveu:

  Veja so o pacote SuppDists, função maxFratio!

 att,
 FH


  2013/9/19 Luiz Henrique da Conceicao Leal richfi...@ig.com.br

 Pessoal, alguém conhece alguma função que realize o teste de Hartley de
 homogeneidade de variâncias
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Teste de Hartley

[FHRB Toledo]

Veja so o pacote SuppDists, função maxFratio!

att,
FH


2013/9/19 Luiz Henrique da Conceicao Leal richfi...@ig.com.br

 Pessoal, alguém conhece alguma função que realize o teste de Hartley de
 homogeneidade de variâncias
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Passar parâmetro para with() dentro de uma função

[FHRB Toledo]

Faça:

planilha - data.frame(x=1:10, v1=rnorm(10), v2=rnorm(10))

with(planilha, v1+v2) # Fora da função funciona

func - function(dados, func) {
 with(dados, func)
}

func(planilha, v1+v2) # Da erro na hora de chamar a função

Talvez seja qq confusão na atribuição de nomes, parâmetros para a função.

att,
FH


2013/8/19 Rodrigo Coster rcos...@gmail.com

 Caros,

 eu tenho uma função que me retorna um data.frame() com várias colunas que
 vou utilizar para outros calculos. Eu queria facilitar a minha vida e
 colocar esses calculos dentro da função, como um parâmetro (dado que não
 existe somente um calculo). Encontrei uma maneira de fazer isso fora da
 função pelo with() ou até mesmo evalq(), só que não consegui fazer isso
 dentro da função. Eu entendo os erros que dão e concordo com eles, mas não
 consegui encontrar uma maneira que funcione.

 Basicamente é isso:

 dados - data.frame(x=1:10, v1=rnorm(10), v2=rnorm(10))

 with(dados, v1+v2) # Fora da função funciona

 func - function(data, func) {
  with(dados, func)
 }
 func(dados, v1+v2) # Da erro na hora de chamar a função
 func(dados, 'v1+v2') # Retorna o string

 Alguma sugestão/solução?

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Atualização R Linux

[FHRB Toledo]

Teodoro,

Veja na página do R www.r-project.org no item instalação no Linux, escolha
a aba ubuntu e veja como alterar o repositório e atualizar para a versão
mais atual!

att,
FH

2013/6/26 Teodoro Calvo teocal...@gmail.com

 Olá, bom dia a todos.

 Sou novo por aqui.

 Atualmente utilizo o RKWard na plataforma ubuntu 13.04, porém o r-base
 esta na versão 2.15.

 Como posso atualiza-lo para a 3.0.1 ?
 Já tentei apt-get pelo terminal, e nada !

 Obrigado desde já.

 Att Téo Calvo
 __**_
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem 
 (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia)
 e forneça código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Selecionar casos - Função

[FHRB Toledo]

row.names a logical value indicating whether the row names of x are to be
written along with x to the file

Colocando row.names = FALSE você resolve seu problema!

att,
FH


2013/6/23 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Fernando

 Não entendi! Qual comando eu uso pra ignorar? Não achei


 Em 23 de junho de 2013 20:04, FHRB Toledo fernandohtol...@gmail.comescreveu:

 André,

 Uma busca no google, página 37:
 http://cran.r-project.org/web/packages/xlsx/xlsx.pdf

 att,
 FH


 2013/6/23 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Resolvi usando a biblioteca xlsx

 for (nm in Nms) write.xlsx(Res[[nm]], paste(nm, 'xls',
 sep='.'),sheetName=paste(nm), col.names=TRUE)

 Porém isso cria nos meus bancos de dados uma variável com o nome em
 branco que é a a posição da linha;
 Ex,:


 codigo sexo idade  1 10001 m 19  2 10001 m 19  3 10001 m 20  4 10001 m
 20
 Como ignorar isso, para que a variável não seja criada no meu banco?

 Abraços


 Em 23 de junho de 2013 17:30, Sérgio Henrique almeida da silva ju 
 sergio.edfis...@gmail.com escreveu:

 Essa biblioteca não está mais disponível! :-(


 Em 23 de junho de 2013 17:28, Raphael Saldanha 
 rfsalda...@outlook.comescreveu:

  Compreendo.

 Tente este:

 library(xlsReadWrite)
 write.xls(mydata, c:/mydata.xls)

 De: http://www.statmethods.net/input/exportingdata.html

 ---

 Atenciosamente,
 Raphael Saldanha

 rfsalda...@outlook.com


 --
 Date: Sun, 23 Jun 2013 16:49:28 -0300
 From: sergio.edfis...@gmail.com

 To: r-br@listas.c3sl.ufpr.br
 Subject: Re: [R-br] Selecionar casos - Função

 o problema é que estes arquivos devem ser salvos dentro de arquivos
 zip com formato xls, a empresa deseja usa diretamente esses arquivos, não
 sei o motivo... para evitar todo trabalho manual eu estou tentando fazer
 isso diretamente pelo R.


 Em 23 de junho de 2013 16:45, Raphael Saldanha rfsalda...@outlook.com
  escreveu:

 Você pode salvar como CSV. O Excel costuma reconhecer e abrir o
 arquivo sem problemas:

 ?write.csv2

 ---

 Atenciosamente,
 Raphael Saldanha

 rfsalda...@outlook.com


 --
 Date: Sun, 23 Jun 2013 16:12:27 -0300
 From: beniltoncarva...@gmail.com
 To: r-br@listas.c3sl.ufpr.br
 Subject: Re: [R-br] Selecionar casos - Função


 Você terá de esperar a ajuda de alguém que use Excel... Salvar como
 texto, como indiquei, não te impede de abrir no Excel...
 On 23 Jun 2013 14:41, Sérgio Henrique almeida da silva ju 
 sergio.edfis...@gmail.com wrote:

 Obrigado amigo

 Mas nesse caso surge um problema em dados como fatores, por exemplo,
 no meu caso ele cria um \m\ ou \f\, como resolver isso?


 Em 23 de junho de 2013 14:35, ctrucios ctruc...@gmail.com escreveu:

 Oi, eu utilizou a função write.table. Aqui um exemplo
 write.table(TC95,TC95.xls, sep=\t,col.names=TRUE, row.names=FALSE,
 quote=TRUE, na=NA)
 Em 23/06/2013 14:31, Sérgio Henrique almeida da silva ju 
 sergio.edfis...@gmail.com escreveu:

 Só mais uma pergunta

 Como posso salvar esses arquivos em excel?
 Eu tentei o comando da biblioteca WriteXLS, mas não rodou

 for (nm in Nms) WriteXLS(Res[[nm]], ExcelFileName = paste(nm, 'xls',
 sep='.'))

 tem algum jeito melhor de salvar?

 Obrigado


 Em 23 de junho de 2013 12:51, Sérgio Henrique almeida da silva ju 
 sergio.edfis...@gmail.com escreveu:

 Obrigado

 Deu tudo certinho!

 Abraços


 Em 23 de junho de 2013 01:28, Benilton Carvalho 
 beniltoncarva...@gmail.com escreveu:

 Para dividir em subconjuntos, use o comando split().

 Res = split(dados,  dados$codigo)

 Daí, combine com um for() loop para gravar (ou mesmo um lapply).

 Nms = names(Res)
 for (nm in Nms) write.table(Res[[nm]], file=paste(nm, 'txt', sep='.'))

 (código não-testado)

 Sobre gravar compactado, vc pode mudar o nome do arquivo para uma
 conexão gzip ou, se bem me lembro, zip. Veja a ajuda para os comandos de
 mesmo nome.

 b
 On 22 Jun 2013 23:23, Sérgio Henrique almeida da silva ju 
 sergio.edfis...@gmail.com wrote:

 Prezados

 Gostaria de fazer uma função que através de um banco formasse diversos
 outros bancos selecionados por uma variável.

 Exemplo:

 codigo - c(rep(10001,10), rep(10005,15),rep(20001,20))
 sexo - c(rep(m,20),rep(f,25))
 idade - rnorm(45,20)

 dados -
 cbind(as.data.frame(codigo),as.data.frame(sexo),as.data.frame(idade))

 Quero quebrar esse banco dados em diversos outros bancos pela variável
 codigo

 como:

 dados1 - dados[which(dados$codigo==10001),]
 .
 .
 .
 dadosn - dados[which(dados$codigo==n),]

 Posso fazer uma table(codigo) e jogar os valores dentro desse comando,
 mas não sei como fazer isso, deixando a função mais automática.

 Outra pergunta tem como eu saber esses bancos comprimidos através do
 R, por exemplo em ZIP ou RAR?

 Obrigado
 --
 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

Re: [R-br] Seleção de plantas ( sorteio)

[FHRB Toledo]

Helio,

Uma sugestão para seu problema...

dados - data.frame(ponto = c(rep(3001,6),rep(3002,6),rep(3003,6)),
planta = c(rep(seq(1:6),3)),
valor = seq(1.5,10,.5))

dados.sub - split(dados, dados$planta)

amostrador - function(x, tamanho) x[sample(1:nrow(x), tamanho, replace =
FALSE), ]

dados2 - do.call(rbind, lapply(dados.sub, amostrador, tamanho = 2))
dados3 - do.call(rbind, lapply(dados.sub, amostrador, tamanho = 3))
dados4 - do.call(rbind, lapply(dados.sub, amostrador, tamanho = 4))
##...

Vale lembrar que pela implementação a amostragem das linhas é sem
reposição, para alterar tal comportamento use o replace = TRUE!

att,
FH


2013/5/24 Rodrigo Coster rcos...@gmail.com

 Helio, dai temos varias possiveis soluções.. ou sortear as 2 colunas no
 fim:

 do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x),
 ...), ], size=1))[,c(1:2, sample(3:4, 1))]

 Ou mexer na função dentro do apply e fazer o sorteio para cada ponto:

 do.call(rbind, by(dados, dados[,1], function(x, ...) { z -
 x[sample(1:nrow(x), ...), c(1:2, sample(3:4, 1))] ; colnames(z) -
 c('ponto', 'planta', 'valor') ; z }  , size=1))


 2013/5/24 Hélio Gallo Rocha heliogalloro...@gmail.com

 Mandei exemplo em que dobrou o numero de plantas, esse seria o exemplo
 correto:

 ponto=c(rep(3001,6),rep(3002,6),rep(3003,6));ponto   ;length(ponto)
 planta=c(rep(seq(1:6),3));planta  ;length(planta)
 valor1=seq(1.25,10,.5);valor1;length(valor1)
 valor2=seq(10,1.25,-.5);valor2;length(valor2)
 dados=cbind(ponto,planta,valor1,valor2);dados


 Grato

 Hélio


 Em 24 de maio de 2013 09:20, Hélio Gallo Rocha heliogalloro...@gmail.com
  escreveu:

 Rodrigo, fiquei em outra situação.

 A função envida funciona, mas se tiver mais uma variável como este
 exemplo:
 ponto=c(rep(3001,12),rep(3002,12),rep(3003,12))
 planta=c(rep(seq(1:6),6))
 valor1=seq(1.25,10,.25)
 valor2=seq(10,1.25,-.25)
 ramo1=c(rep(A,36))
 ramo2=c(rep(B,36))
 dados=cbind(ponto,planta,valor1,valor2);dados

  ponto planta valor1 valor2
  [1,]  3001  1   1.25  10.00
  [2,]  3001  2   1.50   9.75
  [3,]  3001  3   1.75   9.50
  [4,]  3001  4   2.00   9.25
  [5,]  3001  5   2.25   9.00
  [6,]  3001  6   2.50   8.75
  [7,]  3001  1   2.75   8.50
  [8,]  3001  2   3.00   8.25
  [9,]  3001  3   3.25   8.00
 [10,]  3001  4   3.50   7.75
 [11,]  3001  5   3.75   7.50
 [12,]  3001  6   4.00   7.25
 [13,]  3002  1   4.25   7.00
 [14,]  3002  2   4.50   6.75
 [15,]  3002  3   4.75   6.50
 [16,]  3002  4   5.00   6.25
 [17,]  3002  5   5.25   6.00
 [18,]  3002  6   5.50   5.75
 [19,]  3002  1   5.75   5.50
 [20,]  3002  2   6.00   5.25
 [21,]  3002  3   6.25   5.00
 [22,]  3002  4   6.50   4.75
 [23,]  3002  5   6.75   4.50
 [24,]  3002  6   7.00   4.25
 [25,]  3003  1   7.25   4.00
 [26,]  3003  2   7.50   3.75
 [27,]  3003  3   7.75   3.50
 [28,]  3003  4   8.00   3.25
 [29,]  3003  5   8.25   3.00
 [30,]  3003  6   8.50   2.75
 [31,]  3003  1   8.75   2.50
 [32,]  3003  2   9.00   2.25
 [33,]  3003  3   9.25   2.00
 [34,]  3003  4   9.50   1.75
 [35,]  3003  5   9.75   1.50
 [36,]  3003  6  10.00   1.25

 do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x),
 ...), ], size=1))
  ponto planta valor1 valor2
 3001  3   1.75   9.50
 3002  5   6.75   4.50
 3003  4   9.50   1.75

 No caso seria sortear apenas um valor.

 A função que enviou resolve se as colunas valor 1 e valor 2 estiverem
 uma debaixo da outra, no  exemplo seriam 72 linhas.

 Se tiver solução, seria mais fácil,

 Obrigado,
 Hélio


 Em 23 de maio de 2013 21:39, Rodrigo Coster [via R-br] 
 ml-node+s2285057n4659434...@n4.nabble.com escreveu:

 Ops, é sortear plantas dentro de ponto, nao ponto dentro de plantas.
 Troca o dados[,2] por dados[,1]

 do.call(rbind, by(dados, dados[,1], function(x, ...)
 x[sample(1:nrow(x), ...), ], size=1))



 2013/5/23 Rodrigo Coster [hidden 
 email]http://user/SendEmail.jtp?type=nodenode=4659434i=0
 

 Da para fazer com do.call() e by():

 do.call(rbind, by(dados, dados[,2], function(x, ...)
 x[sample(1:nrow(x), ...), ], size=1))

 Depois é só ir mudando o parametro size para outros tamanhos... e se
 quiser amostragem com reposição, adicionar o replace=TRUE depois do size.


 2013/5/23 Hélio Gallo Rocha [hidden 
 email]http://user/SendEmail.jtp?type=nodenode=4659434i=1
 

 Caros da Lista

 Com os dados hipotéticos abaixo:

 ponto=c(rep(3001,6),rep(3002,6),rep(3003,6))
 planta=c(rep(seq(1:6),3))
 valor=seq(1.5,10,.5)
 dados=cbind(ponto,planta,valor);dados

  ponto planta valor
  [1,]  3001  1   1.5
  [2,]  3001  2   2.0
  [3,]  3001  3   2.5
  [4,]  3001  4   3.0
  [5,]  3001  5   3.5
  [6,]  3001  6   4.0
  [7,]  3002  1   4.5
  [8,]  3002  2   5.0
  [9,]  3002  3   5.5
 [10,]  3002  4   6.0
 [11,]  3002  5   6.5
 [12,]  3002  

Re: [R-br] regressão no R

[FHRB Toledo]

Alana,

Leia o guia de postagens!

Não envie arquivos para o email da lista!

Produza um Código mínimo reproduzível!

Sobre o seu problema, note que seu arquivo DadosPanakeia.txt tem a linha
custosPanakeia com os valores numéricos separadas as casas decimais com
virgula ,, enquanto que a linha seguinte hgbPanakeia os valores estão
separados por ponto .!

Padronize, ou tudo com ponto, ou tudo com virgula... tome também o cuidado
de no seu script, na linha do read.table() coloque para o argumento dec o
que você usou como separador de decimal!

Feito isso a regressão é estimada sem problemas.

Valeu destacar que isso em momento algum envolveu um erro ou problema do
programa R, sendo apenas questão de observação, organização e conhecimento
dos próprios dados, atente-se a isso!

Á disposição,
FH


2013/5/15 alanaro...@sapo.pt

 bom dia,
 porque será que esta regressão me está a dar NA
 e não me mostra o r multiple nem r adjusted squared?
 não percebo o erro do gráfico tambem.
 mando em anexo.
 envio tambem o ficheiro excel dos dados originais...
 obrigada
 Ana


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Valor Esperado....

[FHRB Toledo]

André,

O gabarito está certo... ao menos no meu raciocínio, acomanhe!

Na urna havia 2B e 3P e a amostragem é SEM reposição, então:

a Probabilidade de sair 2 bolas brancas é 2/5 * 1/4 = 2/20
a prob de sair uma branca e e preta é 2/5 * 3/4 + 3/5 * 2/4 = 12/20
e a probabilidade de sair duas pretas é 3/5 * 2/4 = 6/20

A soma dessas probabilidades dá 1 (obviamente)

Se X é a variável aleatória de quantas brancas tem X pode ser 0, 1 ou 2,
procedendo  Soma X.p(X), chegamos em 16/20, que é o mesmo que 20/25!

Note que no seu raciocínio você não considerou que o número de bolas pretas
e brancas eram diferentes na urna... !

espero ter ajudado.

att,
FH


2013/5/15 andre...@bol.com.br

 Olá colegas,

 Em primeiro lugar, não é exercicio que estou postando, se trata de uma
 questão de concurso que acredito que não esteja correto o gabarito.
 A questão é a seguinte:

 Numa urna existem 2 bolas brancas e 3 pretas. Serão retiradas da urna 2
 bolas SEM reposição.
 Seja X o número de bolas brancas retiradas. O valor esperado da variável
 aleatória X é:

 a) 1
 b) 25/20
 c) 20/25
 d) 1,5
 e) 0,75

 #

 Então, temos o seguinte:

 Ω = {BB,BP,PP,PB}  # espaço amostral;
 X = {0,1,2}   # valores possíveis que a variável aleatória pode
 assumir.

 em que,
 B = bola branca;
 P = bola preta.

 Assim, teremos a seguinte distribuição de probabilidade:

 ___
 X012
 ___
 p(X)1/42/41/4
 ___

 E(X) = Σxp(x)
= 0.(1/4) + 1.(2/4) + 2.(1/4)
= 1 (letra A da questão)

 Segundo o gabarito da prova, a resposta correta é a letra C), isto é, E(X)
 = 20/25

 Alguém aqui, poderia confrimar ou não esse resultado.

 desde já agradeço!

 *Att.*
 *André*


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Survival

[FHRB Toledo]

Ana Luiza,

google it: [R] Kaplan-Meier

Veja o que sai na primeira entrada:
http://www.stat.ucdavis.edu/~hiwang/teaching/10fall/R_tutorial%201.pdf

Espero ter ajudado!

att,
FH




2013/5/6 Ana Luiza Cassin pesquisaclin...@labcor.com.br

  Prezados,


 Estou começando a usar or programa R e preciso fazer curva de
 sobrevivência (Kaplan-Meier) sobre o acompanhamento de pacientes ao longo
 do tempo.
 Gentileza, se alguém tiver conhecimento ou algum tutorial sobre essa
 análise, o mais detalhado possível, agradeceria.

 Atenciosamente,
 --
   [image: Labcor] http://www.labcor.com.br/
 Ana Luiza Cassin
 Monitora de Pesquisa
 Monitor Research
 Tel.com.: +55(31)3313-1251
 Tel.cel.: +55(31)9308-0072

 Labcor
 www.labcor.com.br

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

email.jpg___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] R 3.0 Ubuntu 13.04

[FHRB Toledo]

Mude o repositório!


2013/4/29 Frederico Steinmetz Alvarez eng.fredericoalva...@gmail.com

 Boa tarde Benilton.

 Estou dando uma olhada nos pacotes disponíveis no servidor da ufpr para o
 ubuntu 13.04 e os mesmos são:


 r-base-core-dbg_3.0.0-2raring2_i386.debhttp://cran-r.c3sl.ufpr.br/bin/linux/ubuntu/raring/r-base-core-dbg_3.0.0-2raring2_i386.deb28-Apr-2013
 07:03 3.5M   [image: [ 
 ]]r-base-core-dbg_3.0.0-2raring_i386.debhttp://cran-r.c3sl.ufpr.br/bin/linux/ubuntu/raring/r-base-core-dbg_3.0.0-2raring_i386.deb
  21-Apr-2013
 21:48 3.5M  [image: [ 
 ]]r-base-core_3.0.0-2raring2_i386.debhttp://cran-r.c3sl.ufpr.br/bin/linux/ubuntu/raring/r-base-core_3.0.0-2raring2_i386.deb
  28-Apr-2013
 07:03 20M  [image: [ 
 ]]r-base-core_3.0.0-2raring_i386.debhttp://cran-r.c3sl.ufpr.br/bin/linux/ubuntu/raring/r-base-core_3.0.0-2raring_i386.deb
  21-Apr-2013
 21:48 20M Não estou encontrando a versão para 64bis.

 Quando coloco apt-get install r-base-core, ele pede a versão 3.0.0.2, mas
 diz que apenas a 2.13.2 está disponível.

 Semana passada fiz esta instalação num debian wheezy e foi sem problemas.

 Saudações.


 Em 29 de abril de 2013 13:50, Benilton Carvalho 
 beniltoncarva...@gmail.com escreveu:

 Para o 13.04, que e' a ultima versao estavel AFAIK (conforme titulo da
 sua msg), tudo deve estar disponivel.

 Para o 13.10, imagino que a forma mais segura de ter o R-3.0.0 e'
 compilando voce mesmo

 b

 Em 29 de abril de 2013 13:37, Frederico Alvarez
 eng.fredericoalva...@gmail.com escreveu:
  Boa tarde a todos.
 
  Estou tentando instalar o R 3.0.0 no ubuntu 13.10 via apt.
  Mas acontece que no repositório da UFPR não há o pacote r-base-core
 para 64
  bit. Porém, há dois pacotes para 32 bits.
 
  Alguém sabe como posso resolver esse problema?
 
  Saudações.
 
  --
  Frederico Steinmetz Alvarez
  Engenheiro de Produção Elétrica
  Mestrando em Energia Eólica
  Laboratório Didático Pedagógico de Fluidos - LDPFLU
  Departamento de Engenharia Mecânica - DEMEC
  Universidade Federal de Pernambuco - UFPE
  (81) 8885-1105 (Oi)
  (48) 8448-5994 (Oi)
  (48) 9143-5141 (Vivo)
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Frederico Steinmetz Alvarez
 Engenheiro de Produção Elétrica
 Mestrando em Energia Eólica
 Laboratório Didático Pedagógico de Fluidos - LDPFLU
 Departamento de Engenharia Mecânica - DEMEC
 Universidade Federal de Pernambuco - UFPE
 (81) 8885-1105 (Oi)
 (48) 9143-5141 (Vivo)

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Matriz X (Y = XB + e)

[FHRB Toledo]

No ajuste do aov, a opção x = TRUE não te serve?

2013/4/23 Tiago Souza Marçal tiagosouzamar...@hotmail.com:
 Prezados,

 gostaria de saber se existe uma função que simplica a obtensão desta matriz?

 m-matrix(c(rep(1,15),rep(1,5),rep(0,10),rep(0,5),rep(1,5),rep(0,5),rep(0,10),rep(1,5)),nr=15)

 Att.

 Tiago.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
 mínimo reproduzível.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] obter media dos coeficientes de modelos

[FHRB Toledo]

Maurício,

Então, se bem entendi o seu problema acho que aí vaí uma solução:

n.trat - 3 # número de tratamentos
n.rep - 2 # número de repetições
n.time - 4 # número de tempos
n.ind - n.trat * n.rep # número de indivíduos
n - n.trat * n.rep * n.time # número total de registros
id - rep(1:n.ind, each = n.time) # identificação do indivíduo
time - rep(1:n.time, times = n.ind) # identificação do tempo
trat - rep(factor(LETTERS[1:n.trat]), each = n.rep*n.time) # tratamentos

dados - data.frame(id, time, trat) # estrutura do delineamento

y - function(x) ifelse(dados$time  2, rbinom(1, 1, .78), rbinom(1, 1,
.35)) ## funcao que simula um ensaio

nsmc - 1000 # numero de simulacoes de monte carlo
varios - replicate(nsmc, cbind(dados, y = y(x)), simplify = FALSE) #
repete o ensaio nsmc vezes

modelo - function( plan ) glm(y ~ trat, data = plan, family =
binomial)$coef # funcao que faz o ajuste

ajustes - do.call(rbind, lapply(varios, modelo)) # varios ajustes

medias - apply(ajustes, 2, mean) # médias dos coeficientes

Primeiro criei uma função que atribui os resultados ao vetor y, segundo
suas premissas. O replicate() repete esse procedimento por nsmc vezes e
armazena cada uma das simulações em uma posição da lista.

Depois disso é só usar a sua função de ajuste (usei a glm normal) em uma
função e aplica com o lapply(). Fiz um artifício de já organizar os
coeficientes do modelo numa matriz, uma linha para cada simulação e uma
coluna por parâmetro.

Espero ter ajudado.

abraço,
FH


2013/4/10 Maurício Lordêlo mslord...@gmail.com

 Caros,
 No CMR abaixo crio uma matriz de delineamento que estah associada ao nro
 de tratamentos, nro de repeticoes de cada tratamento e o tempo que cada
 resposta (neste caso binaria) foi observada.
 Suponha que o vetor y seja gerado levando em consideracao que no tempo=1,
 a probabilidade eh 0.78; nos tempos1, esta probabilidade eh 0.35.

 n.trat=3 # número de tratamentos
 n.rep=2 # número de repetições
 n.time=4  # número de tempos
 n.ind=n.trat*n.rep # número de indivíduos
 n=n.trat*n.rep*n.time # número total de registros
 id=rep(1:n.ind, each = n.time) # identificação do indivíduo
 time=rep(1:n.time, times = n.ind) # identificação do tempo
 trat=rep(factor(LETTERS[1:n.trat]), each = n.rep*n.time) # tratamentos

 y = matrix(data = NA, nrow = n, ncol = 1, byrow = FALSE,dimnames = NULL)
 data1=data.frame(id,time,trat,y)
 for(i in 1:n)
 {
 y[i]=ifelse(data1[i,2]2,rbinom(1,1,.78),rbinom(1,1,.35))
 }
 data2=data.frame(id,time,trat,y)
 data2

 Neste caso tenho formado meu conjunto de dados (data2) já com os valores
 de y gerados.
 Agora vem a minha dificuldade. Preciso rodar um mesmo modelo repetidas
 vezes, como por exemplo:

 modelo1=geeglm(y~trat, id=id, data=data2, family=binomial)

 e depois obter uma media para cada um dos coeficientes (com seus
 erros-padrao) obtidos.
 Ou seja, terei outros conjuntos de dados gerados (com o mesmo nro de
 tratamentos, repeticoes e tempos)
 sendo que a unica diferenca estarah nos valores de y. Preciso obter a
 média dos coeficientes de cada um
 dos modelos rodados a partir destes diferentes y´s.
 Agradeco quem puder auxiliar.
 Mauricio


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Atualização R 2.15.3i

[FHRB Toledo]

Sergio,

Sua distribuicão linux usa um repositório para os programas, inspecione o
arquivo source.list que está em /etc/apt/

Pois bem, quando você instala o R, provavelmente, a versão default que ele
busca nos repositórios padrão da sua distro é outro, que não a mais
recente... deste modo, você deve incluir um repositório da versão mais
recente no seu arquivo source.list e fazer a atualizacão do seu R!

abraco,
FH

2013/3/11 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Oi Fernando

 Obrigado pela ajuda.

 Não entendi muito bem o que você quis dizer...

 Abraços

 Em 11 de março de 2013 21:05, FHRB Toledo fernandohtol...@gmail.comescreveu:

 Use outro repositório, nesse caso o repositório da versão mais atual e
 não o padrão da sua distro!

 2013/3/11 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Prezados

 Estou tentando atualizar meu R para a versão 2.15.3, porém não estou
 conseguindo. Já utilizei os tutorias anteriores aqui da lista mais não deu
 certo.

 Utilizo o Ubuntu 12.10 (quantal) 64 bit e meu R está na versão 2.15.1.

 Alguém tem alguma ideia de como posso atualizar?

 Att

 --
 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

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 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

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Atualização R 2.15.3

[FHRB Toledo]

Use outro repositório, nesse caso o repositório da versão mais atual e não
o padrão da sua distro!

2013/3/11 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Prezados

 Estou tentando atualizar meu R para a versão 2.15.3, porém não estou
 conseguindo. Já utilizei os tutorias anteriores aqui da lista mais não deu
 certo.

 Utilizo o Ubuntu 12.10 (quantal) 64 bit e meu R está na versão 2.15.1.

 Alguém tem alguma ideia de como posso atualizar?

 Att

 --
 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

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] [OFF-TOPIC] Livros de divulgação científica em estatística

[FHRB Toledo]

Uma sugestão é o livro: O livro dos códigos, que conta a história da
criptografia e de Alan Turing, um pouco fora especificamente da estatística
mas muito interessante!

2013/3/10 Simone D. Sartorio sisarto...@yahoo.com.br

 Alexandre, tem um coleção que eu gostei bastante, para iniciar e motivar
 os estudos dos alunos em um formato diferente:

 Guia mangá de Estatística
 http://novatec.com.br/livros/mangaestatistica/

 abraços
 Simone

 *
 ***---***---***---***---***---***---***---***---***---***---***---***---***
 *
 **  Simone Daniela Sartorio*
 *   Professora Adjunta I da UFSCar, Centro de Ciências Agrárias, Campus
 Araras/SP.
 * Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP;
 *   Licenciada em Matemática - UNESP/Rio Claro.
 *   Contatos:
 * tel: (19)3543-2614, ramal 2664
 * cel: (19)8182-0586 - Tim
 *
 ***---***---***---***---***---***---***---***---***---***---***---***---***
 *
 *

 Tenha um bom dia! ;)
 *

   --
 *De:* ASANTOS alexandresanto...@yahoo.com.br
 *Para:* r-br@listas.c3sl.ufpr.br
 *Enviadas:* Domingo, 10 de Março de 2013 12:21
 *Assunto:* [R-br] [OFF-TOPIC] Livros de divulgação científica em
 estatística

 Prezados membros da lista,

 Saberiam indicar bons livros de divulgação científica em
 estatística ao exemplo de Uma senhora toma chá de Salsburg e O andar do
 bêbado de Mlodinow,

 Obrigado,

 -- ==
 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:alexandresanto...@yahoo.com.br
 alexandre.san...@cas.ifmt.edu.br
 ==

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] 'subsetting' de um dataframe, como fazer uma função ou loop?

[FHRB Toledo]

aggregate(), não seria o caso?

On 7 March 2013 09:36, Elias Teixeira Krainski
eliaskrain...@yahoo.com.brwrote:

 sapply(OCTOPUS_S_2003[, -(1:3)],
tapply(OCTOPUS_S_2003[,2:3], mean, na.rm=TRUE))


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Tradução para português do R Reference Card

[FHRB Toledo]

ótima notícia... !

2013/3/6 Mauricio Cardeal mcardeal2...@gmail.com

 Parabéns Conrado pela iniciativa e envolvimento e Walmes pela divulgação.
 Será de grande utilidade!
 Mauricio Cardeal
 UFBA

 Em 5 de março de 2013 18:55, Walmes Zeviani walmeszevi...@gmail.comescreveu:

 Prezados Listeiros,

 Venho divulgar a tradução do R Reference Card que um aluno meu, Conrado
 Oliveira (estatística-ufpr), fez e me passou. Fiquei entusiasmado e
 agradecido com a iniciativa e envolvimento dele que não poderia deixar de
 compartilhar o resultado dessa dedicação. Confesso que não tive tempo de
 revisar o trabalho com detalhe mas quero que se sintam à vontade para usar,
 distribuir e fazer sugestões.

 tradução Conrado Oliveira:
 http://www.leg.ufpr.br/~walmes/cursoR/guia_rapido_R.pdf
 original Tom Short:
 http://cran.r-project.org/doc/contrib/Short-refcard.pdf

 À disposição.
 Walmes.

 ==
 Walmes Marques Zeviani
 LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
 Departamento de Estatística - Universidade Federal do Paraná
 fone: (+55) 41 3361 3573
 VoIP: (3361 3600) 1053 1173
 e-mail: wal...@ufpr.br
 skype: walmeszeviani
 twitter: @walmeszeviani
 homepage: http://www.leg.ufpr.br/~walmes
 linux user number: 531218
 ==

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Mauricio Cardeal
 UFBA

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Quadrado médio de duas características QM(y1+y2)

[FHRB Toledo]

Assim, para o seu CMR, segue o código:

genotipo-gl(10,3)
rep-rep(c(1,2,3),10)
resp1-
c(249.75,227.56,368.18,615.42,756.11,481.54,610.40,534.73,488.87,587.84,627.27,594.91,358.48,341.18,409.65,530.40,644.79,619.42,592.22,588.89,627.72,605.05
,620.72,
711.55,493.93,457.53,366.90,479.98,671.50,471.44)
resp2-
 
c(604.15,527.36,711.87,548.00,59.78,464.04,529.91,436.40,277.88,186.90,478.44,454.60,480.34,238.63,537.39,430.57,352.27,256.40,563.79,862.62,521.22,420.44,
519.95,375.29,425.79,450.88,264.96,611.18,372.21,284.37)

dados-data.frame(genotipo ,rep,resp1,resp2)

ANOVA_resp1 - anova(lm(resp1 ~ genotipo, data = dados)) # SQ da resposta 1
ANOVA_resp2 - anova(lm(resp2 ~ genotipo, data = dados)) # SQ da resposta 2

ANOVA_soma - anova(lm(I(resp1 + resp2) ~ genotipo, data = dados)) # SQ da
soma

Espero ter ajudado.

att,
FH

2013/2/27 Tiago Souza Marçal tiagosouzamar...@hotmail.com

  Fernando agradeço pela ajuda e como você pediu eu estou enviando um CMR
 para um delineamento em dic:

 genotipo-gl(10,3)
 rep-rep(c(1,2,3),10)
 resp1-
 c(249.75,227.56,368.18,615.42,756.11,481.54,610.40,534.73,488.87,587.84,627.27,594.91,358.48,341.18,409.65,530.40,644.79,619.42,592.22,588.89,627.72,605.05
 ,620.72,
 711.55,493.93,457.53,366.90,479.98,671.50,471.44)
 resp2-
  
 c(604.15,527.36,711.87,548.00,59.78,464.04,529.91,436.40,277.88,186.90,478.44,454.60,480.34,238.63,537.39,430.57,352.27,256.40,563.79,862.62,521.22,420.44,
 519.95,375.29,425.79,450.88,264.96,611.18,372.21,284.37)

 dados-data.frame(gen,rep,resp1,resp2)

 Att.

 Tiago.

 --
 Date: Wed, 27 Feb 2013 20:47:01 -0300
 From: fernandohtol...@gmail.com
 To: r-br@listas.c3sl.ufpr.br
 Subject: Re: [R-br] Quadrado médio de duas características QM(y1+y2)


 Tiago,

 Somar as duas variáveis é proibido?

 Se, por exemplo for um DIC e seus dados estão em um data frame chamado
 dados, com a variável resposta chamada Y, faça:

 lm(Y ~ a, data = dados) # soma de quadrados da variável a
 lm(Y ~ b, data = dados) # soma de quadrados da variável b
 lm(Y ~ I(a + b), data = dados) # soma de quadrados da soma

 Vale destacar que sem CMR isso é apenas um chute.

 att,
 FH

 2013/2/27 Tiago Souza Marçal tiagosouzamar...@hotmail.com

  Boa noite pessoal,

 estou querendo estimar a correlação genética entre duas variáveis (y1 e
 y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM
 (a+b).

 Alguém sabe como estimar este quadrado médio para as duas variáveis no R.

 Att.

 Tiago.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___ R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de
 postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo
 reproduz�vel.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Quadrado médio de duas características QM(y1+y2)

[FHRB Toledo]

Tiago,

Somar as duas variáveis é proibido?

Se, por exemplo for um DIC e seus dados estão em um data frame chamado
dados, com a variável resposta chamada Y, faça:

lm(Y ~ a, data = dados) # soma de quadrados da variável a
lm(Y ~ b, data = dados) # soma de quadrados da variável b
lm(Y ~ I(a + b), data = dados) # soma de quadrados da soma

Vale destacar que sem CMR isso é apenas um chute.

att,
FH

2013/2/27 Tiago Souza Marçal tiagosouzamar...@hotmail.com

  Boa noite pessoal,

 estou querendo estimar a correlação genética entre duas variáveis (y1 e
 y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM
 (a+b).

 Alguém sabe como estimar este quadrado médio para as duas variáveis no R.

 Att.

 Tiago.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Pesquisar elementos dentro de uma variável

[FHRB Toledo]

Talvés seja o caso de fazer uso das expressões regulares!

http://stat.ethz.ch/R-manual/R-patched/library/base/html/regex.html

2013/2/26 Manoel Galdino mcz@gmail.com

 Não sei de onde saiu o banco anterior da resposta do Edson. Fiquei
 curioso... Em todo caso, com base apenas no que você postou,

 NEvec - (BA, AL) ## aqui vocÊ cria um vetor com as siglas dos estados
 do NE. Lembrando que AC é Acre, do Norte.

 #Assumindo que seu banco de dados se chame mydf e que a variável seja x',
 então, é preciso criar uma coluna com as siglas.

 mydf$uf - substring(mydf$x, 2,3) # pega as siglas dos estados da variável
 x, que são o segundo e terceito caracteres

 #agora você usa o subset como o Edson sugeriu

 novodf - subset(mydf, uf %in% NEvec)


 abç
 M


 2013/2/26 Edson Lira edinhoes...@yahoo.com.br

 Faça:

 novo-subset(anterior, anterior$estado %in% c('BA','AL','AC'))

 [. ]'s.
 Edson Lira
 Estatístico
 Ma-Am

 Em 26/02/2013, às 12:34, Victor Eduardo victorduc...@gmail.com
 escreveu:

 Pessoal, tem como eu pesquisar dentro de uma variável no R apenas algumas
 siglas ou números? Tenho um banco de dados com mais de 5000 observações e
 estou interessado apenas na região nordeste, ou seja, preciso fazer uma
 busca no meu banco de dados para que pegue somente essa região. A variável
 está dessa maneira:

 [1] EBA292741 EDF530011 EES320531 EMS500271
 EMS500272
   [6] EPI221101 EPR410692 ERJ330455 ERS
 ERS431491
  [11] ESC420541 ESP355030 ETO172101 MAC120021
 MAC120022
  [16] MAC120031 MAC1200340001 MAC120041 MAC1200430001
 MAL270011



 No caso, pegaria BA, AC, AL e os demais estados do nordeste. Como posso
 fazer uma busca dentro dessa variável para me retornar apenas as siglas que
 quero?



 Abraços e ótima semana!

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Manoel Galdino
 https://sites.google.com/site/galdinomcz/

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] grafico de função de autocorrelação

[FHRB Toledo]

Natália,

Meu centavo na discussão, veja a função ellipse do pacote mixtools:
http://www.inside-r.org/packages/cran/mixtools/docs/ellipse

Acho que serve para alguma coisa!

att,
FH

2013/2/22 Robert Iquiapaza rb...@ufmg.br

 Com pequenas alterações

 library(tseries)
 y = log(AirPassengers)
 radf=acf(diff(y))
 li=-sd(diff(y))*1.65
 ss=which(radf$acfli | radf$acf-li)
 library(plotrix)
 for (i in 2:length(ss))

 draw.ellipse(radf$lag[ss[i]],ifelse(radf$acf[ss[i]]0,-li,li),0.02,0.08,border=red)

 Sds

 2013/2/22 Ivan Bezerra Allaman ivanala...@yahoo.com.br:

  library(tseries)
  y = log(AirPassengers)
  acf(diff(y))
 
  library(plotrix)
  draw.ellipse(0.5,0.2,0.5,0.2)
 
 
  (S,f,P)
  Allaman
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.



 --

 Natália da Silva Martins
 Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
 Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP
 Contato: (19) 8306-4743
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] função sumary do R

[FHRB Toledo]

Por que os nomes da variáveis estão começando com números... !

2013/2/20 alanaro...@sapo.pt

 Boa tarde,
 eu estou a trabalhar com as estatisticas descritivas no R. a leitura do
 ficheiro não deu erros mas a função sumary não faz nada.
 código:
 setwd(G:/SIN/Users/Ana Rua/indicadores_ano)

 dados-read.table(**indicadores_08to12.csv,**header=TRUE,sep=;)
 head(dados)
   Clinica X2008 X2009 X2010 X2011 X2012
 1  PBarca   70%   72%   72%   76%   77%
 2   Braga   74%   75%   78%   76%   73%
 3Fafe   74%   75%   74%   72%   76%
 4Maia   70%   75%   71%   76%   80%
 5  VNGaia   80%   83%   84%   84%   81%
 6 SMFeira   73%   75%   74%   78%   75%
 summary(dados)
  Clinica   X2008X2009X2010X2011
  X2012
  Abrantes: 1   74%: 5   75%: 8   74%: 4   78%: 6   81%
  : 7
  Almada  : 1  : 4   76%: 5   79%: 4   76%: 5   76%
  : 4
  Alverca : 1   69%: 3   71%: 4  : 3   75%: 4   80%
  : 4
  Amadora : 1   70%: 3  : 3   72%: 3  : 3   75%
  : 3
  APDP: 1   71%: 3   68%: 2   73%: 3   77%: 3   77%
  : 3
  Barreiro: 1   76%: 3   70%: 2   77%: 3   72%: 2   79%
  : 3
  (Other) :31   (Other):16   (Other):13   (Other):17   (Other):14
 (Other):13
 porque aparece o X nas no nome das variaveis?
 e não percebo porque o sumary apresenta isto.
 Cumprimentos Ana Rocha

 __**_
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem 
 (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia)
 e forneça código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Resultado de prcomp()

[FHRB Toledo]

Veja a posicao loadings do seu ajuste de PCA.

lanche.pca$loadings

att,
FH

2013/2/13 Carlos Andrade prf.canto...@gmail.com

 Amigos da lista,

 Estou praticando análise de componentes principais para reduzir o número
 de variáveis para uma análise de cluster.
 Encontrei um exercício e algumas apostilas do R para análise multivariada.

 Seguem os dados do exercício, que pede para usar 2 componentes principais:

  Lista1MB213
city bread hamburger butter apples tomatoes
 1 Anchorage  70.9 135.6  155.0   63.9100.1
 2   Atlanta  36.4 111.5  144.3   53.9 95.9
 3 Baltimore  28.9 108.8  151.0   47.5104.5
 4Boston  43.2 119.3  142.0   41.1 96.5
 5   Buffalo  34.5 109.9  124.8   35.6 75.9
 6   Chicago  37.1 107.5  145.4   65.1 94.2
 7Cincinnati  37.1 118.1  149.6   45.6 90.8
 8 Cleveland  38.5 107.7  142.7   50.3 93.2
 9Dallas  35.5 116.8  142.5   62.4 90.7
 10  Detroit  40.8 108.8  140.1   39.7 96.1
 11 Honolulu  50.9 131.7  154.4   65.0 93.9
 12  Houston  35.1 102.3  150.3   59.3 84.5
 13   KansasCity  35.1  99.8  162.3   42.6 87.9
 14   LosAngeles  36.9  96.2  140.4   54.7 79.3
 15Milwaukee  33.3 109.1  123.2   57.7 87.7
 16  Minneapolis  32.5 116.7  135.1   48.0 89.1
 17  NewYork  42.7 130.8  148.7   47.6 92.1
 18 Philadelphia  42.9 126.9  153.8   51.9101.5
 19   Pittsburgh  36.9 115.4  138.9   43.8 91.9
 20  StLouis  36.9 109.8  140.0   46.7 79.0
 21 SanDiego  32.5  84.5  145.9   48.5 82.3
 22 SanFrancisco  40.0 104.6  139.1   59.2 81.9
 23  Seattle  32.2 105.4  136.8   54.0 88.6
 24   Washington  31.8 116.7  154.8   57.6 86.6

 dados - Lista1MB213[,2:6]

 standardise - as.data.frame(scale(dados))

 lanche.pca - prcomp(standardise)

  summary(lanche.pca)
 Importance of components:
PC1PC2PC3PC4 PC5
 Standard deviation 1.5591 0.9654 0.9134 0.7252 0.52613
 Proportion of Variance 0.4862 0.1864 0.1669 0.1052 0.05536
 Cumulative Proportion  0.4862 0.6726 0.8394 0.9446 1.0

 A minha dúvida é como distribuir as cinco variáveis no dois componentes
 PC1 e PC2 que juntos explicam 67,26% da
 variância. Como seria a equação de combinação linear?
 Há muito tempo li em um tutorial do SPSS a solução, mas não lembro onde
 encontrá-lo.
 Agradeço antecipadamente qualquer ajuda.


 --
 Atenciosamente,

 Prof. Carlos A. S. de Andrade
 LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de
 Produção
 Universidade Federal de Campina Grande.
  Centro de Humanidades
 Unidade Acadêmica de Economia

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] plot CCA

[FHRB Toledo]

Veja o resultado de

 str(vare.cca)

Provavelmente é uma lista... aí você observa nessa lista, se tem alguma
posicao que se chama sites, dentro dessa posição deve haver uma matriz,
cujos nomes das linhas devem ser esses aí (row1, row2, row3, ...)

Aí é só você mudar o nome das linhas dessa matriz e boa!

Reportenos seus resultados!

Lembrando que esse é apenas um chute, bem alto por sinal, afinal você não
nos forneceu um CMR!

att,
FH

2013/2/6 Renata Carvalho renataecol...@gmail.com

 Bom dia, estou tentando fazer um plot de CCA mas não consigo colocar
 símbolos como referência aos sites. O R renomeou cada site para row1, row2,
 row3, Fui orientada a transformar cada site em número para depois atribuir
 os símbolos mas não estou conseguindo. Até agora tenho utilizado estes
 scripts:

 as.numeric (dados$CODE)

 plot(vare.cca, type=n)

 text(vare.cca, dis=cn)

 text(vare.cca, display = sites, col=blue, cex=0.5)  # aqui os sites
 aparecem no plot mas com o nome de row1, row2...

 Alguém sabe como posso fazer isso?

 obrigada

 Renata

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] plot CCA

[FHRB Toledo]

 ...
   .. ..- attr(*, dimnames)=List of 2
   .. .. ..$ : chr [1:180] row1 row2 row3 row4 ...
   .. .. ..$ : chr [1:4] ricFlor abFlor ricGen abGen
  $ method : chr cca
  $ inertia: chr mean squared contingency coefficient
  $ terms  :Classes 'terms', 'formula' length 3 Y ~ Scale + CODE +
 AWMSI + NumP + dIICM1 + CA
   .. ..- attr(*, variables)= language list(Y, Scale, CODE, AWMSI, NumP,
 dIICM1, CA)
   .. ..- attr(*, factors)= int [1:7, 1:6] 0 1 0 0 0 0 0 0 0 1 ...
   .. .. ..- attr(*, dimnames)=List of 2
   .. .. .. ..$ : chr [1:7] Y Scale CODE AWMSI ...
   .. .. .. ..$ : chr [1:6] Scale CODE AWMSI NumP ...
   .. ..- attr(*, term.labels)= chr [1:6] Scale CODE AWMSI NumP
 ...
   .. ..- attr(*, order)= int [1:6] 1 1 1 1 1 1
   .. ..- attr(*, intercept)= int 1
   .. ..- attr(*, response)= int 1
   .. ..- attr(*, .Environment)=environment: R_GlobalEnv
  $ terminfo   :List of 3
   ..$ terms  :Classes 'terms', 'formula' length 2 ~Scale + CODE + AWMSI +
 NumP + dIICM1 + CA
   .. .. ..- attr(*, variables)= language list(Scale, CODE, AWMSI, NumP,
 dIICM1, CA)
   .. .. ..- attr(*, factors)= int [1:6, 1:6] 1 0 0 0 0 0 0 1 0 0 ...
   .. .. .. ..- attr(*, dimnames)=List of 2
   .. .. .. .. ..$ : chr [1:6] Scale CODE AWMSI NumP ...
   .. .. .. .. ..$ : chr [1:6] Scale CODE AWMSI NumP ...
   .. .. ..- attr(*, term.labels)= chr [1:6] Scale CODE AWMSI
 NumP ...
   .. .. ..- attr(*, order)= int [1:6] 1 1 1 1 1 1
   .. .. ..- attr(*, intercept)= int 1
   .. .. ..- attr(*, response)= num 0
   .. .. ..- attr(*, .Environment)=environment: R_GlobalEnv
   ..$ xlev   :List of 1
   .. ..$ CODE: chr [1:4] CB CO CP FR
   ..$ ordered: Named logi [1:6] FALSE FALSE FALSE FALSE FALSE FALSE
   .. ..- attr(*, names)= chr [1:6] Scale CODE AWMSI NumP ...
  - attr(*, class)= chr cca

 2013/2/6 FHRB Toledo fernandohtol...@gmail.com

 Veja o resultado de

  str(vare.cca)

 Provavelmente é uma lista... aí você observa nessa lista, se tem alguma
 posicao que se chama sites, dentro dessa posição deve haver uma matriz,
 cujos nomes das linhas devem ser esses aí (row1, row2, row3, ...)

 Aí é só você mudar o nome das linhas dessa matriz e boa!

 Reportenos seus resultados!

 Lembrando que esse é apenas um chute, bem alto por sinal, afinal você não
 nos forneceu um CMR!

 att,
 FH

 2013/2/6 Renata Carvalho renataecol...@gmail.com

  Bom dia, estou tentando fazer um plot de CCA mas não consigo colocar
 símbolos como referência aos sites. O R renomeou cada site para row1, row2,
 row3, Fui orientada a transformar cada site em número para depois atribuir
 os símbolos mas não estou conseguindo. Até agora tenho utilizado estes
 scripts:

 as.numeric (dados$CODE)

 plot(vare.cca, type=n)

 text(vare.cca, dis=cn)

 text(vare.cca, display = sites, col=blue, cex=0.5)  # aqui os sites
 aparecem no plot mas com o nome de row1, row2...

 Alguém sabe como posso fazer isso?

 obrigada

 Renata

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] eliminar linhas específicas de um data.frame

[FHRB Toledo]

Aline,

faça:

 dados[which(dados$comprador == 1), ] -  NULL

att,
FH

2013/2/1 Aline Lipsky line_l...@hotmail.com



 Boa tarde,

 Tenho um data.frame similar a esse com mais de 6000 obs.:

data local comprador  29/02/2012RS 0
 29/02/2012   RS 0  29/02/2012   RS 0  29/02/2012
RS 0  29/02/2012   RS 1  29/02/2012RS 1
 29/02/2012   RS 1  27/02/2012   RS 0  27/02/2012   RS
   1  27/02/2012   RS 0  27/02/2012   RS 0  27/02/2012
RS 1  27/02/2012   RS 1  27/01/2012  RS 0
 28/01/2012   RS 0  27/01/2012   RS 0  26/01/2012   RS
   0  26/01/2012   RS 0
 Preciso eliminar as linhas onde a variável comprador é 1. Me interessam
 apenas os dados do comprador 0.

 Agradeço desde já  a atenção,

 Att.,

 Aline

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Comando que mantém a ordem natural dos dados

[FHRB Toledo]

Gisele,

Onde está declarado y?

Leia o guia de postagem (link no rodapé da mensagem).

att,
FH

2013/1/22 Gisele Correia giselem...@gmail.com

 Boa tarde,

 Gostaria de saber como plotar dados mantendo a ordem das variáveis. O R
 está ordenando automaticamente.

 exemplo:
 x - c(2,1,4,2)
 plot(x,y)
 (aqui o x é mostrado na ordem 1,2,3 e 4 no gráfico).

 Desde já agradeço, obrigada!!
 --
 Gisele Correia
 Engenheira de Pesca
 PPG-CIPET/UFAM

  Agradeço por todas as dificuldades que enfrentei (e ainda enfrento)...
 se não fosse por elas, não teria saído do lugar. A facilidade nos impede de
 caminhar. Mesmo as críticas nos auxiliam muito... Chico Xavier

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] encontrar num banco a linha de um nome

[FHRB Toledo]

use which()

2013/1/11 Fátima Lima Paula fatima.lima.pa...@gmail.com

 Consegui usando subset

 Em 11 de janeiro de 2013 12:57, Fátima Lima Paula 
 fatima.lima.pa...@gmail.com escreveu:

 Acho que estou surtando. Não consigo achar o indivíduo pelo nome em um
 data frame.
 Tenho o nome dele e quero achar a linha em que ele está.
 Alguém poderia me ajudar, por faor?
 Obrigada
  Fátima



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] R para Linux

[FHRB Toledo]

Roda!

2012/12/22 Gilbert Queiroz gilbert_quei...@yahoo.com.br

 Pessoal,
 R para Linux, roda legal???
 Abs.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Guardar objetos dentro de uma função

[FHRB Toledo]

Não sei se estou falando besteira!

Mas, acho que seria o caso de definir uma classe para esse caso!

Procede?

att,
FH

2012/12/13 Leandro Marino leandromar...@leandromarino.com.br:
 Caros,

 boa tarde.

 Gostaria de tirar um dúvida. Estou preparando uma função básica e ela tem 3
 objetos distintos de saída sendo um deles o resumo para o usuário e outros a
 informação linha por linha do arquivo de dados de entrada.

 O que gostaria de fazer é colocar estes dois últimos objetos escondidos na
 função de tal forma que quando o usuário rodasse só seria apresentado o
 resumo. Seria algo equivalente ao que o lm() faz. Se rodo o lm() só aparece
 o resumo, mas posso extrair os resíduos com o comando lm()$residuals .

 y - rnorm(100)
 x - rbinom(100,10,prob=.2)
 obj - lm(y~x)
 names(obj)


 []s,
 Leandro

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
 mínimo reproduzível.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Implementação Função

[FHRB Toledo]

Senhores,

Dia ou outro me vejo com a necessidade de aplicar o teste de Bartlett
(função bartllet.test(), do pacote stats).

Em um desses usos regulares notei certo comportamento anômalo na mesma. O
uso que dou a ela era para testar a homogeneidade de variâncias entre
modelos ajustados, ou seja, aplicava-a sobre uma lista de saídas da função
lm(), da forma:

bartlett.test(list(x = lm(...), z = lm(...), w = lm(...), y = lm(...))), De
acordo com a sintaxe apresentada no help() da mesma (
http://stat.ethz.ch/R-manual/R-patched/library/stats/html/bartlett.test.html
).

Para investigar o porque desse comportamento anômalo, além do help()
indicado no link, busquei a citação do help, de autoria próprio Bartlett (
http://www.jstor.org/stable/96803). E a fonte da função disponível no
tar.gz do R 2.15.2, que baixei do CRAN (
http://cran.r-project.org/src/base/R-2/R-2.15.2.tar.gz, acesso hoje 10 de
dezembro de 2012).

Transcrevo no final da mensagem um trecho da fonte, que se encontra em:
~/R-2.15.2/src/library/stats/R/bartlett.test.R.

Pois bem, lendo a código fonte, o artigo do Bartlett e consultando o Livro
Steel, Torrie and Dickey Principles and Procedures of Statistics: A
biometrial Approach, terceira edição, capítulo 19, página 480. Pude
perceber que todo o fundamento do teste aplica-se a variâncias. Entretanto
ao observar a implementação da função (linha destacada) observei que ao
capturar sum(obj$residuals^2), a função estaria atendo-se a soma de
quadrados do resíduo, e não a variância. Obviamente, uma soma de quadrados,
dividido pelo seu respectivo grau de liberdade é uma quadrado médio e por
sua vez uma variância, entretanto, isso não é realizado!

Pelos testes que fiz, ao utilizar a função com diferentes modelos ajustados
com o mesmo grau de liberdade dos resíduos, pela proporcionalidade das
somas de quadrados o resultado do teste é satisfatório, entretanto, para
diferentes modelos, com diferentes graus de liberdades o resultado do teste
não confere com o obtido calculado a mão!

Diante do exposto, gostaria de ouvir a opinião dos senhores quanto a esse
fato, estou certo? Em caso positivo, qual seria um meio elegante de entrar
em contato com o R CORE TEAM no sentido de apresentar esse fato?

Agradeço qualquer explicação, sugestão, etc.

att,
FH

bartlett.test.default -
function(x, g, ...)
{
LM - FALSE
if (is.list(x)) {
if (length(x)  2L)
stop('x' must be a list with at least 2 elements)
DNAME - deparse(substitute(x))
if (all(sapply(x, function(obj) inherits(obj, lm
LM - TRUE
else
x - lapply(x, function(x) x - x[is.finite(x)])
k - length(x)
}
else {
if (length(x) != length(g))
stop('x' and 'g' must have the same length)
DNAME - paste(deparse(substitute(x)), and,
   deparse(substitute(g)))
OK - complete.cases(x, g)
x - x[OK]
g - factor(g[OK])
k - nlevels(g)
if (k  2)
stop(all observations are in the same group)
x - split(x, g)
}

if (LM) {
n - sapply(x, function(obj) obj$df.resid)
v - sapply(x, function(obj) sum(obj$residuals^2))  *
DESTAQUE * DESTAQUE * DESTAQUE *
} else {
n - sapply(x, length) - 1
if (any(n = 0))
stop(there must be at least 2 observations in each group)
v - sapply(x, var)
}

n.total - sum(n)
v.total - sum(n * v) / n.total
STATISTIC - ((n.total * log(v.total) - sum(n * log(v))) /
  (1 + (sum(1 / n) - 1 / n.total) / (3 * (k - 1
PARAMETER - k - 1
PVAL - pchisq(STATISTIC, PARAMETER, lower.tail = FALSE)
names(STATISTIC) - Bartlett's K-squared
names(PARAMETER) - df

RVAL - list(statistic = STATISTIC,
 parameter = PARAMETER,
 p.value = PVAL,
 data.name = DNAME,
 method = Bartlett test of homogeneity of variances)
class(RVAL) - htest
return(RVAL)
}
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Ler um arquivo grande

[FHRB Toledo]

Daniel,

Você pode mandar o link (endereço) do seu blog!

Esse é um assunto que interessa.

Obrigado.

att,
FH

2012/12/6 Daniel Marcelino dmsilva...@gmail.com

 Sim Benilton, acabei relendo o manual do função read e testei isso.
 Deu certo. Escrevi até um post sobre isso no meu blog.

 Valeu.

 Daniel


 2012/12/5 Benilton Carvalho beniltoncarva...@gmail.com:
  se vc souber quais sao as classes das suas colunas, use o truque do
  colClasses...
 
  abuse do seu OS, va' ao terminal e use:
 
  wc -l nome_do_arquivo.txt
 
  anote o numero (que e' o numero de linhas), dai' use:
 
  conteudo = read.table('nome_do_arquivo.txt', nrow=numero de linhas)
 
  vc pode ate' combinar o comando acima com o colClasses p max eficiencia.
 
 
  2012/12/5 Daniel Marcelino dmsilva...@gmail.com
 
  Caros,
  Estou tentando ler um arquivo de texto com 930 mb, mas sem sucesso desde
  ontem.
  Meu computador é multicore Macbook  8 gb memória, mas não consegue
  finalizar de ler o arquivo e congela. Alguém tem alguma ideia de
  como superar isso sem precisar abrir o arquivo e dividir manualmente o
  banco?
 
  Daniel
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
  código mínimo reproduzível.
 
 
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Ajuda com NA

[FHRB Toledo]

!is.na(dados)

2012/11/28 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Olá

 Estou tentando excluir os NAs do meu banco, já fiz de tudo, inclusive
 mandinga, mas não consigo:

  structure(list(afia3 = c(0, 1, 0, 3, 0, 0, 0, 0, 0, 0, 0, 5,
 0, 0, 0, 0, 2, 4, 0, 3, 0, 3, 0, 0, 0, 3, 2, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 5, 0, 1, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0,
 0, 1, 2, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, NA, 0, 0, 0, 0,
 2, 2, 0, 3, 0, 0, 0, 0, 5, 0, 3, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0,
 0, 0, 0, 0, 4, 5, 0, 1, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 4, 4, 0,
 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 3, 0,
 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 2, 1, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0,
 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0,
 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 2, 0, 0, 0, 0, 2, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 2, 0, 3, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 3, 2, 5, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 3, 5, 2, 1, 3, 3, 0, 0, 0, 0, 6, 0, 0, 0, 4,
 0, 0, 0, 0, 0, 3, 4, 4, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 1, 0, 0, 3, 2, 0,
 0, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0,
 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 2, 0, 0, 3, 0, 0, 0,
 0, 3, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 5, 0, 4, 0, 0, 0, 2,
 0, 0, 7, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0,
 0, 0, 0, 3, NA, 0, 0, 0, 5, 0, 0, 4, 0, 0, 0, 0, 0, 0, 5, 6,
 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 2, 3, 0, 3,
 0, 3, 0, 0, 0, 0, 0, 0, 0, 3, 3, 3, 0, 3, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 4, 6, 5, 3, 0, 0, NA, 0, 0, 5, 0, 0, 0, 4,
 0, 0, 0, 6, 0, 0, 0, 0, 3, 0, 0, NA, 0, 0, 0, 3, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 3, 0, 0, 4, 0, 0, 0, NA, 0, 0, 0, 2, 0, 0,
 0, 7, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0,
 0, 0, 5, 3, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 0, 3, 5,
 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0,
 3, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 5, 0, 0,
 3, 0, 0, 4, 0, 0, 2, 3, 0, 0, 0, 0, 2, 0, 0, 0, 6, 0, 0, 0, 0,
 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 3, 4,
 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 2, 0, 0, 0, 0, 3, 0, 0,
 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 5,
 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 4, 0, 0)

 Como posso fazer isso?

 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) 94429486/78101651 id: 123*20942



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Comando install.packages(pacote)

[FHRB Toledo]

Isso deve ser algum conflito do firewall da sua rede!

2012/11/23 Richard Santos jamesrichardsan...@gmail.com

 Pessoal,
 Eu dei o comando install.packages('pacote'), e o R me deu as mensagens:
 --- Please select a CRAN mirror for use in this session ---
 Aviso: unable to access index for repository
 http://cran.fiocruz.br/src/contrib
 Aviso: unable to access index for repository
 http://www.stats.ox.ac.uk/pub/RWin/src/contrib
 Aviso: unable to access index for repository
 http://cran.fiocruz.br/bin/windows/contrib/2.15
 Aviso: unable to access index for repository
 http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15
 Mensagens de aviso perdidas:
 1: In open.connection(con, r) :
   não foi possível conectar a 'cran.r-project.org' na porta 80.
 2: package ‘pacote’ is not available (for R version 2.15.2)

 What the hell?




 Em 23/11/12, Richard Santosjamesrichardsan...@gmail.com escreveu:
  Excelentes dicas, Flávio, obrigado pela atenção. Gostei do que foi
  falado sobre a linguagem, e sobre a instalação de pacotes local. Vou
  testar, realmente, para aprender o R aos poucos, também. Só uma
  confirmação: essas opções do Revolution e dataframe, etc, são para ler
  bases grandes também, né?
  Grato.
  Richard
 
  Em 22/11/12, Flavio Barrosflaviomargar...@gmail.com escreveu:
  - sempre prefiro a programação, mas a programação do R é mais extensa
  do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no
  SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1;
 
  Os programas em R são bem compactos na verdade (em termos de sintaxe),
  mesmo se comparados a outras linguagens como python, ruby ou mesmo o
  matlab. Vide (
 
 http://blog.revolutionanalytics.com/2012/11/which-programming-language-is-the-most-concise.html
  )
 
  Quando a usar X4 = 1 ou dados$x4 = 1 no R é possível das duas formas,
  basta
  'fixar' o conjunto de dados antes com o   attach(dados). Particularmente
  eu
  prefiro escrever x - 4 (padrão do google também!
 
 http://google-styleguide.googlecode.com/svn/trunk/google-r-style.html#assignment
  que diga-se de passagem é um usuário ativo do R.)
 
  - para grandes bases, tem-se algumas opções como Revolution, mas mesmo
  o Revolution parece que apresenta um problema que encontrei similar ao
  Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou
  formato. Eu não queria encontrar nenhuma limitação quanto a ler
  grandes bases;
 
  Na verdade isso não é necessário. No R existem algumas estruturas de
  dados
  como o matrix, o data.frame e o list. De fato, se você quiser ler os
 seus
  dados em uma matriz, todos os dados tem que ser do mesmo tipo, mas em um
  data.frame não! O formato XDF da revolution também não impõe esse tipo
 de
  limitação. Caso as variáveis tenham tamanhos diferentes ainda é possível
  ler como um list de forma que nem o tamanho das variáveis precisa ser
  igual. Mande um e-mail para o suporte da revolution que acredito que o
  problema foi sua implementação. Eles respondem rápido. Já recebei uma
  resposta no domingo depois de 20 minutos!
 
  - para cada procedimento que se quer fazer, pode ser necessário
  instalar um pacote diferente, eu tentei pedir para o administrador do
  meu trabalho que uma pasta minha fiique livre para eu instalar
  pacotes, mas não sei se isso vai ser feito.
 
  Acho que já foi falado aqui, mas não custa lembrar: não é necessário
  privilégio de administrador para instalar pacotes. Caso você não seja
  administrador basta dar um install.packages('pacote') que ele vai criar
  um
  biblioteca local de pacotes. Esse não é um problema.
 
  No mais, é normal para aprender qualquer software ter que gastar
  tempo, mas as questões que citei acima se quiserem podem considerar
  como gosto pessoal mesmo. Mas acho importante eu citar essas questões
  para os empregadores saberem, pelo menos. Pois creio que essas
  questões não permitem ainda dizer para todos que o R é uma opção fácil
  de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não
  que não seja ótimo. Apenas acho que ainda há argumentos para que uma
  empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS,
  quando for necessário.
 
  Na minha opinião existem casos e casos para migrar do SAS ou SPSS para o
  R.
  Uma questão importante é o treinamento: se a empresa já tem
 profissionais
  bem treinados usando SAS, muito código pronto que seria difícil migrar a
  melhor solução é continuar a usar essas ferramentas. Agora no caso de um
  consultor, uma startup ou uma empresa pública que vai optar por um
 pacote
  estatístico acredito que o R é sim uma melhor opção.
 
  O R é muito utilizado no Google por exemplo. Se é bom para o Google dá
  para
  usar em qualquer empresa.
 
  Não advogo todo mundo usar o R, mas a linguagem tem muitas
  características
  interessantes, é bem intuitiva para quem programa em outras linguagens
  (python, C, Java por ex) e tem custo muito menor que o SAS ou SPSS. E é
  open source.
 
  No seu caso 

Re: [R-br] Comando install.packages(pacote)

[FHRB Toledo]

Mas você colocou mesmo 'pacote'? Esse é mesmo uma biblioteca que existe?

Veja no CRAN a respeito da biblioteca que você pretende usar!

attm
FH

On 23 November 2012 11:16, Richard Santos jamesrichardsan...@gmail.comwrote:

 Resolvi, mas aí descobri que package ‘pacote’ is not available (for R
 version 2.15.2), rs.

 2012/11/23, Benilton Carvalho beniltoncarva...@gmail.com:
  setInternet2(TRUE)
 
  2012/11/23 Richard Santos jamesrichardsan...@gmail.com
 
  Pessoal,
  Eu dei o comando install.packages('pacote'), e o R me deu as mensagens:
  --- Please select a CRAN mirror for use in this session ---
  Aviso: unable to access index for repository
  http://cran.fiocruz.br/src/contrib
  Aviso: unable to access index for repository
  http://www.stats.ox.ac.uk/pub/RWin/src/contrib
  Aviso: unable to access index for repository
  http://cran.fiocruz.br/bin/windows/contrib/2.15
  Aviso: unable to access index for repository
  http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15
  Mensagens de aviso perdidas:
  1: In open.connection(con, r) :
não foi possível conectar a 'cran.r-project.org' na porta 80.
  2: package ‘pacote’ is not available (for R version 2.15.2)
 
  What the hell?
 
 
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Seminários dos Usuários de R no Brasil

[FHRB Toledo]

Concordo com o Leandro, acho que seria a hora de preparar um fomulário e
identificar o perfil de nós usuários da lista e também quais expectativas
temos quanto a um evento dessa natureza!

2012/11/23 Leandro Marino leandromar...@leandromarino.com.br

 Considerando a partipação de muitas pessoas aqui da lista, o que vocês
 acham de criarmos um formulário indicando quem poderia ajudar a organizar e
 quem participaria?

 Assim depois poderiamos fazer uma chamada para definir possiveis
 cidades-sedes para um evento nacional, o que acham??

 um abraço,
 Leandro

 Enviado de um dispositivo móvel.
 Sent from a mobile device.
 On Nov 23, 2012 3:02 PM, Simone D. Sartorio sisarto...@yahoo.com.br
 wrote:

 Walmes e Meninos!
 Como sempre, podem contar comigo p o q der e vier!
 ajudo no que precisarem!
 bjaoo
 Si


 *
 ***---***---***---***---***---***---***---***---***---***---***---***---***
 *
 **  Simone Daniela Sartorio*
 *   Professora Adjunta da UNIFEI;
 * Doutora e Mestre em Estatística e Experimentação Agronômica -
 ESALQ/USP;
 *   Licenciada em Matemática - UNESP/Rio Claro.
 *
 ***---***---***---***---***---***---***---***---***---***---***---***---***
 *
 *

 Tenha um bom dia! ;)
 *

   --
 *De:* FHRB Toledo fernandohtol...@gmail.com
 *Para:* r-br@listas.c3sl.ufpr.br; Ari Clecius ari072...@yahoo.com.br
 *Enviadas:* Sexta-feira, 23 de Novembro de 2012 9:22
 *Assunto:* Re: [R-br] Seminários dos Usuários de R no Brasil

 Senhores,

 Parabenizo ao pessoal pela iniciativa, me disponibilizo para qualquer
 apoio que se faça necessário e estiver ao meu alcance!

 att,
 FH

 2012/11/23 Ari Clecius ari072...@yahoo.com.br

 Acho a idéia muito interessante seria um momento de apresentar trabalhos,
 dúvidas e trocar experiências, sou um usuário do R a pouco tempo, em torno
 de um ano, aos poucos vejo que é o melhor sofware estatístico que já usei
 dentre vários, pirouette, Statistica, Minitab, SAS, SPSS.Nem sempre a
 entrada de dados é a mais fácil, mas nos permite estudar e entender o que
 está acontecendo, em outros softwares a análise ocorre de forma
 automática. Sou aluno de Doutorado em Engenharia Civil área de concentração
 Saneamento ambiental, gosto de estudar estatística e me disponibilizo, caso
 seja efetivado a idéia do encontro a tentar trazê-lo para fortaleza, tenho
 certeza que terei o apoio aqui na UFC.

 Att,

 Ari Clecius Alves de Lima
 Engenheiro Químico
 Me. Engenharia Civil
 (085)87687786
 (085)33669042
 *De:* Rodrigo Sant'Ana rodrigo.gri...@gmail.com
 *Para:* r-br@listas.c3sl.ufpr.br
 *Enviadas:* Quinta-feira, 22 de Novembro de 2012 20:39

 *Assunto:* Re: [R-br] Seminários dos Usuários de R no Brasil

 Acho a ideia super pertinente e interessante,
 coloco-me a disposição para ajudar no que precisar.

 E parabenizo-os, Fábio e Walmes, pela iniciativa.
 Com certeza ganharíamos muito com um evento deste.

 ___
 Rodrigo Sant'Ana
 Oceanógrafo


 Em 22 de novembro de 2012 20:36, Fabrício Barth fabricio.ba...@gmail.com
  escreveu:

 Estou utilizando o R a pouco tempo e estou adorando. Gostaria de ver um
 evento aqui no Brasil. :)


 2012/11/22 Leandro Marino leandromar...@leandromarino.com.br

 No que eu puder ajudar estou dentro. :-)
 Tentei fazer no SINAPE 2008 mas mesmo tendo feito a adesão durante o
 evento foi de 4 ou 5 pessoas...

 Enviado de um dispositivo móvel.
 Sent from a mobile device.
 Em 22/11/2012 20:23, Fabio Mathias Corrêa fabio.u...@yahoo.com.br
 escreveu:

  Nobres colegas,

 Em diversos países os usuários do R promovem seminários de 2 ou 3 dias
 para apresentarem seus trabalhos relacionados ao R. No RBras que ocorreu em
 São Carlos realizamos um encontro dos usuários do R bem descontraído de
 forma a conhecermos as pessoas da lista e demais usuários que não
 participam da lista. Não lembro de terem feito este encontro nos eventos
 posteriores.

 Ao conversar com o Walmes hoje, acreditamos que temos potencial para
 expandirmos e promovermos algo mais elaborado. Claro que este tipo de
 evento depende da colaboração de todos!

 Por isso, venho através da lista saber do interesse dos usuários sobre
 este tipo de evento e assim podemos começarmos a planejar.

 Valeu!


 Fábio Mathias Corrêa

Universidade Estadual de Santa Cruz
 Departamento de Ciências Exatas e da Terra - DCET


 Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna
 CEP 45662-900. Ilhéus-Bahia


 Tel.: 73-3680-5076

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Fabrício J. Barth
 http

[R-br] [off topic] Adobe!

[FHRB Toledo]

Senhores,

Esse é um thread totalmente OFF TOPIC, peço a compreensão de vocês quanto a
isso!

Um colega da lista brasileira do DEBIAN (
debian-user-portugu...@lists.debian.org) encabeça essa iniciativa. Acho
extremamente válida, quem for da mesma opinião, por favor assine.

Leiam o email abaixo!

att,
FH

-- Forwarded message --
From: Juliano Cesar cetapsao...@yahoo.com.br
Date: 2012/11/18
Subject: Assinem a petição para a Adobe sempre fazer a versão mais recente
do Adobe Reader para Linux.
To: debian-user-portugu...@lists.debian.org 
debian-user-portugu...@lists.debian.org


Vamos divulgar e assinar
http://www.change.org/pt-BR/peti%C3%A7%C3%B5es/adobe-systems-incorporated-we-want-to-always-have-the-latest-version-of-adobe-reader-on-linux

Se alguém tiver uma ideia de uma melhor descrição em inglês e portarmos a
descrição em vários idiomas, pode ser que isso dê certo.

A meta é de 2 milhões de assinaturas, sabemos que existem alternativas de
leitores de PDF, mas nenhum se compara ao Adobe Reader, é uma pena que a
Adobe parou de atualizar a sua versão para Linux.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] correlação de pearson

[FHRB Toledo]

Vale lembrar que correlações de Pearson podem ser obtidas/estimadas por
meio da função cor() do base do R!

Veja a documentação em ??cor

att,
FH

On 16 November 2012 09:49, Daniel C Bezerra danielcbeze...@gmail.comwrote:

 O pacote só deve estar disponível no CRAN para a versão atual do R
 (2.15.2) e não para a sua 2.15.1.

 Atualize o R e depois tente instalar novamente.

 Abs,

 D

 On Fri, Nov 16, 2012 at 9:44 AM, alanaro...@sapo.pt wrote:

 eu tentei instalar os packages para a correlação e acontece-me o seguinte:

 install.packages(ppcor)

 --- Please select a CRAN mirror for use in this session ---
 Warning: unable to access index for repository
 http://cran.fiocruz.br/bin/**windows/contrib/2.15http://cran.fiocruz.br/bin/windows/contrib/2.15
 Warning: unable to access index for repository
 http://www.stats.ox.ac.uk/pub/**RWin/bin/windows/contrib/2.15http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15
 Warning messages:
 1: In open.connection(con, r) : unable to resolve 'cran.r-project.org'
 2: package ‘ppcor’ is not available (for R version 2.15.1)

 library(ppcor)

 Error in library(ppcor) : there is no package called ‘ppcor’

 require(ppcor)

 Loading required package: ppcor
 Warning message:
 In library(package, lib.loc = lib.loc, character.only = TRUE,
 logical.return = TRUE,  :
   there is no package called ‘ppcor’

  install.packages(ds)

 Warning: unable to access index for repository
 http://cran.fiocruz.br/bin/**windows/contrib/2.15http://cran.fiocruz.br/bin/windows/contrib/2.15
 Warning: unable to access index for repository
 http://www.stats.ox.ac.uk/pub/**RWin/bin/windows/contrib/2.15http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15
 Warning message:
 package ‘ds’ is not available (for R version 2.15.1)

 require(ds)

 Loading required package: ds
 Warning message:
 In library(package, lib.loc = lib.loc, character.only = TRUE,
 logical.return = TRUE,  :
   there is no package called ‘ds’

 obrigada. An Rocha
 __**_
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem 
 (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia)
 e forneça código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Lattice Design

[FHRB Toledo]

Walmes,

Obrigado pela prontidão!

Pois bem, rodei o seu código... eis que: sua implementação confere com os
resultados esperados relativos ao ajuste das médias com recuperação da
informação interbloco (popmeans, do doBy quando o modelos é ajustado pelo
lme()).

No primeiro ajuste,, o do aov, o popmeans retorna as próprias médias
aritméticas dos tratamentos, sem qualquer ajuste!

O teste correto da soma de quadrados de tratamentos ajustados confere com
Cochran  Cox!

Infelizmente ainda não consegui recuperar as médias ajustadas intrabloco,
mas com certeza o que é posto por você ajuda e muito!

Obrigado...

abraço,
FH

2012/11/15 Walmes Zeviani walmeszevi...@gmail.com

 Fernando,

 Já rodei os exemplos do Pimentel e do Ramalho et al. Implementação de
 ambos estão disponíveis nos scripts que levei para o Curso mais recente que
 dei na Embrapa Arroz e Feijão (que terminei ontem). O link para arquivos do
 curso é esse

 http://www.leg.ufpr.br/~walmes/cursoR/cnpaf2/

 Pegue o script lati.R. Eu nunca olhei a implementação do CCox e fiquei
 curioso e imaginando será que os livros divergem? Meus códigos reproduzem o
 Pimentel. Bem, com os seus dados a análise que atualmente faço é essa

 rm(list=ls(all=TRUE)); ls()
 download - read.table('http://dl.dropbox.com/u/38195533/dados_CC.txt',
header = TRUE, # com cabecalho
sep = '\t', # separador de celulas TAB
dec = ',', # separador de decimal , (virgula)
na.string = '.') # indicador de omissao
 str(download)

 ## lendo e transformando variaveis em fator
 dados - transform(download,
trat = factor(trat), # transforma em fator
rep = factor(rep), # idem
bloco = factor(bloco)) # ...

 all(complete.cases(dados)) # completo
 xtabs(~rep+trat, dados)
 xtabs(~rep+bloco, dados)
 xtabs(~bloco+trat+rep, dados)

 str(dados) # estrutura da planilha

 m0 - lm(terms(resp~rep/bloco+trat, keep.order=TRUE), data=dados)
 anova(m0)

 require(doBy)
 popMeans(m0, effect=trat)

 require(nlme)
 dados$blocoin - with(dados, interaction(rep, bloco, drop=TRUE))
 mm0 - lme(resp~rep+trat, random=~1|blocoin, dados)
 anova(mm0)

 popMeans(mm0, effect=trat)

 À disposição.
 Walmes.

 ==
 Walmes Marques Zeviani
 LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
 Departamento de Estatística - Universidade Federal do Paraná
 fone: (+55) 41 3361 3573
 VoIP: (3361 3600) 1053 1173
 e-mail: wal...@ufpr.br
 skype: walmeszeviani
 twitter: @walmeszeviani
 homepage: http://www.leg.ufpr.br/~walmes
 linux user number: 531218
 ==

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] R versus SAS e SPSS

[FHRB Toledo]

Novamente,

Acredito que não só EU, como outros colegas da lista do R-Br ficam se
perguntando qual o REAL objetivo dessas suas manifestações aqui para nós na
lista do R-Br?

Seu primeiro email buscou fazer uma enquete de qual o melhor
programa/linguagem/ambiente para análises estatísticas ou seja lá qual seja
seu propósito... E agora esse em que, sem mais nem menos, se manifesta
dessa forma!

O que espera que respondemos?

Sem mais.

att,
FH

2012/11/16 Richard Santos jamesrichardsan...@gmail.com

 O único software estatístico livre que faz todas as coisas e com bases
 grandes (apenas não pequenas) é o R. Os únicos softwares estatísticos
 não-livres que fazem todas as coisas e com grandes bases (tb não
 apenas com pequenas) é o SPSS e o SAS (acabei encontrando limitações
 no Minitab também. O Stata atende, mas com limitações de opções
 também).

 Agora quando a instituição demonstra dificuldade em obter o SAS ou o
 SPSS (como agora, pela primeira vez encontrei uma assim), optei assim
 por adquirir o SPSS e trazer o notebook para o trabalho. Desisti de
 aprender o R, porque me dei conta de que o R você tem que programar
 para tudo, até para fazer gráficos, e isso até no SAS que antes você
 tinha que programar pois pelo point in click não tinha como fazer um
 gráfico a seu modo já se aboliu. Além do que a programação do R teria
 que aprender tudo agora, e não é a C que já estou acostumado. E como
 lembra a própria mensagem do R, você não tem garantia de nada. E
 para mim que até hoje não sou especialista no R, vai demorar uma
 década ainda para eu fazer tudo o que eu quero rapidamente. Demora
 naturalmente o treinamento, por ser inclusive uma linguagem de
 programação e ampla.

 Então vou ficar por pouco tempo na lista do R, até lá vai ser um
 prazer aprender, mas vou optar em breve por trazer o meu notebook com
 o SPSS comprado. Não vou ter tempo para aprender tudo o que preciso a
 contento (até lá é possível até que tenha trocado por algum motivo de
 instituição - para uma que tem o SPSS).

 Obrigado.
 Abraços,
 Att.
 Richard
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Lista

[FHRB Toledo]

Thor,

Seja bem vindo a lista!

Como você está iniciando no programa sempre leia a documentação da função
que pretende usar, no caso da glm() digite no prompt do seu R help(glm).
Será mostrado para você um arquivo que segue uma estrutura padrão entre
todas as funções do R e é um excelente ponto de partida!

Quanto a um exemplo simples do uso dela minha sugestão é por meio da função
demo().

Digite no seu prompt demo(lm.glm, package = stats) que será mostrado um
bom exemplo de dados simples primeiro ajustados pela lm e depois pela glm...

Qualquer outra dúvida, fico a disposição.

att,
FH

2012/11/16 Narede Golaco naredegol...@hotmail.com

 Olá gostaria de postar uma dúvia, relacionada a função glm:

 Se possivel gostaria de um exemplos simples de analise de variancia usando
 a função glm em conjunto com qualquer pacote de dados do R, pode ser o
 mtcars, airqualiy...

 Gostaria muito mesmo do exemplo, pois estou iniciado no programa e só
 encontro na internet exemplos muito avançados.


 Agradeço desde já a ajuda de todos.

 Se eu puder postar esta dúvida ou não, por favor gostaria de ser avisado


 Obrigado a todos pela condescêndencia.


 Att: Thor Cerqueira
 Graduando em matematica



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Lattice Design

[FHRB Toledo]

Senhores,

Estou com uma dúvida com relação a reprodução de uma análise de um
experimento em delineamento em lattice quadrado balanceado, tipo I, segundo
Pimentel Gomes (cap 10).

Escreve esse email baseado tanto no livro do Pimentel Gomes, como também no
Cochran e Cox, cap 10, pag 274. Para quem interessar segue também um
exemplo (teste) do livro Ramalho et al., Experimentação em Genética ... cap
11, pág 171 (mas não é estelionato, rsrs).

Pois bem, como disse, trata-se da análise de blocos incompletos em lattice
quadrado balanceado. A pergunta é: como obter a soma de quadrados de
tratamentos ajustadas aos efeitos de blocos? E ainda, de que forma obter
também as médias de tratamentos livre dos efeitos de blocos.

Meu ponto de partida é o post do Walmes no wiki do LEG (
http://www.leg.ufpr.br/doku.php/ridiculas#analise_intrabloco_e_interbloco_de_experimento_em_blocos_incompletos_tipo_i).
Todavia, minha implementação não concorda com o que lá é mostrado.

Peço desculpas pela forma que procedi, mas a título didático me prendi
exatamente ao modo como é originalmente apresentado em Cochran e Cox (pag
274) e seguindo o que lá é feito consegui reproduzir ipsis literis a
análise.

Assim, gostaria da ajuda de vocês para aprender um modo menos custoso de
reproduzir o jeito original!

Desde já agradeço qualquer ajuda/sugestão com relação a dúvida! Caso alguém
se interesse disponibilizo outro conjunto de dados, da mesma natureza no
mesmo lik, com o nome de dados_magno.txt.

att,
FH

Segue o CMR:

rm(list = ls(all = TRUE)); ls()

download - read.table('http://dl.dropbox.com/u/38195533/dados_CC.txt',
   header = TRUE, # com cabecalho
   sep = '\t', # separador de celulas TAB
   dec = ',', # separador de decimal , (virgula)
   na.string = '.') # indicador de omissao

## lendo e transformando variaveis em fator
dados - transform(download,
   trat = factor(trat), # transforma em fator
   rep = factor(rep), # idem
   bloco = factor(bloco)) # ...

str(dados) # estrutura da planilha

## fator bloco dentro de repeticoes
dados$blocoDrep - factor(with(dados, paste(rep, bloco, sep = '')))

## niveis de cada fator **apenas auxiliar**
r - nlevels(dados$rep)
k - nlevels(dados$bloco)

## array do experimento **generalizacao da matrix do delineamento**
X - with(dados, tapply(resp, list(rep, blocoDrep, trat), mean, na.rm =
TRUE))

## totais dos fatores
R - apply(X, 1, sum, na.rm = TRUE) # repeticoes
B - apply(X, 2, sum, na.rm = TRUE) # blocos dentro de repeticoes
T - apply(X, 3, sum, na.rm = TRUE) # tratamentos

## somas de quadrados do 'jetinho' Pimentel Gomes
Sx - sum(X, na.rm = TRUE)
Sx2 - sum(X^2, na.rm = TRUE)
C - Sx^2 / sum(!is.na(X))
SQtotal - Sx2 - C
SQr - ( 1 / k^2 ) * sum( R^2 ) - C
SQt - ( 1 / r ) * sum( T^2 ) - C
SQb.unadj - ( 1 / k ) * sum( B^2 ) - C
SQe.DBC - SQtotal - (SQr + SQt)

## totais dos blocos onde os tratamentos ocorrem... ver Cochran e Cox (cap
10, pag 274)
X.tb - apply(X, c(2, 3), function(x) sum(!is.na(x)))
Bt - t(X.tb) %*% B
W - k * T - (k + 1) * Bt + sum(T)

## outras somas de quadrados... blocos ajustados e erros
SQb.adj - sum(W^2) / ( k^3 * ( k + 1 ) )
SQe.intra - SQe.DBC - SQb.adj
Eb - SQb.adj / (r * ( k - 1 ))
Ee - SQe.intra / ( ( k - 1 ) * ( r * k - k - 1 ) )

## coeficiente 'mu' determina se existe ou nao ajuste das medias
if( Eb = Ee ) mu - 0 else mu - ( Eb - Ee ) / ( k * ( r - 1 ) * Eb )

Eel - Ee * ( 1 + ( r * k * mu ) / ( k + 1 ))
EDBC - SQe.DBC / ( ( r - 1 ) * ( k^2 - 1 ) )

## soma de quadrados de tratamentos ajustados **PROBLEMA**
SQt.adj - sum((T + mu * W - mean(T + mu * W))^2) / r

## montando quadro da ANOVA
DofF - c( (r - 1 ), ( k^2 - 1 ), ( k^2 - 1 ), ( r * ( k - 1 ) ), ( ( k - 1
) * ( r * k - k - 1 ) ), ( ( r - 1 ) * ( k^2 - 1 ) ), ( ( k - 1 ) * ( r * k
- k - 1 ) ), ( sum(!is.na(X)) - 1 ))
SSs - c(SQr, SQt, SQt.adj, SQb.adj, SQe.intra, SQe.DBC, Eel * (( k - 1 ) *
( r * k - k - 1 )), SQtotal)
table - cbind(DofF, SSs, SSs / DofF)
dimnames(table) - list(c('Rep', 'Trat. unadj.', '  Adj.', 'Bloco
Adj.', 'Error Intra-bloco', '  DBC', '  Effective', 'Total'),
c('Df', 'Sum Sq', 'Mean Sq'))

## teste F para tratamentos ajustados
F.trat - table[3, 3] / table[7, 3]
p.trat - pf(F.trat, table[3, 1], table[7, 1], lower.tail = FALSE)

## ajuste via aov
ajuste - aov(terms(resp ~ rep/bloco + trat, keep = TRUE),
  contrasts = list(rep = contr.sum, bloco = contr.sum, trat =
contr.sum),
  data = dados)

## comparando resultados

## ANOVAS
table # anova 'made by hands'
c('F value' = F.trat, 'Pr(F)' = p.trat) # teste F 'made by hands'
summary(ajuste) # anova via aov()

## sugestao/implementacao do Walmes (Wiki do LEG)
car::Anova(ajuste, type = 'III')
drop1(ajuste, scope = .~., test = 'F')

## medias de tratamentos ajustadas
data.frame(Unadj. = T / r,
   Adj. = (T + mu * W) / r,
   ## sugestao/implementacao do Walmes (Wiki do LEG)
  

Re: [R-br] Duvida, como exibir help-pages do R no navegador usando Ubunto-linux e emacs?

[FHRB Toledo]

Augusto,

Para você customizar a aparência do seu R veja o post:
http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/

att,
FH

2012/11/1 Augusto Ribas ribas@gmail.com

 Muito obrigados a todos, funcionou perfeitamente :)

 Eu achei diferente, por que no windows o help em html é o default, e no
 linux é diferente.
 Mas criei o Rprofile e coloquei os comandos e agora esta como eu gostaria.

 Em 1 de novembro de 2012 01:03, Fernando Mayer 
 fernandoma...@gmail.comescreveu:

 O .Rprofile não é criado automaticamente. Você deve criá-lo se quiser
 customizar algumas opções do R, como carregar sempre os mesmos
 pacotes, ajustar opções, etc.

 O .Rprofile criado no seu $HOME será global, ou seja, em qualquer
 diretório que vc iniciar o R, ele irá ler o $HOME/.Rprofile.

 Mas vc também pode ter ter um .Rprofile em diretórios específicos, com
 opções mais específicas. Por exemplo para sempre carrgar um pacote
 usado em uma análise que está em um diretório, mas não carregá-lo
 sempre, globalmente. A ordem de leitura será sempre o .Rprofile do
 diretório (se houver), e depois do $HOME.

 O qe vc quer é colocar essa linha que o Walmes mencionou

 options(help_type=html)

 no seu $HOME/.Rprofile para ele ficar global.

 []s,

 ---
 Fernando Mayer
 Universidade Federal de Santa Catarina - UFSC
 Departamento de Ecologia e Zoologia - ECZ/CCB
 URL: http://fernandomayer.github.com
 e-mail: fernandomayer [@] gmail.com


 2012/10/31 Daniel C Bezerra danielcbeze...@gmail.com:
  Outra sugestão é iniciar o help do browser e depois fazer as pesquisas
 por
  lá.
 
  Use:
  help.start()
 
  E depois clique em:  Search Engine  Keywords
 
  Abs,
 
  D
 
 
  2012/10/31 Walmes Zeviani walmeszevi...@gmail.com
 
  help(plot, help_type=html)
 
  Você pode controlar globalmente alterando o options(), algo assim
 
  options(help_type=html)
  help(lm)
 
  À disposição.
  Walmes.
 
 
 ==
  Walmes Marques Zeviani
  LEG (Laboratório de Estatística e Geoinformação, 25.450418 S,
 49.231759 W)
  Departamento de Estatística - Universidade Federal do Paraná
  fone: (+55) 41 3361 3573
  VoIP: (3361 3600) 1053 1173
  e-mail: wal...@ufpr.br
  skype: walmeszeviani
  twitter: @walmeszeviani
  homepage: http://www.leg.ufpr.br/~walmes
  linux user number: 531218
 
 ==
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
  código mínimo reproduzível.
 
 
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 Grato
 Augusto C. A. Ribas

 Site Pessoal: 
 http://recologia.blogspot.com.br/http://augustoribas.heliohost.org
 Lattes: http://lattes.cnpq.br/7355685961127056


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] mclapply

[FHRB Toledo]

Clara,

Dois exemplos:

http://ridiculas.wordpress.com/2011/10/23/paralelizacao-de-processos/

http://ridiculas.wordpress.com/2011/11/28/paralelizacao-de-processos-ii/

À disposição.

att,
FH

2012/10/18 clara luz clara_santa...@yahoo.com.br

 Olá galera do R! Tudo bom?

 Gostaria de saber se alguem conhece ou já trabalhou com a função mclapply,
 a Parallel version of lapply.

 Segundo sua descrição, deve-se inserir como argumento um vetor X, sendo
 que para cada elemento de X será aplicada a função que eu desejo (FUN):
 mclapply(X, FUN, ..., mc.preschedule = TRUE, mc.set.seed = TRUE, mc.silent
 = FALSE, mc.cores = getOption(cores), mc.cleanup = TRUE)

 Porém eu gostaria de aplicar FUN a uma matriz inteira. Ou seja, eu
 precisaria que cada elemento de X fosse uma matriz... X-(matrizA, matrizB,
 matrizC, ...)
 Não sei se me fiz entender, mas será que isso é possível?

 Abraços eRRantes!

 *Clara Luz B. Sant'Anna
 *Bióloga, Bacharel e Licenciada
 Mestranda em Ecologia - Unicamp
 Laboratório de Biogeografia da Conservação - UFG
 (19) 9233.1062

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] ordenar colunas - seleção expandida

[FHRB Toledo]

Marcelo,

Creio que isso deva funcionar!

dados - dados[, order(dados$colunaX)]

att,
FH

2012/10/10 marcelo claro de souza marcelo.claro.so...@gmail.com

 Olá pessoal,
 Preciso reordenar uma coluna de um conjunto de dados, mas gostaria que
 todas as outras fossem reordenadas juntamente a primeira em estilo seleção
 expandida, assim como no excel.
 Tentei o comando sort , mas só consegui ordenar uma coluna de vez.
 Alguma dica?
 Muito obrigado.
 Abraço

 Marcelo

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] ordenar colunas - seleção expandida

[FHRB Toledo]

Paulo Justiniano e Marcelo,

Exatamente, inverti linhas e colunas, peço desculpas... Para um tutorial
completo sobre o tema, consultem:
http://www.ats.ucla.edu/stat/r/faq/sort.htm

att,
FH

2012/10/10 Paulo Justiniano paulo...@leg.ufpr.br

 Fernando, me parece que voce inverteu linhas e colunas.
 seria:

 dados - dados[order(dados$colunaX), ]



 On Wed, 10 Oct 2012, FHRB Toledo wrote:

  Marcelo,

 Creio que isso deva funcionar!

 dados - dados[, order(dados$colunaX)]

 att,FH

 2012/10/10 marcelo claro de souza marcelo.claro.so...@gmail.com**
   Olá pessoal,
   Preciso reordenar uma coluna de um conjunto de dados, mas gostaria
 que todas as outras fossem reordenadas juntamente a primeira
   em estilo seleção expandida, assim como no excel.
   Tentei o comando sort , mas só consegui ordenar uma coluna de vez.
   Alguma dica?
   Muito obrigado.
   Abraço

   Marcelo

   __**_
   R-br mailing list
   R-br@listas.c3sl.ufpr.br
   
 https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
   Leia o guia de postagem 
 (http://www.leg.ufpr.br/r-br-**guiahttp://www.leg.ufpr.br/r-br-guia)
 e forneça código mínimo reproduzível.




 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Evitar gravação

[FHRB Toledo]

Um if() não funciona?

2012/10/9 Carlos Mendonça csaes...@centroin.com.br

 Meus amigos, o meu problema é como fazer com que um arquivo não seja
 gravado se na etapa anterior do programa acontece um erro.

 No exemplo, eu não gostaria que o arquivo .csv fosse gravado pois
 aconteceu um erro na etapa anterior do programa.

 Exemplo:
  arquivo1= merge(arquivo1, arquivo99, by = code6, all.x = TRUE)
 
  arquivo1l$2004t1 = ifelse(is.na(arquivo1$i2004t1), 100,
 arquivo1$i2004t1)
 Erro: constante numérica inesperado em serie_bmovel$2004
 

  write.csv2(arquivo1, \\arquivo1.csv, row.names = FALSE)




 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Somar mais de uma variável

[FHRB Toledo]

Tente algo com o próprio aggregate... da forma:

 aggregate(arquivo[, 2:4], by = list(arquivo[,1]), FUN = sum)

Comigo vem funcionando... reporte qualquer resultado!

att,
FH

2012/10/1 Carlos Mendonça csaes...@centroin.com.br

 Caros,

 existe uma forma mais prática de somar várias variáveis por um determinado
 by, ao invés da maneira que faço atualmente?

 Arquivo de dados (arquivo)

 codigo valor1   valor2 valor3
 a  3   25 30
 a  5   10 50
 a  2   15 20
 b15   85 20
 b  5   15 10
 c  2 5 25
 c  4 7   5
 c  2 5   8
 c  2 3 12


 1) Utilizo o comando arq1 = aggregate(arq1$valor1, list(arquivo$codigo),
 sum) para somar o valor1 por cada codigo.

 2) Utilizo o comando arq2 = aggregate(arq1$valor2, list(arquivo$codigo),
 sum) para somar o valor2 por cada codigo.

 3) Utilizo o comando arq3 = aggregate(arq1$valor3, list(arquivo$codigo),
 sum) para somar o valor3 por cada codigo.

 4) Depois junto o arq1, arq2 e arq3 e fico com um arquivo igual ao abaixo:

 codigo valor1 valor2 valor3
 a 10  50100
 b 20   100 30
 c 10  20 50


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Comando for em cadeia

[FHRB Toledo]

Guilherme,

Nada te impede de aninhar laços for(). O que você chama de dvv?

Qual o resultado numérico dessa sua implementação?

att,
FH

2012/9/17 Guilherme Heiden guilhermehei...@hotmail.com

  Olá, eu talvez não tenha sido claro no assunto do e-mail, mas na minha
 concepção o problema é este: um comando for dentro de outro.

 Eu estou trabalhando com a seguinte matriz:

  [,1] [,2] [,3] [,4] [,5]
 [1,]   205   117   49
 [2,]   189   102   45
 [3,]   11   35   30   157
 [4,]   10374   26

 E a partir desta eu gostaria de gerar outra matriz com a seguinte
 configuração:

  [,1]  [,2] [,3][,4][,5]
 [1,] dv11   dv12   dv13  dv14dv15
 [2,]   -  dv22   dv23   dv24dv25
 [3,]   --   dv33   dv34   dv35
 [4,]   ---   dv44dv45
 [5,]  -- - -   dv55

 *Para achar os devidos valores, eu criei a seguinte função:*
 read.csv2(C:/multivar.csv,header=T)

 multivar-read.csv2(C:/multivar.csv,header=T)

  attach(multivar)



 *##Criar DV11:*

 vetor1-NULL

  for(i in 1:4) {

  vetor1[i]-(multivar[i,1]-multivar[i,1])^2

  }

 dv11-(sum(vetor1)^0.5)

 rm(vetor1)



 *##Criar DV12:*

 vetor1-NULL

 for(i in 1:4) {

 vetor1[i]-(multivar[i,1]-multivar[i,2])^2

 }

 dv12-(sum(vetor1)^0.5)

 rm(vetor1)


 E assim sucessivamente.


 *Minha dúvida: *

 *Há a possibilidade de criar um comando for para os valores de dv serem
 criados automaticamente, ao invés de gerar um de cada vez, manualmente?*

 *
 *

 *
 *

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Selecionar colunas de interesse apartir de uma matriz...

[FHRB Toledo]

André,

Tente:

selecao - function (x) which(x == 7)

so.selecionadas - matriz[, apply(matriz, 2, selecao)]

Sendo matriz a sua matriz apresentada.

att,
FH

2012/9/15 andrebvs andre...@bol.com.br

 Olá colegas!

 Como selecionar apenas as colunas que contenham um determinado número de
 interesse?

 Por exemplo:

 C1   C2C3C4C5C6C7
 1  2  3  4  5  6  7
 0  7  4  9  8  8  1
 0  8  5  6  6  4  1
 3  2  4  4  1  1  8
 5  5  9  8  7  5  2

 Neste exemplo, suponhamos que o valor de interesse seja 7, então, gostaria
 que a saída fosse as seguintes colunas:

 C2C5C7
 2   5  7
 7   8  1
 8   6  1
 2   1  8
 5   7  2

 thanks!

 Att.
 André

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Selecionar colunas de interesse apartir de uma matriz...

[FHRB Toledo]

Sim, sim... ato falho, vamos por partes:

Me precipitei... agora segue código testado, um pouco mais complicado, veja
se compreende a lógica:

matriz - matrix(c(1,2,3,4,5,6,7,0,7,4,9,8,8,1,0,8,5,6,6,4,1,3,2,4,4,
 1,1,8,5,5,9,8,7,5,2),,7,byrow=T)

matriz[,which(apply(matriz, 2, function(x) sum(match(x, 7), na.rm = TRUE))
== 1)]

att,
FH

2012/9/15 andrebvs andre...@bol.com.br

 Obg Fernando, mas aparece a seguinte mensagem quando executo os comandos
 abaixo: *Erro: tipo de subscrito inválido 'list'*

 *CMR:*
 * *
 matriz - matrix(c(1,2,3,4,5,6,7,0,7,4,9,8,8,1,0,8,5,6,6,4,1,3,2,4,4,
  1,1,8,5,5,9,8,7,5,2),,7,byrow=T)

 selecao - function (x) which(x == 7)
 soselecionadas - matriz[,apply(matriz,2,selecao)]



 --
 Em 15/09/2012 19:14, *FHRB Toledo  fernandohtol...@gmail.com * escreveu:
 André,

 Tente:

 selecao

 so.selecionadas

 Sendo matriz a sua matriz apresentada.

 att,
 FH

 2012/9/15 andrebvs 
 andre...@bol.com.brhttp://../../../undefined/compose?to=andre...@bol.com.br
 

 Olá colegas!

 Como selecionar apenas as colunas que contenham um determinado número de
 interesse?

 Por exemplo:

 C1   C2C3C4C5C6C7
 1  2  3  4  5  6  7
 0  7  4  9  8  8  1
 0  8  5  6  6  4  1
 3  2  4  4  1  1  8
 5  5  9  8  7  5  2

 Neste exemplo, suponhamos que o valor de interesse seja 7, então,
 gostaria que a saída fosse as seguintes colunas:

 C2C5C7
 2   5  7
 7   8  1
 8   6  1
 2   1  8
 5   7  2

 thanks!

 Att.
 André

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.brhttp://../../../undefined/compose?to=R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dúvida função no R

[FHRB Toledo]

Felipe,

Para CMR, consulte: https://gist.github.com/1088208

att,
FH

2012/9/13 Felipe Buchbinder felb...@gmail.com

 O k foi definido no início do código, fora da função, como k = 5.

 Desculpe a ignorância, mas o que é um CMR? Como eu posso fazer um para
 enviar à lista?

 Obrigado pela disposição,

 Felipe


 2012/9/13 Leonard de Assis assis.leon...@gmail.com

  Onde foi definido o k?

 []s
 Leonard de Assis
 assis dot leonard at gmail dot com

 Em 13/09/2012 12:56, Felipe Buchbinder escreveu:

 Prezados professores,

  Estou com um problema em uma função no R.
 A função é a seguinte:

   resp - function(subject,item,trait)
 {

   p.vec = NULL

   for(category in 1:k)  # k = 5
   {
 p = ogive ( theta[subject,trait] , a[item] , b[item,category] , D )
# aqui é uma logit (D=1) ou probit (D=1.7)
 p.vec = c(p.vec,p)
   }

   #Generates a random number

   seed-runif(1,0.2,0.8)

   #Selects an answer

   temp = which ( p.vecseed )
   if( length (temp)  0 ) { resp = max (temp) } else { resp = 1 }
   resp
 }



  Esta função está inserida em um loop e deveria me retornar um número
 inteiro entre 1 e 5.
 Da primeira vez, isto acontece.
 Mas a partir da segunda, ela começa a me retornar valores como 20, 16,...

  Alguma idéia de porque isso acontece???

  Desde já, obrigado!

  Felipe Buchbinder




 ___
 R-br mailing 
 listr...@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
 mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] [Dúvida] Programação paralella e distribuida em R.

[FHRB Toledo]

Outra sugestão são as possibilidades com os pacotes snow, foreach e
multicore!

Fiz algumas brincadeiras com o multicore e para mim os resultados foram bem
satisfatórios, veja em:

http://ridiculas.wordpress.com/2011/10/23/paralelizacao-de-processos/
http://ridiculas.wordpress.com/2011/11/28/paralelizacao-de-processos-ii/

att,
FH

2012/8/31 Leonard de Assis assis.leon...@gmail.com

 Quase impossível uma solução interpretada executar com velocidade
 semelhante a uma solução compilada

 É meio lógico isso, pois a solução interpretada exige alguns passos a
 mais na máquina para ser executada

 Mas com a evolução do hardware, cada vez mais essa diferença diminui.
 Para soluções simples, a diferença é ridicula.

 []s
 Leonard de Assis
 assis dot leonard at gmail dot com

 Em 31/08/2012 00:09, Pedro Rafael escreveu:
  O performace chega a se comparar à algumas linguagens compiladas?
  [   ],
  Pedro Rafael Diniz Marinho.
 
 
  Em 30 de agosto de 2012 21:27, diogoferrari [via R-br]
  ml-node+s2285057n465604...@n4.nabble.com escreveu:
  Caro,
 
  De uma olhada no pacote RevoScaleR do Revolutions analytics.
  O processamento paralello funciona bem sim.
 
  Abs
 
 
 
  On 30/08/2012, at 20:14, Pedro Rafael [hidden email] wrote:
 
  Pessoal estou começando a estudar programação paralela e distribuida
  no R. Gostaria de saber se estas tarefas podem ser feitas com o R e se
  os resultados obtidos são interessantes. Estou lendo um pouco sobre
  OpenMP e  OpenMPI, alguem trabalha com essa API utilizando R?
 
  [   ],
  Pedro Rafael Diniz Marinho.
  ___
  R-br mailing list
  [hidden email]
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
  código mínimo reproduzível.
  ___
  R-br mailing list
  [hidden email]
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 
 
  
  If you reply to this email, your message will be added to the discussion
  below:
 
 http://r-br.2285057.n4.nabble.com/R-br-Duvida-Programacao-paralella-e-distribuida-em-R-tp4656043p4656046.html
  To unsubscribe from R-br, click here.
  NAML
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Rstudio n Ubuntu 12.04

[FHRB Toledo]

André,

O que você se refere a ... quando clico para instalar...?

Se você baixou o pacote, extensão *.deb, vá no terminal no diretório onde
está o pacote e como super usuário digite: # dpkg -i nome do pacote
extensão ponto deb!

Reporte mais detalhes do seu problema!

att,
FH

2012/8/30 Andre Oliveira andreolso...@yahoo.com.br

 Baixei o RStudio aqui em meu netbook com Ubuntu 12.04. Quando clico para
 instalar  vai abrir o gerenciador de pacotes,  mas diz que não pode abrir o
 pacote. Já baixei duas vezes e nada.  Alguma dica da possível causa?



 André Oliveira Souza

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Converter list em data.frame

[FHRB Toledo]

Se em todas as posições da sua lista tiver o mesmo tamanho (numero de
linhas você pode usar o do.call()... Veja se desse modo te adianta:

do.call(rbind, resd)

A identificação do índice fica um pouco comprometida :( mas se você souber
se nas várias posições dessa lista haja um padrão, mesmo numero de linhas
pode incluir coisas a mais do tipo, sejam 15 linhas em cada posição!

cbind(rep(1:lenght(resd), each = 15), do.call(rbind, resd)

att,
FH

2012/8/27 Thiago Veloso thi_vel...@yahoo.com.br

  Olá pessoal,

  Uma dúvida básica: possuo uma lista de objetos que contam a frequência de
 valores (count) em classes (value). Nessa lista, cada índice ([[1]], [[2]]
 etc) corresponde a um passo de tempo. Vejam abaixo uma amostra:

  head(resd)
 [[1]]! essa é a primeira unidade de tempo.
   value count
  [1,] 1   267! na classe '1' foram encontradas 267 ocorrências, e
 assim por diante
  [2,] 686
  [3,] 721
  [4,] 8 7
  [5,] 9  3683
  [6,]10   537
  [7,]11   135
  [8,]1217
  [9,]1374
 [10,]NA  2228

 [[2]]! essa é a segunda unidade de tempo.
   value count
  [1,] 1  2257! na classe '1' foram encontradas 2257 ocorrências, e
 assim por diante
  [2,] 257
  [3,] 3   289
  [4,] 4 2
  [5,] 688
  [6,] 713
  [7,] 832
  [8,] 9  1137
  [9,]10   756
 [10,]1199
 [11,]1222
 [12,]1367
 [13,]14 1
 [14,]15 7
 [15,]NA  2228

 [[3]]
   value count
  [1,] 1  2170
  [2,] 261
  [3,] 3   316
  [4,] 4 3
  [5,] 674
  [6,] 7 1
  [7,] 812
  [8,] 9   781
  [9,]10  1100
 [10,]1166
 [11,]12   160
 [12,]1346
 [13,]14 9
 [14,]1528
 [15,]NA  2228

 [[4]]
   value count
  [1,] 1  2174
  [2,] 256
  [3,] 3   300
  [4,] 4 3
  [5,] 662
  [6,] 8 4
  [7,] 9   564
  [8,]10  1317
  [9,]1175
 [10,]12   113
 [11,]1342
 [12,]1485
 [13,]1532
 [14,]NA  2228

 [[5]]
   value count
  [1,] 1  2174
  [2,] 258
  [3,] 3   301
  [4,] 4 5
  [5,] 654
  [6,] 8 4
  [7,] 9   469
  [8,]10  1393
  [9,]1183
 [10,]1276
 [11,]1341
 [12,]14   136
 [13,]1533
 [14,]NA  2228

 [[6]]
   value count
  [1,] 1  2185
  [2,] 252
  [3,] 3   300
  [4,] 4 4
  [5,] 652
  [6,] 8 3
  [7,] 9   397
  [8,]10  1457
  [9,]1183
 [10,]1274
 [11,]1338
 [12,]14   146
 [13,]1536
 [14,]NA  2228

  Como converter essa lista em um data.frame, mantendo o índice, e salvar
 em um arquivo ascii?

  Agradece pelas dicas,
  Thiago.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] OFF - Emacs

[FHRB Toledo]

Sergio,

Digite:

C-x 1 - sua janela será reduzida a uma subdivisão;
C-x 2 - sua janela subdivide em duas verticalmente; e
C-x 3 - sua janela subdivide em duas lado a lado; (é isso que você quer!)

Vale lembrar que para alternar entre buffers digite: C-x o, e caso esteja
em uma dado Buffer (o da esquerda, por exemplo) e você digitar C-x 2 sua
tela ficará subdividida em 3, na esquerda com duas uma sobre a outra e na
direita uma só apenas! Esse raciocínio se estende... faça sua subdivisão
entre janelas ao seu bel prazer, desfrute do EMACS e seja bem vindo!

att,
FH

2012/8/22 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com

 Olá

 Sou iniciante no emacs e segui o tutorial do Walmes (
 http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/), porém
 antes de modificar o meu .emacs as minhas janelas de Buffers ficavam lado a
 lado e agora fica em baixo.

 Como faço para elas ficarem lado a lado, pois prefiro dessa forma? Além
 disso, como faço para mudar a cor da seleção de um texto, pois com o fundo
 escuro a cor atual fica quase imperceptível.

 Obrigado!

 --
 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) 94429486/78101651 id: 123*20942



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] [Dúvida] Comentários multiplos.

[FHRB Toledo]

Grato as diferentes soluções que foram postadas!

2012/8/21 Daniel C Bezerra danielcbeze...@gmail.com

 Uma outra maneira de fazer isso no emacs (que vai funcionar a partir da
 versão 22.1) é usando o modo CUA. É como eu costumo fazer.

 Para ativá-lo clique em Options e depois CUA (ou use M-x cua-mode)

 Além de vc poder usar os atalhos mais comumente usados para copiar (C-c)
 colar (C-v), etc, será possível editar uma coluna inteira com apenas um
 comando.

 Basta fazer assim:
 - coloque o cursor no início do texto que vc quer comentar.
 - Ctrl-ENTER
 - Desça com o cursor até a parte final da região que vc quer comentar
 - Escreva # (sem aspas)
 - Pronto!

 Depois basta teclar Ctrl-ENTER novamente para voltar para o modo de edição
 convencional.

 Se vc quiser que o modo CUA esteja sempre disponível ao iniciar o emacs,
 após iniciá-lo em uma sessão, basta clicar em Options  Save Options
 (adicionando ao seu arquivo .emacs).

 Mais informações sobre o modo CUA estão aqui:
 http://www.emacswiki.org/CuaMode

 Abs,

 Daniel

 2012/8/20 Benilton Carvalho beniltoncarva...@gmail.com

 Eu ainda to no Emacs 23.3... pode ser que haja alguma diferenca entre
 Emacs e XEmacs...

 E e' mesmo uma boa usar os padroes... :) Ha' algum tempo uso o CTRL+D
 e agora e' mesmo forca do habito.

 b

 2012/8/20 Walmes Zeviani walmeszevi...@gmail.com:
  Infelizmente, por alguma razão ligada a peça entre cadeira e monitor,
 isso
  não funcionou pra mim. Frustrado fui a busca e encontrei
 
 
 http://malaybasu.wordpress.com/2006/10/23/comment-out-region-in-r-script-in-ess-mode-in-xemacs/
 
 http://www.quora.com/How-can-I-comment-out-an-entire-block-of-text-in-Emacs/answer/Adam-Hupp
 
  e agora tó usando o atalho nativo M-; para comentar blocos, linhas, etc.
  Obrigado pela dica.
 
  À disposição.
  Walmes.
 
 
 ==
  Walmes Marques Zeviani
  LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759
 W)
  Departamento de Estatística - Universidade Federal do Paraná
  fone: (+55) 41 3361 3573
  VoIP: (3361 3600) 1053 1173
  e-mail: wal...@ufpr.br
  twitter: @walmeszeviani
  homepage: http://www.leg.ufpr.br/~walmes
  linux user number: 531218
 
 ==
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] [Dúvida] Comentários multiplos.

[FHRB Toledo]

Também fico ansioso... !

2012/8/20 Walmes Zeviani walmeszevi...@gmail.com

 Benilton,

 Como usuário de Emacs, confesso que fiquei curioso para saber como
 adicionar o atalho. Atualmente eu .emacs se parece com isso

 http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/

 e sempre que encontro algo útil eu adiciono. Poderia compartilhar comigo
 esse dica para comentar linhas?

 À disposição.
 Walmes.

 ==
 Walmes Marques Zeviani
 LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
 Departamento de Estatística - Universidade Federal do Paraná
 fone: (+55) 41 3361 3573
 VoIP: (3361 3600) 1053 1173
 e-mail: wal...@ufpr.br

 twitter: @walmeszeviani
 homepage: http://www.leg.ufpr.br/~walmes
 linux user number: 531218
 ==

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Eliminar linhas repetidas em um data.frame

[FHRB Toledo]

Uma combinação de ifelse ou all.equal?

2012/8/10 andrebvs andre...@bol.com.br

 Desculpem, realmente, era isso mesmo!
 Obrigado Fernando e Benilton.


 --
 Em 10/08/2012 19:58, *Benilton Carvalho  beniltoncarva...@gmail.com 
 *escreveu:

 diga-se de passagem, a descricao do Fernando faz exatamente o que eu
 entendi da sua explicacao:

 df1 = data.frame(x=c(1, 2, 3, 1, 2), y=c(2, 1, 2, 2, 1))
 unique(df1)
 df1[!duplicated(df1),]


 b

 2012/8/10 Benilton Carvalho :

  Andre,
 
  se vc seguir as recomendacoes da lista e postar um exemplo
  reproduzivel junto com a descricao do problema, tenho certeza de que
  vc tera' ajuda efetiva.
 
  nao funcionou eh relativo qdo nao se tem com precisao o problema a
  ser solucionado.
 
  b
 
  2012/8/10 andrebvs :

  Obrigado Fernando, mas não funcionou. Ambas as função unique e
  df[!duplicated(df), ] é para verificação de elementos repetidos em UM
 vetor,
  no meu caso, estou querendo el iminar vetores iguais.

 
 
  
  Em 10/08/2012 18:45, Fernando Mayer  fernandoma...@gmail.com 
 escreveu:
  Ou mesmo
 
  unique(df)
 
 
  ---
  Fernando Mayer
  Universidade Federal de Santa Catarina - UFSC
  Departamento de Ecologia e Zoologia - ECZ/CCB
  URL: http://fernandomayer.github.com
  e-mail: fernandomayer [@] gmail.com
 
 
  2012/8/10 Fernando Mayer :
  Para um data.frame df
 
  df[!duplicated(df), ]
 
  []s,
 
  ---
  Fernando Mayer
  Universidade Federal de Santa Catarina - UFSC
  Departamento d e Ecologia e Zoologia - ECZ/CCB

  URL: http://fernandomayer.github.com
  e-mail: fernandomayer [@] gmail.com
 
 
  2012/8/10 andrebvs :
  Caros colegas, como eliminar linhas iguais em um data.frame. No
 exemplo
  abaixo, tenho que a linha 3 é igual a linha 10,14 e 85. Gostaria
  que esse vetor, assim como outros que se repitam, aparecessem apenas
 uma
  vez.
 
  1 2 3 4 5 6 8 9 12 15
  2 2 3 4 5 6 8 9 12 13
  3 2 3 4 5 6 8 9 12 13
  4 2 3 4 5 6 8 9 15 16
  5 2 3 4 5 6 8 9 13 15
  6 2 3 4 5 6 9 12 15 16
  7 2 3 4 5 8 9 12 13 15
  8 2 3 4 5 6 8 9 12 15
  9 2 3 4 5 6 8 9 11 12
  10 2 3 4 5 6 8 9 12 13
  11 2 3 4 5 6 8 9 11 15
  12 2 3 4 5 6 8 9 11 15
  13 2 3 4 5 6 8 9 11 15
  14 2 3 4 5 6 8 9 12 13
  15 2 3 4 5 6 8 9 12 15
  .
  .
  .
  84 1 2 4 5 9 10 11 12 15
  ; 85 2 3 4 5 6 8 9 12 13
 
  desde já agradeço!
 
  Att.
  André
 
 
 
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/list 
  info/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
  código
  mínimo reproduzível.
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 
 
  ___
  R-br mailing list
  R-br@listas.c3sl.ufpr.br
  http 
  s://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

  Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código
  mínimo reproduzível.
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dúvidas sobre uma função

[FHRB Toledo]

Olha só... quem é vivo sempre aparece!

Matheus, seu CMR não é muito claro, mas vamos ver se entendi o problema e
essa solução resolve!

# sua matriz...
matriz - matrix(c(1, 2, 4, 4,
   3, 6, 3, 1,
   7, 4, 8, 9,
   1, 2, 2, 9), 4, 4, byrow = TRUE)

# outra matriz um pouco maior!
outra - matrix(sample(1:10, 1000, replace = TRUE), 250, 4)

# a funcao...
funcao - function(x)ifelse(sum(match(x, c(1, 2, 3)), na.rm = TRUE) != 0,
'tem', 'nao')

# aplicando-a a sua e a outra matriz
apply(matriz, 1, funcao)
apply(outra, 1, funcao)

A lógica por tras da função é: faz-se um match nos valore da linha da
matriz com o seus valores alvo (1, 2 e 3), retornando a posição desses
valores na linha ou NA! depois soma-se esses valores, excluindo os NA, se a
soma for diferente de zero pelo menos um dos numeros (1, 2 ou 3) têm na
referida linha!

Resta só agora esperarmos por uma solução muito mais elegante dos colegas
da lista.

abraço,
FH

2012/8/3 Matheus Mendes mendes...@hotmail.com

 Bom dia...

 Estou com um trabalho a ser feito e gostaria de solicitar ajuda.

 Tenho uma planilha com a seguinte estrutura:

 ** **

 Primeiro Segundo Terceiro Quarto

 12  44   

 36  31

 74   8   9

 12   2   9

 ** **

 Preciso identificar quais as linhas que recebem os números 1, 2, 3 e 4,
 não importando a ordem, ou seja, se é primeiro, segundo, ...

 Estou com dúvidas para fazer uma função para este fim...

 Grato

   

 Matheus Henrique Silveira Mendes

 Engenheiro Agrônomo - UFLA

 Mestre em Genética e Melhoramento de Plantas – UFLA

 Doutorando em Genética e Melhoramento de Plantas – UFLA

 (35)91243078 (TIM)

 (35)98871518 (VIVO)

 ** **

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Encontrar caminho de um arquivo no computador

[FHRB Toledo]

A propósito, essas linhas do Benilton é o que se pode chamar de *avec
elgance*!

2012/7/27 Luis Iván Ortiz Valencia liov2...@gmail.com

 Muito obrigado pelas valiosas respostas.

 Em 27 de julho de 2012 13:39, Listeiro 037 
 listeiro_...@yahoo.com.brescreveu:

 Olá.

 O velho DOS, se não me engano, tinha um FIND, mas no momento não há
 como testá-lo aqui.

 Mas há esta página, que talvez possa ser útil:

 http://www.computerhope.com/findhlp.htm

 Quando eu procurava prá saber de um arquivo na hierarquia, eu usava
 DIR /S/A, que era prá exibir subdiretórios abaixo e arquivos de
 qualquer atributo, com o devido nome do arquivo colocado à frente.

 No seu caso, esse FIND pode ser melhor.


 Em Fri, 27 Jul 2012 13:22:06 -0300 (BRT)
 Paulo Justiniano paulo...@leg.ufpr.br escreveu:

  é comando do LINUX
 
 
  On Fri, 27 Jul 2012, Luis Iván Ortiz Valencia wrote:
 
   Prof Paulo
   Não consegui achar a função locate, eu uso OS windows.
  
   abraços
  
   Ivan
  
   Em 27 de julho de 2012 08:47, Paulo Justiniano
   paulo...@leg.ufpr.br escreveu: Complementanto...
  
 $ locate indicadores.pdf
  
 tb pode ser usado
 locate é mais rápido pois consulta uma tabela ao inveś de
   varrer o conteúdo no momento. Entretanto a tabela precisa estar
   atualizada o que pode/é ser feito no sistema regularmente com
   updatedb (root)
  
 Ponrtanto se é um arquivo mais recente find é melhor,
 mas para vários arquivos ou mais antigos vale a pena
   consuiderar locate, eventualmente precedido de updatedb (se de fato
   necessário)
  
  
  
  
 On Thu, 26 Jul 2012, Benilton Carvalho wrote:
  
  
   Ivan,
  
   Não estou bem certo do que vc precisa. Vc quer algo q
   funcione como o 'find . -name indicadores.pdf'?
  
   Se eu tivesse muito q fazer isso, usaria
  
   Fns=list.files('c:/', recursive=TRUE, full=TRUE)
   i=grep(^indicadores.pdf$, basename(Fns))
   path=ifelse(length(i) 0, dirname(Fns[i]), NA)
  
   b
  
   Sent from a mobile device
  
   On Jul 26, 2012 3:16 PM, Luis Iván Ortiz Valencia
   liov2...@gmail.com wrote: Prezados
   Gostaria de saber se é possível no R encontrar o
   caminho de um arquivo específico no computador. Por exemplo, existe
   um arquivo indicadores.pdf
   na pasta C:\ivan\R, existe uma função no R que responda
   com o caminho C:\ivan\R dando como entrada indicadores.pdf?
  
   atte
  
   IVAN
  
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.





 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Usando o R no Ubuntu

[FHRB Toledo]

Uma boa sugestão é o Emacs+ESS, que já foi dita...

Ele entre outras, permite muitas customizações ao seu bel prazer, veja
algumas possibilidades:
http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/

2012/7/24 Leonard de Assis assis.leon...@gmail.com

  to gostando (acho q acostumei com ele) do rstudio

 []s
 Leonard de Assis
 assis dot leonard at gmail dot com

 Em 24/07/2012 08:11, Fabio Mathias Corrêa escreveu:

  Olá Vinícius,

  No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS,
 Gedit+R-Plugin, RStudio e outros.

  Eu particularmente uso o Gvim+Vim-r-plugin.

  Dê uma pesquisada e veja qual lhe agrada mais!

  Valeu!


 Fábio Mathias Corrêa

 Universidade Estadual de Santa Cruz
 Departamento de Ciências Exatas e da Terra - DCET


 Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna
 CEP 45662-900. Ilhéus-Bahia


  Tel.: 73-3680-5076
   --
 *De:* Vinícius Lionel Mateus vinynegre...@gmail.comvinynegre...@gmail.com
 *Para:* R-br@listas.c3sl.ufpr.br
 *Enviadas:* Terça-feira, 24 de Julho de 2012 0:13
 *Assunto:* [R-br] Usando o R no Ubuntu

 Prezados,

 Alguém poderia me indicar algum material para ajudar no uso do R na
 plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar
 usar o GUI também, assim como no Windows.

 --
 --
  Atenciosamente,

 VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665)
 Bacharel em Química - Doutorando em Química Analítica
 Laboratório de Química Atmosférica - Departamento de Química
 Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio)
 Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.:
 22451-900
 Telefone: (+55) (21) 3527-1327
   (+55) (21) 9358-8051
 www.puc-rio.br

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing 
 listr...@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
 mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Função implícita

[FHRB Toledo]

Lagrange?

2012/7/1 Eder David Borges da Silva e...@leg.ufpr.br:
 Davi,
 Sem um CMR fica difícil.Mas segue algumas funções que podem de dar uma
 ajuda:
 ?optim
 ?rootSolve
 Essa são genéricas  e numéricas , a sua função não tem solução analítica?
 Att



 Em 1 de julho de 2012 00:43, Davi davi.butt...@gmail.com escreveu:

 Bom dia pessoal,



 O problema que estou é o seguinte: preciso encontrar uma (ou mais, não
 sei) soluções para uma equação do tipo



 f(x) - c =0



 onde f(x) é uma função implícita e c uma constante qualquer. Tem algum
 método?



 Muito obrigado,



 Davi









 ___

 Davi Butturi-Gomes

 Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP -
 Piracicaba

 Mestre em Biometria - UNESP - Botucatu

 Ecólogo - UNESP - Rio Claro

 MSN: davi_bgo...@hotmail.com

 Skype: davibg87

 (19)9809.5753




 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
 mínimo reproduzível.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.