similar to: type III sums of squares in R

Hi everyone, Is it possible to use tapply(x,y,mean) if not all groups of x by y are of the same length (for example if you have one missing observation)? I tried tapply(x,y,mean,na.omit=T) but it doesn't work! Steffi -- --------------------------------- Stefanie von Felten Doktorandin ETH Z?rich Institut f?r Pflanzenwissenschaften ETH Zentrum, LFW A 2 Telefon: 044 632 85 97 Telefax: 044

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

Hi, I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor? If I type:

Hi, I have done different discriminant function analysis of multivariat data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also need the predicted values to produce a plot of the analysis, as far as I know. Here my code: pcor.lda2<-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr,

Hello, Does anyone know which method from Newcombe (1998)* is implemented in prop.test for comparing two proportions? I would guess it is the method based on the Wilson score (for single proportion), with and without continuity correction for prop.test(..., correct=FALSE) and prop.test(..., correct=TRUE). These methods would correspond to no. 10 and 11 tested in Newcombe, respectively. Can

Hi I have some density estimates obtained from density(). I would like to calculate the sum of squares of these. As the x values of the estimates are not the same, and I would prefer not to restrict the estiomate to a certain range of x values, how can I do the calculation? Lets say: d1 <- density(Data1) d2 <- density(Data2) If the x values would be the same, I would: ssq <- sum(

Hi I am using ggplot2 at the moment and I must say it is definitely better then ggplot - good work. My problem is that I am using facet_grid() in the following way: > p <- ggplot(ssq, aes(x=year, y=-log(ssq))) > p + geom_point() + facet_grid(me*gi~cs*rz) and it works nicely, except that I would like to have, in naddition to the values of me, gi, cs and rz the name of the variable.

Hi, Within my R code, I am using a C++ function as below: overlaps <- matrix(0, nrow=B, ncol=length(N)) overlaps.P <- matrix(0, nrow=B, ncol=length(N)) .C("speedUp", D, S, pD, pS, nrow(D), as.integer(N), length(N), ssq[i], i, as.integer(B), overlaps, overlaps.P, DUP=FALSE) the function "speedUp", is supposed to update matrices overlaps

Thanks! I went through the online posts which supports the power of .Call over .C. But my probably naive question is why does this work for my code with R but not R-devel? And another question is related to using .Call. Based on the manual page, I do not need to change the function parameters when using .Call. So I can run like this: .Call("sppedUp", D, S, pD, pS, nrow(D), as.integer(N),

Hi. I try to install openvpn at my centos 7 server (CentOS Linux release 7.6.1810 (Core)). According to many tutorials at the internet, I was using yum install epel-release -y to enable the repo. When (after installing epel-release) I try to install openvpn using yum install openvpn -y I get this error: Geladene Plugins: fastestmirror Loading mirror speeds from cached hostfile

I have a challenge that I want to share with the group. This is not homework (but I may assign it as such if I teach the appropriate class again) and I have found one solution, so don't need anything urgent. This is more for fun to see if others can find a better solution than I did. The challenge: I want to read a book in a given number of days. I want to read an integer number of

Hello all, I've tried to reproduce the contour plot that appears in the book of Venables and Ripley, at page 255. Is a F-statistic surface and a confidence region for the regression parameters of a non-linear model. It uses the stormer data that are in the MASS package. I haven't been able to reproduce the plot either in R ( version 1.5 ) and S. It makes the axes and it puts the

Hi, I wish to install R on my computer but I do not know how to do this. I have a Windows ME. Please advise on how to install R. Thank you. Stefanie [[alternative HTML version deleted]]

Dear Greg, Thanks for the new release. The decomposition of the SSQ is just what I need! Regards, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle >>> gregory_r_warnes at groton.pfizer.com 10/30/02 07:16PM >>> Version 0.7.3 of the gregmisc package

Hello, I run into a problem: > ftable(table(Fire, Standard, StoAll), col.vars=c("Fire","Standard")) Error in table(Fire, Standard, StoAll) : object 'Fire' not found I do not understand that because when I read the table everything seems correct. > Stocking_all<-read.table("P://Benchmark//analysis//r//stocking_10//stock

Dear R user, I have a question on using R to analyze data with repeated measurements. I have 2 species with several strains (12) per species, each of which has been measured twice with for a given trait. No particular covariance, just two measures. Now I want to analyze the data with an ANOVA (aov) considering these repeated measures to get the MSq and SSq for the species and strain level. I

Hi, I would like to prioritize VOIP traffic when we use the phone, but other times not do traffic shaping at all. Right now I have my openwrt router set up with htb to do shaping. In order to get it to work well I had to set my upload and download speeds much lower than my line speed. With these settings, I get good VOIP reception even while surfing the net and doing a long download. However,

I can't get the y axis to extend the full range that I need, which is -8 to 8 Here's my code. I tried using ylim, but it's still truncating at the extremes in my data. plot(x, y, type ="o", xlab="Panelist", ylab="T value",lwd=1.5,lty=1) xlim=range(1:14),ylim=range(-8:8), las=1) Any suggestions? Thanks! -- View this message in context:

Dear R-list! My question is related to an Anova including within and between subject factors and unequal group sizes. Here is a minimal example of what I did: library(car) within1 <- c(1,2,3,4,5,6,4,5,3,2); within2 <- c(3,4,3,4,3,4,3,4,5,4) values <- data.frame(w1 = within1, w2 = within2) values <- as.matrix(values) between <- factor(c(rep(1,4), rep(2,6))) betweenanova <-

Hi, recently I have started working with R (v. 2.7.2), and I have been using R's internal ARIMA_Like() function (from the "stats" package) to estimate some ARIMA models. In particular, I use ARIMA_Like() in a function "fn()" that I feed to the optim() method; the main goal is to find optimal ARIMA prediction models for some time series. The ARIMA_Like() function returns a