R help es - Feb 2022 - colores de los plots

Estimados
En este script quisieramos poder manejar los colores de los plots y no
vemos como se puede hacer
Por lo que apreciariamos su ayuda
  saludos cordiales.
José


rm(list=ls(all=TRUE))

library(sets)

sets_options("universe", seq(1, 45, 1))

#VARIABLES   no cierra



   variables <- set(
       temp= fuzzy_partition(varnames = c(normal= 36.5,
   moderada= 37.5, alta= 38, muyalta=39), sd = 0.2),

   disnea = fuzzy_partition(varnames = c(no = 0, poca = 20,mucha=40),sd=2.5),

   gusto = fuzzy_partition(varnames = c(no = 0, poca = 30,mucha=40),sd=2.5),

   garganta = fuzzy_partition(varnames = c(no = 0, poca = 25,mucha=40),sd=2.5),

   severidad = fuzzy_partition(varnames = c(baja = 10, media = 20,alta=40),
                               FUN = fuzzy_cone, radius = 10)
   )



   #RULES

   rules <- set(

     fuzzy_rule(temp %is% muyalta && disnea%is% mucha &&
gusto%is%mucha&&
                  garganta%is%mucha, severidad %is% alta),
     fuzzy_rule(temp %is% muyalta&& disnea%is% mucha, severidad %is%
alta),
     fuzzy_rule(temp %is% normal && disnea%is% no &&
gusto%is%no&&
                  garganta%is%no, severidad %is% baja),
     fuzzy_rule(temp %is% normal && disnea%is% no, severidad %is% baja),
     fuzzy_rule(temp %is% moderada && disnea%is% poca &&
gusto%is%poca&&
                  garganta%is%poca, severidad %is% media),
     fuzzy_rule(temp %is% moderada && disnea%is% poca, severidad %is%
media)

   )




model <- fuzzy_system(variables, rules)
print(model)
plot(model)


hacer========================
example.1 <- fuzzy_inference(model, list(temp = 37,
disnea=40,garganta=40,gusto=40))

gset_defuzzify(example.1, "centroid")  # NaN

plot(example.1)



example.2 <- fuzzy_inference(model, list(Edad = 70, IMC=30,Tas=140,
                                          tabaquismo= 40, glicemia=140))
gset_defuzzify(example.2, "centroid")  # NaN

plot(example.2)


example.3 <- fuzzy_inference(model, list(Edad = 75, IMC=30,Tas=140,
                                          tabaquismo= 40, glicemia=150))
gset_defuzzify(example.3, "centroid")  # NaN

plot(example.3)
sets_options("universe", NULL)


--
Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece
Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de
Salud. La persona que envia este correo asume el compromiso de usar el servicio
a tales fines y cumplir con las regulaciones establecidas

Infomed: http://www.sld.cu/