- R-br - DINF UFPR

ajuda xlsReadWrite
by Vinicius Brito Rocha 13 Jul '11

13 Jul '11

Srs. alguem pode me mandar por email o pacote xlsReadWrite 1.54 do http://www.swissr.org/download estou no trabalho e nao consigo baixar. e ainda estou com o problema do rJava e não consigo usar o xlsx 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.

1 0

Ajuda (possivelmente simples) para reestruturar um data.frame
by Jose Claudio Faria 13 Jul '11

13 Jul '11

Pessoal, Ajuda ai que estou travado! #Tenho isso (desconsiderar a numeração das linhas) tra r y 1 A 1 2.70 2 A 2 3.02 3 A 3 3.25 4 A 4 3.21 5 A 5 3.04 6 A 6 3.07 7 B 1 1.16 8 B 2 1.16 9 B 3 0.91 10 B 4 0.93 11 B 5 0.82 12 B 6 0.71 13 C 1 3.01 14 C 2 2.91 15 C 3 3.14 16 C 4 2.96 17 C 5 3.01 18 C 6 2.93 19 D 1 4.88 20 D 2 4.68 21 D 3 5.20 22 D 4 4.94 23 D 5 5.34 24 D 6 5.07 25 E 1 4.89 26 E 2 5.26 27 E 3 5.07 28 E 4 5.14 29 E 5 5.47 30 E 6 4.98 # Para economizar espaço em um relatório (que uso muito) preciso disso (idem): tra r y tra r y tra r y tra r y tra r y 1 A 1 2.70 B 1 1.16 C 1 3.01 D 1 4.88 E 1 4.89 2 A 2 3.02 B 2 1.16 C 2 2.91 D 2 4.68 E 2 5.26 3 A 3 3.25 B 3 0.91 C 3 3.14 D 3 5.20 E 3 5.07 4 A 4 3.21 B 4 0.93 C 4 2.96 D 4 4.94 E 4 5.14 5 A 5 3.04 B 5 0.82 C 5 3.01 D 5 5.34 E 5 5.47 6 A 6 3.07 B 6 0.71 C 6 2.93 D 6 5.07 E 6 4.98 # Fiz na "mão grande" usando "cbind" assim dados.tb <- cbind(dados.q1[1:6,], dados.q1[7:12,], dados.q1[13:18,], dados.q1[19:24,], dados.q1[25:30,]) Contudo, creio que deve ter uma solução mais racional e genérica para quando mudar o número de tratamentos e repetições! Tentei split, pacotes doBy e reshape, mas não fui bem sucedido. Para reprodução de possível solução pode usar 'df' abaixo: df <- data.frame(gl(n=5, k=6, lab=LETTERS[1:5]), rep(1:6, 5), 1:30) names(df) <- c('tra', 'r', 'y') Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 3

Interpolação
by Nelson Betzek 13 Jul '11

13 Jul '11

É possível fazer interpolação pelo inverso da distancia elevado a uma potencia, ou só por Krigagem?

3 2

Modelo linear generalizado misto
by Marcelo Cardoso mello 13 Jul '11

13 Jul '11

Pessoal da lista ... tenho um modelo logístico que segue glm(Doença~Raça+Fazenda+Medicamento) ou glmer(Doença~Raça+Medicamento+(1|Fazenda)) O que vcs acham de colocar fazenda (animais em uma amostra de fazndas) como um fator aleatório no modelo? OU então colocar os próprios indivíduos como um fator aleatório? É uma opção melhor? Posso colocar o intercepto como aleatório também? E a odds ratio como fica? ... posso estimar a odds dos fatores fixos do modelo e dai qual função estima odds com teste de significância a partir de um modelo glmer? Uso logistic.display da epicalc para glm e para glmer?. Na verdade não sei se a glmer é a melhor função para o meu caso. Alguem conhece algum material que possa me ajudar com exemplos no R?Help.Marcelo.

5 5

.net+R
by Vinicius Brito Rocha 12 Jul '11

12 Jul '11

Srs. existe uma maneira de eviar resultados de R diretamente para .net? att. -- 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 2

Dúvida
by Edson Lira 12 Jul '11

12 Jul '11

Boa tarde R-br lista, estou com um pequeno problema. Tenho uma lista de exames que associo ao laboratório respectivo, vejam um exemplo abaixo: Descrição lab Quantidade 1 17 ALFA HIDROXI PROGESTERONA <NA> 1 2 ACIDO URICO Bioquímica 1 3 ACIDO URICO Bioquímica 1 4 ACIDO URICO Bioquímica 1 5 ACIDO URICO Bioquímica 1 6 ACIDO URICO Bioquímica 4 7 ACIDO URICO Bioquímica 9 8 ACIDO URICO Bioquímica 1 9 ACIDO URICO Bioquímica 1480 10 ACIDO URICO Bioquímica 37 11 ACIDO VALPROICO <NA> 1 12 ALANINA AMINOTRANSFERASE (ALT/GPT) IFCC <NA> 1 13 ALBUMINA Bioquímica 205 14 ALBUMINA Bioquímica 17 15 ALFA-FETO-PROTEINA HORMÔNIOS 3 16 ALFA-FETO-PROTEINA HORMÔNIOS 2 17 AMILASE Bioquímica 77 18 AMILASE Bioquímica 4 19 ANDROSTENEDIONA <NA> 1 20 ANDROSTENEDIONA <NA> 1 Para gerar os dados de contagem dos exames por laboratório. Estou implementando uma rotina para otimizar o trabalho de consolidação dos dados. O problema é que alguns exames que já estão associados ao respectivo laboratório quando não são realizados no período, eles não têm registros, mas no script estão relacionados, vejam a rotina abaixo: library(memisc) acj$lab<-recode(acj$Lab, "Bioquímica" <-c("ACIDO URICO","ALBUMINA","BILIRRUBINA DIRETA", "BILIRRUBINA INDIRETA","BILIRRUBINA TOTAL E FRAÇÕES","BILIRRUBINAS TOTAL", "COLESTEROL", "COLESTEROL HDL", "COLESTEROL LDL", "COLESTEROL VLDL", "CREATININA","DHL - DESIDROGENASE LÁCTICA", "FERRITINA", "FERRO SÉRICO", "FERRO SÉRICO LABTEST", "FOSFATASE ALCALINA", "FOSFATASE ÁCIDA FRAÇÃO PROSTÁTICA", "FOSFATASE ÁCIDA PROSTÁTICA (RIE)", "FOSFORO", "GAMA GT- GAMA-GLUTAMIL-TRANSFERASE", "GLICOSE", "GLICOSE COM FLUORETO", "GLICOSE POS-PRANDIAL", "GLICOSE-6-FOSFATO DEHIDROGENASE (GGFD)", "GLOBULINA", "HEMOGLOBINA GLICOSILADA", "LIPIDIOS TOTAIS", "LIPIDOGRAMA", "MAGNESIO", "POTÁSSIO", "PROTEINAS TOTAIS", "SATURAÇÃO DE TRANSFERRINA","SÓDIO","TRANSAMINASE OXALACÉTICA - AST", "TRANSAMINASE PIRÚVICA - ALT", "TRANSAMINASES", "TRANSAMINASES", "TRANSFERRINA", "ALFA-1 GLICOPROTEINA ÁCIDA", "ALT-TRANSAMINASE PIRÚVICA", "AMILASE", "AST-TRANSSAMINASE OXALACÉTICA", "CÁLCIO","CAPACIDADE DE OXIDAÇÃO DO FERRO","COLESTEROL TOTAL","COLESTEROL TOTAL E FRAÇÕES", "ESTUDO DO FERRO", "FOSFATASE ÁCIDA TOTAL", "TESTE DE ROG", "TRIGLICERÍDEOS", "URÉIA", "LÍTIO"), "Citoquímica" <-c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM")) Quando rodo esta rotina o R mostra a mensagem: Erro em factor(y, levels = 1:max(y), labels = newcodes) : invalid labels; length 13 should be 1 or 12 Além disso: Mensagens de aviso perdidas: In recode(acj$Lab, "Bioquímica" <- c("ACIDO URICO", "ALBUMINA", : recoding "Citoquímica" <- c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM") has no consequences O problema é que não existe nenhum exame associado ao laboratório de Citoquímica (FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM). ou seja não foi realizado nenhum exame, portanto, não aparecem no banco, como resolver este problema? []'s Edson Lira Estatístico Manaus-Amazonas

1 0

Equação do plano
by Helber Freitas 12 Jul '11

12 Jul '11

Olá, Pessoal. Fico feliz em poder contar com o grupo. Sou novo por aqui e aproveito para agradecer e parabenizar a iniciativa. Preciso de uma função que forneça a equação do plano dados 3 pontos distintos (x1,y1,z1), (x2,y2,z2) e (x3,y3,z3). Grato, Helber

3 2

RJAVA JRI (problema ao plotar gráficos)
by Emmanuel Arnhold 12 Jul '11

12 Jul '11

Estou tentando fazer o R "conversar" com o JAVA. Estou usando RJAVA e JRI. Tudo certo quando rodo algumas análises. O problema é que os gráficos não estão aparecendo na janela gráfica. Alguma dica?? O que pode estar errado?? Agradeço desde já. Emmanuel

2 2

Re: [R-br] Derivadas Parciais e Ponto Crítico!
by Thiago De paula protásio 11 Jul '11

11 Jul '11

Obrigado Walmes pela dica do pacote... Resolveu o meu problema....a função rms retorna os pontos críticos e os autovalores da equação característica: IΒ −λI = 0 E a análise dos sinais dos autovalores são fundamentais.... Novamente agradeço-lhe!Abraço! 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

Derivadas Parciais e Ponto Crítico!
by Thiago De paula protásio 11 Jul '11

11 Jul '11

Pessoal; Recorro novamente a ajuda de vcs...Estou trabalhando com a seguinte superfície de resposta.... z = 31.13884 -2.23631*x +0.04441*x^2 + 1.25980*y -0.01134*y^2 Contudo, gostaria de determinar no R os pontos críticos e verificar se a natureza do ponto: sela, máximo ou mínimo. A ideia consiste nas derivadas parciais....Dessa forma, gostaria de saber se existe alguma função no R que me auxilie nesse problema. 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

2 1