similar to: how to get accumulate contingency table

Hi, I post the question again. x<-rep(c(2,3,4),20) y<-sample(rep(c(0,1),30)) tr<-rep(c(1:5),6) data1<-data.frame(x,y,tr) data1<-data1[order(data1$tr),] > data1 x y tr 2 1 1 4 0 1 3 1 1 2 1 1 4 1 1 3 1 1 2 0 1 4 0 1 3 0 1 2 0 1 4 0 1 3 1 1 3 0 2 2 0 2 >table(data1) tr = 1 y x 0 1 2 2 2 3 1 3 4 3 1 tr = 2 y x 0 1 2 2 2 3 2 2 4 1 3

Hi all, Suppose: y<-rnorm(100) x1<-rnorm(100) lm.yx<-lm(y~x1) To predict from a new data source, one can use: # works as expected dum<-data.frame(x1=rnorm(200)) predict(lm.yx, newdata=dum) Suppose lm.yx has been run and we have the lm object. And we have a dataframe that has columns that don't correspond by name to the original regressors. I very! naively assumed that doing

Hi, I found the following on Windows 2000/NT R Version 1.9.1 (2004-06-21) (also Version 1.9.0): The S4 group "Math" doesn't work as documented; i.e., "log", "log10", "gamma" and "lgamma" are included in the documentation but don't work. See example code below. Moreover, what about 'genericForPrimitive' which is used in

Hi, For further calculations I need to extract indivdual Variances of different random effects from a fitted model. I found out how to extract the correlations (VarCorr(m1)@reSumry$group1) but I was not able to find a way to extract the other components individually. To extract the Residuals I tried: (ranef(m1)@ stdErr) which unfortunately did not work. Thank you very much for your help!

Colleagues ---------------------------------- System info: R version rw1020 on NT ESS using emacs ver. 20.4 ---------------------------------- I need to get the means for each column of a dataframe in the list created by splitting a data frame. At present, I am getting the mean of all columns in aggregate. The structure of the unsplit data is: > shuttle.tr1[1:10,] juliandate

On Thu, Mar 15, 2018 at 9:09 PM, David Blaikie via llvm-dev <llvm-dev at lists.llvm.org> wrote: > +Chandler who might have some thoughts on this. > > Could you provide an example here of the motivation for the feature you're > missing? Might help motivate the discussion (and/or we'll end up nitpicking > how it could be done differently without that feature... - which

Thanks. The motivating example can be seen in this review: https://reviews.llvm.org/D44382. In that review, I am unit testing .debug_line parsing, specifically, the behaviour when the parser is fed a malformed section. Most of the code under test goes through some slight variations in the code path, depending on a) the DWARF version (interesting cases are 3, 4 and 5), and b) whether the DWARF is

I've been doing some experiments with the Boost Graph Library (BGL) and interfacing to llvm's GraphWriter system. I'm using GBL for some custom work and needed a way to hook BGL's write_graphviz into GraphWriter to use the existing display infrastructure. To do this I added a HAVE_BOOST configure item and #ifdef'd the Boost support code in GraphWriter. Is this useful to

Hello All! I am trying to plot the frequency of species coocurrance. If given a data set similar like this...(V1="species A", V2="species B", V3="frequency of cooccurance") > data V1 V2 V3 1 A B 0 2 A C 2 3 A D 5 4 B C 0 5 B D 1 6 C D 0 > data1<-as.data.frame(lapply(data,function(x)(rep(x,data$V3)))) > as.data.frame(data1[-1]) >

Got problems with plotCI (plotrix) I only want to plot the upper part of the error bar in my barplot I had the exact same commands working two months ago Now I wanted to change the legend and when I ran it again plotCI stopped working. To me it seems like there must some bug in R library(plotrix) library (graphics) x = matrix(c( 13.75516276, 3.944604404,

The following example, extracted from pp. 38-42 of Chambers (1998) Programming with Data (Springer), works for me in S-Plus 6.2 but not R 1.9.0pat, which returns the indicated error messages: > setClass('track', representation(x="numeric", y="numeric")) [1] "track" > tr1 <- new("track", x=1:3, y=4:6) > setMethod("plot", +

Hello all, I have an ISDN termination box (TR1) that converts ISDN(Bri) to 2 normal analogue lines. The same number is assigned to these lines. These lines are connected to 2 spa3k registered to my asterisk box. When calls arrive, TR1 try to pass call to the first spa. If spa not takes the call immediately then try to pass to the other spa. The only configuration I found works is to put the

In C++ code for use in a R-3.1.0 package, my specific problem is that I would like to use <unordered_map> if it is available, or <tr1/unordered_map> if not, or <map> if all else fails. I (think I) can accomplish this with configure.ac as AC_INIT("DESCRIPTION") CXX=`"${R_HOME}/bin/R" CMD config CXX` CXXFLAGS=`"${R_HOME}/bin/R" CMD config

2008/12/27 Mark Kromis <greybird at mac.com> > Is this something your planning as putting in the tree, >> thus require pulling in changes from google (license allowing), or does user >> need to have the libraries/headers pre-installed? >> > > Including it in the tree is the most reasonable thing to do. No point in > inconveniencing the user over tiny libraries

Greetings R Users, I have a system of equations for which I would like to output all the R-squares. Assume there are four equations in my system, the only way I found to output all the R-squares is by calling them out one by one as this: summary(fitSUR$eq[[1]])$r.squared summary(fitSUR$eq[[2]])$r.squared summary(fitSUR$eq[[3]])$r.squared summary(fitSUR$eq[[4]])$r.squared But isn't there a

On Dec 27, 2008, at 11:01 PM, Misha Brukman wrote: > 2008/12/27 Mark Kromis <greybird at mac.com> >> Is this something your planning as putting in the tree, thus >> require pulling in changes from google (license allowing), or does >> user need to have the libraries/headers pre-installed? >> >> Including it in the tree is the most reasonable thing to do.

Hi all, I asked something like this earlier but decided that a proper minimal example might be helpfull ;0) Why does this work with regards to the expression (substitution): require(stats) plot(cars) text(5,120,labels=substitute(i^{z+phantom()}*"("*a*" AMU)",list(i="yx",z=2,a=0))) text(c(5,5),c(115,110),labels=c("One","Two")) But adding this

This shortens runtime of piglit test fp-long-alu to ~22s No piglit regressions observed on nvc0! Signed-off-by: Tobias Klausmann <tobias.johannes.klausmann at mni.thm.de> --- src/gallium/drivers/nouveau/codegen/nv50_ir.cpp | 6 +++--- src/gallium/drivers/nouveau/codegen/nv50_ir.h | 7 ++++--- src/gallium/drivers/nouveau/codegen/nv50_ir_peephole.cpp | 2 +-

I am using the tm package to do text miniing: I have a huge list of stopwords (2000+) that are in a csv file. I read it as follows: stopwordlist <- read.csv("stopwords to be Removed 10042011.csv") myStopwords <- as.character(stopwordlist$stopwords) When try removing the stopwords using tr1=tm_map(tr1,removeWords,myStopwords) I am getting the following error: Error in

Dear all, Can anyone help me, my R software can not run a nested linear regression by using the lme funcion. The message that appears is  Error: could not find function "lme" I already downloaded and loaded the package, please see below. Thank you in advance for any help! Nadia. > data<-read.csv("/Users/nadiasan1/Desktop/MOE and MOR.csv")> attach(data)>