s******a 发帖数: 184 | 1 假如一个变量 X 是 另外几个变量的函数, X=f(Y1, Y2,...,Yn), 假设Y1,。。YN
的分布是知道的,如何得出X的 分布
是不是 用Monte Carlo simulation , 但Y1,...YN之间是相关的,这样的情况如何
处理? |
|
s******a 发帖数: 184 | 2 假如一个变量 X 是 另外几个变量的函数, X=f(Y1, Y2,...,Yn), 假设Y1,。。YN
的分布是知道的,如何得出X的 分布
是不是 用Monte Carlo simulation , 但Y1,...YN之间是相关的,这样的情况如何
处理? |
|
s******a 发帖数: 184 | 3 假如一个变量 X 是 另外几个变量的函数, X=f(Y1, Y2,...,Yn), 假设Y1,。。YN
的分布是知道的,如何得出X的 分布
是不是 用Monte Carlo simulation , 但Y1,...YN之间是相关的,这样的情况如何
处理? |
|
d*****n 发帖数: 754 | 4 Let x2=y2... xn=yn, use pdf(y) / det(jocabian(y, x) ) as the pdf (x)
[在 surfreta (retababy) 的大作中提到:]
:假如一个变量 X 是 另外几个变量的函数, X=f(Y1, Y2,...,Yn), 假设Y1,。。YN
的分布是知道的,如何得出X的 分布
:是不是 用Monte Carlo simulation , 但Y1,...YN之间是相关的,这样的情况如
何处理? |
|
n***y 发帖数: 2730 | 5 你要是没有Y1,...,YN的joint distribution(联合分布),这是没有唯一答案的。一
般在这种情况下是做一些model假设,比如One factor Gaussian Copula,然后你要有
些data去calibrate一下。总之你这样general的问题是没有唯一解的。 |
|
n***y 发帖数: 2730 | 6 你要是没有Y1,...,YN的joint distribution(联合分布),这是没有唯一答案的。一
般在这种情况下是做一些model假设,比如One factor Gaussian Copula,然后你要有
些data去calibrate一下。总之你这样general的问题是没有唯一解的。 |
|
W**o 发帖数: 547 | 7 Do you know the correlation between Y1,...YN?
YN |
|
d*****r 发帖数: 39446 | 8 【此篇文章是由自动发信系统所张贴】
y1 已经退出本俱乐部, 特此通知. |
|
d*****r 发帖数: 39446 | 9 【此篇文章是由自动发信系统所张贴】
y1 申请加入本俱乐部的请求已经通过审批, 成为本俱乐部的正式成员, 特此通知. |
|
m**********r 发帖数: 507 | 10 1. No. e.g., n=2 set Y1 = Y2 = X = Ber(1/2)
-> Y1|X is independent of Y2|X
(1) P(Y1=1, Y2=1|X=1) = 1 = P(Y1=1|X=1) P(Y2=1|X=1)
(2) P(Y1=1, Y2=0|X=1) = 0 = P(Y1=1|X=1) P(Y2=0|X=1)
(3) P(Y1=0, Y2=1|X=1) = 0 = P(Y1=0|X=1) P(Y2=1|X=1)
(4) P(Y1=0, Y2=0|X=1) = 0 = P(Y1=0|X=1) P(Y2=0|X=1)
similarly check ALL combinations of Y1 and Y2 conditioning on the event X=0
-> marginally, both Y1 and Y2 are Ber(1/2)
Thus P(Y1=1) * P(Y2=1) = (1/2) * (1/2) = 1/4
-> But Y1 is NOT independent of Y2, because
P(Y1=1, Y... 阅读全帖 |
|
x********g 发帖数: 595 | 11 Seriously,
r -> (r2^2 - X2^2 + Xc^2 - Y2^2 + Yc^2 + (Y2*(4*r1^2*r2^2*Y1 - 4*r1*r2^3*Y1
- 4*r2^2*X1^2*Y1 + 4*r2^2*X1*X2*Y1 - 4*r1^2*X2^2*Y1 + 4*r1*r2*X2^2*Y1 +
4*X1^2*X2^2*Y1 - 4*X1*X2^3*Y1 + 4*r2^2*X1*Xc*Y1 + 8*r1^2*X2*Xc*Y1 -
8*r1*r2*X2*Xc*Y1 - 4*r2^2*X2*Xc*Y1 - 8*X1^2*X2*Xc*Y1 + 4*X1*X2^2*Xc*Y1 + 4*
X2^3*Xc*Y1 -
4*r1^2*Xc^2*Y1 + 4*r1*r2*Xc^2*Y1 + 4*X1^2*Xc^2*Y1 + 4*X1*X2*Xc^2*Y1
- 8*X2^2*Xc^2*Y1 - 4*X1*Xc^3*Y1 + 4*X2*Xc^3*Y1 - 4*r2^2*Y1^3 + 4*X2^2*Y1^3 -
8*X2*Xc*Y1^3 + |
|
s***e 发帖数: 911 | 12
y{\dot}a=d is one equation:
a1*y1+a2*y2+a3*y3=d;
y{\cross}a are three equations:
a3*y2-a2*y3=c1
-a3*y1+a1*y3=c2
a2*y1-a1*y2=c3.
The matrix A={{0,a3,-a2},{-a3,0,a1},{a2,-a1,0}}is singhular(the rank is 2.
So you have to solve this equation accampany with a1*y1+a2*y2+a3*y3=d.
Express y2,y3 by y1:
y2=(a2/a1)*y1-c3/a1;
y3=(a3/a1)*y1+c2/a1.
Then
a1*y1+a2*y2+a3*y3={y1*a1+[(a2^2/a1)*y1-c3*a2/a1]+[(a3^2/a1)*y1+c2*a3/a1]}
=d
so,
y1=[a1*d-(c \cross a)_1]/|a|^2;
Then u can get y2,y3 by sim |
|
w******a 发帖数: 25 | 13 Here is an R example to impute one or two missing data in each record:
The data will look like
col1 col2 col3
x
x x x
x x
x x
x x x
x
x x x
...
library(Rlab)
alp = 1
K_delta = 2
len_Y1 = 200
#Sample setting:
#Measurment N_
patient Percent
# 1 12
0.18
# 1 2 4
0.05... 阅读全帖 |
|
p***l 发帖数: 134 | 14 Use a conditional-unconditional formula, the idea is crystal clear and easy.
It is a classic technique for martingale/stochastic process problem.
Let X=# of tosses until you observe two heads in a row, Y1=the result of the
first toss. Also denote x=E(X) for consistency with your equation.
x=E(X)=E[E(X|Y1)]= E(X|Y1=T)/2 + E(X|Y1=H)/2
Note that Y1 is independent of the remaining tosses. Provided Y1=T, the "
head counter" is reset. Hence E(X|Y1=T)=E(1+X)=1+x. Provided Y1=H, you have
.5:.5 condition... 阅读全帖 |
|
r********n 发帖数: 7441 | 15 发信人: conti (conti), 信区: TsinghuaCent
标 题: zz pk饶毅进入第二轮的张旭的背景
发信站: 水木社区 (Sun Aug 28 15:38:09 2011), 站内
如果神经方向就只能上一个候选人,那么张旭入选还真是一点也不出格
发信人: WeiLiao (WeiLiao), 信区: Biology
标 题: 张旭的publications
发信站: BBS 未名空间站 (Sun Aug 28 01:39:20 2011, 美东)
去查了一下,张旭总共发表了90多篇英文文章(不算中文)。回国之后作为联系作者也
有快20篇了。其中有cell, neuron, Journal of neuroscience, PNAS 好多篇。
而且他
1995年就回国了,那个时候国内的环境有多差不用说了吧。这也快20年了。比饶益
2007
年才全职回去多了12年。比起培养的中国学生也多多了,算国内贡献肯定不会小。如果
他一直留在国外,成就未必就比饶差。谁把饶益的文章列一列。国内多少篇,总共多少
篇,培养了多少中国学生?
1. Li KC, Wang F, ... 阅读全帖 |
|
k***t 发帖数: 276 | 16 interviewstreet上的。brute-force很straightforward。
有没有妙招降复杂度?
===================================================
Vertical Sticks
Given array of integers Y=y1,...,yn, we have n line segments such that
endpoints of segment i are (i, 0) and (i, yi). Imagine that from the top of
each segment a horizontal ray is shot to the left, and this ray stops when
it touches another segment or it hits the y-axis. We construct an array of n
integers, v1, ..., vn, where vi is equal to length of ray shot from the top... 阅读全帖 |
|
W*****o 发帖数: 1780 | 17 去查了一下,张旭总共发表了90多篇英文文章(不算中文)。回国之后作为联系作者也
有快20篇了。其中有cell, neuron, Journal of neuroscience, PNAS 好多篇。而且他
1995年就回国了,那个时候国内的环境有多差不用说了吧。这也快20年了。比饶益2007
年才全职回去多了12年。比起培养的中国学生也多多了,算国内贡献肯定不会小。如果
他一直留在国外,成就未必就比饶差。谁把饶益的文章列一列。国内多少篇,总共多少
篇,培养了多少中国学生?
1. Li KC, Wang F, Zhong YQ, Lu YJ, Wang Q, Zhang FX, Xiao HS, Bao L, Zhang X
. (2011) Reduction of follistatin-like 1 in primary afferent neurons
contributes to neuropathic pain hypersensitivity. Cell Res, 21: 697-699.
2. Li KC, Zhang FX, Li CL, Wang F, Yu MY, Z... 阅读全帖 |
|
c*******y 发帖数: 98 | 18 DP好难,我来贴个无脑的,不知道题目理解对了没有。。
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
int findMaxAllOneTopLeftSquare(vector>& matrix)
{
int i, j;
int m = matrix.size(); if (!m) return 0;
int n = matrix[0].size(); if (!n) return 0;
int max = min(m, n);
for (i=0; i
for (j=0; j<=i; j++){
if (matri... 阅读全帖 |
|
H****h 发帖数: 1037 | 19 椭圆上点(x1,y1)的外法向量方向是(x1/a^2,y1/b^2).
(x1,y1)和(x0,y0)的联线方向是(x0-x1,y0-y1).
你需要两者方向一致,于是有
(x1/a^2)/(x0-x1)=(y1/b^2)/(y0-y1)>0
所以有:(x1/a^2)(y0-y1)=(y1/b^2)(x0-x1).
得PQ |
|
k**y 发帖数: 920 | 20 Corr(X1+X2,Y1+Y2)可以等于 [-1,1]上任何数,对吗?
用老周的方法,画图,X1,Y1是夹角小于arccos0.9的两个向量,长度任意
X2,Y2是夹角小于arccos0.9的两个向量,长度任意
问向量X1+X2和Y1+Y2的夹角。
画个图就可以发现,这夹角可以[0,180]
比如假设corr(x1,y1)=corr(x2,y2)
那么当x2=y1, y2=x1的时候,corr(x1+x2,y1+y2)=1
当x2=-y1, y2=-x1的时候,corr(x1+x2,y1+y2)=-1 |
|
k**y 发帖数: 920 | 21 Corr(X1+X2,Y1+Y2)可以等于 [-1,1]上任何数,对吗?
用老周的方法,画图,X1,Y1是夹角小于arccos0.9的两个向量,长度任意
X2,Y2是夹角小于arccos0.9的两个向量,长度任意
问向量X1+X2和Y1+Y2的夹角。
画个图就可以发现,这夹角可以[0,180]
比如假设corr(x1,y1)=corr(x2,y2)
那么当x2=y1, y2=x1的时候,corr(x1+x2,y1+y2)=1
当x2=-y1, y2=-x1的时候,corr(x1+x2,y1+y2)=-1 |
|
s***e 发帖数: 911 | 22
首先这两个线段不一定有交点. 我来作个完全的推导:
线段一定义了一个矢量:
L1=((x2-x1),(y2-y1))
线段二定义另外一个适量:
L2=((x2'-x1'),(y2'-y1'))
线段一上的任何一个点R(p)可以这样确定:
R1(p)=(x1,y1)+p*L1
其中p是个参量,在[0,1]间变化. p=0就对应端点(x1,y1);p=1对应(x2,y2)
同理线段二上任何一点被确定为:
R2(q)=(x1',y1')+q*L2
对某个特殊的值(p,q), R1(p)=R2(q). 我们就来求这对(p,q)存在的条件. 如果(p,q)
存在,并且都在[0,1]内, 这两线段就相交,反之就不交. 于是就得两个方程:
(x1-x2)+p*(x2-x1)-q*(x2'-x1')=0
(y1-y2)+p*(y2-y1)-q*(y2'-y1')=0
这就是个二维线性方程组, 你解(p,q)就可以乐... |
|
S*****0 发帖数: 15 | 23 My solution is:
data t1;
input x y;
cards;
1 5
2 6
3 6
7 9
8 10
11 12
13 15
14 17
;
run;
data t2(rename=(m=x n=y));
set t1 end=eof;
retain x1 y1 m n;
if _N_=1 then
do;
x1=x;
y1=y;
end;
else
do;
if y1
do;
m=x1; n=y1;
output;
x1=x; y1=y;
end;
else
... 阅读全帖 |
|
l******r 发帖数: 18699 | 24 answer is no change 给你个proof
let x1,x2,...,xn be iid rv, its ordered variables are y1<=y2<=..<=yn
let P_Y be pdf of yi's and P_X be the pdf of xi's, let pdf of xi be F
note that P_Y(y1,y2,...,yn)=P(Y1<=y1,..,Yn<=yn)=P(Y1<=y1)...P(Yn<=yn)=F(y1).
..F(yn)
P_X(x1,...,xn)=F(x1)...F(xn)
the information-theoretic entropy of yi's is
-sum F(yi) ln F(yi)=-sum F(xi) ln F(xi)= the entropy of xi's
so the information-theoretic entropy is the same after reordering the rv |
|
b********2 发帖数: 5 | 25 我现在对如下的时滞微分方程组写了个小code:
y1'=-a*y1^2-pai*S*(c-y1-y2)+k*y2*y2(t-\tau);
y2'=ip*y2*(c-y1-y2)-k*y2*y2(t-\tau)
(这里\tau是时滞。 y1,y2是变量, a, pai, k, c, ip是参数。
\tau=2, t in [0, 5].)
这个程序运行起来没有问题,但是我不知道我的code是不是正确的描述了这个方程组。
1。 特别是我的ylag1(2)的表达
2。 还有在figure(2) 和figure(4)中,我想分别画 时间t与 y1, y1与y2的图形,
不知
道我的code的表达对不对? 我用了 sol.y(1,:), sol.y(2,:).这样写是正确的吗?
谢了
。
function v=manuf(t, Z, y)
v=zeros(2, 1);
a=2; pai=3; c=20; k=3; ip=2;
ylag1=Z(:,1);
v(1)=-a*y(1)*y(1)-pai*y(1)*(c-y(1)-y(2))+k*y(2)*ylag1(2);
v(2)=ip*... 阅读全帖 |
|
b********2 发帖数: 5 | 26 我现在对如下的时滞微分方程组写了个小code:
y1'=-a*y1^2-pai*S*(c-y1-y2)+k*y2*y2(t-\tau);
y2'=ip*y2*(c-y1-y2)-k*y2*y2(t-\tau)
(这里\tau是时滞。 y1,y2是变量, a, pai, k, c, ip是参数。
\tau=2, t in [0, 5].)
这个程序运行起来没有问题,但是我不知道我的code是不是正确的描述了这个方程组。
1。 特别是我的ylag1(2)的表达
2。 还有在figure(2) 和figure(4)中,我想分别画 时间t与 y1, y1与y2的图形,
不知
道我的code的表达对不对? 我用了 sol.y(1,:), sol.y(2,:).这样写是正确的吗?
谢了
。
function v=manuf(t, Z, y)
v=zeros(2, 1);
a=2; pai=3; c=20; k=3; ip=2;
ylag1=Z(:,1);
v(1)=-a*y(1)*y(1)-pai*y(1)*(c-y(1)-y(2))+k*y(2)*ylag1(2);
v(2)=ip*... 阅读全帖 |
|
f********t 发帖数: 14 | 27 一个貌似简单的问题,觉得对但又不敢确定,请高手指教。
U(x,y) is a continuously differentiable function of (x,y).
For three fixed values y1, y2 and x0,
a*U(x0,y2)+(1-a)*U(x0,y1)>U(x0,a*y2+(1-a)*y1),
where 0=
Add a small term k*a*(1-a) to x0, for k small enough, I claim the inequality
is still true, as follows,
a*U(x0+k*a*(1-a),y2)+(1-a)*U(x0+k*a*(1-a),y1)>U(x0,a*y2+(1-a)*y1).
Let f(a)=a*U(x0+k*a*(1-a),y2)+(1-a)*U(x0+k*a*(1-a),y1)-U(x0,a*y2+(1-a)*y1),
the diff |
|
f*******a 发帖数: 663 | 28 开始忘贴代码了,有朋友要求,就把修改后的代码贴在这里。改动不多,可以部分提升
效率。原来的也没删,注释掉了。供参考。
=========================================================================
clear;
close all
disp('The only input needed is a distance matrix file')
disp('The format of this file should be: ')
disp('Column 1: id of element i')
disp('Column 2: id of element j')
disp('Column 3: dist(i,j)')
if(0)
% mdist=input('name of the distance matrix file (with single quotes)?\n'
);
mdist = 'example_distances.dat';
disp('Reading input dis... 阅读全帖 |
|
w****x 发帖数: 2483 | 29 sqrt 就不做了,二分比较熟悉.
second largest:
//Second largest integer
void _inner_2nd(int a[], int n, int& n1, int& n2)
{
if (1 == n)
{
n1 = a[0];
n2 = INT_MIN;
return;
}
int x1,x2;
_inner_2nd(a, n/2, x1, x2);
int y1,y2;
_inner_2nd(a+n/2, n-n/2, y1, y2);
if (x1 > y1)
{
n1 = x1;
if (y1 > x2)
n2 = y1;
else n2 = x2;
}
else
{
n1 = y1;
if (x1 > y2)
n2 = x1;
else n2 = y2... 阅读全帖 |
|
j********g 发帖数: 80 | 30 int twoRobs(vector > &MP){
int M = MP.size();
int N = MP[0].size();
int P_LEN = M + N -1;
vector > dp1 = vector >(N, vector(N, 0));
vector > dp2 = vector >(N, vector(N, 0));
dp1[0][0] = MP[0][0];
dp2[0][0] = MP[0][0];
for (int p = 2; p <= P_LEN; p++){//p is the length of the path
for (int x1 = 0; x1 < N; x1++){
for (int x2 = x1; x2 < N; x2++){
int y1 = p - 1 -... 阅读全帖 |
|
f*******t 发帖数: 7549 | 31 暴力做法蛮简单的吧,所有矩形面积之和减掉矩形相交面积之和
func CalcTotalArea(rectangles []Rectangle) int {
area := 0
for r, idx := range rectangles {
area += (r.x2 - r.x1) * (r.y2 - r.y1)
for i := idx + 1; i < len(rectangles); i++ {
r2 := rectangles[i]
x1 := max(r.x1, r2.x1)
y1 := max(r.y1, r2.y1)
x2 := min(r.x2, r2.x2)
y2 := min(r.y2, r2.y2)
if x1 < x2 && y1 < y2 {
area -= (x2 - x1) * (y2 - y1)
}
}
}
return area
} |
|
a*******g 发帖数: 1221 | 32 噢,对,我想起来了,我当时也是这样的。我之前提到,你可以算一下斜率,如果接近
竖直,就循环Y值。当然这种方法也不太好。我想到能不能这样。
我先计算一个单位长的方向向量t=(dx,dy)=(x1-x2,y1-y2)/sqrt((x1-x2)^2+(y1-y2)^2)
然后从某一点开始,利用这个单位长的方向向量一直延伸,把剩下的点都画出来
譬如(x1+dx,y1+dy), (x1+2dx,y1+2dy), ..., (x1+ndx,y1+ndy)。画直线的话就要从(
x1,y1)向(dx,dy)延伸,也同时向(-dx,-dy)延伸。
你说的情况我都考虑到了。 |
|
b***y 发帖数: 2799 | 33 ☆─────────────────────────────────────☆
observer (笑看人生) 于 (Tue Jun 12 13:04:21 2007) 提到:
N * N的 点阵,点有黑有白,找一个矩形,四角点全黑,面积最大
好象还没人给解法,继续攒人品。
1. Scan 矩阵,找黑点,记录坐标(xi,yi),存到array, O(N^2)
得到 M 个黑点
2.a 土办法是scan 黑点array, 对任意2个点(x1, y1), (x2, y2)
如果(x1, y2), (x2, y1)也是黑点,(直接查原点阵)
那么是一个合格矩形,
找出max area O(M^2)
2.b 另外一个办法,根据x,y坐标,sort 黑点array,
得到2个sorted array, AX, AY, O(MlgM)
scan AX, 对2个x座标相同的黑点,(x0, y1), (x0, y2)
查 AY找y1, y2, 在y1, y2的黑点中找x座标相同的点,
(x3,y1), (x3,y2), x3 != x0, 那么就是一个合格的矩形
找到max area
一般情 |
|
r****y 发帖数: 1437 | 34 【 以下文字转载自 Mathematics 讨论区 】
【 原文由 rossby 所发表 】
假设测量一个物理量随时间变化的情况
y = y(t) t=1, 2, 3, 4, ...,
每一时刻测量的误差也是知道的,为err(t),每一次的误差都是独立随即变量
如果此物理量随时间的变化可以表示成
y0(t)= at + b + y1(t)
where y1(t+T)= y1(t)
如果先用线性拟合估计a, b及其error bar,
y1(t)的数值估计 mean{(y(t)-at-b)} for t=i, i+T, i+2T,..., i+NT
如果我想给出y1(t)的数值估计的error bar,该怎么给出呢?
如果没有线性函数at+b,误差err(t)=err=Costant, 那么y1(t)的
error bar 是 err/ sqrt(N). 加了这个线性项, 我就不知道怎么做了。 |
|
s****h 发帖数: 921 | 35 2个曲线都是有N个采样点组成.
y1(k) y2(k)
In MATLAB:
y1,y2 are columen vector.
I use the index=
sqrt((y1-y2)'*(y1-y2))
__________________________
sqrt(y1'*y1)
Is it reasonable? |
|
b*******n 发帖数: 5065 | 36
here?
[r,x0,y0]=
[[r = 1/2*(r1^2*X2-Y1^2*X2+Y1^2*Xc-X1*Yc^2-X2*X1^2+X1^2*Xc-X1*Xc^2+X2*Xc^2
+Yc^2*X2-r2^2*X1+X2^2*X1+Y2^2*X1+r2^2*Xc-X2^2*Xc-Y2^2*Xc-2*Y1*RootOf(-2*r2^2
*Xc^3*X1+2*r1^3*r2*Xc^2-2*r1^2*r2^2*Y1^2+2*X1^4*Xc*X2-2*r1^2*r2^2*Xc^2+2*r2^
2*Yc^2*Xc^2+2*r1^2*X2^2*X1^2-2*r1^2*Xc^3*X1+2*r1^3*r2*X2^2-2*r1^3*r2*Yc^2-2*
Y1^2*X1*Yc^2*X2+2*r1^3*r2*Y2^2-2*X1^3*Xc^2*X2-2*r1^2*X2*Xc^3+2*r1^2*X2^2*Xc^
2+2*r1^2*Y1^2*X2^2+4*r1^2*r2^2*Yc^2-2*r1^2*Y2^2*r2^2-2*r1*r2*Yc^4+r2^2*Xc^4+
r2^2*Yc^4+r2^2 |
|
o*******w 发帖数: 349 | 37 能不能进一步指点一下. Chebshev极限定理的英文名称是什么? Chebeshev limit
theorem?
我感觉这个问题可规结为这么个问题:
我们可以把函数的取值看成在一个随机抽样.
k
Y1= SUM pi*yi
1
N
Y2 = SUM p'i * y'i
1
yi 和 y'i 在同一个样本空间中去值. 我们知道 Y1 的方差比Y1的大.
但是我需要证明 max(Y1) > max(Y1) Or y 是一个恒常函数 即 Y1 = O(Y2) |
|
l******i 发帖数: 1404 | 38 我只会最笨的解法,
Step1: Given Y1, Y2 are iid from N(0,1), you can write the pdf of jointly
distribution of (Y1, Y2), denoted by f(y1, y2).
Step2: Change of variable:
X1 = Y1;
X2 = \rho*Y1+((1-\rho^2)^(1/2))*Y2;
Then X1, X2 are two random variable ~ N(0,1) with correlation \rho.
Now the pdf of (X1, X2) should be
f(y1(x1,x2), y2(x1,x2))/abs(Jacobi determinant of the change of variable)
Step3:
cdf of (X1, X2) can be calculate from step2, denoted by F(x1, x2).
Thus F(m,m) is the cdf of max(X1,X2).
Thus F'(... 阅读全帖 |
|
w******a 发帖数: 25 | 39 Here is an R example to impute one missing data in each record,half of the code is to make data sample, you probably only need second half,but including them here helps you understand what is going on:
The data will look like
col1 col2
x
x x
x
x x
x x
...
library(Rlab)
alp = 1
Prob_R1 = 0.5
Prob_R0 = 1 - Prob_R1
len_Y1 = 200
K_delta = 2
Y1 = rnorm(len_Y1,mean=0,sd=1)
R1 = rbinom(n=len_Y1, size=1, prob=Prob_R1)
Y2 = rnorm(n=len_Y1,... 阅读全帖 |
|
z**o 发帖数: 149 | 40 For data a1:
data a1;
input x$ y1 y2 y3 y4;
datalines;
1 2 1 2 3
2 4 2 5 8
3 6 9 7 8
4 8 4 6 5
;
want to check if y1-yx =1 or not:
for _n_=1, I want to check if y1=1;
for _n_=2, I want to check if y1-y2=1;
for _n_=3, I want to check if y1-y3=1;
is there a way to put a value of one variable into another variable name?
like: array y(*) y1-y&x.; ?
Thanks a lot for BIG NIU!! |
|
h*********g 发帖数: 18 | 41 I have 2 groups of data. Group1: x1,x2,x3,x4. Group2:y1,y2,y3,y4. I want to
get the ratio of y over x (Y/X) and also needs to know the SD or SEM of the
ratio. How to calculate the SD or SEM of Y/X? In my data, I can use y1/x1,
or y1/x2 or y1/x3 or y1/x4. x1 is not only related to y1, but can also
related y2 or y3 or y4, (what i mean is I can also use y2/x1 or y3/x1 or y4/
x1)
Please help.
Thanks! |
|
r*******t 发帖数: 8550 | 42 开始时,我说估计寄时会有x个货(手中有x-y,也告诉神医了)。但是去了几家,很多
家都卖光了,只拿到y1个。上星期二又跑了几家,拿到y2个,跟神医说了我又拿到y1+
y2个(星期三早上神医发email之前就通知他了)。
然后神医星期三早上说他"already paid",那我就drop le。现在才发现他星期三"
already paid"只cover x-y个货, not x-y + y1 + y2个货. 查的钱是神医y1 + y2个
货的钱,加上 x-y + y1 + y2 达到免shipping label的最低数。 |
|
k***s 发帖数: 206 | 43 X=
Y=
There's a function F(X,Y), and if you expand that function, you get
1+x1*y1 + x2*y2 + x1*x2*y1*y2
Find F(X,Y).
(Hint: dot product of X * Y is (x1*y1+x2*y2) )
这个我想 F(X,Y) 至少因该是 X*Y + 1, 但是这个 x1*x2*y1*y2 项怎么处理? |
|
l********n 发帖数: 9 | 44 存在30对随机变量(x1,y1),(x2,y2),......,(x30,y30),x1与y1,x2与y2,......,x30与y30
之间相关,但相关关系不完全相同.(x1,y1)的分布已知,同样其它29对随机变量的分布也已
知,按它们的分布,随机产生25组样本,每组样本中有100个个体,每个个体包含(x1,y1),(x2
,y2),......,(x30,y30)共30对数据.对每个个体,把y1,y2,......,y30相加得YS,但只把x1
,x3,x4,x7,x8,x13,x15,x17,x23,x25,x27,x30(即只把其中若干个相加)相加得XS,这样得
到
新的25组样本,每组样本中包含100对(XS,YS).对新的每组样本,按XS从小到大排序,取出排
在95位的XS及其对应的YS(不管YS的排序),这样得到25对(XS,YS).请问得到的25个YS服从
什
么分布?这里25组是随便定,如果需要可以取多于25组.要知道它的分布,是否还需要其它条
件?这涉及到统计中的什么问题?
望各位大侠赐教. |
|
r****y 发帖数: 1437 | 45 是比较难讲清楚,这么说吧
x = 0:1:10
y1 = x
y2 = exp(-x)
理想的figure是这样的
single curve (x, y1), not two curves
y1-left axis, label at 0, 1, 2,.., 10
y2-right axis, label at 1, 0.9, 0.8, ... (not same as y1)
The point is to show corresponding between y1 and y2 on the figure
in addtion to the curve.
This could be done by manipulating the properties of right axis,
e.g. set(ax(2), 'YTick', ****, 'YTickLabel', ****)
|
|
d****i 发帖数: 77 | 46 饶来饶去一晚上了,被饶的晕头转向
还请大虾们帮忙看看
X1>=X2>=...>=Xm >0
Y1>=Y2>=....>=Ym >0
X1+X2...+Xm = Y1+Y2+...Ym
and X1< 2*Y1
X1<2*Xm
Y1<2*Ym
Xm>Ym
prove: X1+...+Xk <= 2*(Y1+...Yk)
for any k<=m |
|
s*******a 发帖数: 921 | 47 可以用二次规划来解.将原问题转化成:
min f=y1^2+y2^2+y1+y2
s.t. x1<=y1
-x1<=y1
x2<=y2
-x2<=y2
y1>=0, y2>=0
matlab里面可以用quadprog来求解. |
|
A****s 发帖数: 129 | 48 yes..and health's proof is quite clear.
I just had another similar proof. for any x1, x2, and x3 as (x1+x2)/2
y3 as the point on the boundary of the set where d(x3,y3)=f(x3)
then I could draw two lines from (x1,y1) (x2,y2) on the 2-dimension plane
(x1,x2,y3) which parallel to (x3, y3), on the same side of (x1,x2),
and make their length as f(x1), f(x2). So y1 y2 are in this convex
set and (y1+y2)/2 is,too. Moreover, (y1+y2)/2 lies on (x3,y3) while y3 is on
the boundary. so d(x3,(y1+y2)/2) is smal |
|
m******s 发帖数: 1 | 49 来自主题: Mathematics版 - 问个小问题 现在假定有两个函数S(x1,x2)和E(y1,y2)
分别满足方程
\partial_{x1}\partial_{x2}S=0
和\partial_{y1}\partial_{y2}E=0.
现在给出一个statement 1:
\partial_{x1}\partial_{x2}S=0
<==> \partial_{y1}\partial_{y2}E=0.
现在假定S和E是通过线性变换联系起来:
S(x1,x2)=\int dy1 dy2 K(x1,x2,y1,y2) E(y1,y2)
问题是,要求statement 1 成立,对K的限制是什么?
或者,K的可能表达式是什么?
不是数学系的问题,不需要太严格,说个大概就可以了。谢谢!~! |
|
G******i 发帖数: 163 | 50 Let x0=0, y1=0,y2=1.
Consider U(x,y)=(1-x)f(y),
where f(0)=0, f(1)=1.
Your condition
a*U(x0,y2)+(1-a)*U(x0,y1)>U(x0,a*y2+(1-a)*y1) for 0
<=> f(a)
Your claim
a*U(x0+k*a*(1-a),y2)+(1-a)*U(x0+k*a*(1-a),y1)
>U(x0,a*y2+(1-a)*y1) for 0
<=> f(a)< a -k*a^2*(1-a) for 0
So, essentially, you are saying:
if f(0)=0, f(1)=1, f(a)
then there exists a small k0>0 such that
f(a)< a -k*a^2*(1-a) for 0
But this is incorrect, considering the case f(a)=a-a^3( |
|