- R-br - DINF UFPR

Modelo de Ising em duas dimensões para imagens
by Mauro Sznelwar 03 Set '11

03 Set '11

Existe em algum lugar, ou alguém tem implementado este no modelo no R para imagens?

3 2

Correlação Canônica! Cargas Padronizadas!
by Thiago De paula protásio 02 Set '11

02 Set '11

Pessoal; Gostaria da valiosa ajuda de vcs... Alguém conhece algum comando no R para que eu possa obter as cargas canônicas e as cargas canônicas cruzadas padronizadas??? Muito obrigado!!!! 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

1 0

sqlSave - pacote RODBC
by Vinicius Brito Rocha 02 Set '11

02 Set '11

Senhores, estou tentando, sem sucesso, dar um apend em uma tabela criada em um banco Teradata, usei o própiro sqlSave para criar a tabela, porém quando uso o sqlSave para inputar a tabela com parâmeto apend=T e inputar os arquivos tenho a seguinte mensagem de erro. sqlwrite returned [RODBC] Failed exec in Update 23000 -2802 [NCR][ODBC Teradata Driver][Teradata Database] Duplicate row error in SB_MERCADO.ech_amostra_4. alguem sabe o que pode ter ocorrido? -- 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.

2 1

movimentar legenda eixo Y
by Marcelo Claro de Souza 01 Set '11

01 Set '11

Gostaria de saber como faço para alterar a distância das legendas nos eixos x e y, pois quando escrevo cm2 g-1 no eixo Y os expontes aparecem cortados quando salvo a figura em pdf ou jpg. Muito obrigado.

5 4

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