Enhancing SAP PA Summary for Custom R

customsapparenhance

	#Create your own S class, pass your results as object, use your own summary function
	Summary_Printer <- function(obj1) UseMethod("Summary_Printer")
	Summary_Printer.default <- function(regression,sandwich)
	{
	  printer <- NULL
	  printer$regression <- regression
	  printer$sandwich <- sandwich
	  class(printer) <- "Summary_Printer"
	  printer
	}

	#Your own summary function of the Summary_Printer Class
	#This will be output as summary in PA
	summary.Summary_Printer <- function(printer)
	{
	  summary(printer$regression)
	  print("**** SANDWICH ESTIMATOR ***\n")
	  print(printer$sandwich)
	}

	REG_SANDWICH <- function(data,x,y)
	{
	  formattedString<-paste(x,collapse='+')
	  formatedString <- as.formula(paste(paste(y," ~ ",collapse=""),formattedString,collapse=""))
	  results <- lm(formatedString,data)
	  df <- data.frame(results$coefficients)
	  df <- cbind(row.names(df),df)
	  
	  #Visualisation
	  par(mfrow=c(2,2))
	  plot(results)
	  plot(results)
	  plot(results)
	  plot(results)
	  
	  sandwich_test <- bptest(results,studentize=FALSE)
		
	  results1 <- Summary_Printer(results)

	  return(list(out=df, model=results1))
	}