similar to: variance component analysis for nested model

Suppose I have a set of binned data, counts exceeding a series of arbitrary thresholds, a total N, a minimum and maximum, those sorts of things. Is there a "standard" method for estimating arbitrary quantiles from this? My initial thought is that the counts and min/max give me solutions at various points along the empirical cdf. As the data are roughly log-normal, I thought maybe I

I am trying to use rsync for making snapshots, sharing hardlinks at the destination using the --link-dest option. In conjunction with this, I would like for --dry-run to report some reliable-ish number for how much space it is going to consume at the destination so that I can arrange to free such space prior the launching rsync "for effect". Right now, that isn't working: #

Hi everyone I'm using Tukey HSD as post-hoc test following a lme analysis. I'm measuring hemicelluloses in different species treated with three different CO2 concentrations (l=low, m=medium, h=high). The whole experiment is a split-plot design and the Tukey-function from the package multcomp is suitable for lme-analysis with random factors. The analysis works fine but I get a non

Hi all I hope you might help me with some aspects of producing a graph in lattice. There are three things I have struggling with and that is: 1. to separate the horizontal box rows from each other; 2. to change the colour of the horizontal and vertical strips to white; and 3. to place the axes labels on the left y axes and on the bottom x axes. I would really appreciate some help. I have put the

Hi all I was wondering if it is possible to get rid of the horizontal strips and produce each barchart with a left y axes and lower x axes only. Also can you specify an exact size of graph ie 88mm wide with a font size of 'x'. library(lattice) library(latticeExtra) n=as.factor(c(1:5,1:5))

Hello List, I am new to R and statistics in general. My two colleges and I would like to get some training in the proper usage of R. A quick internet search provided no results for the region of Basel in Switzerland. Does some one on this list knows of any trainers or schools in the region in question or where to look next? Regards, Dominik Riva Universitätsspital Basel Med.

Dear R-list I am (forced) to change from Linux to Mac and am now looking for a new editor for R. I would like one that features a split window (console + editor) as well as syntax highlighting. Can anyone help? Especially the split-window feature does not seem to be easily available in the editors desribed on the R-help site, except Emacs, which I am reluctant to start using. I am on a MacBook Air

Dear list, My apologies if a solution / explanation to this already exists on the list, but it is difficult to assign it to a certain keyword. test<-c(1:3) testfct <- function(x) {test[1]<-100} test [1] 1 2 3 testfct(1) [1] 1 2 3 Basically, I would like to write data into an external table that the function does not know. Why is this not working / what alternatives exist? Thanks,

hi, this may be trivial, but we can't seem to find anything adequate, (although there is a work around with match() ). We are looking for something along the lines of plot(table1[table1$var2==or("a","b","c","d"),"var1"]) would be handy, with the potential or() function leading to what plot(table1[table1$var2=="a" |

BaselR - The new R meeting We are pleased to announce the new R meeting to be held in Basel, Switzerland. BaselR will be held from 6:30-9:30pm on Tues, Apr 27 at TransBARent: http://transbarent.business.sv-group.ch Doors open at 6:30,pm with the presentations starting at 7:00pm Introduction: What is Basel R? Andreas Krause:... Graphing Pharma Data Yann Abraham: Graphics Charles Roosen: Web

Dear R developers, I want to advertise a small project for the educational purpose to show people how to link native routines. In R programming, we need to link to native routines in C, C++, or Fortran from the R environment. In most cases, the linking works within one package, namely the R code in a package calls the native code in the same package. This is not necessarily always the case.

Dear R-list member, I'm using the function pmnorm() (-->library(mnormt)) within a for-loop. Certain parameter values leads to an error message: "(In sqrt(diag(S)) : NaNs produced, In sqrt(1/diag(V)) : NaNs produced, In cov2cor(S) : diag(.) had 0 or NA entries; non-finite result is doubtful)" obviously because "NaNs" were produced. Is it possible to tell R that it

Dear all I imported a Stata .dta file with the read.dta-function from the foreign-package. The dataframe's dimensions are > dim(d.apc) [1] 15806 1300 Importing needs up to 15 min and calculations with these data are rather slow (although I subset the data before starting analyses). My questions are: 1. Has someone experiences importing Stata files (alternatives to read.dta) ? 2.

Hello Marc, samba is 4.1.22. I let the dcpromote choose the DC. DNS is internal Yes the windows DC has also global catalog Regards Olivier On Fri, Jan 8, 2016 at 4:38 PM, Marc Muehlfeld <mmuehlfeld at samba.org> wrote: > Hello Olivier, > > Am 08.01.2016 um 12:02 schrieb Olivier Weinstoerffer: > > I have a domain composed by 3 linux samba 4 AD servers. > > - Which

Hi everyone I have the following problem: My dataframe has 3 variables: ID, Year and and an outcome variable. The dataframe contains repeated measurements because the subjects filled out a questionnaire every year. The time span covers 2 years. Now I want to check if there is a significant change in the outcome over the 2 years with a paired wilcox.test. The problem: Not every subject

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

Dear all, I need to build a PBX based on Asterisk for a call center. I have worked with raw Asterisk but it's hard to work for big implementations think. Also I have worked with Trixbox CE for a small bussines and it was prette good, but I have not have many features like ACD. I know there is another version called Trixbox PRO -specially Call Center edition- that's not free but has got

Hi, For a n-level factor, I'd like to specify the first contrast and have the remaining n-2 constructed automatically so that the set is orthogonal. I then test the contrasts with summary.lm(anova-object). In S-Plus, the following works: >y.anova <- aov( y ~ C(CO2,c(1,0,-1)) ) >summary.lm(y.anova) In R, it fails with the following error: >levels(CO2) [1]

Dear all I'm not sure if I did the right analysis for my specific split splot design. We are studying biomass increase with different CO2 concentrations with four different functional plant groups (e.g. grasses, herbs, broad-leafed trees and conifers). Of each functional plant group we have four species. The design is orthogonal. The design is: Blocks: 2 (climate chambers, called

Hello A very short question: Using multiple plots with par(mfrow=c(3,3)), how can I get R to indicate one common y- and x-axis label? I tried to use text() in par, but this is then overwritten when I plot the graphs. -- Sebastian Leuzinger Institute of Botany, University of Basel Schönbeinstr. 6 CH-4056 Basel Ph. 0041 (0) 61 267 3511 fax 0041 (0) 61 2673504 email: Sebastian.Leuzinger@unibas.ch