Fevereiro 2013 - R-br

Comparação entre duas populações com distribuição de Poisson
by Alexandre Santos 12 Fev '13

12 Fev '13

Bom dia pessoal, Tenho duas populações com distribuição de Poisson e gostaria de saber qual a abordagem mais correta para comprar essa duas populações, se empregando análise de variância sobre um GLM ajustado ou teste de Chi quadrado sendo que os dois testes fornecem o mesmo resultado, sendo o CRM: # Dados y<- c(mapply(rpois, lambda=c(5,5), MoreArgs=list(n=15)))##Criação da variável resposta trat <- gl(2,15) ##Criação da variável resposta #Usando GLM m.1 <- glm(y~trat,family="poisson") anova.m.1 <- anova(m.1, test="Chi") anova.m.1 # #Usando teste Chi^2 chisq.test(y,correct=TRUE) # Obrigado, Alexandre

2 2

BLUP
by Ze Henrique 10 Fev '13

10 Fev '13

Bom dia a todos. Gostaria de saber se alguem conhece como trabalhar com BLUP no R, pacote ou algo que faça este tipo de estatística. Desde já agradeço a ajuda. Obrigado. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2

Método de Newton Rhapson...
by andrebvs＠bol.com.br 09 Fev '13

09 Fev '13

2 9

Join my network on LinkedIn
by Juliana Passos via LinkedIn 09 Fev '13

09 Fev '13

LinkedIn ------------ Juliana Passos requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Juliana Passos http://www.linkedin.com/e/-33g73l-hcytrjre-4i/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Juliana Passos http://www.linkedin.com/e/-33g73l-hcytrjre-4i/rso/230788258/wFJc/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hcytrjre-4i/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Programando com C e R
by ctrucios 08 Fev '13

08 Fev '13

Boa tarde pessoal, Estou tentando utilizar C e R juntos, fiz 3 funções em C e chamo elas do R. Mas estive testando e testando e deu este erro. Na verdade as vezes da este erro e as vezes trava de vez. Não sei se é pela forma em que estou programando em C (na verdade estou bem enferrujado e estou utilizando ponteiros) deixo aqui o código quem puder de uma olhadinha. [image: Imagem inline 2] Na verdade é um mesmo código que faço três vezes (assim que descobri o erro). O problema aparece quando rodo mais de dois vezes (por isso que aqui copiei o codigo 3 vezes) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) as funções em C seguem no seguinte link. https://www.dropbox.com/sh/5f8bpx6pt2cgam6/fcuvUtFili/Codigos *Carlos Trucíos Maza ** : ctrucios(a)gmail.com

3 3

ESTIMATIVA DA CONCENTRAÇÃO LETAL MÉDIA (CL 50)
by Tiago Souza Marçal 07 Fev '13

07 Fev '13

Ola pessoal alguém conhece alguma rotina para a estimamação da CL 50 no R? Att. Tiago.

3 3

como contar elementos em um data.frame
by Aline Lipsky 07 Fev '13

07 Fev '13

> Boa tarde, > > Eu tenho um data.frame como o abaixo: > > Data Embarcação Municipio > Jan A RG > Jan B PEL > Jan J PEL > Dez J RG > Mar A RG > Set F PEL > Set G RG > Dez G RG > > > Preciso contar as embarcações por mês e Município e obter uma tabela assim: > > Mês RG PEL > Jan 5 16 > Fev 26 31 > Mar 7 2 > > Grata pela atenção, > > Aline

2 4

Duvida sobre como alterar o contraste no options de "contr.poly" para "contr.sum" e entender o que acontece.
by Augusto Ribas 07 Fev '13

07 Fev '13

Ola caros colegas da lista, eu estava tentando estabelecer a matriz de contrastes para usar na função lm, mas não estou conseguindo entender como as coisas funcionam no R. Eu li um tutorial antigo, esse aqui (http://wwwuser.gwdg.de/~cscherb1/content/Statistics%20Course%20files/Workin…) Eu não entendo se as coisas mudaram desde que esse tutorial foi escrito, mas com matriz de soma como no exemplo acima tudo é simples, mas eu não consigo mais replicar o exemplo no R. Mais precisamente o que eu estava tentando era o seguinte: #Dados de exemplo: variavel<-factor(rep(letters[1:3],each=10)) resposta<-rnorm(30,rep(c(1,1,2),each=10)) #Ai eu ajusto um modelo linear #modelo usando opções default modelo<-lm(resposta~variavel) summary(modelo) #Eu imagino que a matriz que deveria fazer os contrastes como acima seria essa no caso de usar soma como no texto cbind(c(-1,1,0),c(-1,0,1)) Mas o default usa outro esquema de matriz de contrastes, logo so colocar isso nos contrastes não funciona #O default não usa esse sistema certo? options("contrasts") #Eu imaginei que essa opção que fosse a que eu precisaria mudar options("contrasts")$contrasts[2] #Mas tentando mudar de forma simples assim eu obtenho um erro #mas não é uma mudança simples assim, quando eu tento isso: options("contrasts")$contrasts[2]<-"contr.sum" #eu tenho essa mensagem de Erro Erro em options("contrasts")$contrasts[2] <- "contr.sum" : destino da atribuição expande para um objeto fora da linguagem Mas como eu faço ao invéz de escrever isso, contrasts(variavel) #Escrever isso, ou algo assim cbind(c(-1,1,0),c(-1,0,1)) #E obter o mesmo resultado que esta no summary(modelo) la em cima, se não me perdi em algum momento mais para cima. #Porque mesmo a opção contr.sum parece que ele trabalha de outra forma do que estou imaginando. contr.sum(3,contrasts=T,sparse=F) Bem é isso, é apenas uma duvida sobre como o R funciona que não estou conseguir entender, se alguém puder dar uma luz eu serei grato. Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 5

Contar sequência de números...
by andrebvs＠bol.com.br 07 Fev '13

07 Fev '13

3 2

percentage de error da predicao de GLM
by Humberto Hazin 07 Fev '13

07 Fev '13

Bom dia Pessoal estou tentando calcular o percentual de predição de erro via boot, porem o resultado final não me parece o que quero (percentual de predição de erro ) alguém poderia me ajudar? train.pos.setbsh data set com ~27000 linhas (https://www.dropbox.com/s/9flxc9tpttq1mru/train.pos.setbsh.csv) test.setbsh data set com ~1433 linhas (https://www.dropbox.com/s/61fbts7bbizb7cr/test.setbsh.csv) bt <- function(formula, data, indices) { d <- data[indices,] fit <- glm(formula,family=gaussian,data=d) pred<-predict(fit,test.setbsh,type="response") error<- (test.setbsh$lgcpue-pred)/test.setbsh$lgcpue return(error) } # bootstrapping with 300 replications results <- boot(data=train.pos.setbsh, statistic=bt, R=300, formula=lgcpue~Y+S+Q+A) boxplot(as.vector(results$t)) abline(0,0) desde já agradeço Humberto Enviado do Email do Windows

1 0