[R] ANOVA MS residuals estimation: two models

Hello I have an experiment with two factors MANAGEMENT and REGION. And I have three MANAGEMENT types and three REGIONS. Like this. REGION: A, B and C MANAGEMENT: N, P, ST the independent variable is called PHt I have set the aov in R in two possible ways. The first model I understand, but I

[R] Duncan test: 2-way ANOVA without repetition, but with multiple subjects

Hello I have a question on how to perform a Duncan test after I set a model like this (see below). The data consist of a dependent variable (PHt) and two dependent variables (REGION (3 levels) AND MANAGEMENT (3 levels)). MANAGEMENT N PS REGION A 196 196

[R] p value for mu: anova()

Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run it in R, I do not get the pvalue for mu. Can some one help me

[R] resample from data frame: unlinked columns

Hello I am trying to resample from this data set (see below). The function I am using so far is doing it by considering A and B columns as linked. I used this function. NUCh_rep-replicate(500,data[sample(1:nrow(data),replace=T),],simplify=F) What I need is to resample (500 times) with

[R] spatial analysis

Hello I am using library(spatstat). I managed to test of the random distribution of a collection of sampled trees in the forest. But what I want now is to test if their diameter is randomly distributed across the forest, for example the bigger trees are not clustered together. So far I have

[R] Multiple regression Categorical data

Hello It is possible to set up an lm() model where none of the categories of the categorical independent variable need to be used as references, I mean use the total mean instead. Thanks /R __ R-help@r-project.org mailing list

[R] box.cox

Hello I am using box.cox() and I get this error message: Warning message: 'box.cox' is deprecated. Use 'bcPower' instead. See help(Deprecated) and help(car-deprecated). I went to help but I did not understand the explanation, I am still wondering what is really happening. Thanks /R

[R] random effects

Hello I have a data set with fixed and random effects, therefore I am using the lme function: lm(y ~ xfixed, random=~1|xrandom, data) After this I want to get the F-values for both the fixed and random predictors. I can easily get the F-value and df for the xfixed predictors (anova()), but

[R] random effects

Hello I have a data set with fixed and random effects, therefore I am using the lme function: lm(y ~ xfixed, random=~1|xrandom, data) After this I want to get the F-values for both the fixed and random predictors. I can easily get the F-value and df for the xfixed predictors (anova()), but

[R] greyPalette() in legend plot

Hello I have a barplot where I get a grey gradient colors for each bar, I want to append a legend with the names of each bar and next to each name I need a box with the corresponding grey color. I managed to get the legend without boxes, but never got the wanted grey gradient in the legend:

[R] prcomp

Hello I am trying to run a PCA on the attached file, but I get this error message: pc-prcomp(data[,-(1:2)],scale=T)$x Error in svd(x, nu = 0) : infinite or missing values in 'x' Thanks in advance /R x y x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 1 25.49 45.62 125 156 165 130 179 152 82 165

[R] persp()

Hello I am trying to draw a basic black and white map of two European countries. After searching some key words in google and reading many pages I arrived to the conclusion that persp() could be used to draw that map. I have prepared three small example files, which are supposed to be the

[R] export/import matrix

Hello I have a problem on keeping the format when I export a matrix file with the write.table() function. When I import the data volcano from rgl package it looks like this in R: data[1:5,] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 100 100 101

[R] Factor Analysis with orthogonal and oblique rotation

Hello I seem to find only two types of rotation for the factanal function in R, the Varimax and Promax, but is it possible to run a orthogonal and oblique rotations in R? Thanks in advance Rosario __ R-help@r-project.org mailing list

[R] R in MAC add many extra ´s

Hello I have a annoying problem with R (which I am running in MAC). Every time I enter a function and I press enter to run it, I get many additional ´s, it is really bothering me. I have tried to google to find out how to fix the problem but I did not succeed. Anyone with similar experience

[R] EMA package

Hello I want to use the test.nested.model, in order to compare two nested models, which is under package EMA, however, I did not find it in the list of packages to install in R. How is this possible? Thanks for help in advance. Rosario __

[R] multiple variables Y and X

Hello I have a model with several hundred Y variables, and also several 1000 X variables. The model is linear lm(Y ~ X). My questions are: 1.- how to avoid writing all Xs variables? is list() the right function? 2.- about the multiple Ys with dependence among some of them, how to incorporate

[R] Large number of Y and X variables

Hello I have two issues 1. I have a linear model with several thousands of Y and X variables, so if I do not want to write them down one by one in the lm(Y ~ X) model how should I go about it. 2. Also, if I suspect that some of the Y variables could not be independent, what function should I

[R] matrix inverstion

Hello I have this matrix which I am trying to invert. I get a message about reciprocal condition number, what that does mean? XT_X [,1] [,2] [,3] [,4] [,5] [1,]30021 [2,]02011 [3,]00211 [4,]21140 [5,]1

[R] extract MSS from a lme model

Hello I am using a mixed model (lme Model- height ~ site + family + site*family) and I would like to estimate the additive variance (Vf) for which I need the MSS. When I carry out the Anova I do not get the MSS. Could someone explain me why? Thanks in advance. /R

[R] linear regression and t-distribution

Hello I have a data set with outlier and it is not normally distributed. I would instead like to use a more robust distribution like t-distribution. My question is if the coefficients of the regression are different from zero, but assuming a t-distribution. Could someone hint me what package

[R] multiple regression plane

Hello I have a multiple linear regression with two cofactors, I would like to represent a plane but I could not find any help which worked out. Any suggestions. Regards and thanks in advance. Rosario __ R-help@r-project.org mailing list

[R] Error in ANOVA for model comparison

Hello I am trying to compare two models using anova(), however I get a message error (see below). In the net I only found some information on certain library(car) for which one should use anova with A capital letter (Anova instead of anova), but I could not find car library as it says it does

[R] set Gamma parameters in glm model

Hello I hope my question makes sense. It is possible to specify the shape paramenters in a glm model with family Gamma? Thanks in advance Rosario [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] model values for alpha and beta in a glm gamma inverse distribution

Hello How could I include in the model a value for lambda in a glm family=poisson model? Or alpha and beta in a glm family=gamma(link=inverse) model? Thanks Rosario __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] initial values for alpha and beta in Gamma inverse

Hello I am trying to fit a model using glm function with family=Gamma(link=inverse). Is it possible to give initial values, and how? I tried to search for info but not managed to get any answer. Kind regards and thanks in advance. Rosario [[alternative HTML version deleted]]

[R] initial values for alpha and beta in Gamma inverse

Hello I am trying to fit a model using glm function with family=Gamma(link=inverse). Is it possible to give initial values, and how? I tried to search for info but not managed to get any answer. Kind regards and thanks in advance. Rosario [[alternative HTML version deleted]]

[R] extracting P values from lm model

Hello I am trying to get out of an lm model the fstatistics, however after I run the model I write names(Model) and the fstatistic does not appear only these. names(Model) [1] coefficients residuals effects rank fitted.values [6] assignqrdf.residual

[R] multivariate analysis

Hi I have 1800 response variables to regress on two factors (latitude and age), what is the script to run all response variables at once instead of writing 1800 models? Thanks R. Lähettäjä: r-help-boun...@r-project.org [r-help-boun...@r-project.org]

Re: [R] multivariate analysis

Hi I have 1800 response variables to regress on two factors (latitude and age), what is the script to run all response variables at once instead of writing 1800 models? Thanks R. __ R-help@r-project.org mailing list

[R] Export lmList

Hello I am trying to export and lmList by write.table(x, x.csv) and it complains. I have tried to convert the lmList into data.frame but also complains. Any suggestion? Thanks Rosario __ R-help@r-project.org mailing list

[R] individual intercept and slope

Hello I would like to extract the estimates for the intercept and slope by individual for growth from a lm fit. Any advice? Individual Time point Height 1 1 10 1 2 11 1 3 23 1

[R] plot by cathegories within a factor

Hello I have a data set summarized like this: File name= Height Group Ind Age Trait 1 1 1 20 1 1 2 21 1 2 1 22 1 2 2 21 1 3 1 24 1 3 2 45 1 4 1 23 1 4 2 26 2

[R] Error: could not find function lmList

Dear colleges, I would like to use the lmList function. I have installed the lme4 library but when I try to use the lmList function I get this error message. Error: could not find function lmList Here you can see what kind of messages I am getting when installing lme4 package. Does anyone

[R] FitARp

Dear R users, I have a data set with time series as continuous (time(day of measurement) = 1,7,14,21...), for each time I have measured height in a total of 100 individuals. I would like to correlate height with chlorophyll content (single measurement per individual), and for that I would

[R] (no subject)

Dear R users, I am trying to draw error bars in a bar plot, I use this code (tried many others which did not work): install.packages() library(gplots) y -c(39.02, 46.42) se - c(7.57,7.35) plot - barplot(y, beside=TRUE, col=0, ylim=c(0,47), axis.lty=1, main=far-red, xlab=latitude,

[R] restricting coefficient value

Dear colleges, I am traying to restrict the value of one coefficient in my glm model, see example: medparmodel - glm(los ~ white * (-0.21) + type2 + type3 + hmo, family=poisson(), data=medpar) but this does not work, and searching in other help webpages I did not manage to get any good hint.

[R] getting p values

Dear colleges I need to get the p values for a table with 15000 entries of t values. Does any of you know how to do it? I can, of course, get one by one but that is not sensible. Thanks Rosario __ R-help@r-project.org mailing list