- R-br - DINF UFPR

Factanal - Erro!
by Lincoln Thiago Neves 05 Fev '12

05 Fev '12

Amigos, Gostaria de ma ajuda com o procedimento de Análise Fatorial no R. Estou tentando rodar um programa com o factanal, mas o R apresenta o seguinte erro: Erro em solve.default(cv) : sistema é computacionalmente singular: condição recíproca número = 2.65204e-20 Gostaria de entender sobre esse erro. Quem puder me dar essa ajuda serei grato! O programa que estou rodando é o seguinte: escola_2<-read.csv2('Input_Q_escola2.csv') outputaf_2<-factanal(escola_2,factors= 9,scores='regression') rot_escola_2<-varimax(loadings(outputaf_2), normalize = FALSE) matrizL_2<-outputaf$loadings escores_2<-outputaf$score write.table(file='matrizL_2.csv', matrizL_2, sep=';', dec=',') write.table(file='escores_2.csv', escores_2, sep=';', dec=',') O erro acontece abaixo do "outputaf_2". A base "Input_Q_escola2.csv" segue anexa. Obrigado! -- Lincoln Neves MSN: lincolnthiago(a)hotmail.com E-mail: lincolnthiago.neves(a)gmail.com Web: http://ltnestsof.blogspot.com/ Cel: 8883-3632

2 1

Ponto máximo equações quadráticas
by Alexandre Santos 04 Fev '12

04 Fev '12

Boa tarde pessoal, Tenho várias equações quadráticas como por exemplo: fx <- function(x) -6.685+8.227*x-0.841*x^2, gostaria de calcular as coordenadas dos pontos máximos/mínimos para as equações e gostaria de saber se existe uma função para tal ou tenho que encontrar a derivada, depois encontrar as raízes, etc. Muito obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

6 6

cluster hierarquico ward
by Vinicius Brito Rocha 03 Fev '12

03 Fev '12

Sr. estou com uma dúvida sobre o método hierarquico de clusterização de ward. Não em como utiliza-lo no R e sim como funciona o seu algoritmo? Estou com uma bibliografia, mas honestamente não estou entendendo. Alguém poderia me dar uma luz? 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.

2 2

mudar cor grafico - Plot(x,y)
by jose coelho 01 Fev '12

01 Fev '12

Vivam todos, Alguém me diz com mudar a cor num gráfico feito com Plot (xy), mas mudar a penas a cor da área de desenho. tentei com bg="gray" mas altera toda a área! Muito obrigado -- jmmcoelho

3 4

Factor para Dummy
by Vinicius Brito Rocha 01 Fev '12

01 Fev '12

Sr., Acredito que essa seja fácil. Qual o pacote trasforma uma variável tipo factor em uma matrix de variáveis Dummy? 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.

2 2

ajuste com nls
by Gilson Sanchez 31 Jan '12

31 Jan '12

Pessoal mais uma duvida, eu fiz um ajuste com o nls, mas eu quero saber como utilizar a e b ou seja se tem algum comando para que eu possa utilizar os valores que foram estimadas pela equação Segue o CMR que estou trabalhando # Planta 1 planta1=foto[1:90,-3:-6] names(planta1) [1] "especie" "esp" "PAR" "Tleaf" "W" "Pn" "E" [8] "C" "int_CO2" colnames(planta1)<-c("especiea", "espa", "PARa","Tleafa", "Wa", "Pna", "Ea", + "Ca", "int_CO2a") # Obtendo PAR planta1.luz =tapply(planta1$PARa, INDEX=planta1$PARa, FUN=mean) planta1.luz 0 25 50 75 100 500 1000 1500 2000 0 25 50 75 100 500 1000 1500 2000 planta1.me.Pna=tapply(planta1$Pna, INDEX=planta1$PARa, FUN=mean) planta1.me.Pna 0 25 50 75 100 500 1000 1500 2000 -0.428 0.144 0.411 1.275 2.052 4.349 5.082 4.078 3.794 planta1.ajuste=nls((planta1.me.Pna~(a+0.428)*(1-exp(-b*planta1.luz/(a+0.428)))-0.428), + start=list(a=200, b=0.5)) planta1.ajuste Nonlinear regression model model: planta1.me.Pna ~ (a + 0.428) * (1 - exp(-b * planta1.luz/(a + 0.428))) - 0.428 data: parent.frame() a b 4.37564 0.02943 residual sum-of-squares: 1.258 Number of iterations to convergence: 5 Achieved convergence tolerance: 3.175e-06 max=c(planta1.ajuste$a) = ????? alfa=c(plnat1.ajuste$b) = ??????? Desde já muito obrigado. -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

3 4

(sem assunto)
by Paulo Ricardo Silva Moreira 31 Jan '12

31 Jan '12

Boa tarde Galera, Gostaria de saber como faço para postar várias séries no mesmo gráfico, e como faço para colocar uma cor diferente para cada série, de forma a diferenciá-las. Att, Paulo Ricardo.

3 2

Inserir legenda em Boxplot
by Marcelo Claro de Souza 31 Jan '12

31 Jan '12

Tentei inserir uma legenda dentro de um boxplot utilizando o comando abaixo e sempre recebo a informção de erro inesperado. legend("topleft","(x,y)",pch=1,tittle=" A",inset=0.02) Teria alguma outra opção? Abraço e muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

4 3

Dúvida sobre teste
by Katia Alves Campos 31 Jan '12

31 Jan '12

2 1

Importar Microdados - arquivos .txt sem separação entre as colunas
by Roney Fraga 31 Jan '12

31 Jan '12

Caros, Diante da necessidade de importar microdados, principalmente do IBGE, por meio do pacote IBGEPesq()* *que vem apresentando alguns problemas, pergunto se algum usuário tem alguma sugestão para importar microdados alternativa ao pacote IBGEPesq(). Microdados consiste geralmente em grandes pesquisas que são disponibilizadas em arquivos .txt apenas com números sem nenhuma separação entre as colunas, sendo o arquivo dicionário disponibilizado a parte em outro arquivo .txt, nesse arquivo dicionário consta em que coluna começa e termina cada variável. No link http://www.datafilehost.com/download-2565fe79.html segue exemplo do arquivo .txt com os dados da PNAD 2009 banco domicílios. E no link http://www.datafilehost.com/download-fe9c31a8.html segue exemplos de dicionários para o SAS, disponibilizado pelo IBGE, e para o Stata que fiz com base no dicionário anterior. Se necessário o mesmo dicionário pode ser encontrado no R através dos comandos install.packages('dicionariosIBGE')library(dicionariosIBGE) data(dicPNAD2009) dicdom2009 No Stata uma vez que tenha o dicionário, na extensão .dct, e o arquivo .txt com os dados a serem lidos basta executar o comando infile using “/Caminho/Arquivo/dicionario.dct” no R não sei se é possível seguir o mesmo raciocínio para importar dados em arquivos .txt, caso algum usuário possa contribuir será uma grande ajuda. Os arquivos em .txt geralmente são grandes, por exemplo, a PNAD 2009 arquivo pessoas tem 320MB. Uma vez ultrapassada a dificuldade de importar os dados, será mais fácil utilizar o R para trabalhar com os dados públicos das seguintes pesquisas: PNAD, POF, PME, Enade, Censo Escolar, Censo da Educação Superior, Censo dos Profissionais do Magistério, Saeb, Enem, Provão, PNERA, Prova Brasil, Pesquisa de Ações Discriminatórias no Âmbito Escolar, etc... Atenciosamente Roney

6 6