Wallee

Wallee · 2015年3月5日

老师我远程登录的时候显示 unable to connect to service . 登录 unbuntu 查看 rstudio-server verify-installation 一切正常：

rstudio-server verify-installation

rstudio-server stop/waiting

rstudio-server start/running, process 28470

然后选择rstudio-server active-sessions

rstudio-server active-sessions

PID TIME COMMAND

什么session都没有 rstudio 重启还是打不开新的rsession 也不报错这个怎么解决啊已经关闭 apparmor了

Wallee · 2014年12月19日

没有配置R到系统变量……

Wallee · 2014年12月12日

如何使用提取一个数组中每个元素和之后10个元素之和组成新元素。

例如 a 为1:100 的数组有 100个元素为1,2,3,4,5,6...100

需要生成数组b b 的每一个元素是 a 中每一个元素和他后面10个数字之和

b 为 (1+2+3...+10),(2+3+4+5...+11),...,(91+92+93+...100) 这样一个1*100的数组。
for循环可以解决能否用lapply 解决呢？因为没办法确定带入FUN的a 元素的下标啊如果用which的话其中如果出现了两个相同元素怎么搞定呢？

Wallee · 2014年9月28日

Java 调用R的挖掘结果问题。

R：代码:

a<-rpart(Species~.,iris)"

a

显示结果：

n= 150

node), split, n, loss, yval, (yprob)

* denotes terminal node

1) root 150 100 setosa (0.33333333 0.33333333 0.33333333)

2) Petal.Length< 2.45 50 0 setosa (1.00000000 0.00000000 0.00000000) *

3) Petal.Length>=2.45 100 50 versicolor (0.00000000 0.50000000 0.50000000)

6) Petal.Width< 1.75 54 5 versicolor (0.00000000 0.90740741 0.09259259) *

7) Petal.Width>=1.75 46 1 virginica (0.00000000 0.02173913 0.97826087) *

java调用这个R 用这个包代码：

r.eval("library(rpart)");

r.eval("a<-rpart(Species~.,iris)");

REXP trees=r.eval("a");

REXP no1=r.eval("labels(a)");

System.out.println(trees);

System.out.println(no1);

结果：

[VECTOR ([VECTOR ([FACTOR {levels=("<leaf>","Petal.Length","Petal.Width"),ids=(1,0,2,0,0)}], [INT* (150, 50, 100, 54, 46)], [REAL* (150.0, 50.0, 100.0, 54.0, 46.0)], [REAL* (100.0, 0.0, 50.0, 5.0, 1.0)], [REAL* (1.0, 1.0, 2.0, 2.0, 3.0)], [REAL* (0.5, 0.01, 0.44, 0.0, 0.01)], [INT* (3, 0, 3, 0, 0)], [INT* (3, 0, 3, 0, 0)], [REAL* (1.0, 1.0, 2.0, 2.0, 3.0, 50.0, 50.0, 0.0, 0.0, 0.0, 50.0, 0.0, 50.0, 49.0, 1.0, 50.0, 0.0, 50.0, 5.0, 45.0, 0.3333333333333333, 1.0, 0.0, 0.0, 0.0, 0.3333333333333333, 0.0, 0.5, 0.9074074074074074, 0.021739130434782608, 0.3333333333333333, 0.0, 0.5, 0.0925925925925926, 0.9782608695652174, 1.0, 0.3333333333333333, 0.6666666666666666, 0.36, 0.30666666666666664)])], [INT* (2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, ... (50 more values follow))], [NULL ], [NULL ], [REAL* (0.5, 0.44, 0.01, 0.0, 1.0, 2.0, 1.0, 0.5, 0.06, 1.21, 0.8200000000000001, 0.11, 0.04836665518033404, 0.060969937947592937, 0.03192700006786314)], [STRING "class"], [VECTOR ([REAL* (0.3333333333333333, 0.3333333333333333, 0.3333333333333333)], [REAL* (0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 1.0, 1.0, 0.0)], [REAL* (1.0)])], [VECTOR ([INT* (20)], [REAL* (7.0)], [REAL* (0.01)], [INT* (4)], [INT* (5)], [INT* (2)], [INT* (0)], [INT* (30)], [INT* (10)])], [VECTOR ([NULL ], [NULL ], [NULL ])], [INT* (5)], [REAL* (150.0, 150.0, 150.0, 150.0, 0.0, 0.0, 0.0, 100.0, 100.0, 100.0, 100.0, 0.0, 0.0, 0.0, -1.0, -1.0, -1.0, 1.0, -1.0, -1.0, 1.0, -1.0, -1.0, -1.0, -1.0, -1.0, -1.0, -1.0, 50.0, 50.0, 34.16405023547881, 19.038507534082754, 1.0, 0.92, 0.8333333333333334, 38.969404186795494, 37.35353535353535, 10.686868686868685, 3.555555555555557, 0.91, 0.73, 0.67, 2.45, 0.8, 5.45, 3.3499999999999996, 0.8, 5.45, 3.3499999999999996, 1.75, 4.75, 6.15, 2.45, 4.75, 6.15, 2.95, 0.0, 0.0, 0.0, 0.0, 1.0, 0.76, 0.5, 0.0, 0.0, 0.0, 0.0, 0.8043478260869565, 0.41304347826086957, 0.2826086956521739)], [REAL* (88.96940418679549, 81.3449555415529, 54.0960582510677, 36.013092487572635)], [INT* (1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, ... (50 more values follow))], [BOOLi* ])]

我知道结果不一样的原因是 a 在R里面存着是以对象方式存储他的构成为：

names(a)

[1] "frame" "where" "call"

[4] "terms" "cptable" "method"

[7] "parms" "control" "functions"

[10] "numresp" "splits" "variable.importance"

[13] "y" "ordered"

由这些结果构成a在R中的存储。但是现在问题时我如何将a 在R的显示结果就是这个“树形图“

> a

n= 150

node), split, n, loss, yval, (yprob)

* denotes terminal node

1) root 150 100 setosa (0.33333333 0.33333333 0.33333333)

2) Petal.Length< 2.45 50 0 setosa (1.00000000 0.00000000 0.00000000) *

3) Petal.Length>=2.45 100 50 versicolor (0.00000000 0.50000000 0.50000000)

6) Petal.Width< 1.75 54 5 versicolor (0.00000000 0.90740741 0.09259259) *

7) Petal.Width>=1.75 46 1 virginica (0.00000000 0.02173913 0.97826087) *

传递到Java当中呢？实在找不到办法。使用print（a）传递过来，还是一样，不变。

Wallee · 2014年9月19日

我发现了一个一个办法…………

自己发现了

http://cos.name/cn/topic/102766

Wallee · 2014年9月19日

数据结构和往常的不太一样是这样子的：

A B quantity

1 339 1

1 308 2

3 339 1

5 225 1

其中A列是订单编号 B列是产品编号， quantity 是B列产品编号的购买数量。想进行关联分析就是选择了339的人有多大概率同一时间选择225 这个是典型的关联分析问题。但是问题出在数据格式上。

标准的关联分析需要把数据变成对象。同一个订单的物品信息放在一起格式应该是：

{339,308}

{339}

{225}

但是现在按照订单编号，物品编号单独列出来了。应该怎么用这个函数呢？

谢谢大神了

发现一个解决办法正在尝试:

http://cos.name/cn/topic/102766

Wallee · 2013年2月22日

<br />
> summary(fit2011)</p>
<p>Call:<br />
lm(formula = y ~ ns(x, 7))</p>
<p>Residuals:<br />
ALL 8 residuals are 0: no residual degrees of freedom!</p>
<p>Coefficients:<br />
              Estimate Std. Error t value Pr(>|t|)<br />
(Intercept) -3.925e-17         NA      NA       NA<br />
ns(x, 7)1    1.509e-01         NA      NA       NA<br />
ns(x, 7)2    4.656e-01         NA      NA       NA<br />
ns(x, 7)3    5.942e-01         NA      NA       NA<br />
ns(x, 7)4    8.618e-01         NA      NA       NA<br />
ns(x, 7)5    8.878e-01         NA      NA       NA<br />
ns(x, 7)6    1.031e+00         NA      NA       NA<br />
ns(x, 7)7    9.608e-01         NA      NA       NA</p>
<p>Residual standard error: NaN on 0 degrees of freedom<br />
Multiple R-squared:     1,      Adjusted R-squared:   NaN<br />
F-statistic:   NaN on 7 and 0 DF,  p-value: NA<br />

</p>

是指的是高阶多项式吗？谢谢啦

然后如何去预测值呢？用什么方法？就是说知道x 求出y

Wallee · 2013年2月21日

例如一串数组，一共 1万个数

我要统计最大10% 的均值，次大10%均值。中间20% 均值这样子有嘛命令么有呢？

Wallee · 2013年2月20日

回复第6楼的 superabe：我在5楼发现了点问题你看下~ 谢谢啦~~~~~

Wallee · 2013年2月20日

回复第4楼的 superabe：嗯我现在试了下这么写但是后来在算的时候的出来的数组是 “NULL” 疯了

</p>
<p>taxPerMonNew<-function(sinRawData){<br />
y<-sinRawData<br />
if (y<=3500) {0<br />
}else{<br />
  y<-sinRawData-3500}</p>
<p>if (y>0 && y<=1500)  y*0.03<br />
if (1500<y & y<=4500)  y*0.1-105<br />
if (4500<y & y<=9000) y*0.2-555<br />
if (9000<y & y<=35000) y*0.25-1005<br />
if (35000<y & y<=55000) y*0.3-2755<br />
if (55000<y & y<=80000)  y*0.35-5505<br />
if (80000<y)  y*0.45-13505}</p>
<p>taxFunPerMonNew<-function(rawData,taxFun){<br />
   lapply(rawData,taxFun)</p>
<p>}</p>
<p>taxDataPerMonNew2011<-taxFunPerMonNew(rawData2011,taxPerMonNew)</p>
<p>

其中 rawData2011 是很多0到10000之间的数字。应用玩方程后 taxDataPerMonNew2011 里面的数据全部都是NULL 咋回事呢？
</p>

Wallee · 2013年2月20日

回复第2楼的 superabe：能不能用 lapply呢？

Wallee · 2013年2月20日

扣除额为3500元

税率表：y=x-3500（x为扣除三险一金后的月收入）

IF 0<y≤1500，应缴纳个人所得税=y*3%；

IF 1500<y≤4500，应缴纳个人所得税=y*10%-105；

IF 4500<y≤9000，应缴纳个人所得税=y*20%-555；

IF 9000<y≤35000，应缴纳个人所得税=y*25%-1005；

IF 35000<y≤55000，应缴纳个人所得税=y*30%-2755；

IF 55000<y≤80000，应缴纳个人所得税=y*35%-5505；

IF 80000<y，应缴纳个人所得税=y*45%-13505；

rawData 是一个个人收入样本大概2万个左右我要计算样本中的每个数字的所得税

y就是所得税

然后存储到一个数组中。

y=rawData-3500

if ( 0<y≤1500)

(y*3%)

else if (1500<y≤4500)

(y*10%-105)

else if (4500<y≤9000)

(y*20%-555)

else if (9000<y≤35000)

(y*25%-1005)

else if (35000<y≤55000)

(y*30%-2755)

else if (55000<y≤80000)

(y*35%-5505)

else (80000<y)

(y*45%-13505)

我的语句首先总是报告 else if 不对

其次报告 if 不能应用在向量

Wallee · 2013年2月20日

厉害~回复第2楼的 ypchen：厉害~~好厉害~~~~~~

Wallee · 2013年2月20日

a<-runif(10000,0,1)

将a 中的每一个值代入公式：

y=0.7263x^2 +0.9695x+6.2322

如何写命令啊~~~我疯了~~~~

菜鸟一枚急求助