24小时热门版块排行榜    

查看: 951  |  回复: 18
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 寒山听雪 的 6 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会

寒山听雪

铁虫 (小有名气)


[交流] R 包 qkerntool已经更新到 1.19版本,欢迎使用

说明

这是一个关于以核方法为主的r包介绍,这个包包含了许多用于分类,聚类降维和可视化的非线性机器学习工具。该包使用了条件负定核(cnd kernel),共集成了26个核,包含q核版本和cnd版本,其中包含dbscan,谱聚类,广义判别分析 (gda),主成分分析(pca),isomap,线性判别分析(lle),sammon mapping和tsne。

1.在r中直接输入代码在cran上查找安装

install.packages("qkerntool" )


2 直接下载好的qkerntool_1.19.tar.gz

install.packages("e:/qkerntool_1.18.tar.gz" )


3 rstudio中在菜单栏tools-install packages中查找安装
install from下拉菜单可以选择是从cran还是本地安装包安装

cnd kernel
由于用于数据分析的经典核方法通常被认为仅限于正定核,因此,在qkerntools中提出了新的距离算法,新核方法计算高维投影空间中的距离公式如下
在这里插入图片描述
qkernel使用两个向量参数之间的内核函数值,即通常依赖于由定义的内核矩阵h之类的表达式
在这里插入图片描述
更改核的参数会更改数据的投影,从而更改投影空间中的距离。
为了通过r实现负定核表达,构建了一个s4类r包,其中表示为类“qkernel”或“cndkernel”的对象。这些函数可以作为参数传递给泛型函数,这些函数可以评估更有用的内核表达式,如核矩阵

qkernmatrix 由核方法构成的核矩阵

> library("qkerntool" )
> set.seed(20181110)
> x <- matrix(rnorm(300), ncol = 100)
> qkfunc <- rbfbase(sigma=10,q=0.8)
> qkfunc
qbase qkernel function.
hyperparameter : sigma =  10  q =  0.8
> qkfunc(x[1, ], x[2, ])
         [,1]
[1,] 4.976754

> qkernmatrix(qkfunc,x)
an object of class "qkernmatrix"
         [,1]     [,2]     [,3]
[1,] 0.000000 4.976754 4.969811
[2,] 4.976754 0.000000 4.960555
[3,] 4.969811 4.960555 0.000000

   
创造自己的核
自己输入核函数

> rbf1 <- function(x, y){
+   sdist <- crossprod(x) + crossprod(y) - 2*crossprod(x,y)
+   sdist <- (1 - exp(log(0.8)*sdist/10))/(1-0.8)
+   return(sdist)
>+ }
> class(rbf1) <- "qkernel"
> rbf1(x[1, ], x[2, ])

   
qkpca
参数名         内容
x         n x d矩阵(n个样本,d个特征)或继承cndkernmatrix ,qkernmatrix的核矩阵
kernel         核类型
qpar         核参数
features         选择的成分个数
th         忽略主成分的特征值的值

以iris数据集为例,随机提取二十个样本为测试集,对剩余样本使用“rbfbase”核进行降维以及对应的预测,例子使用了iris数据集,代码及结果如下

data(iris)
test <- sample(1:150,20)
qkpc <- qkpca(~.,data=iris[-test,-5],kernel="rbfbase",
              qpar=list(sigma=50,q=0.8),features=2)

# print the principal component vectors
pcv(qkpc)
#plot the data projection on the components
plot(rotated(qkpc),col=as.integer(iris[-test,5]),
     xlab="1st principal component",ylab="2nd principal component" )

# embed remaining points
emb <- predict(qkpc,iris[test,-5])
points(emb,col=as.integer(iris[test,5]))

[ last edited by jjdg on 2019-4-26 at 20:55 ]

[ Last edited by jjdg on 2019-4-26 at 20:55 ]
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

deephill

铁杆木虫 (职业作家)



寒山听雪(金币+1): 谢谢参与
看上去很强大啊
4楼2019-04-26 22:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

寒山听雪

铁虫 (小有名气)


欢迎各位继续使用,目前发布一年多
18楼2020-04-20 15:43:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2019-04-26 09:26   回复  
寒山听雪(金币+1): 谢谢参与
2019-04-26 10:00   回复  
寒山听雪(金币+1): 谢谢参与
发自小木虫IOS客户端
2019-04-26 23:23   回复  
寒山听雪(金币+1): 谢谢参与
tzynew6楼
2019-04-26 23:52   回复  
寒山听雪(金币+1): 谢谢参与
syhorchid7楼
2019-04-27 00:08   回复  
寒山听雪(金币+1): 谢谢参与
xhmaohan8楼
2019-04-27 03:49   回复  
寒山听雪(金币+1): 谢谢参与
发自小木虫Android客户端
erjianxin9楼
2019-04-27 07:32   回复  
寒山听雪(金币+1): 谢谢参与
祝福 发自小木虫Android客户端
youngen10楼
2019-04-28 18:14   回复  
寒山听雪(金币+1): 谢谢参与
wxw082511楼
2019-04-28 21:11   回复  
寒山听雪(金币+1): 谢谢参与
发自小木虫Android客户端
liubq0812楼
2019-04-28 21:46   回复  
寒山听雪(金币+1): 谢谢参与
1 发自小木虫IOS客户端
nono200913楼
2019-05-02 10:24   回复  
寒山听雪(金币+1): 谢谢参与
·
Biotage14楼
2019-05-02 22:07   回复  
寒山听雪(金币+1): 谢谢参与
2019-05-02 22:23   回复  
寒山听雪(金币+1): 谢谢参与
96021316楼
2019-05-21 11:15   回复  
wxw082517楼
2019-09-02 07:41   回复  
hi_world_119楼
2020-05-03 17:49   回复  
相关版块跳转 我要订阅楼主 寒山听雪 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 多所高校已采取教师退出机制或教师末位淘汰机制。大家怎么看? +6 LNP@mRNA 2024-12-27 8/400 2024-12-27 11:12 by LNP@mRNA
[论文投稿] 学位论文漏引了一篇文献,怎么办 20+5 青椒1234 2024-12-25 8/400 2024-12-27 10:31 by LIU_V
[硕博家园] 2025年秋季博士招生——氢能与燃料电池课题组 +3 dawnhf 2024-12-25 3/150 2024-12-27 10:27 by Kevin-梁
[考博] 2025博士申请 +3 Apple@888 2024-12-24 8/400 2024-12-27 09:38 by Apple@888
[基金申请] 现在报基金要打招呼吗 +5 akslis2024 2024-12-25 5/250 2024-12-27 06:52 by llhljsy
[考博] 想找一位年轻的985博导 +10 艾弗森迪尔 2024-12-21 14/700 2024-12-27 06:48 by 艾弗森迪尔
[考博] 要不要加博导微信 +7 姑且随风而去 2024-12-24 16/800 2024-12-26 18:08 by 姑且随风而去
[考博] 2025博士招生 +3 生生不息厉害啊 2024-12-23 6/300 2024-12-26 10:33 by 生生不息厉害啊
[考博] 定向就业,毕业证学位证由工作单位领取? +6 河西夜郎 2024-12-22 7/350 2024-12-26 06:54 by 柳成阴
[考研] 西安电子科技大学 图情 +4 五官在线 2024-12-25 5/250 2024-12-25 23:17 by r50@zifb
[硕博家园] 四线城市公办高校无编or二线城市民办高校 +17 joshhartford 2024-12-20 37/1850 2024-12-25 20:52 by joshhartford
[教师之家] 有两个省级人才帽子跳槽找工作,难度怎么样 +9 akslis2024 2024-12-23 9/450 2024-12-25 13:01 by qsx04152006
[考博] 2025届化工专业博士申请 5+3 Hkeyan 2024-12-23 5/250 2024-12-25 08:47 by highxixi
[论文投稿] 求助,要不要根据审稿人的要求引用他们的文献? 20+3 angel8533 2024-12-24 7/350 2024-12-24 22:32 by nono2009
[教师之家] 咨询一下,是不是教授之间比较忌讳谈论文造假,因为可能很多人都是造假 +11 akslis2024 2024-12-21 12/600 2024-12-24 14:10 by 六两废铜
[考研] 26考研 +9 青云要冲沪深 2024-12-21 9/450 2024-12-24 11:09 by zjhzclm
[硕博家园] 南京大学金融科技研究中心诚招会计审计财务分析方面的顾问专家(可远程兼职) +4 cowhorce 2024-12-23 5/250 2024-12-24 08:34 by cowhorce
[论文投稿] 请问ICASSP必须要去参会吗? 5+3 火炎焱烤猪 2024-12-21 4/200 2024-12-23 16:39 by bluestork
[考博] 2025申博求助 +7 125814 2024-12-22 7/350 2024-12-22 17:07 by 毕生所学
[硕博家园] 脂质氧化 +3 song794 2024-12-22 4/200 2024-12-22 14:12 by song794
信息提示
请填处理意见