[R-br] Ajuda com a função biplot.

ivanalaman ivanalaman em yahoo.com.br
Quarta Março 23 12:22:39 BRT 2011


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-tp3399760p3399760.html
Sent from the R-br mailing list archive at Nabble.com.


Mais detalhes sobre a lista de discussão R-br