thr3ads.net - R help es - [R-es] Matriz funcion de pertenencia [Jan 2013]

If this information is useful, please help other people find it:
Share via:

Ivan Santos

2013-Jan-16 10:41 UTC

[R-es] Matriz funcion de pertenencia

Buenos dias

Estoy modelando un sistema de control difuso, y tengo un problema a la hora de
obtener la matriz con las reglas difusas. Me explico.

El sistema modela la velocidad de un ventilador, en funcion de la temperatura y
la humedad

1. He definido las funciones de pertenencia de temperatura, humedad y velocidad

#Función de pertenencia temperatura 

Tbaja<-vector(length=length(T), mode="numeric")
Tbaja[T>40 & T<=60] <- 1
Tbaja[T>60 & T<=80] <- (80-T[T>60 & T<=80])/(80-60)
Tbaja[T>80] <- 0

Tmedia<-vector(length=length(T), mode="numeric")
Tmedia[T<60] <- 0
Tmedia[T>60 & T<=80] <- (T[T>60 & T<=80]-60)/(80-60)
Tmedia[T>80 & T<=100] <- (100-T[T>80 & T<=100])/(100-80)
Tmedia[T>100] <- 0

Talta<-vector(length=length(T), mode="numeric")
Talta[T<80] <- 0
Talta[T>=80 & T<=100] <- (T[T>=80 & T<=100]-80)/(100-80)
Talta[T>100] <- 1

#Funcion pertenencia humedad

Hbaja<-vector(length=length(H), mode="numeric")
Hbaja[H<25] <- 1
Hbaja[H>=25 & H<=50] <- (50-H[H>=25 & H<=50])/(50-25)
Hbaja[H>50] <- 0

Hmedia<-vector(length=length(H), mode="numeric")
Hmedia[H<25] <- 0
Hmedia[H>=25 & H<50] <- (H[H>=25 & H<50]-25)/(50-25)
Hmedia[H>=50 & H<=75] <- (75-H[H>=50 & H<=75])/(75-50)
Hmedia[H>75] <- 0

Halta<-vector(length=length(H), mode="numeric")
Halta[H<50] <- 0
Halta[H>=50 & H<=75] <- (H[H>=50 & H<=75]-50)/(75-50)
Halta[H>75] <- 1

#Funcion pertenencia velocidad

Wbaja<-vector(length=length(W), mode="numeric")
Wbaja[W<250] <- 1
Wbaja[W>250 & W<=500] <- (500-W[W>250 &
W<=500])/(500-250)
Wbaja[W>500] <- 0

Wmedia<-vector(length=length(W), mode="numeric")
Wmedia[W<250] <- 0
Wmedia[W>=250 & W<500] <- (W[W>=250 &
W<500]-250)/(500-250)
Wmedia[W>=500 & W<=750] <- (750-W[W>=500 &
W<=750])/(750-500)
Wmedia[W>750] <- 0

Walta<-vector(length=length(W), mode="numeric")
Walta[W<500] <- 0
Walta[W>=500 & W<=750] <- (W[W>=500 &
W<=750]-500)/(750-500)
Walta[W>750] <- 1

2. Una vez obtenidas las funciones de pertenencia, obtenemos las implicaciones.
#Construimos las reglas de implicacion
#velocidad baja, temperatura baja, velocidad baja
if (Tbaja && Hbaja) Wbaja

#velocidad media, temperatura baja, velocidad media
if (Tmedia && Hbaja) Wmedia

#velocidad alta, temperatura baja, velocidad baja
if (Talta && Hbaja) Wmedia

3. Con las reglas deberia construir una matriz, para obtener las funciones de
pertenencia de las reglas. Pero aqui me encuentro con el problema de como meter
dichas reglas, en una matriz 3x3. He probado a asignar directamente las reglas a
una posicion concreta pero me da error.

4. Como podria obtener la superficie de control.

Muchas gracias
 		 	   		  
	[[alternative HTML version deleted]]

daniel

2013-Jan-16 12:34 UTC

head link

[R-es] Matriz funcion de pertenencia

Si entiendo bien, quieres asignar a una matriz objetos que no son vectores
atómicos. Francamente, nunca lo intenté, pero ese puede ser el problema.
¿No te sirven las listas?

Daniel Merino


El 16 de enero de 2013 07:41, Ivan Santos <lyaran2003@hotmail.com>
escribió:
> Buenos dias
>
> Estoy modelando un sistema de control difuso, y tengo un problema a la
> hora de obtener la matriz con las reglas difusas. Me explico.
>
> El sistema modela la velocidad de un ventilador, en funcion de la
> temperatura y la humedad
>
> 1. He definido las funciones de pertenencia de temperatura, humedad y
> velocidad
>
> #Función de pertenencia temperatura
>
> Tbaja<-vector(length=length(T), mode="numeric")
> Tbaja[T>40 & T<=60] <- 1
> Tbaja[T>60 & T<=80] <- (80-T[T>60 & T<=80])/(80-60)
> Tbaja[T>80] <- 0
>
> Tmedia<-vector(length=length(T), mode="numeric")
> Tmedia[T<60] <- 0
> Tmedia[T>60 & T<=80] <- (T[T>60 & T<=80]-60)/(80-60)
> Tmedia[T>80 & T<=100] <- (100-T[T>80 &
T<=100])/(100-80)
> Tmedia[T>100] <- 0
>
> Talta<-vector(length=length(T), mode="numeric")
> Talta[T<80] <- 0
> Talta[T>=80 & T<=100] <- (T[T>=80 &
T<=100]-80)/(100-80)
> Talta[T>100] <- 1
>
> #Funcion pertenencia humedad
>
> Hbaja<-vector(length=length(H), mode="numeric")
> Hbaja[H<25] <- 1
> Hbaja[H>=25 & H<=50] <- (50-H[H>=25 &
H<=50])/(50-25)
> Hbaja[H>50] <- 0
>
> Hmedia<-vector(length=length(H), mode="numeric")
> Hmedia[H<25] <- 0
> Hmedia[H>=25 & H<50] <- (H[H>=25 & H<50]-25)/(50-25)
> Hmedia[H>=50 & H<=75] <- (75-H[H>=50 &
H<=75])/(75-50)
> Hmedia[H>75] <- 0
>
> Halta<-vector(length=length(H), mode="numeric")
> Halta[H<50] <- 0
> Halta[H>=50 & H<=75] <- (H[H>=50 &
H<=75]-50)/(75-50)
> Halta[H>75] <- 1
>
> #Funcion pertenencia velocidad
>
> Wbaja<-vector(length=length(W), mode="numeric")
> Wbaja[W<250] <- 1
> Wbaja[W>250 & W<=500] <- (500-W[W>250 &
W<=500])/(500-250)
> Wbaja[W>500] <- 0
>
> Wmedia<-vector(length=length(W), mode="numeric")
> Wmedia[W<250] <- 0
> Wmedia[W>=250 & W<500] <- (W[W>=250 &
W<500]-250)/(500-250)
> Wmedia[W>=500 & W<=750] <- (750-W[W>=500 &
W<=750])/(750-500)
> Wmedia[W>750] <- 0
>
> Walta<-vector(length=length(W), mode="numeric")
> Walta[W<500] <- 0
> Walta[W>=500 & W<=750] <- (W[W>=500 &
W<=750]-500)/(750-500)
> Walta[W>750] <- 1
>
> 2. Una vez obtenidas las funciones de pertenencia, obtenemos las
> implicaciones.
> #Construimos las reglas de implicacion
> #velocidad baja, temperatura baja, velocidad baja
> if (Tbaja && Hbaja) Wbaja
>
> #velocidad media, temperatura baja, velocidad media
> if (Tmedia && Hbaja) Wmedia
>
> #velocidad alta, temperatura baja, velocidad baja
> if (Talta && Hbaja) Wmedia
>
> 3. Con las reglas deberia construir una matriz, para obtener las funciones
> de pertenencia de las reglas. Pero aqui me encuentro con el problema de
> como meter dichas reglas, en una matriz 3x3. He probado a asignar
> directamente las reglas a una posicion concreta pero me da error.
>
> 4. Como podria obtener la superficie de control.
>
> Muchas gracias
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Daniel

	[[alternative HTML version deleted]]

Ivan Santos

2013-Jan-16 12:49 UTC

head link

[R-es] Matriz funcion de pertenencia

Es mi primera incursion en R, y no se muy bien como hacerlo, por eso preguntaba.
¿Como podria hacerlo con una lista?

Muchas gracias

Date: Wed, 16 Jan 2013 09:34:27 -0300
Subject: Re: [R-es] Matriz funcion de pertenencia
From: daniel319@gmail.com
To: lyaran2003@hotmail.com
CC: r-help-es@r-project.org

Si entiendo bien, quieres asignar a una matriz objetos que no son vectores
atómicos. Francamente, nunca lo intenté, pero ese puede ser el problema. ¿No te
sirven las listas?
Daniel Merino


El 16 de enero de 2013 07:41, Ivan Santos <lyaran2003@hotmail.com>
escribió:

Buenos dias



Estoy modelando un sistema de control difuso, y tengo un problema a la hora de
obtener la matriz con las reglas difusas. Me explico.



El sistema modela la velocidad de un ventilador, en funcion de la temperatura y
la humedad



1. He definido las funciones de pertenencia de temperatura, humedad y velocidad



#Función de pertenencia temperatura



Tbaja<-vector(length=length(T), mode="numeric")

Tbaja[T>40 & T<=60] <- 1

Tbaja[T>60 & T<=80] <- (80-T[T>60 & T<=80])/(80-60)

Tbaja[T>80] <- 0



Tmedia<-vector(length=length(T), mode="numeric")

Tmedia[T<60] <- 0

Tmedia[T>60 & T<=80] <- (T[T>60 & T<=80]-60)/(80-60)

Tmedia[T>80 & T<=100] <- (100-T[T>80 & T<=100])/(100-80)

Tmedia[T>100] <- 0



Talta<-vector(length=length(T), mode="numeric")

Talta[T<80] <- 0

Talta[T>=80 & T<=100] <- (T[T>=80 & T<=100]-80)/(100-80)

Talta[T>100] <- 1



#Funcion pertenencia humedad



Hbaja<-vector(length=length(H), mode="numeric")

Hbaja[H<25] <- 1

Hbaja[H>=25 & H<=50] <- (50-H[H>=25 & H<=50])/(50-25)

Hbaja[H>50] <- 0



Hmedia<-vector(length=length(H), mode="numeric")

Hmedia[H<25] <- 0

Hmedia[H>=25 & H<50] <- (H[H>=25 & H<50]-25)/(50-25)

Hmedia[H>=50 & H<=75] <- (75-H[H>=50 & H<=75])/(75-50)

Hmedia[H>75] <- 0



Halta<-vector(length=length(H), mode="numeric")

Halta[H<50] <- 0

Halta[H>=50 & H<=75] <- (H[H>=50 & H<=75]-50)/(75-50)

Halta[H>75] <- 1



#Funcion pertenencia velocidad



Wbaja<-vector(length=length(W), mode="numeric")

Wbaja[W<250] <- 1

Wbaja[W>250 & W<=500] <- (500-W[W>250 &
W<=500])/(500-250)

Wbaja[W>500] <- 0



Wmedia<-vector(length=length(W), mode="numeric")

Wmedia[W<250] <- 0

Wmedia[W>=250 & W<500] <- (W[W>=250 &
W<500]-250)/(500-250)

Wmedia[W>=500 & W<=750] <- (750-W[W>=500 &
W<=750])/(750-500)

Wmedia[W>750] <- 0



Walta<-vector(length=length(W), mode="numeric")

Walta[W<500] <- 0

Walta[W>=500 & W<=750] <- (W[W>=500 &
W<=750]-500)/(750-500)

Walta[W>750] <- 1



2. Una vez obtenidas las funciones de pertenencia, obtenemos las implicaciones.

#Construimos las reglas de implicacion

#velocidad baja, temperatura baja, velocidad baja

if (Tbaja && Hbaja) Wbaja



#velocidad media, temperatura baja, velocidad media

if (Tmedia && Hbaja) Wmedia



#velocidad alta, temperatura baja, velocidad baja

if (Talta && Hbaja) Wmedia



3. Con las reglas deberia construir una matriz, para obtener las funciones de
pertenencia de las reglas. Pero aqui me encuentro con el problema de como meter
dichas reglas, en una matriz 3x3. He probado a asignar directamente las reglas a
una posicion concreta pero me da error.




4. Como podria obtener la superficie de control.



Muchas gracias



        [[alternative HTML version deleted]]




_______________________________________________

R-help-es mailing list

R-help-es@r-project.org

https://stat.ethz.ch/mailman/listinfo/r-help-es




-- 
Daniel
 		 	   		  
	[[alternative HTML version deleted]]

Leonardo Hernández Pérez

2013-Jan-16 12:59 UTC

head link

[R-es] Interpretación del chisq.test

Hola amigos:

Estoy haciendo un curso y estamos estudiando el chisq.test este que les 
dejo debajo es un ejemplo que tiene rstudio.

Pueden decirme como interpretar el resultado?

Un saludo

Leonardo

x <- trunc(5 * runif(100))
chisq.test(table(x))


X-squared = 3.9, df = 4, p-value = 0.4197
------------ próxima parte ------------
An embedded and charset-unspecified text was scrubbed...
Name: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130116/0fcfcf7e/attachment.pl>

Javier Marcuzzi

2013-Jan-16 13:15 UTC

head link

[R-es] Interpretación del chisq.test

Estimado Leonardo

Si copia y pega lo siguiente se dara cuenta
x <- trunc(5 * runif(100))
x
table(x)
chisq.test(table(x))



El 16 de enero de 2013 09:59, Leonardo Hernández Pérez <
leonardo.hernandez@etecsa.cu> escribió:
>
> Hola amigos:
>
> Estoy haciendo un curso y estamos estudiando el chisq.test este que les
> dejo debajo es un ejemplo que tiene rstudio.
>
> Pueden decirme como interpretar el resultado?
>
> Un saludo
>
> Leonardo
>
> x <- trunc(5 * runif(100))
> chisq.test(table(x))
>
>
> X-squared = 3.9, df = 4, p-value = 0.4197
> ---
> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE
> running at host imx3.etecsa.cu
> Visit our web-site: <http://www.kaspersky.com>,
<http://www.viruslist.com>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
	[[alternative HTML version deleted]]

Leonardo Hernández Pérez

2013-Jan-16 14:46 UTC

head link

[R-es] Interpretación del chisq.test

Hola Javier:

Lo que necesito es que me digan cual es la interpretación literal del 
resultados, por ejemplo para es resultado "X-squared = 3.9, df = 4, 
p-value = 0.4197" el valor de p me dice que no existen diferencias 
estadísticamente significativas pero en table(x) lo que tengo es la 
frecuencia de ocurrencia de 0 , 1, 2, 3 y 4 entonces supongo que entre 
estos valores (la frecuencia para cada uno de estos números) no existen 
diferencias estadísticamente significativas.

Es así?

Un saludo

Leonardo


On 16/01/13 08:15, Javier Marcuzzi wrote:> Estimado Leonardo
>
> Si copia y pega lo siguiente se dara cuenta
> x <- trunc(5 * runif(100))
> x
> table(x)
> chisq.test(table(x))
>
>
>
> El 16 de enero de 2013 09:59, Leonardo Hernández Pérez 
> <leonardo.hernandez en etecsa.cu <mailto:leonardo.hernandez en
etecsa.cu>>
> escribió:
>
>
>     Hola amigos:
>
>     Estoy haciendo un curso y estamos estudiando el chisq.test este
>     que les dejo debajo es un ejemplo que tiene rstudio.
>
>     Pueden decirme como interpretar el resultado?
>
>     Un saludo
>
>     Leonardo
>
>     x <- trunc(5 * runif(100))
>     chisq.test(table(x))
>
>
>     X-squared = 3.9, df = 4, p-value = 0.4197
>     ---
>     This message was processed by Kaspersky Mail Gateway
>     5.6.28/RELEASE running at host imx3.etecsa.cu
<http://imx3.etecsa.cu>
>     Visit our web-site: <http://www.kaspersky.com>,
>     <http://www.viruslist.com>
>
>     _______________________________________________
>     R-help-es mailing list
>     R-help-es en r-project.org <mailto:R-help-es en r-project.org>
>     https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>
>
> ---
> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running
at host imx2.etecsa.cu
> Visit our web-site: <http://www.kaspersky.com>,
<http://www.viruslist.com>
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130116/4b676929/attachment-0001.html>
------------ próxima parte ------------
An embedded and charset-unspecified text was scrubbed...
Name: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130116/4b676929/attachment-0001.pl>

Marcuzzi, Javier Rubén

2013-Jan-16 14:52 UTC

head link

[R-es] Interpretación del chisq.test

Estimado Leonardo

No es costumbre de esta lista responder dudas estadísticas, si de R o R
y algo de estadística pero no tan general.

Como usted esta realizando un curso y estudiando el test, sería bueno
que consulte a sus profesores.

Javier

El mié, 16-01-2013 a las 09:46 -0500, Leonardo Hernández Pérez
escribió:> Hola Javier: 
> 
> Lo que necesito es que me digan cual es la interpretación literal del
> resultados, por ejemplo para es resultado "X-squared = 3.9, df = 4,
> p-value = 0.4197" el valor de p me dice que no existen diferencias
> estadísticamente significativas pero en table(x) lo que tengo es la
> frecuencia de ocurrencia de 0 , 1, 2, 3 y 4 entonces supongo que entre
> estos valores (la frecuencia para cada uno de estos números) no
> existen diferencias estadísticamente significativas. 
> 
> Es así?
> 
> Un saludo
> 
> Leonardo
> 
> 
> On 16/01/13 08:15, Javier Marcuzzi wrote:
> 
> > Estimado Leonardo
> > 
> > 
> > Si copia y pega lo siguiente se dara cuenta
> > x <- trunc(5 * runif(100))
> > x
> > table(x)
> > chisq.test(table(x))
> > 
> > 
> > 
> > 
> > El 16 de enero de 2013 09:59, Leonardo Hernández Pérez
> > <leonardo.hernandez en etecsa.cu> escribió:
> >         
> >         Hola amigos:
> >         
> >         Estoy haciendo un curso y estamos estudiando el chisq.test
> >         este que les dejo debajo es un ejemplo que tiene rstudio.
> >         
> >         Pueden decirme como interpretar el resultado?
> >         
> >         Un saludo
> >         
> >         Leonardo
> >         
> >         x <- trunc(5 * runif(100))
> >         chisq.test(table(x))
> >         
> >         
> >         X-squared = 3.9, df = 4, p-value = 0.4197
> >         ---
> >         This message was processed by Kaspersky Mail Gateway
> >         5.6.28/RELEASE running at host imx3.etecsa.cu
> >         Visit our web-site: <http://www.kaspersky.com>,
> >         <http://www.viruslist.com>
> >         
> >         _______________________________________________
> >         R-help-es mailing list
> >         R-help-es en r-project.org
> >         https://stat.ethz.ch/mailman/listinfo/r-help-es
> >         
> > 
> > 
> > 
> > 
> > --- 
> > This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE
running at host imx2.etecsa.cu
> > Visit our web-site: <http://www.kaspersky.com>,
<http://www.viruslist.com>
> 
> --- 
> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running
at host imx3.etecsa.cu
> Visit our web-site: <http://www.kaspersky.com>,
<http://www.viruslist.com>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

Carlos Ortega

2013-Jan-16 15:09 UTC

head link

[R-es] Interpretación del chisq.test

Hola,

Mira como indica Javier en otras listas.
Por ejemplo aquí:

http://stats.stackexchange.com/search?q=chisq.test

Saludos,
Carlos Ortega
www.qualityexcellence.es


El 16 de enero de 2013 13:59, Leonardo Hernández Pérez <
leonardo.hernandez@etecsa.cu> escribió:
>
> Hola amigos:
>
> Estoy haciendo un curso y estamos estudiando el chisq.test este que les
> dejo debajo es un ejemplo que tiene rstudio.
>
> Pueden decirme como interpretar el resultado?
>
> Un saludo
>
> Leonardo
>
> x <- trunc(5 * runif(100))
> chisq.test(table(x))
>
>
> X-squared = 3.9, df = 4, p-value = 0.4197
> ---
> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE
> running at host imx3.etecsa.cu
> Visit our web-site: <http://www.kaspersky.com>,
<http://www.viruslist.com>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

Carlos Ortega

2013-Jan-16 15:12 UTC

head link

[R-es] Matriz funcion de pertenencia

Hola,

Además de poder leer alguno de los varios documentos gratuitos de cómo usar
R en sus funciones básicas (entre ellas "list") puedes mirar aquí
ejemplos
de  cómo enfocar tu problema:

http://stackoverflow.com/questions/9195112/reshape-matrix-into-a-list-of-lists

Saludos,
Carlos Ortega
www.qualityexcellence.es


El 16 de enero de 2013 13:49, Ivan Santos <lyaran2003@hotmail.com>
escribió:
> Es mi primera incursion en R, y no se muy bien como hacerlo, por eso
> preguntaba. ¿Como podria hacerlo con una lista?
>
> Muchas gracias
>
> Date: Wed, 16 Jan 2013 09:34:27 -0300
> Subject: Re: [R-es] Matriz funcion de pertenencia
> From: daniel319@gmail.com
> To: lyaran2003@hotmail.com
> CC: r-help-es@r-project.org
>
> Si entiendo bien, quieres asignar a una matriz objetos que no son vectores
> atómicos. Francamente, nunca lo intenté, pero ese puede ser el problema.
> ¿No te sirven las listas?
> Daniel Merino
>
>
> El 16 de enero de 2013 07:41, Ivan Santos <lyaran2003@hotmail.com>
> escribió:
>
> Buenos dias
>
>
>
> Estoy modelando un sistema de control difuso, y tengo un problema a la
> hora de obtener la matriz con las reglas difusas. Me explico.
>
>
>
> El sistema modela la velocidad de un ventilador, en funcion de la
> temperatura y la humedad
>
>
>
> 1. He definido las funciones de pertenencia de temperatura, humedad y
> velocidad
>
>
>
> #Función de pertenencia temperatura
>
>
>
> Tbaja<-vector(length=length(T), mode="numeric")
>
> Tbaja[T>40 & T<=60] <- 1
>
> Tbaja[T>60 & T<=80] <- (80-T[T>60 & T<=80])/(80-60)
>
> Tbaja[T>80] <- 0
>
>
>
> Tmedia<-vector(length=length(T), mode="numeric")
>
> Tmedia[T<60] <- 0
>
> Tmedia[T>60 & T<=80] <- (T[T>60 & T<=80]-60)/(80-60)
>
> Tmedia[T>80 & T<=100] <- (100-T[T>80 &
T<=100])/(100-80)
>
> Tmedia[T>100] <- 0
>
>
>
> Talta<-vector(length=length(T), mode="numeric")
>
> Talta[T<80] <- 0
>
> Talta[T>=80 & T<=100] <- (T[T>=80 &
T<=100]-80)/(100-80)
>
> Talta[T>100] <- 1
>
>
>
> #Funcion pertenencia humedad
>
>
>
> Hbaja<-vector(length=length(H), mode="numeric")
>
> Hbaja[H<25] <- 1
>
> Hbaja[H>=25 & H<=50] <- (50-H[H>=25 &
H<=50])/(50-25)
>
> Hbaja[H>50] <- 0
>
>
>
> Hmedia<-vector(length=length(H), mode="numeric")
>
> Hmedia[H<25] <- 0
>
> Hmedia[H>=25 & H<50] <- (H[H>=25 & H<50]-25)/(50-25)
>
> Hmedia[H>=50 & H<=75] <- (75-H[H>=50 &
H<=75])/(75-50)
>
> Hmedia[H>75] <- 0
>
>
>
> Halta<-vector(length=length(H), mode="numeric")
>
> Halta[H<50] <- 0
>
> Halta[H>=50 & H<=75] <- (H[H>=50 &
H<=75]-50)/(75-50)
>
> Halta[H>75] <- 1
>
>
>
> #Funcion pertenencia velocidad
>
>
>
> Wbaja<-vector(length=length(W), mode="numeric")
>
> Wbaja[W<250] <- 1
>
> Wbaja[W>250 & W<=500] <- (500-W[W>250 &
W<=500])/(500-250)
>
> Wbaja[W>500] <- 0
>
>
>
> Wmedia<-vector(length=length(W), mode="numeric")
>
> Wmedia[W<250] <- 0
>
> Wmedia[W>=250 & W<500] <- (W[W>=250 &
W<500]-250)/(500-250)
>
> Wmedia[W>=500 & W<=750] <- (750-W[W>=500 &
W<=750])/(750-500)
>
> Wmedia[W>750] <- 0
>
>
>
> Walta<-vector(length=length(W), mode="numeric")
>
> Walta[W<500] <- 0
>
> Walta[W>=500 & W<=750] <- (W[W>=500 &
W<=750]-500)/(750-500)
>
> Walta[W>750] <- 1
>
>
>
> 2. Una vez obtenidas las funciones de pertenencia, obtenemos las
> implicaciones.
>
> #Construimos las reglas de implicacion
>
> #velocidad baja, temperatura baja, velocidad baja
>
> if (Tbaja && Hbaja) Wbaja
>
>
>
> #velocidad media, temperatura baja, velocidad media
>
> if (Tmedia && Hbaja) Wmedia
>
>
>
> #velocidad alta, temperatura baja, velocidad baja
>
> if (Talta && Hbaja) Wmedia
>
>
>
> 3. Con las reglas deberia construir una matriz, para obtener las funciones
> de pertenencia de las reglas. Pero aqui me encuentro con el problema de
> como meter dichas reglas, en una matriz 3x3. He probado a asignar
> directamente las reglas a una posicion concreta pero me da error.
>
>
>
>
> 4. Como podria obtener la superficie de control.
>
>
>
> Muchas gracias
>
>
>
>         [[alternative HTML version deleted]]
>
>
>
>
> _______________________________________________
>
> R-help-es mailing list
>
> R-help-es@r-project.org
>
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>
>
> --
> Daniel
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

miguel.angel.rodriguez.muinos en sergas.es

2013-Jan-17 10:44 UTC

head link

[R-es] BioStatFLOSS y EpiLinux

Hola.

Por si es de interés...

Hemos publicado la versión 1.1 de BioStatFLOSS
http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=62658&idioma=es

En ella se incluye la versión 2.15.2 de R junto con RStudio (0.97.248),
RCommander (1.9-2), Deducer, Red-R, Tinn-R, ... y una serie de librerías para
bioestadística y epidemiología. Todo ello PORTABLE (no necesita instalación)
para Windows.

Para los que quieran algo homólogo en Linux (aquí con el Sistema Operativo
completo) os recomiendo que le echéis un ojo a EpiLinux 3
http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=50178&idioma=es


Un Saludo,
_____________________________
Miguel Ángel Rodríguez Muíños
Coordinador de los Proyectos EpiLinux y BioStatFLOSS
Dirección Xeral de Innovación e Xestión da Saúde Pública
Consellería de Sanidade
Xunta de Galicia

http://dxsp.sergas.es












Nota: A información contida nesta mensaxe e os seus posibles documentos adxuntos
é privada e confidencial e está dirixida únicamente ó seu destinatario/a. Se
vostede non é o/a destinatario/a orixinal desta mensaxe, por favor elimínea. A
distribución ou copia desta mensaxe non está autorizada.

Nota: La información contenida en este mensaje y sus posibles documentos
adjuntos es privada y confidencial y está dirigida únicamente a su
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje,
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.

See more languages: http://www.sergas.es/aviso_confidencialidad.htm

Griera

2013-Jan-17 13:58 UTC

head link

[R-es] BioStatFLOSS y EpiLinux

Muchas gracias por esta inciativa de un dvd-linux con estos programas.

Saludos.  

On Thu, 17 Jan 2013 11:44:10 +0100
<miguel.angel.rodriguez.muinos en sergas.es> wrote:
> Hola.
> 
> Por si es de interés...
> 
> Hemos publicado la versión 1.1 de BioStatFLOSS
>
http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=62658&idioma=es
> 
> En ella se incluye la versión 2.15.2 de R junto con RStudio (0.97.248),
RCommander (1.9-2), Deducer, Red-R, Tinn-R, ... y una serie de librerías para
bioestadística y epidemiología. Todo ello PORTABLE (no necesita instalación)
para Windows.
> 
> Para los que quieran algo homólogo en Linux (aquí con el Sistema Operativo
completo) os recomiendo que le echéis un ojo a EpiLinux 3
>
http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=50178&idioma=es
> 
> 
> Un Saludo,
> _____________________________
> Miguel Ángel Rodríguez Muíños
> Coordinador de los Proyectos EpiLinux y BioStatFLOSS
> Dirección Xeral de Innovación e Xestión da Saúde Pública
> Consellería de Sanidade
> Xunta de Galicia
> 
> http://dxsp.sergas.es
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> Nota: A información contida nesta mensaxe e os seus posibles documentos
adxuntos é privada e confidencial e está dirixida únicamente ó seu
destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe, por
favor elimínea. A distribución ou copia desta mensaxe non está autorizada.
> 
> Nota: La información contenida en este mensaje y sus posibles documentos
adjuntos es privada y confidencial y está dirigida únicamente a su
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje,
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.
> 
> See more languages: http://www.sergas.es/aviso_confidencialidad.htm
> 
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

miguel.angel.rodriguez.muinos en sergas.es

2013-Jan-18 08:14 UTC

head link

[R-es] BioStatFLOSS y EpiLinux

Considerando el no instalarla.... Es mucha mejor opción el USB....
:-)

Un Saludo,
Miguel.



-----Mensaje original-----
De: r-help-es-bounces en r-project.org [mailto:r-help-es-bounces en
r-project.org] En nombre de Griera
Enviado el: jueves, 17 de enero de 2013 14:58
Para: r-help-es en r-project.org
Asunto: Re: [R-es] BioStatFLOSS y EpiLinux


Muchas gracias por esta inciativa de un dvd-linux con estos programas.

Saludos.

On Thu, 17 Jan 2013 11:44:10 +0100
<miguel.angel.rodriguez.muinos en sergas.es> wrote:
> Hola.
>
> Por si es de interés...
>
> Hemos publicado la versión 1.1 de BioStatFLOSS
> http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=62658&idioma
> =es
>
> En ella se incluye la versión 2.15.2 de R junto con RStudio (0.97.248),
RCommander (1.9-2), Deducer, Red-R, Tinn-R, ... y una serie de librerías para
bioestadística y epidemiología. Todo ello PORTABLE (no necesita instalación)
para Windows.
>
> Para los que quieran algo homólogo en Linux (aquí con el Sistema
> Operativo completo) os recomiendo que le echéis un ojo a EpiLinux 3
> http://www.sergas.es/MostrarContidos_N3_T01.aspx?IdPaxina=50178&idioma
> =es
>
>
> Un Saludo,
> _____________________________
> Miguel Ángel Rodríguez Muíños
> Coordinador de los Proyectos EpiLinux y BioStatFLOSS Dirección Xeral
> de Innovación e Xestión da Saúde Pública Consellería de Sanidade Xunta
> de Galicia
>
> http://dxsp.sergas.es
>
>
>
>
>
>
>
>
>
>
>
>
> Nota: A información contida nesta mensaxe e os seus posibles documentos
adxuntos é privada e confidencial e está dirixida únicamente ó seu
destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe, por
favor elimínea. A distribución ou copia desta mensaxe non está autorizada.
>
> Nota: La información contenida en este mensaje y sus posibles documentos
adjuntos es privada y confidencial y está dirigida únicamente a su
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje,
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.
>
> See more languages: http://www.sergas.es/aviso_confidencialidad.htm
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Nota: A información contida nesta mensaxe e os seus posibles documentos adxuntos
é privada e confidencial e está dirixida únicamente ó seu destinatario/a. Se
vostede non é o/a destinatario/a orixinal desta mensaxe, por favor elimínea. A
distribución ou copia desta mensaxe non está autorizada.

Nota: La información contenida en este mensaje y sus posibles documentos
adjuntos es privada y confidencial y está dirigida únicamente a su
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje,
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.

See more languages: http://www.sergas.es/aviso_confidencialidad.htm

R help es - Jan 2013 - Matriz funcion de pertenencia

[R-es] Matriz funcion de pertenencia

[R-es] Matriz funcion de pertenencia

[R-es] Matriz funcion de pertenencia

[R-es] Interpretación del chisq.test

[R-es] Interpretación del chisq.test

[R-es] Interpretación del chisq.test

[R-es] Interpretación del chisq.test

[R-es] Interpretación del chisq.test

[R-es] Matriz funcion de pertenencia

[R-es] BioStatFLOSS y EpiLinux

[R-es] BioStatFLOSS y EpiLinux

[R-es] BioStatFLOSS y EpiLinux