<?xml version="1.0" encoding="utf-8" ?><rss version="2.0"><channel><title><![CDATA[vicky428的博客]]></title><description><![CDATA[]]></description><link>https://blog.youkuaiyun.com/vicky428</link><language>zh-cn</language><generator>https://blog.youkuaiyun.com/</generator><copyright><![CDATA[Copyright &copy; vicky428]]></copyright><item><title><![CDATA[assignment1.R-20170914]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100237479</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100237479</guid><author>vicky428</author><pubDate>Tue, 03 Sep 2019 12:43:21 +0800</pubDate><description><![CDATA[#数据预处理
data=data.frame(XY_train)
colnames(data)=c('y',paste("x",1:34,sep=""))
#离散因素数据处理
data[,'y']=factor(data[,'y'])
data[,'x33']=factor(data[,'x33'])
data[,'x34']=factor(data[,'x34'])
#转换因子水平为数值
le...]]></description><category></category></item><item><title><![CDATA[logit2.R-20170913]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100236798</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100236798</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:46:42 +0800</pubDate><description><![CDATA[library(caret) 
library(MASS)
library(glmnet)
data=as.matrix(oadata2)
colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))

c=10 #验证c次
set.seed(5000)
r3=matrix(0,c,1)
a=matrix(0,c,2)
merror=m...]]></description><category></category></item><item><title><![CDATA[nnls.R-20170908]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100236313</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100236313</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:45:01 +0800</pubDate><description><![CDATA[cdata=condata2
cdata=as.matrix(cdata)
colnames(cdata)=c('y','x1','x2','x7^2','x7','x17',paste("x3",0:1,sep="="),paste("x4",0:6,sep="="),
                     paste("x5","=",0:1,sep=""),paste("x6","="...]]></description><category></category></item><item><title><![CDATA[selectvar.R-20170907]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100236054</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100236054</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:44:29 +0800</pubDate><description><![CDATA[# datap=cbind(condata[,1],pc)
# colnames(datap)=c('y',paste("x",1:9,sep=""))
# l=lm(y~.,data=data.frame(datap))
# summary(l)

condata2=condata3
condata2=as.matrix(condata2)
colnames(condata2)=c('y','...]]></description><category></category></item><item><title><![CDATA[lm.R-20170905]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100235815</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100235815</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:43:52 +0800</pubDate><description><![CDATA[#colnames(condata)=c('y','x1','x2','x17',paste("x",3:16,sep=""))
condata=as.matrix(condata)
colnames(condata)=c('y','x1','x2','x17',paste("x3",0:1,sep="="),paste("x4",0:6,sep="="),
                  ...]]></description><category></category></item><item><title><![CDATA[dividedata2.R-20170905]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100235477</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100235477</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:43:17 +0800</pubDate><description><![CDATA[library(caret) 
library("kernlab")
data=as.matrix(oadata)
colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))

c=1000 #验证c次
set.seed(1000)
r1=matrix(0,c,1)
r2=matrix(0,c,1)
r4=matrix(0,c,1)
...]]></description><category></category></item><item><title><![CDATA[fit.R-20170905]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100235278</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100235278</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:42:29 +0800</pubDate><description><![CDATA[yinterval=as.matrix(yinterval)
f=as.matrix(yinterval[,3:18])
xint=seq(5000,155000,10000)
par(mfrow=c(2,2))
plot(xint,f[4,])

m=matrix(0,nrow(f),1)
v=matrix(0,nrow(f),1)
xmu=matrix(0,nrow(f),1)
xsigma...]]></description><category></category></item><item><title><![CDATA[bootstrap.R-20170901]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100234998</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100234998</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:41:47 +0800</pubDate><description><![CDATA[# #bootstrap验证
# #返回回归系数向量的自编函数
# bs &lt;- function(formula, data, indices){
#  d &lt;- data[indices,]
#  fit &lt;- lm(formula,data=d)
#  return(coef(fit))
#  }
# #自助抽样
# set.seed(1234)
# a = lm(y~0+...]]></description><category></category></item><item><title><![CDATA[logit.R-20170901]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100234782</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100234782</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:41:08 +0800</pubDate><description><![CDATA[library(caret) 
data=as.matrix(oadata2)
colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))

c=5000 #验证c次
set.seed(5000)
r1=matrix(0,c,1)
r2=matrix(0,c,1)
p1=matrix(0,12,c)

for (k in 1:c) {...]]></description><category></category></item><item><title><![CDATA[dividedata.R-20170830]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100234516</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100234516</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:40:27 +0800</pubDate><description><![CDATA[#data=cbind(ypart,xcovpart2)
#sample(2, nrow(data), replace=T, prob=c(0.7,0.3))
library(caret) 

data=as.matrix(oadata)
colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))

c=5000 #验证c次
set....]]></description><category></category></item><item><title><![CDATA[covdesign2.R-20170829]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100234169</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100234169</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:39:18 +0800</pubDate><description><![CDATA[data=as.matrix(oadata)
colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))

# library(DoE.base)
# oa.design(ID=L9.3.4)
# oa.design(nfactors=6, nlevels=3)

#leave-one-out验证
library(caret) 
#r...]]></description><category></category></item><item><title><![CDATA[covdesign-20170826]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100233793</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100233793</guid><author>vicky428</author><pubDate>Mon, 02 Sep 2019 12:38:03 +0800</pubDate><description><![CDATA[x1&lt;-xcovpart[,1] 
x2&lt;-xcovpart[,2] 
x3&lt;-xcovpart[,3] 
x4&lt;-xcovpart[,4] 
x5&lt;-xcovpart[,5] 
x6&lt;-xcovpart[,6] 
x7&lt;-xcovpart[,7] 
x8&lt;-xcovpart[,8] 
x9&lt;-xcovpart[,9] 
x10&lt;-xc...]]></description><category></category></item><item><title><![CDATA[svm.R-20170824]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100173140</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100173140</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:10:28 +0800</pubDate><description><![CDATA[library(e1071)
m1&lt;-svm(xcovpart2,ypart)
summary(m1)
ypart.pre3&lt;-predict(m1,xcovpart2)
resi3=abs(ypart-ypart.pre3)/ypart
mean(resi3) #误差
yall.pre3=predict(m1,newdata=data.frame(xcovall2))

m2&lt...]]></description><category></category></item><item><title><![CDATA[savedata.R-20170704]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/100173106</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/100173106</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:09:30 +0800</pubDate><description><![CDATA[write.table(des1, file = "/Users/vicky/Documents/code/R/des1.txt")
write.table(des2, file = "/Users/vicky/Documents/code/R/des2.txt")
write.table(des3, file = "/Users/vicky/Documents/code/R/des3.txt...]]></description><category></category></item><item><title><![CDATA[hello.py-20190220]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95882879</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95882879</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:06:03 +0800</pubDate><description><![CDATA[#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Feb 19 08:40:55 2019

@author: vicky
"""

from flask import Flask
app = Flask(__name__)

@app.route('/')
def hello_world():
    retur...]]></description><category></category></item><item><title><![CDATA[app.py-20190229]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95882988</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95882988</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:05:58 +0800</pubDate><description><![CDATA[from flask import Flask
from flask import request

app = Flask(__name__)

#GET /：首页，返回Home；
@app.route('/', methods=['GET', 'POST'])
def home():
    return '&lt;h1&gt;Home&lt;/h1&gt;'

#GET /signin：登...]]></description><category></category></item><item><title><![CDATA[hello_jinjia.py hello_celery.py -20190219]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95883059</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95883059</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:05:52 +0800</pubDate><description><![CDATA[from jinja2 import Template
template = Template('Hello {{ name }}!')
template.render(name='John Doe')




#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Feb 19 11:38:12 2019

@aut...]]></description><category></category></item><item><title><![CDATA[jm3.py-20180919]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95883261</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95883261</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:05:47 +0800</pubDate><description><![CDATA[#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Sep 17 03:01:06 2018

@author: vicky
"""
import pandas as pd 
import numpy as np
import matplotlib.pyplot as plt 
import seaborn as s...]]></description><category></category></item><item><title><![CDATA[jm4.py-20180918]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95883717</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95883717</guid><author>vicky428</author><pubDate>Sat, 31 Aug 2019 15:05:39 +0800</pubDate><description><![CDATA[#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Sep 17 12:46:07 2018

@author: vicky
"""

import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.patches as mpatches
i...]]></description><category></category></item><item><title><![CDATA[jm2.py-20180919]]></title><link>https://blog.youkuaiyun.com/vicky428/article/details/95883229</link><guid>https://blog.youkuaiyun.com/vicky428/article/details/95883229</guid><author>vicky428</author><pubDate>Wed, 24 Jul 2019 12:45:48 +0800</pubDate><description><![CDATA[#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Sat Sep 15 23:00:00 2018

@author: vicky
"""

import  pandas as pd
import numpy as np  
import matplotlib.pyplot as plt
from sklearn.clus...]]></description><category></category></item></channel></rss>