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 do
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
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 to
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 replac
Hello
I have done resampling from a matrix data with two variables (A,B). After
resampling a get a set of lists (see below example for resample size 2). In
both lists I have the same two variables (A,B).
I want to export this data to an excel file where I have only two columns, one
per variabl
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 no
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
https://stat.ethz.ch/ma
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
_
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
ho
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
ho
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:
>
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
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 files
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 1
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
https://stat.ethz.
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
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-help@r-pro
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 t
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
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
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
__
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 t
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
https://stat.e
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 n
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
https://stat.ethz.ch/m
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
PL
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]]
___
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]]
___
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] "assign""qr"
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
https://stat.ethz.ch/mailman/lis
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]
kä
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
https://stat.ethz.ch/mailman
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
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
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 k
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 like
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", names.a
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.
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
https://stat.ethz.c
39 matches
Mail list logo