Ajuda com a função biplot.
Boa tarde senhores!! Estou com dois problemas com a função biplot. O primeiro é o seguinte: Quando tento colorir os pontos referentes a cada grupo, só sai uma cor em todos os grupos. Vejam. data(iris)#Já está na base do R pca <- prcomp(~Petal.Length + Sepal.Length + Petal.Width + Sepal.Width, data=iris, scale=TRUE) biplot(pca, xlabs=rep(c("1","2","3"),c(50,50,50)),col=c("red","blue","green"))#só é plotada a cor "red". Um outro problema, é que quando tento fazer a mesma coisa na minha base de dados, a função solta o seguinte erro. head(dmof) LIN IDADE BL REP TEMPO ESTA CP CC AC A1 A2 A3 D1 D2 L1 1 Verm 85 1 1 1 outono 8.3 25.9 26.4 32.1 28.1 11.8 38.2 11.9 16.5 2 Verm 85 1 2 1 outono 10.5 31.1 31.3 42.2 37.9 14.3 45.2 18.4 19.0 3 Verm 85 1 3 1 outono 9.2 28.1 28.4 36.7 30.1 11.8 42.5 16.2 17.1 4 Verm 85 1 4 1 outono 9.8 28.2 27.4 38.1 30.7 12.4 39.7 16.3 16.8 5 Verm 85 1 5 1 outono 9.1 28.4 29.2 36.2 31.2 12.4 42.7 18.2 18.4 6 Verm 85 1 6 1 outono 8.5 25.3 26.1 34.1 29.8 11.7 35.5 20.4 17.2 L2 L3 1 9.2 4.6 2 10.8 6.0 3 9.3 4.6 4 9.2 4.5 5 9.6 5.0 6 10.2 4.5 dmof$LIN Levels: Com Tai UFLA Verm pca1 <- prcomp(~CP+CC+AC+A1+A2+A3+D1+D2+L1+L2+L3, data=dmof, scale=TRUE) biplot(pca1, xlabs=rep(c("1","2","3","4"),c(200,200,200,200)), xlim=c(-0.1,0.1)) Erro em dimnames(x) <- list(xlabs, dimnames(x)[[2L]]) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Calls: biplot -> biplot.prcomp -> biplot.default Não entendo o erro, pois estou obedecendo as dimensões da minha base de dados. dim(dmof) 800 17 Se alguém puder me ajudar eu agradeço. Allaman (S,f,P) -- View this message in context: http://r-br.2285057.n4.nabble.com/Ajuda-com-a-funcao-biplot-tp3399760p339976... Sent from the R-br mailing list archive at Nabble.com.
O Ivan, Utiliza o pacote bpca. Ele é bem mais flexível que a função biplot! Valeu!!! Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155 ----- Mensagem original ---- De: ivanalaman <ivanalaman@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Quarta-feira, 23 de Março de 2011 12:22:39 Assunto: [R-br] Ajuda com a função biplot. Boa tarde senhores!! Estou com dois problemas com a função biplot. O primeiro é o seguinte: Quando tento colorir os pontos referentes a cada grupo, só sai uma cor em todos os grupos. Vejam. data(iris)#Já está na base do R pca <- prcomp(~Petal.Length + Sepal.Length + Petal.Width + Sepal.Width, data=iris, scale=TRUE) biplot(pca, xlabs=rep(c("1","2","3"),c(50,50,50)),col=c("red","blue","green"))#só é plotada a cor "red". Um outro problema, é que quando tento fazer a mesma coisa na minha base de dados, a função solta o seguinte erro. head(dmof) LIN IDADE BL REP TEMPO ESTA CP CC AC A1 A2 A3 D1 D2 L1 1 Verm 85 1 1 1 outono 8.3 25.9 26.4 32.1 28.1 11.8 38.2 11.9 16.5 2 Verm 85 1 2 1 outono 10.5 31.1 31.3 42.2 37.9 14.3 45.2 18.4 19.0 3 Verm 85 1 3 1 outono 9.2 28.1 28.4 36.7 30.1 11.8 42.5 16.2 17.1 4 Verm 85 1 4 1 outono 9.8 28.2 27.4 38.1 30.7 12.4 39.7 16.3 16.8 5 Verm 85 1 5 1 outono 9.1 28.4 29.2 36.2 31.2 12.4 42.7 18.2 18.4 6 Verm 85 1 6 1 outono 8.5 25.3 26.1 34.1 29.8 11.7 35.5 20.4 17.2 L2 L3 1 9.2 4.6 2 10.8 6.0 3 9.3 4.6 4 9.2 4.5 5 9.6 5.0 6 10.2 4.5 dmof$LIN Levels: Com Tai UFLA Verm pca1 <- prcomp(~CP+CC+AC+A1+A2+A3+D1+D2+L1+L2+L3, data=dmof, scale=TRUE) biplot(pca1, xlabs=rep(c("1","2","3","4"),c(200,200,200,200)), xlim=c(-0.1,0.1)) Erro em dimnames(x) <- list(xlabs, dimnames(x)[[2L]]) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Calls: biplot -> biplot.prcomp -> biplot.default Não entendo o erro, pois estou obedecendo as dimensões da minha base de dados. dim(dmof) 800 17 Se alguém puder me ajudar eu agradeço. Allaman (S,f,P) -- View this message in context: http://r-br.2285057.n4.nabble.com/Ajuda-com-a-funcao-biplot-tp3399760p339976... Sent from the R-br mailing list archive at Nabble.com. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (2)
-
Fabio Mathias Corrêa -
ivanalaman