楼主: liqiaoshan
705 4

[统计软件] 求一个R抽样代码 [推广有奖]

  • 0关注
  • 0粉丝

等待验证会员

高中生

12%

还不是VIP/贵宾

-

威望
0
论坛币
2079 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
142 点
帖子
25
精华
0
在线时间
22 小时
注册时间
2015-6-7
最后登录
2020-6-4

楼主
liqiaoshan 发表于 2018-12-6 14:55:29 |AI写论文
1000论坛币
琢磨了好久,实在解决不了,向各位大神求教,问题如下:我有一个调查样圆,面积是1000平米(半径约17.8米),调查了样圆里所有的树木,记录了每个树木到圆心的距离和偏角(树木到圆心的连线与正北方向的夹角),然后我要在这个1000平米样圆里随机抽10个100平米的样圆,返回每个100平米样圆里的所有树木信息,请问在R里如何实现,谢谢!!非常感谢!!
sampledata.xlsx (13.55 KB)

最佳答案

atijinmao 查看完整内容

randsample.fn=function(data,area1,area2,N){ #area1 is the area of observed plot, area2 is the area of expected plot r1=sqrt(area1/pi) r2=sqrt(area2/pi) expplot=data.frame() for (j in 1:N){ R=runif(1,0,r1-r2) theta=runif(1,0,360) xa=R*sin(theta*pi/180) ya=R*cos(theta*pi/180) expdata=data[sqrt((data$x-xa)^2+(data$y-ya)^2)
关键词:如何实现 距离和

回帖推荐

atijinmao 发表于4楼  查看完整内容

randsample.fn=function(data,area1,area2,N){ #area1 is the area of observed plot, area2 is the area of expected plot r1=sqrt(area1/pi) r2=sqrt(area2/pi) expplot=data.frame() for (j in 1:N){ R=runif(1,0,r1-r2) theta=runif(1,0,360) xa=R*sin(theta*pi/180) ya=R*cos(theta*pi/180) expdata=data[sqrt((data$x-xa)^2+(data$y-ya)^2)

沙发
atijinmao 发表于 2018-12-11 10:31:19
randsample.fn=function(data,area1,area2,N){   #area1 is the area of observed plot, area2 is the area of expected plot
  
  r1=sqrt(area1/pi)
r2=sqrt(area2/pi)
expplot=data.frame()

  for (j in 1:N){

R=runif(1,0,r1-r2)
    theta=runif(1,0,360)
    xa=R*sin(theta*pi/180)
ya=R*cos(theta*pi/180)
expdata=data[sqrt((data$x-xa)^2+(data$y-ya)^2)<r2,]
    lable=paste("plot",j,sep="")
    expdata[,6]=lable
    colnames(expdata)[6]="plot"
    expplot=rbind(expplot,expdata)
  }
  return(expplot)
}

已有 1 人评分论坛币 收起 理由
yangyuzhou + 100 热心帮助其他会员

总评分: 论坛币 + 100   查看全部评分

藤椅
liqiaoshan 发表于 2018-12-6 15:13:55
或者大家有没有什么函数或者思路可以帮助确定大圆里每个树的位置{:2_26:}

板凳
liqiaoshan 发表于 2018-12-6 22:01:37
走过路过求分享相关经验啊{:2_35:}

报纸
liqiaoshan 发表于 2018-12-17 13:09:21
atijinmao 发表于 2018-12-6 14:55
randsample.fn=function(data,area1,area2,N){   #area1 is the area of observed plot, area2 is the area ...
太棒了!!!感激涕零

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
jg-xs1
拉您进交流群
GMT+8, 2025-12-22 21:39