search for: oneway

Hello everyone: I'm a new member of this group. I have a question about "oneway.test". When I use "anova(lm(....))" to analysis the ANOVA, I can get the information about Sum Sq and Mean Sq. (The R code and the results are as follows.) > anova(lm(BackCalac~factor(Assay),da...

Hi all, I want to write a wrapper for an analysis of variance and I face a curious problem. Here are two different wrappers: fun.1 <- function(formula) { summary(aov(formula)) } fun.2 <- function(formula) { oneway.test(formula) } values <- c(15, 8, 17, 7, 26, 12, 8, 11, 16, 9, 16, 24, 20, 19, 9, 17, 11, 8, 15, 6, 14) group <- rep(1:3, each=7) # While the first wrapper works just fine: fun.1(values ~ group) # the second throws an error: fun.2(values ~ group) Error in substitute(formula)[[...

...;m looking for an R function to fit a one-way ANOVA with one factor containing 10 levels. The factor levels have different numbers of observations (varying between 20 to 40). For most of the dependent variables i'm testing there are unequal variances among the factor levels. I see the function oneway.test: oneway.test(variable ~ factor, data=dataset) which by default does not assume equal variances. I also see the basic ANOVA function (aov), which, I think, is OK for unbalanced designs. I get a different F ratio, however, when using these two functions on the same data, which seems to indicat...

Hi all, I did not get this error in Asterisk 1.2.5 release. I am testing on Asterisk SVN-trunk-r15187 to avail the PARKEDAT variable. - I park the call using ParkAndAnnounce - plays moh. - accept the call using ParkedCall The following errors are coming on the console and there is oneway audio - no audio after Music-On-Hold at caller's side. Please advice. I am testing using cisco 7902 phones and using cisco 2800 router. Codec is g711ulaw regards, -- Executing ParkedCall("SIP/192.168.50.2-09cbd610", "366") -- Channel SIP/192.168.50.2-09cbd610 conne...

...hoc, I use library(userfriendlyscience) and I obtain Games-Howell output but it's impossible for me to obtain a letter superscript comparison between virus groups as it is obtained through library(agricolae). The code used it was the following: library(userfriendlyscience) data(sweetpotato) oneway<-oneway(sweetpotato$virus, y=sweetpotato$yield, posthoc = 'games-howell') oneway I try with cld() importing previously library(multcompView) but doesn't work. Can somebody could helps me? Thanks in advance, David Bars.

Hi all, I apologize for what might be a silly question. I am interested in doing a one way anova. This is not too hard in and of itself, either with anova, aov or oneway.test . However, I need to 1) get pvalues, 2) do a posthoc analysis with Tukey HSD, 3) and have (sometimes) an unbalanced design. I just can't seem to put all the pieces together. Any suggestions? Thanks in advance, Dan. -- Daniel C. Jupiter, Ph.D. Postdoctoral Research Associate Departme...

Hi, I have the Cisco PIX 506 firewall right in front of the asterisk and I am getting a one-way audio. I need your help/guidance to resolve this problem. I have the "fixups" disabled for SIP in the Cisco PIX 506. Any help rendered by you in this subject is greatly appreciated. I have been breaking my head trying to resolve this problem for more than one month. I have included the

Hi, I am encountering a difficulty I don't understand. Be patient, I'm very new to analysis of variance. If I load this data: example12_7=read.table("http://msemac.redwoods.edu/~darnold/math15/data/chapter12/example12_7.dat",header=TRUE) The run the oneway.test: oneway.test(time~drug,data=example12_7,var.equal=TRUE) I get these results: data: time and drug F = 4.1881, num df = 2, denom df = 16, p-value = 0.03445 Now, I've done the problem by hand and this result agrees with my calculations. Now I try aov and get these results: res1 <- a...

Hi helpeRs. When I do two-ways tables dimnames in columns are almost left-adjusted: <- factor(rbinom(30,1,.3),labels=c("first","second")) > b <- a > table(a,b) b a first second first 21 0 second 0 9 If I do a oneway table however, dimnames are right-adjusted: table(a) a first second 21 9 Actually I must admit I dont know if its the dimnames or/and the cell counts that are adjusted ;-) Is there a way that I can get left-adjusted dimnames i oneway tables. Thans in advance Steen Steen Ladelund,...

...#39;m running R 2.7.2 on WinXP cheers *>data("Gasoline",package="Ecdat") >Gasoline_plm<-plm.data(Gasoline,c("country","year")) >pdim(Gasoline_plm) **Balanced Panel: n=18, T=19, N=342 * *>summary(plm(lgaspcar~lincomep,data=Gasoline_plm**)) **Oneway (individual) effect Within Model Call: plm(formula = lgaspcar ~ lincomep, data = Gasoline_plm) Balanced Panel: n=18, T=19, N=342 Residuals : Min. 1st Qu. Median 3rd Qu. Max. -0.40100 -0.08410 -0.00858 0.08770 0.73400 Coefficients : Estimate Std. Error t-value Pr(>|t|)...

Hi I am attempting Anova analysis to compare results from four groups (Samp1-4) which are lists of intensities from the experiment. I am doing this by first creating a structured list of the data and then conducting the ANOVA (Script provided below). Im an R beginner so am not sure if I am using this correctly. Two major questions I have are: 1) Is using the code (zzz.aov <- aov(Intensity ~

Does oneway.test do what you want? Hope this helps, Matt Wiener -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Sven Hartenstein Sent: Monday, June 21, 2004 3:23 PM To: r-help at stat.math.ethz.ch Subject: [R] Welch-JM-Test or Br...

Dear list members, I am looking for an alternative function for a two-way ANOVA in the case of variance heterogeneity. For one-way ANOVA, I found oneway.test(), but I didn't find anything alike for two-way ANOVA. Does anyone have a suggestion? Thank you! Maike Luhmann Freie Universit?t Berlin

My system is: Redhat 9.0 + Asterisk + ISDN4Linux + Teles 16.3 ISA Passive card I haven't sound card. Comunication between two SIP Clients is OK Comunication between PSTN and SIP Client is OneWay (i cant recive dtmf and voice from pstn) is it needed sound card ?

Hello Using R 2.2.1 on a Windows machine. Has anyone programmed the Welch test for equality of variances? I tried RSiteSearch, but this gave references to t test and oneway.test, which are not quite what I need.....I need the Welch test itself, for use in a meta-analysis (to determine if variances are equal). TIA Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development and Research I...

...data, model="pooling") fixed.mod1<- plm(Y~X1+X2,data = Pdata, model="within") rand.mod1<- plm(Y~X1+X2,data = Pdata, model="random") summary(pool11) summary(fixed.mod1) summary(fixef(fixed.mod1)) summary(rand.mod1) ##################output#################### Oneway (individual) effect Pooling Model Call: plm(formula = Y ~ X1 + X2, data = Pdata, model = "pooling") Unbalanced Panel: n=6, T=4-6, N=34 Residuals : Min. 1st Qu. Median 3rd Qu. Max. -19.400 -9.810 -0.648 8.490 23.900 Coefficients : Estimate Std. Error t-value P...

...calling asterisk always tries switch simple_bridge to native_rtp Bridge 0422bfa0-9d22-4bba-9108-a3f14d7d1cab: switching from simple_bridge technology to native_rtp in endpoints table ?direct_media? sets to ?no? on all endpoints but it doesn?t help. if native_rtp not work for some reason I have oneway audio. how can I fix this? if I add mix_monitor it works, but it?s not a right way to fix this issues. Asterisk 13.2.0 -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.digium.com/pipermail/asterisk-users/attachments/20150318/81c272b4/attachment.html&...

...X1+X2,data = Pdata, model="within") > rand.mod1<- plm(Y~X1+X2,data = Pdata, model="random") > > > summary(pool11) > summary(fixed.mod1) > summary(fixef(fixed.mod1)) > summary(rand.mod1) > > > ##################output#################### > > Oneway (individual) effect Pooling Model > > Call: > plm(formula = Y ~ X1 + X2, data = Pdata, model = "pooling") > > Unbalanced Panel: n=6, T=4-6, N=34 > > Residuals : > Min. 1st Qu. Median 3rd Qu. Max. > -19.400 -9.810 -0.648 8.490 23.900 > > Coeffi...

Does anyone know how can I calculate power in wilcoxon.test, kruskal. test and oneway.test in an F? I have already found power.t.test for calculating power in t.test but it seems that there is nothing for the others. Any answer could be useful. Thanks for your interest. --------------------------------- [[alternative HTML version deleted]]

...ata and I am getting an error. Please let me know if I am missing anything from my code. coag<- matrix( scan("//Samba3/nair/R/blood.dat", sep=","), 24, 3, byrow=TRUE) colnames(coag) <- c("time","diet","order") coag <- as.data.frame(coag) oneway.test(time ~ diet, data=coag, var.eq=TRUE) coag.mod <- aov(time ~ diet, data=coag) options(scipen=5) options(digits=5) TukeyHSD( coag.mod ) I get the following error when I run the above code Read 72 items Error in rep.int(n, length(means)) : Unimplemented feature in rep In addition: Warning m...