Upload retail precision marketing sqlr demo

samples/features/r-services/Retail Precision Marketing/Data/README.md

@@ -0,0 +1,11 @@
+**List of data sets** 
+
+|  Data Set Name | Link to the Full Data Set   | Full Data Set Size (MB)  | Link to Report |
+| ---:| ---: | ---: | ---: |
+| CDNOW_master.csv | [link](https://github.com/ZhouFang928/sql-server-samples/blob/master/samples/features/r-services/Retail%20Precision%20Marketing/Data/CDNOW_master.csv) | 1.55MB | N/A|
+
+**Description of data sets**
+
+* The CDNOW data contains the entire purchase history up to the end of June 1998 of the cohort of 23,570 individuals who made their first-ever purchase at CDNOW in the first quarter of 1997. This CDNOW dataset was first used by Fader and Hardie (2001). Each record in this file, 69,659 in total, comprises four fields: the customer's ID, the date of the transaction, the number of CDs purchased, and the dollar value of the transaction.
+
+

samples/features/r-services/Retail Precision Marketing/R/CRMonCDNOWInMemory.R

@@ -0,0 +1,237 @@
+#############################################################
+# Title : CRM Demo in-memory
+# Author: Microsoft
+# Date: Dec, 2015
+#############################################################
+
+# Install package
+
+install.packages("rmarkdown")
+install.packages("fpc")
+
+# Set directory
+
+wd <- getwd()
+data.path <- file.path(wd, "Data", "CDNOW_master.csv")
+
+# Connect to SQL database using ODBC and read data from SQL via Open Source R
+
+library(RODBC)
+getSqlTypeInfo()
+
+# Connect from local PC
+
+channel <- odbcDriverConnect("driver={SQL Server Native Client 11.0};
+                              server=tcp:sqlserver2012-81yms1ai.cloudapp.net,57500;
+                              database=RREDemoSql;
+                              uid=******;
+                              pwd=******;")
+
+df <- sqlFetch(channel, 'CDNOW')
+df <- sqlQuery(channel, paste("select * from dbo.CDNOW"))
+
+df$Date<-as.Date(df$Date)
+
+str(df)
+head(df)
+
+# Remove the rows with the duplicated IDs to see how many customers in total
+
+uid <- df[!duplicated(df[,"ID"]), ]
+
+dim(uid)
+
+# Step 1: RFM analysis
+
+# Call RFM source code 
+
+source(wd, "R", "RFM_Analysis_R_Source_Codes_V1.3.R")
+
+# Set the startDate and endDate, we will only analysis the records in this date range
+
+startDate <- as.Date("19970101","%Y%m%d")
+endDate <- as.Date("19980701","%Y%m%d")
+
+# Calculate RFM value
+
+df <- getDataFrame(df, startDate, endDate, tIDColName="ID", tDateColName="Date", tAmountColName="Amount")
+
+head(df)
+
+# Obtain independent RFM score
+
+df1 <-getIndependentScore(df)
+
+head(df1)
+
+# Draw the histograms in the R, F, and M dimensions 
+
+drawHistograms(df1)
+
+S500 <- df1[df1$Total_Score > 500, ]
+dim(S500)
+
+S400 <- df1[df1$Total_Score > 400, ]
+dim(S400)
+
+# Obtain RFM score with breaks
+# Take a look at the distribution of R, F, M
+
+par(mfrow = c(1,3))
+hist(df$Recency)
+hist(df$Frequency)
+hist(df$Monetary)
+
+# Set the Recency ranges as 0-120 days, 120-240 days, 240-450 days, 450-500 days, and more than 500 days.
+
+r <- c(120, 240, 450, 500)
+
+# Set the Frequency ranges as 0-2 times, 2-5 times, 5-8 times, 8-10 times, and more than 10 times.
+
+f <- c(2, 5, 8, 10)
+
+# Set the Monetary ranges as 0-10 dollars, 10-20 dollars, and so on.
+
+m <-c(10,20,30,100)
+
+# Calculate RFM score with breaks
+
+df2 <- getScoreWithBreaks(df, r, f, m)
+drawHistograms(df2)
+
+S500 <- df2[df2$Total_Score > 500, ]
+dim(S500)
+
+S400 <- df2[df2$Total_Score > 400, ]
+dim(S400)
+
+target <- df2[df2$Total_Score >= 441,]
+dim(target)
+
+# Obtain RFM scores with quantiles as breaks
+
+r <-c(cutpoint(df$Recency))
+f <-c(cutpoint(df$Frequency))
+m <-c(cutpoint(df$Monetary))
+
+df3 <- getScoreWithBreaks(df, r, f, m)
+
+str(df3)
+head(df3)
+tail(df3)
+
+RFM_Result <- subset(df3,
+                     select=c("ID", "Recency", "Frequency", "Monetary",
+                              "R_Score", "F_Score", "M_Score", "Total_Score"))
+
+colnames(RFM_Result) <- c("ID", "R", "F", "M", "R_Score", "F_Score", "M_Score", "Total_Score")
+
+head(RFM_Result)
+
+time <- system.time({
+
+sqlSave(channel,
+        RFM_Result,
+        rownames=FALSE,
+        append=FALSE,
+        varTypes=list(numeric="float",
+                      integer="int"))
+})
+
+sqlUpdate(channel, df)
+odbcClose(channel)
+
+# Clustering using RFM
+
+library(fpc)
+library(cluster)
+
+# Kmeans clustering with number of cluster equal to 8
+
+cl.fit1 <- kmeans(RFM_Result[, 2:8],
+                  centers=8,
+                  iter.max=10,
+                  nstart=1)
+
+cl.fit2 <- kmeans(RFM_Result[, 2:8],
+                  centers=8,
+                  iter.max=20,
+                  nstart=200)
+
+summary(cl.fit1)
+
+cluster<-cl.fit1$cluster
+centers<-cl.fit1$centers
+size<-cl.fit1$size
+
+plot(RFM_Result[, 2:4], col=cl.fit1$cluster)
+title(main="K-means",line=3)	
+
+# Classification using RFM
+
+# Create IsVIP variable
+
+IsVIP <- ifelse(RFM_Result[,'Total_Score'] >= 441, 1, 0)
+Cluster <- cl.fit1$cluster
+RFMVIPCluster <- cbind(RFM_Result, IsVIP, Cluster)
+
+# Create training/testing data set
+
+RD <- sample(1:10, dim(RFMVIPCluster)[1], replace=TRUE)
+
+str(RD)
+table(RD)
+
+RFMVIPCluster$RD <- RD;
+
+urv <- factor(ifelse(RD <= 8,'TRAIN','TEST'))
+TrainTest <- cbind(RFMVIPCluster, urv)
+Train <- TrainTest[which(TrainTest$urv == "TRAIN"), ]
+Test <- TrainTest[which(TrainTest$urv == "TEST"), ]
+
+# Logistic model
+# Build our Logistic Regression Model with IsVIP as response
+
+r1 <- glm(IsVIP~R+F+M, data=Train, family = binomial)
+summary(r1)
+
+p1 <- predict.glm(r1, data=Test, type="response")
+head(p1)
+tail(p1)
+
+# Decision tree
+# Grow tree 
+
+fit <- rpart(Cluster~R+F+M,
+           	 method="class", 
+             data=Train)
+
+# Display the results 
+
+printcp(fit)
+
+# Visualize cross-validation results 
+
+plotcp(fit) 
+
+# Detailed summary of splits
+
+summary(fit) 
+
+# Plot tree 
+
+library(rpart)
+
+plot(fit, uniform=TRUE, main="Classification Tree for CDNOW")
+text(fit, use.n=TRUE, all=TRUE, cex=.8)
+
+# Prune the tree 
+
+pfit <- prune(fit, cp=fit$cptable[which.min(fit$cptable[,"xerror"]), "CP"])
+
+# Plot the pruned tree 
+
+plot(pfit, uniform=TRUE, 
+  	 main="Pruned Classification Tree for CDNOW")
+text(pfit, use.n=TRUE, all=TRUE, cex=.8)
+