similar to: summary stats on continuous data

I have a list of data frames like this: a<- data.frame(x=runif(10), y = runif(10), Acc = 1) b<- data.frame(x=runif(10), y = runif(10), Acc = 0) ls<- list(a,b) and I want to remove the data frames from ls that have Acc values other than 1. How do I do that? Thanks for any help! Matthew

I wonder if someone could explain the behavior of the anova() and lm() functions in the following situation: I have a standard 3x2 factorial design, factorA has 3 levels, factorB has 2 levels, they are fully crossed. I have a dependent variable DV. Of course I can do the following to get the usual anova table: > anova(lm(DV~factorA+factorB+factorA:factorB)) Analysis of Variance Table

Dear All, I'm struggling a little with the behaviour of R with GLM interactions. In particular, I have a dataset with two factors - call them factor A and factor B, where I would like to fit a GLM that is factor A + (grouped factor A):factor B. To try to isolate this, I've ignored the original "factor A" part, as that I have this as a separate column in my data. So, it

Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct. Consider an analysis consisting of two factors A and B. When all factor level combinations are present all appears fine: library(sciplot) #all data

How to generate a new factor variable by two other factor variables? For example, if I have two factor variables, factorA and factorB, factorA factorB 0 0 0 0 1 0 0 1 1 1 Is there a simple way to generate a new 4-levels factor variable as factorC factorA factorB 0 0 0 0 0 0 1 1 0 2 0 1

Hi, I have a lattice xyplot that contains panels according to FactorA, and curves for the 2 levels of Factor B within a panel. I try to add text in the panels of a lattice graph. I suppose I have to write a custom function (panel.txt). What I really would like is to adapt the text in the panel according to the levels of FactorA. In the manuals, I find examples for the strips using which.given

Hello, I have the following data.frame structure(list(Study = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,

Dear list, i have a quick and (hopefully) straightforward question regarding the plot-function after running aov. if i plot an equation like this: plot(dataSubjects~factorA, data=mydata) R gives me the boxplots for this particular factor A. my model, however contains several factors. is there a straightforward way to plot barplots for a specific factor with the constraint that those values

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Hi, I'm doing robust regression with the following command rlm(dip~ind1+ind2-1,method="M",psi=psi,maxit=1000,acc=1e-15) now when I ask for a summary summary(rlm(dip~ind1+ind2-1,method="M",psi=psi,maxit=1000,acc=1e-15)) I get Coefficients: Value Std. Error t value ind1 -0.0377 0.0000 -24203.1415 ind2 1.0370 0.0000 668735.7195 taht is

Dear expeRts, I am trying to wrap up a package "utilities" (for my internal use). After adding a function datNAtreat that uses mapply, R CMD check gives WARNINGs for "S3 generic/method consistency", "checking replacement functions" and?"checking foreign function calls", all of which are accompanied by the following error message: Error in .try_quietly

Dear all, A question on modelling proportional data in R. I have a test experiment that was designed in a particular way, and which I can analyse "by hand" to an extent. I am really struggling to get R to give me sensible results in modelling it "properly", so must be doing something wrong here. As background, I conduct a series of experiments and count the

Hi guys, I am working on some code for automatically translating recursive functions into looping functions to implemented tail-recursion optimisations. See https://github.com/mailund/tailr As a toy-example, consider the factorial function factorial <- function(n, acc = 1) { if (n <= 1) acc else factorial(n - 1, acc * n) } I can automatically translate this into the loop-version

Interestingly, the <<- operator is also a lot faster than using a namespace explicitly, and only slightly slower than using <- with local variables, see below. But, surely, both must at some point insert values in a given environment ? either the local one, for <-, or an enclosing one, for <<- ? so I guess I am asking if there is a more low-level assignment operation I can get my

Following up on this attempt of implementing the tail-recursion optimisation ? now that I?ve finally had the chance to look at it again ? I find that non-local return implemented with callCC doesn?t actually incur much overhead once I do it more sensibly. I haven?t found a good way to handle parallel assignments that isn?t vastly slower than simply introducing extra variables, so I am going with

No clue, but see ?assign perhaps if you have not done so already. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Feb 27, 2018 at 6:51 AM, Thomas Mailund <thomas.mailund at gmail.com> wrote: > Interestingly, the

I have load / store instructions that require accumulator. So a store looks like.. mov 3, acc st acc, addr I have specified "acc" as a separate register class containing only one register which is the "acc". The instr patterns are then splitted into: set imm:$src, ACCClass:$dst (generating the "mov" above) set ACCClass:$src, mem:$dst (generating the

Hi I never seen this before, and dont understand where too look for Please share some light on this. the host output is adding a extra domain. example bellow, its showing right IPs followed by NXDOMAIN !!! [root at 171-SYSLOG ~]# host 171-dc-a.xxxx.acc 171-dc-a.test.acc has address 10.254.228.226 171-dc-a.test.acc has address 10.254.225.45 Host 171-dc-a.test.acc not found: 3(NXDOMAIN) Host

Dear all, I installed "dovecot dovecoot-mysql postfix and postfix-mysql" from debian repository 7. I start them with "/etc/init.d/postfix start" and "/etc/init.d/dovecot start" but When i use nmap localhost I see the following output: root at sito:/etc/dovecot# nmap localhost Starting Nmap 6.00 ( http://nmap.org ) at 2013-09-12 07:35 IRDT Nmap scan report for

On Jan 17, 2008, at 4:12 AM, Sanjiv Gupta wrote: > I have load / store instructions that require accumulator. > So a store looks like.. > > mov 3, acc > st acc, addr > > I have specified "acc" as a separate register class containing only > one register which is the "acc". > The instr patterns are then splitted into: > > set imm:$src,