Gráficos de contorno (Superfície de resposta)

Ver o tópico anterior Ver o tópico seguinte Ir em baixo

Gráficos de contorno (Superfície de resposta)

Mensagem  marciel lelis em Qua Ago 08, 2012 9:25 pm

Estou tentando criar um gráfico de contorno, o gráfico da certo mas quando tento mudar o nome do eixo y aparece a seguinte mensagem:

Erro em localPlotWindow(xlim, ylim, ...) :
argumento formal "ylab" corresponde a múltiplos argumentos especificados

estou usando os seguintes comandos

> contour(sup20, ~ Doses_de_K + Doses_de_S, date = TRUE, xlab="Doses de K (mg/dm³)", ylab="Doses de S (mg/dm³)")

É uma analise de superfície de resposta e abaixo segue uma parte dos dados e os comandos que utilizei nas análises.

Doses_de_K Doses_de_S Bloco H DC MSPA MSR MST
1 0 0 1 8.50 2.29 0.49 0.23 0.72
2 0 20 1 12.20 1.69 0.51 0.22 0.73
3 0 40 1 17.00 3.35 1.04 0.38 1.42
4 0 60 1 44.80 5.22 4.32 1.68 6.00
5 0 80 1 41.00 5.63 3.32 1.14 4.46
6 50 0 1 14.00 3.65 0.24 0.13 0.37
7 50 20 1 51.63 6.30 5.18 2.04 7.23
8 50 40 1 38.50 5.08 2.90 0.97 3.87
9 50 60 1 52.20 5.54 1.26 1.95 3.21
10 50 80 1 45.50 6.24 4.33 1.28 5.61
11 100 0 1 32.20 5.47 3.25 1.40 4.65
12 100 20 1 18.50 3.19 0.65 0.40 1.05
13 100 40 1 15.00 5.04 2.64 1.50 4.14
14 100 60 1 12.00 2.49 0.36 0.10 0.46
15 100 80 1 6.30 1.46 0.12 0.03 0.15
16 150 0 1 11.50 2.68 0.49 0.29 0.78
17 150 20 1 23.00 3.14 1.08 0.36 1.44
18 150 40 1 27.50 4.65 3.08 1.09 4.17
19 150 60 1 23.80 2.71 1.09 0.23 1.32
20 150 80 1 41.23 6.47 5.58 2.78 8.36
21 200 0 1 18.00 3.37 1.03 0.37 1.40
22 200 20 1 51.00 6.64 5.40 2.65 8.05
23 200 40 1 24.97 4.44 2.08 0.75 2.84
24 200 60 1 10.00 2.97 5.83 3.11 8.94
25 200 80 1 19.00 2.47 0.61 0.22 0.83
26 250 0 1 13.00 2.11 0.40 0.18 0.58
27 250 20 1 30.50 5.19 3.33 0.50 3.83
28 250 40 1 8.00 3.15 0.75 0.44 1.19
29 250 60 1 29.50 4.69 2.04 0.99 3.03
30 250 80 1 14.00 2.63 0.58 0.17 0.75
31 300 0 1 15.50 3.12 0.75 0.28 1.03
32 300 20 1 20.00 3.91 1.48 0.52 2.00
33 300 40 1 22.50 4.46 1.66 0.57 2.23
34 300 60 1 13.00 2.89 0.56 0.22 0.78
35 300 80 1 35.50 6.31 2.77 1.08 3.85


dados<-read.table("dados.txt",head=TRUE)
dados
dados <- data.frame(dados)
summary(dados)
variaveis <- dados[,-c(1:3)]

require(rsm)
sup20 <- rsm(H ~ SO(Doses_de_K,Doses_de_S), data = dados)
sup20
summary(sup20)
sup21 <- rsm( DC ~ SO(Doses_de_K,Doses_de_S), data = dados)
sup21
summary(sup21)
sup22 <- rsm( MSPA ~ SO(Doses_de_K,Doses_de_S), data = dados)
sup22
summary(sup22)
sup23 <- rsm( MSR~ SO(Doses_de_K,Doses_de_S), data = dados)
sup23
summary(sup23)
sup24 <- rsm( MST~ SO(Doses_de_K,Doses_de_S), data = dados)
sup24
summary(sup24)





marciel lelis

Mensagens : 4
Data de inscrição : 27/06/2012

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Ver o tópico anterior Ver o tópico seguinte Voltar ao Topo

- Tópicos similares

 
Permissão deste fórum:
Você não pode responder aos tópicos neste fórum