similar to: Histograms per coding variable

Dear all Is there any way of doing a Tukey HSD post-hoc test after a Kruskal- Wallis or Friedman rank sum test (in the ctest package)? Thanks in advance, Albertus Dr. Albertus J. Smit Department of Botany University of Cape Town Private Bag Rondebosch 7700 South Africa Tel. +27 21 689 3032

Dear all I have made some lattice graphs (bwplot(dead ~ treat | group,...). 'group' is a factor with three levels (Artemia larvae, Abalone larvae and Abalone spat), and the result is a graph with three panels corresponding to the levels in 'group', ordered in alphabetical order from bottom to top, as expected. How does one re-order the order in which the levels of

Dear all I need to obtain an estimate of the 50% lethal dose (LD50) from a logistic regression model obtained by applying the glm procedure to some binomial data. The model appears to fit the data very well. I used dope.p from MASS to try and find LD50. The following output appears: > dose.p(iso.glm.logit, cf = c(1,3), p = 1:3/4) Error in dose.p(iso.glm.logit, cf = c(1, 3), p = 1:3/4) :

Dear all How does one get a levelplot (lattice library) to plot NAs in a different colour to that specified in the default colorkey? Thanks in advance, Albertus -- Dr. Albertus J. Smit Department of Botany University of Cape Town PO Box Rondebosch 7700 SOUTH AFRICA ------------------------------------------------- This mail sent through IMP: http://horde.org/imp/

Dear all I have just had the question from a colleague. I know that it is not directly related to R (I will probably use R to do the analysis), but I hope someone can give us some insight: Thanks, AJ Smit I sampled populations of a seaweed in the intertidal in order to estimate the standing biomass of that seaweed at that site. Due to clumped distribution patterns, I chose a stratified

Dear all I have a question regarding repeating a series of nested ANOVAs. The model is: fit <- aov(response ~ treatment + plot/rep) Treatment has three levels, there are three plots and each plot has four replicates. Now, the experiment was also repeated six times during the year. Ideally I would have liked to include sampling time (month) as an additional factor in the analysis, but

Dear all I have the following data (a shortened extract shown; some replictates of time deleted) to which I fitted the linear model given below: time group mass 11 control 0.019 11 control 0.014 14 control 0.0306 14 control 0.0289 14 control 0.0236 17 control 0.0469 17 control 0.0709 11 five 0.0077 11 five

Dear all I was looking at the aov documentation page and came across the following which seems like a contradiction to me: " This provides a wrapper to |lm| for fitting linear models to balanced or unbalanced experimental designs." (I presume 'This' refers to aov) and "|aov| is designed for balanced designs, and the results can be hard to interpret without

Dear R users I recently migrated from Statistica/SigmaPlot (Windows) to R (Linux), so please excuse if this may sound 'basic'. When running a nonlinear regression (V = Vmax * conc / (Ks + conc), i.e. Michaelis-Menten) on SigmaPlot, I get the output listed below: >>>Begin SigmaPlot Output<<< R = 0.94860969 Rsqr = 0.89986035 Adj Rsqr = 0.89458984 Standard Error of

Good day, Please am new on CentOS, may you help me with the upgrade from 5.8 to 6.2 using? Thanks a lot -- -- You Truly Eric Kom System Administrator - Metropolitan College _________________________________________ / You are scrupulously honest, frank, and \ | straightforward. Therefore you have few | \ friends. / -----------------------------------------

Hello I want to make some variables with the ifelse-function, but i don't know how to do it. I want to make these five variables; b2$PRRSvac <- ifelse(b2$status=='A' | b2$status=='Aa',1,0) b2$PRRSdk <- ifelse(b2$status=='B' | b2$status=='Bb',1,0) b2$sanVac <- ifelse(b2$status=='C' | b2$status=='sanAa',1,0) b2$sanDk <-

Hi all, This morning I received this notification from mdadm: This is an automatically generated mail message from mdadm running on server-mail.mydomain.kom A Fail event had been detected on md device /dev/md1. Faithfully yours, etc. In /proc/mdstat I see this: Personalities : [raid1] md1 : active raid1 sdb2[2](F) sda2[0] 77842880 blocks [2/1] [U_] md0 : active raid1 sdb1[1] sda1[0]

Hi... I have tumpled over quite a big problem on my machine. I'm running Debian Potato with OpenSSH-1.2.3 protocol 1.5 installed from the debian packages. When I ssh from any computer (including my own) then I can't get any programs to display X. I have turned ForwardAgent and ForwardX11 on in my ssh_config. When I turn on debugging this is what I get: 17:14 simpsons at lada% ssh -v

I am doing some testing on the adjtime syscall, and I want to dump out the value of the kernel variable "timedelta", and I am failing at it. Some of my attempts got "Module is no longer loaded" and some got "Symbol table entry is not a data object" I found some info that dtrace can print out kernel global variable, so what am I missing here? -- blu "Remember

Hello, I want to instal a normal windows-programm from http://www.tess-kom.de:8080/ with wine, but its not working. Error message: > > the program calcs could not execute > regards manu

Hi,I try to calcualte AIC or Loglik to GARCH model,But the Packege Tseries do not deal with them.How can I calculate AIC or Loglike to GARCH Model By Packege Tseries? Thanks. ____________________________________________________ Free Internet Access NOW! In Alexandria, Ismaileya, Suez, Portsaid, Hurgadha, Sharm Banha, Shebin El-Kom, Damietta, Tanta, Zagazig, Mansoura, Damanhour, Assyout, Qena

W dniu 18.04.2024 o?18:11, David Mulder via samba pisze: > On 4/18/24 1:03 AM, Jaros?aw K?opotek - INTERDUO via samba wrote: >> Hi all, >> >> I run cmd: >> samba-tool gpo manage scripts startup add \ >> {31B2F340-016D-11D2-945F-00C04FB984F9} \ >> /var/lib/samba/sysvol/fartest.local/scripts/startup.bat >> >> with result: >> [cut] >>

Hi, I am running the following loop, but it takes hours to run as n is big. Is there any way "apply" can be used? Thanks. ### Start nclass <- dim(data)[[2]] - 1 z <- matrix(0, ncol = nclass, nrow = nclass) n <- dim(data)[[1]] x <- c(1:nclass) # loop starts for(loop in 1:n) { r <- data[loop, 1:nclass] classified <- x[r == max(r)]

Hi Louis, You could try this: # find the index of the maximum value in each row of _data_, # disregarding the last column classified <- apply(data[,-(nclass+1)],1,which.max) ## or, if the maximum may be repeated: classified <- apply(data[,-(nclass+1)], 1, FUN = function(x) which(x == max(x))) # the variable _truth_ is just the last column of _data_ ? truth <- data[,nclass + 1] ?table

Hello everybody, This is a bug involving functions in core R package: graphics::hist.default, grDevices::nclass.FD, and base::pretty.default. It is not yet on Bugzilla. I cannot submit it myself, as I do not have an account. Could somebody else add it for me, perhaps? That would be much appreciated. Kind regards, Sietse Sietse Brouwer Summary ------- Floating point errors can cause a data