similar to: Fwd: Re: extracting index list when using tapply()

Hello, The quick version of my question is how can I extract a matrix instead of a vector using tapply()? I would like to be able to access both the results of tapply() and also the index variables. In case further explanation would help: I am analyzing a large (3million rows x 9 columns) spatial/temporal dataset and am attempting to calculate the number of unique years containing any data

Hi R users Is there anything wrong in "is" function? (R 2.7.2) I believe that everyone will agree that "7" is an integer, right? but why R shows 7 is not an integer > is.integer(7) [1] FALSE > is(7,"integer") [1] FALSE > is(as.integer(7), "integer") [1] TRUE Thank you very much in advance Chunhao

Hi R users Is there anything wrong in "is" function? (R 2.7.2) I believe that everyone will agree that "7" is an integer, right? but why R shows 7 is not an integer > is.integer(7) [1] FALSE > is(7,"integer") [1] FALSE > is(as.integer(7), "integer") [1] TRUE Thank you very much in advance Chunhao

Hi R users, I have one simple question but I really don't know why I can't get it work. The data was adopted from Efron's An introduction to the bootstrap book. > nlaw LSAT GPA [1,] 576 3.39 [2,] 635 3.30 [3,] 558 2.81 [4,] 578 3.03 [5,] 666 3.44 [6,] 580 3.07 [7,] 555 3.00 [8,] 661 3.43 [9,] 651 3.36 [10,] 605 3.13 [11,] 653 3.12 [12,] 575 2.74

Hi R users, Is is possible for me to use the try function with boot? I would to do the bootstraping with a nonlinear model(it works well when R < 1000). But it does not work very well (when R is large) thus I try to use "try" to resolve. I put the try function in two cases: case1: put the try in front of the boot > c1.try<-try(boot(c1data, statistic = c1.fun,

First thanks for Jinsong's suggestions I would like to do a bootstrap in a nonlinear model. But it fails to converge in most of time. (it did converge if I just use nls without boot). Thus, I use "try" function to resolve my problem. This following code is from Jinsong's suggestion. h1a.nls<-nls(density~nmf(time, alpha, delta, psi, tau, gamma),data=h1a,

Hi Rers, I have a silly question. I don't know how to express the loglikelihood function of 1/(x!) where x=x1,x2,....xn in R. Could anyone give me a hint? Thank you in advance. Chunhao Tu

Hi R Users, I was accidentally save R file when I quit R program. I try to delete but I am not sure the R saved file was successfuly delete or not. Now everyt time when I open the R, it always appears [R.app GUI 1.34 (5589) i386-apple-darwin9.8.0] [Workspace restored from /Users/ctu/.RData] I have read google, R help and Nabble but I still can't solve this problem. Could someone please

Dear R Helpers, I try to write a small package that based on nlme however my code does not work. R always appears this message: Error in eval(expr, envir, enclos) : object "y" not found where y is the response variable. Please help me out! This is my code: require(nlme) AMPmixed<-function(y, x, S1=c("asymptotic","logistic"), random,data,

Thanks Daniel for the info. Siddharth On Thu, Sep 21, 2017 at 10:37 PM, Dániel Krupp <daniel.krupp at ericsson.com> wrote: > Hi Siddharth, > > > > since the CTU functionality is not yet part of the mainline clang (review > ongoing), you will need to clone our patched version of clang5 > > https://github.com/Ericsson/clang > > (ctu-clang5 branch) > > >

Hi R users, I google the website and I found that there are three ways to perform repeated measure ANOVA: aov, lme and lmer. http://www.mail-archive.com/r-help at stat.math.ethz.ch/msg58502.html But the questions are which one is good to use and how to do post-hoc test? I use the example that is provided in the above link and I try > tt<-aov(p.pa~group*time+Error(subject/time),data=P.PA)

Hello, CTU Student union - Silicon Hill club in cooperation with FEE CTU, Department of Computer Science, Prague, Czech rep., invites you (as a speaker) to 4th international and academic conference OpenWeekend 2004 taking place in Prague, October 16-17th 2004. It's topic is "open technologies and free software for broad range of users" with focus on use in academic environment. Its

I am having trouble getting an install of RCurl to work properly on a Unix server. The steps I have taken are: 1. installed cUrl from source without difficulty 2. installed RCurl from source using the command ~/R_HOME/R-devel/bin/R CMD INSTALL -l ~/R_HOME/R-devel/library ~/RCurl_0.8-0.tar.gz I received no errors during this install 3. when I go back to R and require(RCurl), I get >

I am having trouble getting an install of RCurl to work properly on a Unix server. The steps I have taken are: 1. installed cUrl from source without difficulty 2. installed RCurl from source using the command ~/R_HOME/R-devel/bin/R CMD INSTALL -l ~/R_HOME/R-devel/library ~/RCurl_0.8-0.tar.gz I received no errors during this install 3. when I go back to R and require(RCurl), I get >

R developers: I have some parallel C code, written with the MPI library, that I would like to call from R, but I get the error message below. cocite.mat.true.parallel.so compiles without complaint and I have MPI in my path, but R isn't recognizing one of the MPI symbols. I have a feeling that this is going to be tricky, but I need to get this to work. Helpful suggestions? I'd be

I am experiencing an OOLOCK error when attempting to update packages to my site-library. I have an identical setup on my local Linux machine, but this problem on the remote Linux machine (PPC cluster) started about 3 weeks ago and is present whether I used R-devel or R-2.6.0 or R-2.5.1. Also, of course I manually remove the 00LOCK folder after each failed attempt. Below is my error output,

Dear Josh, Thanks for your help! Does your answer mean, that you agree the two methods should do the same, and what I was guessing, despite the small differences? What I prefer about ci.pd is, that the help clearly says which method is implemented, which is not the case for prop.test. But I do not know who has programmed the function. Best wishes Steffi Stefanie von Felten, PhD Statistician

Hi, Works fine with Win2K or XP, but when I go to a Win98 system. "You must supply a password to make this connection: Resource: \\bigred\ipc$ password: i put in all the password, but none of them work. Any advice would be greatly appreciated. -Paul Hong

Dear R users, I have a question how to use 2 "ifelse" to sort my data. Such as from 11 to 20 assign to A; 6 to 10 assign to B, and the rest of them assign to C a<-1:20 tt<-ifelse(a>10, "A",no=ifelse( 5< a <=10, "B", "C")) Many Thanks Chunhao -- View this message in context:

Hi R users, I have a question. How can I use for loop to do pair comparisons. For example, > x<-c(1,2,3) > result<-matrix(data=NA, nrow=choose(3,2), ncol=1) > for(i in 1: length(x)) + { result[i,]<-ifelse(x[i] > x[i+1], yes="Big", no="Small") + result} > result [,1] [1,] "Small" [2,] "Small" [3,] NA