similar to: Multiple comparisons when interacction

I am returning this to the R-help list. Please keep followups on the list. Yes, it can be done. It is not currently easy because multcomp doesn't have the syntax yet. Making this easy is on Torsten's to-do list for the multcomp package. See the MMC.WoodEnergy example in the HH package. The current version on CRAN is HH_1.17. Please see the discussion of this example in R-help:

Hello, With this program I try to repeat analysis for different years. The results of the analysis are not printed when in the loop, except for the year sequence. What is wrong? Thanks a lot. for (i in 92:99){ cat("\n", "=============================================================\n", "YEAR =",i,"\n",

Hello, I send the message again with the data file as txt because it seems not to be accepted as csv in the R-help list. Data comes from a multiyear field experiment in which 4 levels of a treatment (2, 3, 4, 6) are compared to see the effect on yield. It is a randomized complete block design. The SAS code follows: options ls=95; data uno; infile 'data.txt'

I’m trying to obtain a plot of given function, say the density of a normal, in which I could indicate the area under the curve between two points (a,b). I.e. the P(a<X<b). Could anyone indicate me which is the instruction for obtain such a shadow area under a curve? I guess the solution is not difficult, but I can not find the procedure. Thanks Albert Sorribas Departament de Ciències

This is a generic request concerning were to look for finding appropriate information on a precise procedure in R. I’m using R for teaching introductory statistics and my students are learning how to deal with it. However, I find it difficult to locate some of the procedures. For instance, for basic crosstabulation, it is obvious that basic functions as table, ftable, and prop.table can be used.

Hi, A day or so after starting samba, some daemons (diferent forks) begin to hang. Then, the WinXP clients hang too completely. When I try to figure out what is happen, I see that smbd daemons hangs always in a fcntl64() call: # strace -p 6414 Process 6414 attached - interrupt to quit fcntl64(14, F_SETLKW64, {type=F_WRLCK, whence=SEEK_SET, start=3684, len=1} <unfinished ...> Process 6414

Hi, I would like to plot the following equation: dF(x)/dx=(k1+k2F(x))(1-F(x)) where k1 and k2 are parameters that I have estimated already. How can I plot the curve in R? Thanks! Montserrat Rue Universitat de Lleida (Spain) [[alternative HTML version deleted]]

Hi, I would like to compare the hazard functions of two samples using the Cox proportional hazards model. For sample 1 I have individual time-to- event data. For sample 2 I don't have individual data, but grouped data that allows to obtain a hazard function. I am wondering if there is an R function that allows to obtain a hazard ratio of the two hazard funtions (under the

Library gmodels include a function CrossTable that is useful for crosstabulation. In the help, it is indicated that one can call this function as CrossTable(data), were data is a matrix. However, when I try to use this option, it doesn't help. Any idea? Is there a bug? Thanks for your help. Prof. Albert Sorribas Grup de Biomatem?tica i Bioestad?stica Departament de Ci?ncies M?diques B?siques

This question is about including code and R results in Latex trough Sweave. It works fine and is of great help. Now may question is how to get the code and the output computed by R within a frame in Latex. Does any one have an idea on that. I tried \fbox, etc. but it seems to have problems with the code chunk generation of Sweave and I ger errors. I would appreciate any suggestion. I'll

I'm running samba-3.0.2-7.FC1 on a Fedora Core 1 box. I see this error message in the log file for each user on my network; [2004/03/26 10:55:35, 0] rpc_server/srv_util.c:get_domain_user_groups(372) get_domain_user_groups: primary gid of user [glen] is not a Domain group ! get_domain_user_groups: You should fix it, NT doesn't like that What does it mean? Martin Stacey IT Support

I've written down a simple yet a bit obfuscated perl script to put samba logs on a MySQL database, inspired on a little howto i found some time ago googleing around. It works by reading logs from a fifo instead of parsing actual log files, but allows to create sparse files too. Log lines are then parsed by matching to a regexp to extract desired values then inserting the record on mysql. It

I am sorry for the long *involved* post, however we are at wits end and advice would be greatly appreciated. Background: We have three servers in this example. Two of the servers (we will call Server A & B) each have a 2.7T raid 5 array that is mirrored with DRBD and monitored for failover with Heartbeat. NFS is running and exporting the array on Heartbeats virtual IP. "Server C"

I have several devices, mostly amd64 and armhf architectures, all of which updated samba to 4.17.6 without any problems. But one which is armel architecture failed to update; instead I get the following when I try to update: root at goflex:~# samba -V Version 4.17.5-Debian root at goflex:~# apt update Hit:1 http://ftp.uk.debian.org/debian bullseye InRelease Get:2 http://ftp.uk.debian.org/debian

This code below produces empty XXX.png files - if I use plot(), it works and if I enter the commands png(filename=fn) bubble(positions, do.sqrt=FALSE, main=q) dev.off() manually, it works as well. I am lost - any help appreciated. The weird thing is that it worked before I made some changes... (I don't have a copy of the working version...). R version: platform i686-pc-linux-gnu arch

Hello all, My statistical analysis training up until this point has been entirely done in SAS. The code I frequently used was: *Yield Champagin; data yield; set stress; if field='YV' then delete; if field='HB' then delete; if barcode='16187DD4015' then delete; if barcode='16187DD6002' then delete; if barcode='16187DD2007' then delete; if

Hi- I'm wondering if anyone can help me with my code. I'm coming up dry when I try to get a p-value from the following code. If I make a histogram of my resampled distribution, I find the difference between by groups to be significant. I've ranked the data since I have outliers in one of my groups. mange= c(35, 60, 81, 158, 89, 130, 90, 38, 119, 137, 52, 30, 27,

Hello, I am working on tests to compare figures. I have been using ImageMagick, which creates a figure signature, and I can compare a "test" figure signature against a saved "reference" figure signature. It seems to work pretty well. However, it is slow as it requires reading from the file system. Are there any options to compare figures on memory? For example, if I

Regarding point 3, as a moderator I have been helping Andrew get this post out to the list over the past week. His previous attempts were encoded in some way that the listserv rejected. He sent me the post via his gmail account and viewing the source I saw it had at least both plain test and HTML an I said it was worth a try to post it. Certainly on my mail client his post displays acceptably

Hi, I'm analyzing a longitudinal data set with 387 cows were observed in 63 days divided into 6 groups, and every 30 days was found to produce milk. Does not aim to model the time using regression. Only compare the groups differ in terms of milk production. There are many missing observations. Because the data are correlated I used the SAS program: proc mixed data=univar