- R-br - DINF UFPR

Contrastes e função
by Leandro Marino 18 Ago '11

18 Ago '11

#Caros, existe alguma forma de melhorar este exemplo: vetor <- rep('0',108) z=1 for(i in 0:5){ for(j in 0:2){ for(k in 1:2){ for (l in 0:2){ vetor[z] <- paste(i,j,k,l,sep='') z = z + 1 } } } } vetor #Aproveitando, se quero usar um tapply extraindo média, #desvio-padrão e uma função personalizada, alguém saberia com fazer as #três funções de uma vez só? #Gostaria de executar o comando abaixo de uma única vez em um tapply. #Alguma sugestão?! set.seed(1) x <- rnorm(1000,500,100) y <- trunc(runif(1000,1,10)) fper <- function(x){mean(x-500)/100} dataframe <- t(rbind(tapply(x,y,mean),tapply(x,y,sd),tapply(x,y,fper))) colnames(dataframe) <- c('media','desvp','fper') dataframe Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 R (21) 2103-9600 R.:236 ( leandro(a)cesgranrio.org.br  Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

3 3

Material Novo no CRAN
by Marcos Silva 17 Ago '11

17 Ago '11

Só para avisar que existe um novo material sobre o R em portugues no CRAN. http://cran.r-project.org/doc/contrib/Provete-Estatistica_aplicada.pdf Valeu. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

2 1

Malha sobre poligono
by Eder David Borges da Silva 17 Ago '11

17 Ago '11

#Caros colegas, #minha duvida é uma implementação para gerar uma malha da seguinte forma: # 1) Imaginem que tenho um poligono como esse? x <- c(-0.512320711, -0.503592775, -0.669423554, -0.841072957, -0.459953097, -0.186477777, -0.221389520, -0.000281815, 0.319742495, 0.569943319, 0.503029145, 0.209188642, 0.022992680, -0.169021906, -0.241754703, -0.305759565, -0.352308556, -0.390129611, -0.422132042, -0.422132042,-0.512320711) y <- c( 0.73217213, 0.38635022, 0.08640264, -0.43585901, -0.80638249, -0.44997419, -0.03357638, 0.12169059,-0.13238264, 0.07228747, 0.48162768, 0.81333441 , 1.10622439, 1.13092596, 1.14504114, 1.13445475, 1.12739716, 1.12033957, 1.06035006, 1.00388934, 0.73217213) plot(x,y,type='l') # 2) dentro deste poligono quero fazer quadrados irregulares como este px1 <- c(-0.6131326, -0.4633484, -0.3762646, -0.5034070) py1 <- c(-0.6037305,-0.7381249, -0.5774359, -0.3992173) polygon(px1,py1) # 3) quero fazer isso com o mouser clicando 4 vezes e guardando as coordenadas e cada coordenada ganha um ID, ficando assim: data.frame(id=1:length(px1),px1,py1) # 4) este quatros pontos forman um elemento (quadrado irregular), onde o elemento guarda o ID os pontos clicados. elemento <- c(1,1,2,3,4)#Idelemento - 4 ids dos pontos # 5) Esse procedimento é repetido n vezes cobrindo toda a região com poligonos irregulares, como resultado final: # um data frame com id e x, y do ponto # um data frame com id do elemento # restrições para o segundo elemento eu devo clicar por exemplo no ponto 2 e 3 e fazer mais o 5 e 6 para fechar ele, desta forma, # os elementos tem de ficar plotados para aproveitar os pontos já feitos. #OBS: para os deterministicos de plantão isso é uma malha para modelos numericos de mecanica de fluidos. #OBS1: Sugestoes de softwares livres que faz isso são bem vindo. # Diante do problema alguem tem sugestoes de funçoes que posso utilizar, na verdade é um mix de locator + identifique, pois não # consigui fechar todo o raciocinio de como fazer isso. # Obrigado pessoal #Att

2 2

Modelo multionomial
by Rodrigo Pereira 17 Ago '11

17 Ago '11

Prezados colegas, bom dia. Estou tentando ajustar um modelo logístico multinomial no R, mas creio que não estou tendo sucesso. A variável resposta tem três níveis (0, 1 e 2) e com 3 variáveis resposta, sendo duas na forma dicotômica e uma contínua. Tentei ajustar um modelo utilizando a função 'glm', mas para a opção de família não há explícita a multinomial. Devo utilizar a binomial e automaticamente a função entenderá que deverá ocorrer uma modificação? Há algum pacote que forneça uma análise completa para uma regressão multinomial? Além disso, preciso de algumas outras estatísticas, tais como odds ratio ou risco relativo (se possível numa comparação duas a duas), intervalos de confiança para os parâmetros, avaliação de efeitos, etc. Desde já agradeço a atenção dispensada. Atenciosamente, Rodrigo Pereira

6 9

Ajuda para simulações
by sznelwar 17 Ago '11

17 Ago '11

4 5

[OFF-Topic] Mendelev
by FHRB Toledo 17 Ago '11

17 Ago '11

Pessoal da R-Br, olá Não sei se todos estão altamente "antenados", mas parece que estamos vivendo uma revolução... Surgem grandes oportundiades p/ o pessoal envolvido em pesquisa ao usar o Mendelev! Pois bem, estou cursando uma disciplina e o professor faz parte de algum "circulo maior" da organização do mendelev! Para incentivar o seu uso ele agendou um webnário sobre o ambiente... Seguem as informações: Atendendo ao foi solicitado, o pessoal do Mendeley vai apresentar um seminário em português especial para os alunos da disciplina e estou tomando a liberdade de extender p/ turma da lista! Data e horário: Terça-feira, dia 16 de Agosto, 11:00 da manhã (hora de Brasília). Local: internet (é um "webnário"). Vocês devem se cadastrar antes no link https://mendeley.webex.com/mendeley/onstage/g.php?t=a&d=663141002 e na hora do seminário basta entrar na mesma página. **IMPORTANTE:** uma vez que isso será feito porque nós pedimos, será bem chato se poucos participarem, então por favor não deixem de se inscrever e entrar na página no horário correto. Estou esperando pelo menos uns 30 participantes. Abraços, FH

9 14

Ajuda com modelo não linear
by Eduardo Lucas Terra Peixoto 17 Ago '11

17 Ago '11

Boa noite a todosEstou fazendo doutorado na UEL no curso de Ciência animal. Ainda sou novato no uso do R e também não tenho base estatística aprofundada em modelos não lineares.E preciso de um help, estou adaptando uma fórmula de modelo não linear do proc do SAS para o R, mas li o help, busquei dicas no site http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase30.html#x32-20100030 e fui seguindo o modelo explicado por esse autor, mas tem algumas coisas que não entendi e não estou conseguindo rodar o modelo.Primeiramente explicarei resumidamente o que se trata o modelo. Esse modelo é para poder determinar cinco outras variáveis, através da produção de gás pela fermentação de bacteriais do líquido ruminal in vitro, simulando assim o rumen do animal. As variáveis que tenho em mãos são o volume de gás produzido e o tempo da leitura, a partir deles preciso determinar:VcnF = volume de gás produzido pelo carboidratos não fibrosos (mL de gás)kdcnf = taxa de degradação dos carboidratos não fibrosos (%/hora)VcF = volume de gás produzido pelo carboidratos fibrosos (mL de gás)kdcf = taxa de degradação dos carboidratos fibrosos (%/hora)L= tempo de latência (h), este é o tempo que as bactérias gastam para começar a degradar o substrato.Mas ao rodar o modelo abaixo aparece o seguinte erro: fit30=nls(Y~ (VcnF/(1+exp(2-4*kdcnf*I(L-Tempo)))+ VcF/I(1+exp(2-4*kdcf*(L-Tempo)))), + data=cra30, + start=list(VcnF = 10,kdcnf = 0.01,L =0.1,VcF = 50,kdcf =0.01, Tempo=1), + alpha= (0.05))Error in qr.qty(QR, resid) : NA/NaN/Inf in foreign function call (arg 5)In addition: Warning message:In Ops.factor(lhs, rhs) : - not meaningful for factors Aproveitando o ensejo gostaria de perguntar como faço para limitar os parâmetros pois no proc do SAS está assim:VcnF = 10 to 200 by 50 kdcnf = 0.01 to 0.8 by 0.05 L = 0.1 to 5 by 0.2 VcF = 50 to 150 by 30 kdcf = 0.01 to 0.8 by 0.05; bounds L >=0; Caso alguém tenha materiais relativos a modelos não linear para que possa aprofundar melhor, também contribuirá bastante. Agradeço a todos, e espero que alguém possa me ajudar! Eduardo Lucas Terra PeixotoDiscente de doutorado - Ciência AnimalUniversidade Estadual de Londrina Cel. (43) 9653-5574

2 1

Fwd: [Computing] package parallel in R-2.14.0
by Leonard Assis 16 Ago '11

16 Ago '11

Boa notícia lmassis <at> yahoo <dot> com <dot> br assis.leonard <at> gmail <dot> com ---------- Forwarded message ---------- From: Jan de Leeuw <deleeuw(a)stat.ucla.edu> Date: Tue, Aug 16, 2011 at 1:40 PM Subject: [Computing] package parallel in R-2.14.0 To: statcompute <computing(a)stat.ucla.edu> parallel is a recommended package in R-2.14.0 -- I attach the vignette. === Jan de Leeuw; Distinguished Professor and Chair, UCLA Department of Statistics; Editor: Journal of Multivariate Analysis, Journal of Statistical Software; US mail: 8125 Math Sciences Bldg, Box 951554, Los Angeles, CA 90095-1554 phone (310)-825-9550; fax (310)-206-5658; email: deleeuw(a)stat.ucla.edu .mac: jdeleeuw ++++++ aim: deleeuwjan ++++++ skype: j_deleeuw homepages: http://gifi.stat.ucla.edu ++++++ http://www.cuddyvalley.org ------------------------------------------------------------------------------------------------- No matter where you go, there you are. --- Buckaroo Banzai http://gifi.stat.ucla.edu/sounds/nomatter.au ------------------------------------------------------------------------------------------------- _______________________________________________ Computing mailing list Computing(a)lists.stat.ucla.edu http://lists.stat.ucla.edu/mailman/listinfo/computing

1 0

dúvida em glm
by Vinicius Brito Rocha 16 Ago '11

16 Ago '11

Senhores, ajustei um glm com família binommial e link logit, onde minha variável resposta é a probabilidade do evento (chamadas telefônicas repetidas em 24hs / total de chamadas telefônicas), obtidos a aprtir de um aggregate em um banco de dados com informação de telefones. A minha varíavel independente é um fator (UF), Não estou conseguindo interpretar os coeficientes do meu modelo. tenho 27 coef, devido ao meu fator ser UF +1 intercepto. Alguém pode me dar uma força? Abs. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 6

Re: [R-br] Análise Envoltória de Dados (DEA)
by Roney Fraga Souza 16 Ago '11

16 Ago '11

Daniel, bom dia! Ainda não trabalhei com DEA SBM sem orientação, contudo estou aqui pensando se vc utilizar no comando dea() as opções ORIENTATION='in-out' com DIRECT diferente de NULL será que chega perto do seu objetivo? Outra dúvida, em um tópico anterior também sobre DEA, no link, http://r-br.2285057.n4.nabble.com/R-br-Modelo-linear-generalizado-misto-td3… vc postou que, Pacote Benchmarking: Vantagem - Te retorna os pesos(contribuição) de cada DMU onde a soma é igual a 1. Desvantagem - Não te fornece os slacks (quanto a DMU deve aumentar ou diminuir de produto ou insumo para se tornar eficiente) com relação a não retornar os slacks, quando utilizo em ccr.in <- dea(insumos, produtos, RTS='crs', ORIENTATION='in', SLACK=TRUE) o comando summary(ccr.in) já aparece se existe folga em x ou em y. E quando chamo names(ccr.in) cbind(ccr.in$sx,ccr.in$sy,ccr.in$slack) tenho as folgas de x e y. Minha pergunta é, minha análise está equivocada, ou realmente é possível obter as folgas pelo pacote Benchmarking? Atenciosamente Roney

1 0