similar to: nested ANCOVA

Hi. I have a factor and I want to extract just those elements that appear exactly once. How to do this? Toy example follows. > a <- as.factor(c(rep("oak",5) ,rep("ash",1),rep("elm",1),rep ("beech",4))) > a [1] oak oak oak oak oak ash elm beech beech beech beech Levels: ash beech elm oak > table(a) a ash beech elm oak

Dear list members, I'd like to perform a glm analysis with a hierarchically nested design. In particular, I have one fixed factor ("Land Use Classes") with three levels and a random factor ("quadrat") nested within Land Use Classes with different levels per classes (class artificial = 1 quadrat; class crops = 67 quadrats; and class seminatural = 30 quadrats). I have four

Dear all, I am trying to validate a model by comparing simulated output values against observed values. I have produced a simple X-y scatter plot with a 1:1 line, so that the closer the points fall to this line, the better the 'fit' between the modelled data and the observation data. I am now attempting to quantify the strength of this fit by using a statistical test in R. I am no

Hi list, I have been trying to pileup dataframes with different column names using rbind. I am getting the following error message. I was wondering if there is a way to go around this minor problem? Error in match.names(clabs, names(xi)) : names don't match previous names: G Thanks indeed for your information Regards Mahdi -- ----------------------------------- Mahdi Osman (PhD)

Hi together, I always get an error message with using ylim in termplot(), like this: > x<-(1:10) > y<-(10:1) > l<-lm(y~x) > termplot(l,ylim=c(1,2)) Is this a bug, or is there another possibility to do that? Especially, I would like to use term.plot() for gamlss objects. Thanks for your help! Andreas -- Echte DSL-Flatrate dauerhaft f?r 0,- Euro*!

Hi, MixMonitor saves partial conversation when non-standard voice packet size is set (Asterisk 1.4.18.1). For example, if SIP-peer has alaw:30 then saved file would contain only 67% of total conversation. With alaw:20 MixMonitor saves 100% of conversation. It seems that MixMonitor has hardcoded "packets per second" or "samples per packet" values. I did a lot of googling, but

Have you set port 2 as 'NT' in the mISDN config file (not the Asterisk one)? Also, you will probably need to set it to ptmp. You need to configure them in misdn.conf (the Asterisk one this time). Here's the tail of my misdn.conf (4 x BRI): [trunks] ports = 1,2 ; physical port numbers (as defined in mISDN.conf) context = inbound ; context for incoming calls in

Yes it sounds silly, but its true. I made a menu with some submenu with the simple menu system. After pressing 23 times return to get into the submenu or back to the mainmenu it freeze. I think this is a bug of syslinux. I tried it on 3 different PC/Servers always the same effect. Here is my default and one of my submenu files: default: #Projekt PXE-Boot von Franz Lang www.franz-lang.com

Dear readers, Is it possible to specify a model y=X %*% beta + Z %*% b ; b=(b_1,..,b_k) and b_i~N(0,v^2) for i=1,..,k that is, a model where the random slopes for different covariates are i.i.d., in lmer() and how? In lme() one needs a constant grouping factor (e.g.: all=rep(1,n)) and would then specify: lme(fixed= y~X, random= list(all=pdIdent(~Z-1)) ) , that?s how it's done in the

Hello, I have a win 2000 domain with 5 XP SP2 clients logging on to that domain.Let`s say Domain A In my network there is another domain, a samba 2.2.7a Suse, running.Let`s say Domain B. There are different Shares on this Domain B which belong to the Domain Users A (the same username and password on Domain A and B). Now first after booting the System and logging in the Domain A the Xp Clients

Hello, I want to ask if CentOS will support P3Scan (http://p3scan.sourceforge.net/) as RPM delivered with the distribution? For explanation: P3Scan is an transparent mailproxy. It uses ClamAV and SpamAssassin to distinguish between Spam- Virus- and normal Mail. Its more usefull to get an update via the repos instead of compiling it by myself. I hope you understand my englisch. ;) Greets Jens

Hi all, i'm having a rather strange problem with rsync (v2.6.8). Rsync is run to synchronize two identical, rather old Solaris 2.6 servers. All filesystems rsync fine except for one. Rsync bails out on this particular FS with the following messages: # /var/tmp/bin/rsync -avz -e "ssh -i $keyfile" --exclude-from=/tmp/rsync_data.excl server:/data/ /data receiving file list ... done ./

Hello, when testing FreeBSD-7.1-BETA i discovered that the floppy disk controller doesn't work correctly. Trying to format a floppy (perhaps with bad blocks) i get: Processing fdformat: ioctl(FD_FORM): Device not configured instead of the normal E letter. I then checked the same problem is present on FreeBSD-6.3 and it has been reported by Beech Rintoul (*) in 2006! Of course the floppy disk

Dear Wolfgang, Thanks for your response. No, I haven't tried doing the complete reinstallation as you suggest. If I recall correctly, this has happened on previous occasions but like I said in my email it's not all that disruptive to my workflow. It's more like a 'first world problem'. I was encouraged to submit the query because I had correspondence with someone involved

Hi all, I am having memory allocation problem with my R 2.2.1 for Mac OS. The following is the error message that I get. I do not get this message if I break down the large dataset in to sub datasets. I think breaking up the dataset is not a sustainable solution in the long run. The data that I am analysing is essentially big, and it would be reasonable to do the analyis on the whole dataset

Hi I''m very very new to tc iproute etc and have read the LARTC howto. What I want to do is create some "master" classes of bandwidth limit and below that per ip address which "inherits" from this master class. Example: one queue for 128Kbps other queue for 256Kbps What I want now is that for example in "class" 128Kbps the ip 10.0.0.5, 10.0.0.8 etc. goes

Hello, I am unable to figure out how to code a new column in a data frame based on an existing column that matches a column in a reference data frame, in a relational-db fashion. I would like this to maintain a minimum set of reference tables that may be reused over several similar datasets. Specifically, I have two data frames as listed below, 'Bos' and 'tree.' For each case in

Hello, I have used R in the past to conduct multiple comparisons on standard linear models, but am a bit confused as to how to go about doing it with a mixed effects model. I am conducting a bioindication study using carabid beetles in which I have four treatment types (forest harvest types with varying levels of canopy structure retention), and am using canopy closure percent as a covariate in

Dear colleagues, I wish to report a problem I encounter when trying to save a graphic to file. When I produce a graphic and try to save it R becomes unresponsive and I must force quit, and then restart R. The problem occurs when I try to overwrite an existing graphic: for example when I made changes to the graphic and want to save the graphic using the original file name. It only happens when I

Greetings, you kind and wonderful people. I've spent the past few hours or so (read: all day) attempting to get a basic PDC up and running, with no success. I've even gone to the extent of copying John Herpstra's config, which contains all sorts of weird and wonderful stuff, in the hope that he had something I missed. He didn't. I've got one box running debian testing, with