小编use*_*089的帖子

使用R中的ggmap在get_map中出错

大约90%的时间,在使用get_mapggmap,我收到以下错误.有人可以告诉我为什么吗?

map <- get_map(location = 'Australia', zoom = 4)
Run Code Online (Sandbox Code Playgroud)

在download.file错误(!URL,destfile = destfile,安静=短信,模式= "WB"):无法打开URL" http://maps.googleapis.com/maps/api/staticmap?center=Australia&zoom=4&size= %20640x640&scale =%202&maptype = terrain&sensor = false '另外:警告消息:在download.file(url,destfile = destfile,quiet =!messaging,mode ="wb"):无法打开:HTTP状态为'403 Forbidden'

r ggmap

26
推荐指数
2
解决办法
2万
查看次数

使用dendextend在R中绘制tanglegrams子图

我正在使用dendextend在R中绘制Tanglegrams.我想知道是否有可能使用多个子图par(mfrow = c(2,2))

我似乎无法弄明白.

谢谢

library(dendextend)
dend15 <- c(1:5) %>% dist %>% hclust(method = "average") %>% as.dendrogram
dend15 <- dend15 %>% set("labels_to_char")
dend51 <- dend15 %>% set("labels", as.character(5:1)) %>% match_order_by_labels(dend15)
dends_15_51 <- dendlist(dend15, dend51)

par(mfrow = c(2,2))
tanglegram(dends_15_51)
tanglegram(dends_15_51)
tanglegram(dends_15_51)
tanglegram(dends_15_51)
Run Code Online (Sandbox Code Playgroud)

plot visualization r dendrogram dendextend

10
推荐指数
1
解决办法
905
查看次数

在条形图中设置条形图的透明度,并将y轴设置为对数比例 - 但两者似乎都不能在MATLAB中工作

在MATLAB中,我想将条形图中条形图的透明度设置为0.3,并将y轴设置为对数刻度 - 但两者似乎都不起作用......

subplot('Position',[0.15 0.7 0.45 0.15]);
data = [1 2 5 4 7 4 1];
B = bar(data,'g');
ch = get(B,'child');
set(ch,'facea',.3)
Run Code Online (Sandbox Code Playgroud)

这工作正常,但然后添加:

set(gca,'YScale','log');
Run Code Online (Sandbox Code Playgroud)

并且透明度设置不起作用.有任何想法吗?谢谢!

matlab matlab-figure

7
推荐指数
1
解决办法
916
查看次数

设置图例的Alpha值(透明度)以匹配MatLab中的绘图中的alpha

我之前看过这个问题,但似乎没有解决方案,所以我只是想知道它是否完全可能.

我在MatLab中有一个条形图并设置了透明度:

B = bar(x,y,'stacked');
set(B(1),'facecolor',[0 0.3906 0]) 
set(B(2),'facecolor',[0.5625 0.9297 0.5625])
ch1 = get(B(1),'child');
set(ch1,'facea',.5)
ch2 = get(B(2),'child');
set(ch2,'facea',.5)
Run Code Online (Sandbox Code Playgroud)

我希望情节中的透明度能够反映在传奇中:

BL = legend ((B([1 2])),{'data1','data2'},'fontsize',10);
Run Code Online (Sandbox Code Playgroud)

但是,图例中的alpha值似乎为1.

有任何想法吗?谢谢.

matlab

6
推荐指数
2
解决办法
1万
查看次数

在MATLAB中有一个更好/更快的方法随机改组矩阵吗?

在MATLAB中,我使用shake.m函数(http://www.mathworks.com/matlabcentral/fileexchange/10067-shake)随机随机播放每列.例如:

a = [1 2 3; 4 5 6; 7 8 9]
a =

     1     2     3
     4     5     6
     7     8     9

b = shake(a)
b =

     7     8     6
     1     5     9
     4     2     3
Run Code Online (Sandbox Code Playgroud)

这个函数正是我想要的,但是我的列长(> 10,000,000),所以这需要很长时间才能运行.有谁知道更快的实现方法?我试过分别摇动每个列向量,但这并不快.谢谢!

arrays matlab matrix

5
推荐指数
3
解决办法
1855
查看次数

与 R 中的链接面对面绘制系统发育树

我想使用该ape包在 R 中绘制两个相对的系统发育图。一棵树有 40 个节点,一棵树有 26 个节点:

library(ape)
tree1 <- rtree(40)
tree2 <- rtree(26)
Run Code Online (Sandbox Code Playgroud)

cophyloplot函数使用指定的链接面对面地绘制这些图。

我在指定链接时遇到问题。

请注意,在我的实际nexus树文件中,提示标签是文本(如果需要,我不确定如何将它们更改为数字......)。

链接应如下所示:

如果,在tree1Nexus文件中,序列的尖端标签是1-40。在tree2Nexus 文件中,提示标签为 1-26。那么链接应该是:

a <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40)
b <- c(14,1,4,1,9,12,2,10,6,3,13,5,14,15,18,19,19,7,14,9,10,11,25,22,21,16,23,24,26,17,1,12,12,21,15,16,21,8,20,21) 
association <- cbind(a, b)
Run Code Online (Sandbox Code Playgroud)

(即 中的序列 1 与 中tree1的序列 14 连接tree2

所以,我用这样的东西来绘制树木:

cophyloplot(tree1, tree2, assoc=association,length.line=4, space=28, gap=10, rotate=TRUE)
Run Code Online (Sandbox Code Playgroud)

并计算距离矩阵:

dist.topo(tree1, tree2, method = "PH85")
Run Code Online (Sandbox Code Playgroud)

我不太确定我哪里出错了。任何帮助,将不胜感激!

r phylogeny

5
推荐指数
1
解决办法
2392
查看次数

在matlab中,如何使用imagesc在2D热图中"平滑"像素

使用MATLAB,我有一个矩阵(data),我正在绘制imagesc(data)用于生成热图:

data = [1 1 1 1 1 1 1 1 1 1; 1 1.04 1.04 1.04 1.03 1 1.01 1.01 1.03 1.01; 1.36 1.3 1.25 1.2 1.15 1.1 1.2 1.13 1.07 1.11; 3.65 3.16 2.94 2.68 2.39 2.22 2.17 1.95 1.79 1.81; 5.91 5.75 5.47 5.3 4.98 4.79 4.62 4.55 4.38 4.19; 6 6 5.99 5.83 5.49 5.33 5.14 4.94 4.77 4.74];
imagesc(data)
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

有没有办法"平滑"像素,以产生更像这样的东西: 在此输入图像描述

matlab image image-processing

4
推荐指数
2
解决办法
6469
查看次数

R geom_point和ggmap

我想在澳大利亚的R中绘制一个ggmap,并在其数据点上放置一层标记,这些标记对应于以下数据指定的大小:

sitename   lat    lon    sitenum
Sydney     -34    151    1
Melbourne  -37    144    4
Adelaide   -34    138    7
Run Code Online (Sandbox Code Playgroud)

这是我的代码,但是不起作用...

library(ggmap)
map <- get_map(location = 'Australia', zoom = 4)
mapPoints <- ggmap(map) + geom_point(aes(x = lon, y = lat, size = sitenum), alpha = .5)
Run Code Online (Sandbox Code Playgroud)

r ggmap

3
推荐指数
1
解决办法
6820
查看次数

带散点图数据的 Matlab“热图”

抱歉,如果这个问题已在其他地方得到解答,但我无法找到解决方案......

我有scatterplot数据、xy、 和颜色数据c(所有这些都是向量)。

x = rand(100,1);
y = rand(100,1);
c = rand(100,1)*50;
scatter(x,y,c,'fill')
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

但是,我想将其显示为“堆图”,这样c它就不再是大小,而是表示强度,看起来像这样: 在此输入图像描述

我不想使用surf并且希望成为这样的 2D。如果您能提供帮助,谢谢。

matlab

2
推荐指数
1
解决办法
5145
查看次数

Matlab未定义变量

说我有以下代码:

[w1, d1]=fit1; 
    for i = 1:10, 
      p = w1*d1+i;
      C(i,:) = p;
    end

[w2, d2]=fit2; 
    for i = 1:10, 
      q = w2*d2+i;
      D(i,:) = q;
    end
Run Code Online (Sandbox Code Playgroud)

功能fit1:

function[w1, d1] = fit1
w1 = rand(1);
d1 = rand(1);
Run Code Online (Sandbox Code Playgroud)

和功能fit2:

function[w2, d2] = fit2
w2 = w1+0.2;
d2 = d1-0.1;
Run Code Online (Sandbox Code Playgroud)

我收到错误: Undefined function or variable 'w1'.

但是w1之前定义的参数是fit2被调用的,为什么这不起作用呢?

谢谢!

matlab

1
推荐指数
1
解决办法
168
查看次数

在matlab中,找到矩阵中唯一行出现的频率

在Matlab中,我说有以下矩阵,代表10个人口:

pop = [0 0 0 0 0; 1 1 1 0 0; 1 1 1 1 1; 1 1 1 0 0; 0 0 0 0 0; 0 0 0 0 0; 1 0 0 0 0; 1 1 1 1 1; 0 0 0 0 0; 0 0 0 0 0];
Run Code Online (Sandbox Code Playgroud)

oneszeros定义6个不同"类型"的个体.

a = [0 0 0 0 0];
b = [1 0 0 0 0];
c = [1 1 0 0 0];
d …
Run Code Online (Sandbox Code Playgroud)

arrays matlab matrix

1
推荐指数
1
解决办法
164
查看次数