Agosto 2011 - R-br

gerar dados
by Ricardo Pereira 04 Out '11

04 Out '11

Outra duvida. Estou tentando simular dados, tal como um exemplo em outro programa... Exemplo abaixo: proc iml; do i=1 to 1000; seed=456; var=150; w = RANNOR(seed)*SQRT(var); end; *tentei:* set.seed(456) var <- 150 for (i in 1:1000) { w [i] <- rnorm * sqrt(varg) } *ERRO:* Erro em W[j] <- 1 : objeto 'W' não encontrado

10 15

Erro em função
by Carlos Mendonça 23 Set '11

23 Set '11

Caros, estou utilizando a função abaixo que calcula os quartis ponderados que funcionou muito bem até agora. tmp = subset(arq1, select=c(codigo, qtd, val)) equantileByCounts <- function(x, counts, qs){ tot <- sum(counts) i <- order(x) x <- x[i] counts <- counts[i] partial <- cumsum(counts) qsObs <- qs*tot start <- floor(qsObs) end <- ceiling(qsObs) i <- length(qs) aux1 <- sapply(1:i, function(.x, start, end, x){ coord <- c(start[.x], end[.x]) bins <- c(1, partial) grps <- cut(coord, bins, labels=FALSE, include.lowest=T) mean(x[grps]) }, start=start, end=end, x=x) names(aux1) <- names(qs) return(aux1) } aux1 <- do.call(rbind, lapply(split(tmp, tmp$codigo), function(y) equantileByCounts(y$qtd, y$val, c(.25, .75)))) O meu processamento é por trimestres e em apenas um deles deu a seguinte mensagem de erro: "Erro em cut.default(coord, bins, labels = FALSE, include.lowest = T) : 'breaks' não são únicos" Alguém saberia me dizer o que posso fazer? Um abraço, Carlos Mendonça.

2 14

Re: [R-br] [Dúvida] Invocar funções C++ no R
by Pedro Rafael 09 Set '11

09 Set '11

Quando eu falo conectado quero dizer que ele não está sendo reconhecido no DOS. Tava lendo um tutorial que ele manda testar se as variáveis ambientes adicionadas foram colocadas de forma correta. Para ver se o gcc está sendo reconhecido pelo DOS o tutorial manda digitar gcc --help. Para ver se o R está podendo ser chamado pelo DOS sem referenciar o caminho visto que o patch foi add corretamente ele manda fazer R.exe. Já para ver se a variável ambiente do Rtools está adicionada de forma correta ele manda digitar o comando grep. Contudo consegui corrigir o problema a poucos minutos atrás. Muito obrigado Benilton pelas dicas. Qualquer dúvida a mais que tiver eu irei aperriar (risos). Em 26 de agosto de 2011 08:55, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Quando eu falo conectado quero dizer que ele não está sendo reconhecido no > DOS. Tava lendo um tutorial q > > Em 26 de agosto de 2011 08:41, beniltoncarvalho [via R-br] < > ml-node+3770614-1267126705-223914(a)n4.nabble.com> escreveu: > > O que exatamente vc quer dizer com 'RTools nao esta conectado'? RTools >> e' apenas um conjunto de ferramentas (MinGW, Perl, etc) e essas voce >> acessa pelo terminal do DOS. Note que, com o RTools, vc nao precisa do >> gcc. >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3770614&i=0> >> 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-Invocar-funcoes-C-no-R-tp3768… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho. > > > -- Saudações, Pedro Rafael Diniz Marinho.

6 23

R: Script x Programa
by Lucas Ferreira 01 Set '11

01 Set '11

Desculpem a ingenuidade...Gostaria de saber se um Script feito para o R pode ser considerado um programa.Por exemplo, vejo algumas pessoas com o currículo possuindo : Programa X sem registro.Um dos meus scripts é usado por pessoas em universidades diferentes. sera que eu poderia declara-lo no no lattes como um "programa sem registro " ?existe algum impedimento legal em fazer isso? Não tenho nem ideia. obrigado e desculpem ------------------------------------------------------------------ -Lucas Ferreira da Silva Mestrando em Biotecnologia - UFESLaboratório de Biotecnologia Aplicada ao AgronegócioNúcleo de Biotecnologia/CCS/UFESAv. Marechal Campos, 146829040-090, Vitória, ES, BrasilTel/Fax. 55.27. 3335-7348 ----------------------------------------------------------------- ------------------------------------------------------------------

7 11

Ajuda com Função poder
by Victor Eduardo 31 Ago '11

31 Ago '11

Gente, alguém poderia me dar uma luz de como construir a função poder para um teste de hipóteses desse jeito: lançando uma moeda 10 vezes e quero testar isso: H0: p=0.5 H1: p diferente de 0.5 A questão fornece a regiao critica e pede para construir a função poder. Resolvendo tudo, cheguei : P(0.6 - p)/sqrt(p*(1-p)/10) - P (0.4 - p)/sqrt(p*(1-p)/10) Como poderia começar a construir a função poder nesse caso? Abraços!

1 0

equivalente ao "select case" do VB em R
by Antonio Silva 31 Ago '11

31 Ago '11

Olá, Estou passando para R uma rotina em visual basic. No R, qual o comando mais similar ao "select case" do VB? Por exemplo: select case int(AG) case 0, 5, 8, 10 VAR = 37 case 1, 4 VAR = 70 case 3, 6, 9 VAR = 90 case 2, 7 VAR = 112 end select Uma saída seria utilizar "if else" de forma encadeada, mas isto não seria muito elegante. Alguém saberia informar outra opção? Obrigado! Antonio Olinto

3 3

MS-VAR
by julio cesar araujo 31 Ago '11

31 Ago '11

Boa noite Alguem saberia me dizer se o R tem alguma rotina para estimar MS-VAR (Markov Switching-Autoregressions Vector)? Atc Julio

2 2

Dúvida em função de maximização (maxLik)
by Orleans Martins 31 Ago '11

31 Ago '11

Prezados, estou tentando maximizar uma função de verossimilhança para estimar 5 parâmetros. Contudo, meus resultados estão dando NULL. Sou iniciante na manipulação do R, por isso posso ter cometido erros primários. Favor alguém me ajude. Meus dados são: B <-c(1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,1,0,1,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,1,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0) S <- c(0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,1,0,1,1,1,0,1,1,1,0,1,1,1,1,1,1,1,0,1,0,0,0,0,0,1,1,1,1,1,1,1,1,1) A distribuição é de Poisson: B <- rpois(50, 1) S <- rpois(50, 1) A equação que utilizo é: loglik <- function(param) { alpha <- param[1] mu <- param[2] Eb <- param[3] Es <- param[4] delta <- param[5] ll <- (1-alpha)*exp(-Eb)*((Eb^B)/factorial(B))*exp(-Es)*((Es^S)/factorial(S)) + alpha*delta*exp(-Eb)*((Eb^B)/factorial(B))*exp(-(mu+Es))*((mu+Es)^S/factorial(S)) + alpha*(1-delta)*exp(-(mu+Eb))*((mu+Eb)^B/factorial(B))*exp(-Es)*((Es^S)/factorial(S)) } res <- maxLik(loglik, grad=NULL, hess=NULL, start=c(0,1), constraints=NULL) print(res) coef(res) Acho que estou errando nos argumentos dessa equação. Alguém tem alguma sugestão?

2 2

Correlação Canônica!!
by Thiago De paula protásio 31 Ago '11

31 Ago '11

Pessoal, novamente recorro a valiosa ajuda de vcs... Estou querendo utilizar as variáveis canônicas para um agrupamento preliminar de clones de Eucalyptus..... Dessa forma, estou utilizando o pacote CCA e função cc(). Eu queria plotar os escores do primeiro par de variáveis canônicas e analisar..... Contudo, a função plt.indiv() desse mesmo pacote me retorna um gráfico que não estou compreendendo. Me parece que não está sendo plotado os escores.... Talvez, o R esteja considerando as variáveis canônicas padronizadas, mas não tenho certeza.Olhei no help da função, mas as informações são escassas. Gostaria de perguntar se alguém já trabalhou com essa função. Muito obrigado...Todas as contribuições são muito bem vindas.... Olhei no help, mas não tem muitas informações.... Gostaria de valiosa ajuda do senhor para resolver esse problema. Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

4 3

Erro relacionado a algum loop?
by . . 31 Ago '11

31 Ago '11

Olá pessoal, estou recebendo o seguinte erro, Error: evaluation nested too deeply: infinite recursion / options(expressions=)? mas nem sequer há um loop explícito no código. O que pode ser? Obrigado. func <- Vectorize(function(x, a, sad, trunc=0, ...) { result <- function(x) { f1 <- function(n) { dcom <- paste("d", deparse(substitute(sad)), sep="") dots <- c(as.name("n"), list(...)) f <- do.call(dcom, dots) g <- dpois(x,a*n) f * g } integrate(f1, 0, Inf)$value } return(result(x) / (1 - result(trunc))) }, "x") func(1, 0.1, exp)

1 0