///加载自己写的R语言算法库
 public List<double> GetZTFB(double[] data)
        {
            List<double> par = new List<double>();
            try
            {
                //调用R语言算法
                REngine.SetEnvironmentVariables();
                REngine engine = REngine.GetInstance(null, true, null, null);
                NumericVector x = engine.CreateNumericVector(data);
                engine.SetSymbol("x", x);
                string path = System.Windows.Forms.Application.StartupPath + "\\R_File\\ztfb.R";
                path = "source(\"" + path + "\")";
                path = path.Replace("\\", "\\\\");//这一步至关重要，否则会报错
                engine.Evaluate(path);
                engine.Evaluate("a<- ztfb(x)");

                //                engine.Evaluate(@"normal <- function(theta){
                //                                              mu <- theta[1]
                //                                              sigma <- theta[2]
                //                                              n <- length(x)
                //                                              logL <- -0.5*n*log(2*pi) - n*log(sigma) - sum(0.5*(x - mu)^2/sigma^2)
                //                                              return (-logL)
                //                                            } 
                //                a<- nlminb(c(0,1),normal)");
                NumericVector a = engine.Evaluate("a$par[1]").AsNumeric();
                NumericVector b = engine.Evaluate("a$par[2]").AsNumeric();
                par.Add(Convert.ToDouble(a[0]));
                par.Add(Convert.ToDouble(b[0]));
            }
            catch (Exception ex)
            {
                throw new Exception(ex.Message);
            }
            return par;
        }

///加载R语言下载库
 public List<double> GetZTFB(double[] data)
        {
            List<double> par = new List<double>();
            try
            {
                //调用R语言算法
                REngine.SetEnvironmentVariables();
                REngine engine = REngine.GetInstance(null, true, null, null);
				///加载R语言库
				engine.Evaluate("library(\"bit\", lib.loc=\"C:/Program Files/R/R-3.3.2/library\")")
                NumericVector x = engine.CreateNumericVector(data);
                engine.SetSymbol("x", x);
                string path = System.Windows.Forms.Application.StartupPath + "\\R_File\\ztfb.R";
                path = "source(\"" + path + "\")";
                path = path.Replace("\\", "\\\\");//这一步至关重要，否则会报错
                engine.Evaluate(path);
                engine.Evaluate("a<- ztfb(x)");

                //                engine.Evaluate(@"normal <- function(theta){
                //                                              mu <- theta[1]
                //                                              sigma <- theta[2]
                //                                              n <- length(x)
                //                                              logL <- -0.5*n*log(2*pi) - n*log(sigma) - sum(0.5*(x - mu)^2/sigma^2)
                //                                              return (-logL)
                //                                            } 
                //                a<- nlminb(c(0,1),normal)");
                NumericVector a = engine.Evaluate("a$par[1]").AsNumeric();
                NumericVector b = engine.Evaluate("a$par[2]").AsNumeric();
                par.Add(Convert.ToDouble(a[0]));
                par.Add(Convert.ToDouble(b[0]));
            }
            catch (Exception ex)
            {
                throw new Exception(ex.Message);
            }
            return par;
        }