similar to: Error con package agricolae

Hello, My boss wants me to do a Duncan's test, which is under the agricolae package. Unfortunately I am not versed enough in R to run my data. I have 7 subspecies of deer mouse for which I have 23 measurements which are my variables of interest. I have run an ANOVA for each of the set of subspecies and variables, my data look like this: subspecies WMF 1 rowleyi 2.50 2 rowleyi 2.30

Hi there, I am trying to run LSD.test(model) I used the following commands: attach(model) m1<- glm(ttl.m ~ site+year, family=quasipoisson, data= model) df<-df.residual(m1) MSerror<-deviance(m1)/df The following command did not work: comparison<- LSD.test (ttl.m, site+year, df, MSerror, alpha = 0.05, group=FALSE) I get an error message: Error: could not find function

I have used the correlation analysis (pearson) in the agricolae package to analyse my data and got unexpectedly low p-values (therefore making many more highly significant correlations in my data than I had expected). I am wondering if the p-values given should be subtracted from 1 to give the real p-value, because for each variable compared against itself has a p-value of 1 and I thought it

Dear list, I'd like to use path.analysis in the package agricolae in batch format on many files, retrieving the path coefficients for each run and appending them to a table. I don't see any posts in the help files about this package or the path.analysis package. I've tried creating an object out of the call to path.analysis, but no matter what I try, the function automatically prints

Dear all, W2k, R 2.5.1 I am trying out the AMMI function in the agricolae package, to analyse the dependence of environment for a certain cultivar. The function responds to four basic variables: ENV Environment GEN Genotype REP Replication Y Response My question concerns the 'REP' term. When I normally do an analysis of variance, the replication would refer to repeated observations

Hi, I did the analysis of variance of a split-plot and the effect of treatment was significant. I would like compare treatment means using Tukey. I can't extract the mean square to apply HSD.test to use in agricolae package. anava = aov(ganhos ~ Blocos + Trat*Supl + Error(Blocos/Trat)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

Hi there, I am doing multiple comparisons for data that is not normally distributed. For this purpose I tried both functions kruskal{agricolae} and kruskalmc{pgirmess}. It confuses me that these functions do not yield the same results although they are doing the same thing, don't they? Can anyone tell my why this happens and which function I can trust? kruskalmc() tells me that there are no

Cada vez que intento cargar el paquete "agricolae" me da el siguiente mensaje: Error: package or namespace load failed for 'agricolae' in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]): there is no package called 'spData' He intentado todo actualizar los paquetes, volverlos a instalar pero no hay forma. Si sabeis la solución ?? Un saludo. Juan

Estimados usuarios de R: escribo porque tengo un problema con el paquete rgl que no puedo resolver desde hace un tiempo mi sessionInfo es. sessionInfo() R version 2.15.3 (2013-03-01) Platform: i686-pc-linux-gnu (32-bit) el problema específico es el siguiente: library(rgl) > demo(rgl) demo(rgl) ---- ~~~ Type <Return> to start : >

muy buenas! acabo de empezar con linux y me he embarcau en un proyecto que = me queda un poco grande y que tengo que entregar dentro de poco... el caso es que quiero hacer funcionar el cbq.init y no puedo (el cbq.init es un script para controlar el ancho de banda de tu red). el caso es que tengo un debian woody 3.0 con el kernel 2.4, puestos todos los modulos de QoS e instalado el iproute.

Estimados: Estoy entrampado con el siguiente problema. Poseo un dataframe (f), con datos x, y, z, los cuales no son coordenadas ni alturas, sino niveles de tres variables cualquiera. Lo que busco es presentar en un gráfico es: a) el scatterplot de las variables ''x'' e ''y'', y además b) las isolineas de ''z'' para los distintos niveles de

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Dear all, I am trying to repeat an example from Sokal and Rohlfs "Biometry" -- Box 11.4, example of a randomized-complete-blocks experiment. The data is fairly simple: series genotype weight 1 pp 0.958 1 pb 0.985 1 bb 0.925 2 pp 0.971 2 pb 1.051 2 bb 0.952 3 pp 0.927 3 pb 0.891 3 bb 0.892 4

Hello, Can someone point me to what I should do to install the missing files? I am trying to generate self-signed certificates using mkcert.sh but I get the following error: $ /usr/local/share/dovecot/mkcert.sh error on line -1 of ./dovecot-openssl.cnf} 6213:error:02001002:system library:fopen:No such file or

Newbie needs advice, please. Sorry for the length of posting: I use rsync because it is such a great tool to reduce my download time for ISO files of different distros on remote servers that I want to update on my local disk. I've read the man file, the FAQ and the Wiki, and all have helped me understand much better my use of rsync for what I'm trying to do. However, I'm stuck on

Hi, I have a question about using "lme" and "aov" for the following dataset. If I understand correctly, using "aov" with Error term in the formula is equivalent to using "lme" with default settings, i.e. both assume compound symmetry correlation structure. And I have found that equivalency in the past. However, with the follwing dataset, I got different

Hi, I have a question about using "lme" and "aov" for the following dataset. If I understand correctly, using "aov" with an Error term in the formula is equivalent to using "lme" with default settings, i.e. both assume compound symmetry correlation structure. And I have found that equivalency in the past. However, with the follwing dataset, I got different

Hi all, Sent the following question yesterday, but haven't got any suggestions yet. So just trying again, can anyone comment on the problem that I have? Thank you! ------------- Hi, I have a question about using "lme" and "aov" for the following dataset. If I understand correctly, using "aov" with an Error term in the formula is equivalent to using

Hi Bert, Thanks for the suggestions. I tried lme with different control parameters, and also tried using "ML", instaed of "REML", but still got the same answers. Yes, I hope some gurus on this list could give me some hints. Thanks --- "Gunter, Bert" <bert_gunter at merck.com> wrote: > But they are close. This is almost certainly a > numeric issue --

Hola, buenas estimados. Acabo de instalar, después de algún tiempo,  R versión 4.3  y Rcommander en Manjaro Linux, 6.3.5  Y me he dado cuenta que no existe la opción de crear un Nuevo conjunto de datos. ¿eso es así ahora? Muchas gracias por su tiempo. Saludos. -- Ps. Félix Pérez Villoria Centro Kumon & Ayudantía Pedagógica. https://www.ayudantiapedagogica.cl General Ordóñez 30 oficina