similar to: Introduction of NA:s

Hi R community, These days, I am writing some functions to work with 2-ways frequency tables ; you know all this tuff about measures of association: Chisq and derived (phi, cramer's v), tau b, tau c, somer's d and so on. So I consider all those functions could be gathered in a single R file, as they are coherent and dealing with the same problem / objective (analysis of crosstables).

Hello, I wonder if someone could send me suggestions on how to solve the following problem: I have a vector of an arbitrary size (ex. data<-c(10,10,11,10,12,11,10,12,11,11,10,11)) and use the table function, which gives the following result 10 11 12 5 5 2 that''s fine, but what I would like to do now is: construct new classes based on the number of classes from table, 10

Dear R Help, I am trying to create a boolean vector that is TRUE whenever a particular value occurs in a numeric vector, and FALSE otherwise. For example, suppose that > y <- c(5, 2, 4, 3, 1) > y [1] 5 2 4 3 1 and suppose that I want to find where 3 occurs in y. Then, the following yields the solution: > y == 3 [1] FALSE FALSE FALSE TRUE FALSE My problem arises when the

Hello again, I was too quick before. What I was looking for was a function that constructs the design (or incidence) matrix (X in a linear model) from a factor. Uwe Ligges suggested using model.matrix and this does almost what I want, but it is first necessary to construct a data variable. It also asigns ones to all rows of the first column (because this is set to be the contrast, not really what

Hello, are there any good methods in R that will test if two correlation matrices (obtained in different ways) are equal? Something better than the Mantel test would be preferable. Regards, Patrik Waldmann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I wonder if someone could help me adjusting the following code to parallelized snow code: #Creating a data set (not needed to be parallel) n<-100 p<-100 x<-matrix(rnorm(n*p),p) y<-rnorm(n) # Bootstrapping nboot<-1000 alpha<-0.05 rhoboot <- array(0, dim=c(p,nboot)) bootranks <- array(0, dim=c(p,nboot)) bootsamples <- array( floor(runif(n*nboot)*n+1), dim=c(n,nboot)) for

Hello, are there any functions in R for estimation of the mode of a MCMC-chain? Best, Patrik Waldmann [[alternate HTML version deleted]]

Hello, I have an unbalanced data set with some factors and several variables. The missing values are spread unequally over the variables. Because of the unbalanced structure have i tried NLME, but I cannot understand how I should set up the na.omit. I would not like to use na.omit before setting up the model, as that would exclude too many observations (on the whole structure). Lets say I have

I haven't found a function that directly calculates the matrix inverse, does it exist? Otherwise what would be the fastest way to do it? Patrik Waldmann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Hello, I would like to try Tcl/Tk library on Windows. I installed ActiveTCL 8.3.4.3 and add the environment variable TCL_LIBRARY linking to the directory of installation. When I load the Tcl/Tk package, I have the following error: +--------------------------------------------------------------+ Error in firstlib(which.lib.loc, package) : Can't find a usable init.tcl in the

Does this help on our backtest reporting? Regards, Gordon Morrison Global Head of Quantitative Research > * + 44 20 7653 7642 > Mob: + 44 7867 801951 > fax: + 44 20 7645 7442 > * mailto:gordon.morrison at commerzbankib.com > web: http://www.cbksec.com/research/quant > * Commerzbank Securities > 60 Gracechurch Street > London EC3V 0HR, U.K. >

Hello R communauty, I have submitted to CRAN the package R2HTML. The purpose of this package is twice: - provide base HTML exportation functions. Basically, you will use them exactly as a print function (or a cat(...file=)). Thus, making dynamic HTML reports is possible. - making live HTML output with frames and command, which is very usefull for teaching purpose, as the student can keep a

Hi (again), Please find here a correction of the function I send previously --------------------------------------------------------- tableChisqMH=function(x) { n=sum(x) G2=(n-1)*(tablepearson(x)$estimate^2) dll=1 p.value=1-pchisq(G2,dll) out=list(estimate=G2,dll=dll,p.value=p.value,dim=dim(x),name="Mantel-Haenszel Chi-square") return(out) }

Hello UseRs: Excuses for my english. I have a dataset with 65000 records and I'd like to make a summary where I can view all the values (with the number of times that it repeats) that there are each column of my dataset. I tried with summary( ), str( ), but nothing gives me the result that I am loking for. Thank you very much.

Hello again, let say we have following data: Dat1 <- structure(list(factor.sample.LETTERS.1.3...6..replace...T.. = structure(c(1L, 3L, 2L, 1L, 3L, 3L), .Label = c("A", "B", "C"), class = "factor"), factor.sample.letters.1.2...6..replace...T.. = structure(c(2L, 2L, 1L, 1L, 2L, 1L), .Label = c("a", "b"), class =

Hi R-listers, I am trying to better undertand what we would call "functional paradigm" use of S/R to better map my programming activities in other languages. This little function is aimed to create an object (at the end end, it would have it's own class): -- myObject =function(){ list( a=1, foo=function(b) { cat("b:",b)

Dear listers, I'm having a bad R day. I just can't think of the vectorized equivalent of: for (ii in 1:n) aa[ii] = bb[ii,cc[ii]] Any suggestion received with embarrassment and gratitude Simon Gatehouse CSIRO Exploration and Mining, Newbigin Close off Julius Ave North Ryde, NSW Mail: PO Box 136, North Ryde NSW 1670, Australia

https://bugzilla.mindrot.org/show_bug.cgi?id=2331 Bug ID: 2331 Summary: ssh-copy-id -i id_new.pub fails to copy that id, also: wrong error msg Product: Portable OpenSSH Version: -current Hardware: All OS: Linux Status: NEW Severity: major Priority: P5 Component:

Hi, is there a way to reference to a data object without copying it? For example I have a huge matrix called dist and I want two objects obj1 and obj2 to have a memeber dist that points to the matrix, but I don't want, for memory reasons, to copy the matrix twice. As far as I understand the following code will generate three copies of my data: dist <- some_code_that_generates_data

Hi all, Is it possible to add a permil (or per mille) symbol to an R plot (I couldn't find this symbol under demo(Hershey) or the plotmath information). In some ascii tables it is symbol no. 137. cheers, Matt. [[alternative HTML version deleted]]