[Ruler-commits] r26 - pkg
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Tue Jul 31 14:54:10 CEST 2012
Author: merysionek
Date: 2012-07-31 14:54:10 +0200 (Tue, 31 Jul 2012)
New Revision: 26
Added:
pkg/ruleR.R
Modified:
pkg/New.R
Log:
.
Modified: pkg/New.R
===================================================================
--- pkg/New.R 2012-07-31 10:06:58 UTC (rev 25)
+++ pkg/New.R 2012-07-31 12:54:10 UTC (rev 26)
@@ -163,7 +163,7 @@
setClass("DivDoubleRule",contains="DoubleRule",S3methods=TRUE)
-setMethod("calculateDoubleSpecific",
+setMethod("calculateSpecific",
signature(x="DivDoubleRule", y="numeric", z="numeric"),
function(x,y,z){
return(y%/%z)
@@ -173,7 +173,7 @@
setClass("ModuloDoubleRule",contains="DoubleRule",S3methods=TRUE)
-setMethod("calculateDoubleSpecific",
+setMethod("calculateSpecific",
signature(x="ModuloDoubleRule", y="numeric", z="numeric"),
function(x,y,z){
return(y%%z)
@@ -182,7 +182,7 @@
#[6] EXPONENTIAL FUNCTION (Philipp)
setClass("ExpDoubleRule", contains="DoubleRule",S3methods=TRUE)
-setMethod("calculateDoubleSpecific",
+setMethod("calculateSpecific",
signature(x="ExpDoubleRule", y="numeric", z="numeric"),
function(x,y,z){
return(y^z)
@@ -287,7 +287,7 @@
#a list of double rules
-doubleRules<-list("AddDoubleRule","MultDoubleRule","SubsDoubleRule")
+doubleRules<-list("AddDoubleRule","MultDoubleRule","SubsDoubleRule","DivDoubleRule","ModuloDoubleRule","ExpDoubleRule")
Added: pkg/ruleR.R
===================================================================
--- pkg/ruleR.R (rev 0)
+++ pkg/ruleR.R 2012-07-31 12:54:10 UTC (rev 26)
@@ -0,0 +1,468 @@
+#-------------------------------------------------------------------------------------------
+#------------------------------SingleRules--------------------------------------------------
+#-------------------------------------------------------------------------------------------
+
+#VIRTUAL CLASS FOR RULES OPERATING ON SINGLE ARGUMENTS
+
+setClass("SingleRule",
+ representation = representation(previousRule="SingleRule"),
+ S3methods=TRUE)
+
+
+calculateSpecific <- function(x,y,z=NULL){
+ return(y)
+}
+
+
+setMethod("calculateSpecific",signature(x="SingleRule", y="numeric"),
+ function(x,y){
+ return(y)
+ })
+
+
+
+#[1] RULE 1 - ADDING A CONSTANT
+
+setClass("AddConstSingleRule",
+ contains="SingleRule",
+ representation(constantVal="numeric"),
+ S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="AddConstSingleRule", y="numeric"),
+ function(x,y){
+ return(x at constantVal+y)
+ })
+
+
+#[2] RULE 2 - MULTIPLYING BY A CONSTANT
+
+setClass("MultConstSingleRule",
+ contains="SingleRule",
+ representation(constantVal="numeric"),
+ S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="MultConstSingleRule", y="numeric"),
+ function(x,y){
+ return(x at constantVal*y)
+ })
+
+#[3] SUBSTRACTING A CONSTANT
+
+setClass("SubsConstSingleRule",
+ contains="SingleRule",
+ representation(constantVal="numeric"),
+ S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="SubsConstSingleRule",y="numeric"),
+ function(x,y){
+ return(y-x at constantVal)
+ })
+
+#[4] DIGITSUM
+digits <- function(x) {
+ if(length(x) > 1 ) {
+ lapply(x, digits)
+ } else {
+ n <- nchar(x)
+ rev( x %/% 10^seq(0, length.out=n) %% 10 )
+ }
+}
+
+
+
+setClass("DigSumSingleRule", contains="SingleRule",S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="DigSumSingleRule",y="numeric"),
+ function(x,y){
+ return(sum(digits(y)))
+ })
+
+#[5] NEGATIVE
+
+setClass("NegativeSingleRule", contains="SingleRule",S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="NegativeSingleRule",y="numeric"),
+ function(x,y){
+ return(-y)
+ })
+
+
+#[4] IDENTICAL FUNCTION (input=output) used in random sequence generation
+setClass("IdenSingleRule",contains="SingleRule",S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="IdenSingleRule",y="numeric"),
+ function(x,y){
+ return(y)
+ })
+
+
+
+
+#EXECUTING RULES REFERING TO SINGLE ARGUMENT
+
+
+calculate <- function(x,y,z=NULL){
+ return(y)
+}
+
+
+setMethod("calculate",signature(x="SingleRule", y="numeric"), #both [1] and [2] inherit from class 'SingleRule'
+ function(x, y){
+ result<-y
+ if(!is.null(x at previousRule)){ # if there are some rules nested inside 'x'
+ result <- calculate(x at previousRule,result)
+ }
+ return(calculateSpecific(x,result)) # if there are no more nested functions, execute
+ })
+
+#-------------------------------------------------------------------------------------------
+#------------------------------------------------DoubleRules--------------------------------
+#-------------------------------------------------------------------------------------------
+#VIRTUAL CLASS FOR RULES OPERATING ON TWO ARGUMENTS
+
+# firstRule - operations to be executed at the first element of the numeric sequence (it is an object of class 'SingleRule')
+# secondRule - operations to be executed at the second element of the numeric sequence (it is an object of class 'SingleRule')
+# nextSingle - operation to be executed at the result of function DoubleRule
+
+
+setClass("DoubleRule", representation = representation(firstRule="SingleRule", secondRule="SingleRule",nextSingle="SingleRule"),
+ S3methods=TRUE)
+
+
+
+#[1] ADD TWO PREVIOUS EXPRESSIONS
+
+setClass("AddDoubleRule", contains="DoubleRule",S3methods=TRUE)
+
+setMethod("calculateSpecific",signature(x="AddDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y+z)
+ })
+
+#[2] MULTIPLY TWO PREVIOUS EXPRESSIONS
+
+setClass("MultDoubleRule",contains="DoubleRule",S3methods=TRUE)
+
+
+setMethod("calculateSpecific",signature(x="MultDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y*z)
+ })
+
+#[3] SUBSTRACT TWO PREVIOUS EXPRESSIONS
+
+setClass("SubsDoubleRule",contains="DoubleRule",S3methods=TRUE)
+
+
+setMethod("calculateSpecific",signature(x="SubsDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y-z)
+ })
+
+#[4] DIVIDING TWO NUMBERS (Philipp)
+
+setClass("DivDoubleRule",contains="DoubleRule",S3methods=TRUE)
+
+setMethod("calculateDoubleSpecific",
+ signature(x="DivDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y%/%z)
+ })
+
+#[5] MODULO (Philipp)
+
+setClass("ModuloDoubleRule",contains="DoubleRule",S3methods=TRUE)
+
+setMethod("calculateDoubleSpecific",
+ signature(x="ModuloDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y%%z)
+ })
+
+#[6] EXPONENTIAL FUNCTION (Philipp)
+setClass("ExpDoubleRule", contains="DoubleRule",S3methods=TRUE)
+
+setMethod("calculateDoubleSpecific",
+ signature(x="ExpDoubleRule", y="numeric", z="numeric"),
+ function(x,y,z){
+ return(y^z)
+ })
+
+
+
+#EXECUTING RULES OPERATING ON TWO ARGUMENTS
+
+#calculate <- function(x,y,z){stop ("No function to execute this.")} #throw a mistake, because you should execute just single functions contained by cladd DoubleRule
+
+
+setMethod("calculate",signature(x="DoubleRule", y="numeric", z="numeric"),
+ function(x, y, z){
+ firstArg <- y #first element of the sequence
+ secondArg <-z #second element of the sequence
+
+ if(!is.null(x at firstRule)){ #if there are some rules nested inside
+ firstArg <- calculate(x at firstRule,firstArg) #execute first single-argument rule
+ }
+
+ if(!is.null(x at secondRule)){
+ secondArg <- calculate(x at secondRule,secondArg) #execute second single-argument rule
+ }
+
+ result<-calculateSpecific(x,firstArg, secondArg) #if there are no more nested rules, execute
+
+ if(!is.null(x at nextSingle)){
+ result<-calculate(x at nextSingle, result)
+ }
+ return(result)
+
+ })
+
+
+
+#-------------------------------------------------------------------------------------------
+#----------------------------------generating sequences-------------------------------------
+#-------------------------------------------------------------------------------------------
+
+#a list of single rules
+singleRules<-list("IdenSingleRule","AddConstSingleRule","MultConstSingleRule","SubsConstSingleRule","DigSumSingleRule","NegativeSingleRule")
+
+
+#a list of double rules
+doubleRules<-list("AddDoubleRule","MultDoubleRule","SubsDoubleRule","DivDoubleRule","ModuloDoubleRule","ExpDoubleRule")
+
+
+
+#A FUNCTION TO CREATE SINGLE RULES
+# 'a1' is an index from table SingleRule (default a=NULL) //if 'a' is NULL it will be generated
+# 'n' how many rules are to be nested (default=NULL menas that I want to generate it automatically let's say from c(0,1,2)
+# n=0 would mean that I want to create just one rule with nothing nested inside)
+# 'cv1' is a constant value
+# '...' if I would like to add some rules nested I can provide their parameters cv must be always supplied #9even if the function doesn't require that
+
+createSR<-function(a1=NULL,cv1=NULL,n=NULL,...){
+ p<-list(...)#arguments for nesting other functions
+ p<-unlist(p)
+
+ #if(!is.null(n) && length(p)!=2*n) stop (paste("parameters of functions to be nested do not match n=",n))
+ #for(i in seq(1,length(p),by=2)){if(k[i]>length(p)) stop (paste("List of rules is shorter than ",k[i]))}
+
+ if(is.null(a1)) {a1<-sample(2:length(singleRules),1)} #generate 'a' if no is supplied (we don't want to generate a=1 because it is identical function)
+ if(is.null(cv1)) {cv1<-sample(1:100,1)} # generate a constant value if no is supplied
+ if(is.null(n)){n<-sample(c(0,1,2),1,prob=c(3/6,2/6,1/6)) #nesting more than two rules would be impossible to guess
+ p<-as.vector(matrix(replicate(n,c(sample(1:length(singleRules),1),sample(1:100,1))),1,2*n))
+ } # generate 'n' if it is set as null with different probabilities
+
+
+ if("constantVal"%in%slotNames(singleRules[[a1]])){m<-new(singleRules[[a1]],constantVal=cv1)} else{m<-new(singleRules[[a1]])}
+
+ if(n!=0) {createSR(p[[1]],p[[2]],n-1,p[-c(1,2)]); m at previousRule<-k
+ }else{return(m)}
+
+ return(m)
+ }
+
+
+
+# A FUNCTION TO COMBINE DOUBLE RULES - it generates all parameters automatically
+
+# 'a' is index from a list of DoubleRules
+#'fr' firstRule argument of an object of class doubleRule
+#'sr' secondRule argument of an object of class doubleRule
+#'ns' nextSingle argument of an object of class doubleRule
+#'
+createDR<-function(a=NULL,fr=NULL,sr=NULL,ns=NULL){
+ if(!is.null(a) && a>length(doubleRules)) stop (paste("The list of doublrRules is shoreter than ",a, ".Please specify 'a' value, which is smaller than or equal to",length(doubleRules)))
+ if(!inherits(fr,"singleRule") && !is.null(fr))stop(paste("'fr' argument must inherit from class singleRule"))
+ if(!inherits(sr,"singleRule") && !is.null(sr))stop(paste("'sr' argument must inherit from class singleRule"))
+ if(!inherits(ns,"singleRule") && !is.null(ns))stop(paste("'ns' argument must inherit from class singleRule"))
+
+ if(is.null(a)) a<-sample(1:length(doubleRules),1) #generate an index of a doubleRule from the list of doubleRules
+ a<-doubleRules[[a]]
+ #print(a)
+
+ if(is.null(fr)) fr<-sample(c(k=createSR(),k=new("IdenSingleRule")),1,prob=c(0.5,0.5))# firstRule is chosen from an automatically generated SingleRule or identical rule returning the input
+
+ if(is.null(sr)) sr<-sample(c(k=createSR(),k=new("IdenSingleRule")),1,prob=c(0.3,0.7)) #because adding more and more rules makes the rule very difficult I would generate identical function with greater probability
+
+ if(is.null(ns)) ns<-sample(c(k=createSR(),k=new("IdenSingleRule")),1,prob=c(0.3,0.7))
+
+ p<-new(a,firstRule=fr$k, secondRule=sr$k,nextSingle=ns$k)
+ return(p)
+
+ }
+
+
+#A FUNCTION TO GENERATE NUMERIC SEQUENCE OF DECLARED LENGTH
+# 'n' is the length of the numeric sequence (default value is 6)
+# 'x1', 'x2' are the first elements of the numeric sequence (you don't always need 'x2')
+sequence<-function(x1,x2=NULL,rule,n=6){
+ if(inherits(rule,"DoubleRule") && is.null(x2)) stop (" If you want to use a DoubleRoule you need to specify x2")
+ if(class(x1)!="numeric" ||(class(x2)!="NULL" && class(x2)!="numeric")) stop ("arguments 'x1', 'x2' must be of type 'numeric'.")
+ if(!inherits(rule,"SingleRule") && !inherits(rule,"DoubleRule")) stop ("'rule' argument must inherit from 'SingleRule' or 'DoubleRule' class")
+ if(n<3) stop("sequence must be longer than 3")
+
+ k<-list()
+ k[1]=x1
+
+
+
+ if(inherits(rule,"SingleRule")){
+ for(i in 2:n){
+ k[i]<-calculate(x=rule,y=k[[i-1]])
+ }
+
+ }else{
+ k[2]=x2
+ for(i in 3:n){
+ k[i]<-calculate(x=rule,y=k[[i-2]],z=k[[i-1]])
+ }
+
+ }
+ return(list(k,rule))
+
+ }
+
+
+
+
+
+# checking if a vector is in any row of the matrix
+#'mx' is a matrix in which I am searching
+#'vec' is a vector which is being checked
+# result TRUE means that there is already such vector in the matrix
+
+duplicate<-function(mx,vec){
+ return(any(apply(mx, 1, function(x, want) isTRUE(all.equal(x, want)), vec)))
+ }
+
+
+
+
+
+#CHECKING IF THE SEQUENCE IS NOT CONSTANT
+# it returns '0' when teh sequence is constant and '1' when the sequence is not constant
+# a function examines three last elements of a sequence, so even sequences like 27,9,9,9,9 ... are excluded
+
+conCheck<-function(seq){
+
+ if(class(seq)!="list") stop("sequence must be of type 'list'")
+
+ m<-length(seq)
+
+ if(identical(seq[m],seq[m-1]) && identical(seq[m],seq[m-2]) ) {return(0)} else {return(1)}
+ }
+
+
+
+# checking whether the sequence is not constant, numbers are not greater than 1000 or no lesser than -1000
+# type=1 is totally automatic / type=2 generates rules from 'MyRules' list
+check<-function(seqlen,items,type){
+ x1<-as.numeric(sample(1:100,1)) #generate the first element of numeric sequence
+ x2<-as.numeric(sample(1:100,1)) # generate the second element of numeric sequence
+
+ if(type==1){ # type=1 means automatic tests
+
+ m<-sample(c(1,2),1) #if m=1 I will create a singleRule, if m=2 rule will be a combination of singleRules, if m=3 rule is a doubleRule
+ if(m==1){rule<-createSR()} else{rule<-createDR()} } else
+ {z<-sample(2:length(MyRules),1); rule<-MyRules[[z]] }# if m=1 create singleRule else create doubleRulr
+
+ result<-sequence(x1,x2,rule,n=seqlen)[[1]]
+
+ fun<-sequence(x1,x2,rule,n=seqlen)[[2]]
+
+
+ if(conCheck(result)==0 || result[length(result)]>1000 || result[length(result)]< -1000||duplicate(mx=items,vec=unlist(result[[1]])))
+ {check(seqlen,items,type)} else{return(list(result,fun))}
+
+ }
+
+
+
+
+
+
+# AUTOMATIC TEST GENERATION
+# random
+# 'seqlen' specyfies how long should a single sequence be
+# 'type' = 1 everything generated automatically/ type=2 functions generated from the MyRules matrix
+# 'testlen' specyfies how many sequences (item positions) are there to be in a test
+automaticTest<-function(testlen,type=1,seqlen=6){
+ if(class(type)!="numeric") stop ("argument 'type' must be of class 'numeric'")
+ items<-matrix(NA,testlen,seqlen) #I will store generated items in a matrix
+ rules<- list() # I will keep the rules on a list
+ noise<-matrix(NA,testlen,5) # 5 noise answers will be stored in this matrix
+
+ for(i in 1:testlen){
+ b<-check(seqlen,items,type)
+ items[i,]<-unlist(b[[1]]) # I treat last element as a correct answer
+ rules[i]<-b[2]
+ noise[i,] <-sample(c((items[i,seqlen]-6):(items[i,seqlen]+6))[-(seqlen+1)],5) # generate a sample of 5 elements from a set (correct answer-6: correct answer +6) excluding the correct answer
+ }
+
+ return(list(items=items,noise=noise,rules=rules))
+ }
+
+
+
+
+#--------------------------------------------------------------------------------------------------------------------------------
+#---------------------------------------- rules defined by a user----------------------------------------------------------------
+#--------------------------------------------------------------------------------------------------------------------------------
+
+#LIST OF USER'S RULES
+
+
+# Verification table - in which I will store results of the rules
+# starting values are always 10 and 20, constantVal=14, seqlen=6
+# I check whether sequence generated for those values are identical for a new rule - if so - I won't add a new rule to the list
+
+VerifTable<-matrix(NA,1,6)
+
+MyRules<-list(NA)
+
+AddRule<-function(rule){
+ if(!inherits(rule,"SingleRule") && !inherits(rule,"DoubleRule")) stop (paste("new Rule must be an object of class 'SingleRule' od 'DoubleRule'. It cannot be of class",class(rule)))
+
+ s<-sequence(x1=10,x2=20,rule,n=6)
+
+ if(duplicate(VerifTable,unlist(s[[1]]))==TRUE){stop ("There already is a rule that gives the same result.")} else{
+ VerifTable<<-rbind(unlist(s[[1]]),VerifTable)
+ MyRules[length(MyRules)+1]<<-rule
+ print("Your rule succesfully added tp 'MyRule' list")
+ }
+ }
+
+
+#---------------------------------------------------------------------------
+#------function to print objects of class SingleRule and DoubleRule---------
+#---------------------------------------------------------------------------
+
+print.SingleRule<-function(x){
+ pr<-function(x){
+ cat(paste("\nname:", class(x)[1]))
+ if("constantVal"%in%slotNames(x)) {cat(paste(", constant value: ", x at constantVal))}
+ }
+ pr(x)
+
+ if(!class(x at previousRule)=="SingleRule"){x<-x at previousRule; pr(x)}
+
+ }
+
+
+
+print.DoubleRule<-function(x){
+ x1=x
+
+
+ cat(class(x1)[1])
+
+ if("firstRule"%in%slotNames(x)){cat("\n FIRST RULE:"); x<-x at firstRule;print.SingleRule(x)}
+ x=x1
+ if("secondRule"%in%slotNames(x)){cat("\n SECOND RULE:"); x<-x at secondRule;print.SingleRule(x)}
+ x=x1
+ if("nextSingle"%in%slotNames(x)){cat("\n NEXT SINGLE:"); x<-x at nextSingle;print.SingleRule(x)}
+
+ }
+
+
+
More information about the Ruler-commits
mailing list