大约90%的时间,在使用get_map时ggmap,我收到以下错误.有人可以告诉我为什么吗?
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'
我正在使用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) 在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中有一个条形图并设置了透明度:
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中,我使用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),所以这需要很长时间才能运行.有谁知道更快的实现方法?我试过分别摇动每个列向量,但这并不快.谢谢!
我想使用该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)
我不太确定我哪里出错了。任何帮助,将不胜感激!
使用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)

有没有办法"平滑"像素,以产生更像这样的东西:

我想在澳大利亚的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) 抱歉,如果这个问题已在其他地方得到解答,但我无法找到解决方案......
我有scatterplot数据、x和y、 和颜色数据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。如果您能提供帮助,谢谢。
说我有以下代码:
[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中,我说有以下矩阵,代表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)
行ones和zeros定义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) matlab ×7
r ×4
arrays ×2
ggmap ×2
matrix ×2
dendextend ×1
dendrogram ×1
image ×1
phylogeny ×1
plot ×1