daiwenyan

  •  
  • 2010年1月14日
  • 注册于 2009年1月3日
  • 恩。呵呵,把代码修改了,给忘了。谢谢。
  • 师兄用了4个月就全部学会了,呵呵,还是比较容易上手的。
  • 呵呵,虽然看不到,但有这么多人关注。
  • 不知道怎么改了一下,就使得打开R就加载一些包。



    载入需要的程辑包:tcltk

    载入Tcl/Tk接口... 完成

    载入需要的程辑包:Hmisc



    载入程辑包:'Hmisc'





            The following object(s) are masked from package:base :



             format.pval,

             round.POSIXt,

             trunc.POSIXt,

             units



    载入需要的程辑包:R2HTML

    载入需要的程辑包:svMisc



    怎么设置能想以前一样自己设置。
  • R中如何做图形的叠加(用的是同样的坐标,levelplot想和points做出来的图放到一张图上)
  • [quote]引用第13楼大树2009-01-03 11:23发表的“”:

    错误于plot.xy(xy.coords(x, y), type = type, ...) : 还没有调用plot.new[/quote]

    我也在调,还没有解决。我觉得是因为循环中连续使用points的原因,第一次可以执行,以后每次都没有plot,而直接调用points的原因。
  • 仿照师兄两次的灌水帖整理:
    <br />
    library(tuneR)<br />
    d = 20000<br />
    do=sine(512,bit = 16, duration = d)<br />
    re=sine(576,bit = 16, duration = d)<br />
    mi=sine(640,bit = 16, duration = d)<br />
    fa=sine(683,bit = 16, duration = d)<br />
    so=sine(768,bit = 16, duration = d)<br />
    la=sine(852,bit = 16, duration = d)<br />
    si=sine(960,bit = 16, duration = d)<br />
    dos=sine(512,bit = 16, duration = d/2)<br />
    res=sine(576,bit = 16, duration = d/2)<br />
    mis=sine(640,bit = 16, duration = d/2)<br />
    fas=sine(683,bit = 16, duration = d/2)<br />
    sos=sine(768,bit = 16, duration = d/2)<br />
    las=sine(852,bit = 16, duration = d/2)<br />
    dol=sine(512,bit = 16, duration = d*2)<br />
    sol=sine(768,bit = 16, duration = d*2)<br />
    music = bind(do,re,mi,fa,so,so,sos,fas,mi,fa,fa,fas,mis,re,do,mi,sol,do,re,mi,fa,so,so,sos,fas,mi,fa,fa,fas,mis,re,do,mi,dol,<br />
            la,la,las,sos,fa,so,so,sos,fas,mi,fa,fa,fas,mis,re,do,mi,sol,la,la,las,sos,fa,so,so,sos,fas,mi,fa,fa,fas,mis,re,do,mi,dol)<br />
        play(music)<br />
        y = c(11, 8, 2, 11, 10, 8, 5, 3, 11, 9, 8, 5, 1, 12, <br />
                11, 8, 7, 6, 5, 4, 3, 2, 1, 11, 9, 8, 5, 11, 10, 8, 5, 3, <br />
                1, 8, 2, 11, 10, 9, 8, 7, 6, 5, 4, 3, 11, 8, 11, 8, 7, 6, <br />
                5, 4, 3, 2, 1, 12, 8, 8, 50, 8, 4, 9, 4, 12, 11, 10, 7, 6, <br />
                5, 4, 10, 7, 4, 10, 7, 4, 10, 7, 4, 10, 9, 8, 7, 6, 5, 4, <br />
                3, 2, 1, 10, 7, 4, 10, 7, 4, 10, 7, 4, 10, 7, 4, 10, 4, 50, <br />
                9, 8, 7, 50, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 9, 1, <br />
                8, 6, 2, 10, 6, 3, 10, 6, 4, 12, 11, 10, 9, 8, 7, 6, 5, 10, <br />
                6, 4, 10, 6, 3, 10, 9, 8, 7, 6, 2, 6, 1, 50, 50, 7, 2, 11, <br />
                10, 9, 8, 7, 3, 11, 7, 5, 4, 11, 7, 11, 7, 1, 11, 9, 8, 7, <br />
                6, 5, 4, 3, 2, 1, 12, 7, 12, 7, 5, 12, 7, 4, 7, 3, 2, 7, <br />
                50, 50, 50, 50, 10, 9, 8, 7, 6, 5, 4, 2, 50, 50, 50, 50)<br />
        r = c(3, 5, 5, 10, 4, 6, 2, 9, 2, 9, 2, 1, 1, 2, 2, <br />
                7, 3, 3, 3, 10, 3, 3, 3, 3, 2, 1, 3, 1, 12, 2, 3, 3, 3, 8, <br />
                3, 3, 6, 2, 1, 1, 2, 6, 4, 2, 3, 10, 2, 3, 3, 3, 1, 1, 1, <br />
                1, 1, 8, 1, 1, 1, 1)<br />
        x = rep(1:60, r)<br />
        plot(x, y,type="n",xlim = c(0, 60), ylim = c(-20, 30))<br />
        for (i in 1:10000) {<br />
            points(x, y, col=rainbow(length(x))[1:length(x)])<br />
            Sys.sleep(0.2)<br />
            points(x, y,col="white")<br />
            x = ifelse(x == 1, x * 60, x - 1)<br />
        }    <br />
    


    祝大家新年快乐。
  • 果然啊!哈哈,确实强!只是N能怎么能变成大写的啊?
  • 我试了一下,确实不行,期待高手。
    <br />
    > paste("//N", "a")<br />
    [1] "//N a"<br />
    > paste("/N", "a")<br />
    [1] "/N a"<br />
    > paste("\N", "a")<br />
    [1] "N a"<br />
    Warning messages:<br />
    1: 不认识字符串里'\N'这样的逸出号 <br />
    2: 从"\N"里删除了不认识的逸出号 <br />
    > paste("\\N", "a")<br />
    [1] "\\N a"<br />
    
  • > z<- read.csv("a.txt",header=FALSE)

    > z

    V1 V2 V3

    1 我 是 辉

    2 我 是 辉

    >



    我不明白作者的意思。有帮助吗?
  • > seedrates.lm1<-lm(grain~rate,data=seedrates)

    > seedrates.lm2<-lm(grain~rate+I(rate^2),data=seedrates)

    > anova(seedrates.lm2,seedrates.lm1)

    Analysis of Variance Table

    Model 1: grain ~ rate + I(rate^2)

    Model 2: grain ~ rate

    Res.Df   Res.Sum     Sq Df     Sum Sq     F value       Pr(>F)

    1 2       0.026286

    2 3         0.187000     -1   -0.160714     12.228     0.07294

    什么意思没明白,是将两次拟合结果做对比,看有无显著性差异吗?还是什么别的?