Displaying 20 results from an estimated 900 matches similar to: "Anova"
2012 Oct 04
4
Creating vegetation distance groups from one column
Hi R listers,
I am trying to group distances of nests to the vegetation into classes that
are define by (0-5m, 6-10m, 11-15m, 16-20m, 21-25m, 26-30m, 31-35m, 36-40m,
41-45m, 46-50m, 51-55m, 56-60m). Each row is a nest and all the distances to
the vegetation is in a column.
In plyr, I have tried - below script but I think I am going about this the
wrong way and am not successful.
#Veg index
2012 May 19
3
Q - scatterplot, plot function & trellis linear regressions???
Hi R-listers,
Q1) What is the difference between the scatterplot and plot function?
Q2) I am able to make a graph with the scatterplot function:
scatterplot(DevelopIndex ~ Veg,
+ data = Turtle,
+ xlab = "Vegetation border (m)",
+ ylab = "Embryonic development index")
And have been successful. But I do not know if the lines are for:
2012 Oct 08
2
aov() usage
Hi R-listers,
I am wondering if the function aov() in plyr is appropriate for two
different types of tests:
1) > summary(aov(EDI ~ VegIndex, data=data.to.analyze))
AND
2) > summary(aov(HSuccess ~ VegIndex + Aeventexhumed +
VegIndex:Aeventexhumed, data=data.to.analyze))
the later inclusive of an interaction of the two explanatory variables.
This is for obtaining P value and F
2012 Jan 25
6
How do I compare 47 GLM models with 1 to 5 interactions and unique combinations?
Hi R-listers,
I have developed 47 GLM models with different combinations of interactions
from 1 variable to 5 variables. I have manually made each model separately
and put them into individual tables (organized by the number of variables)
showing the AIC score. I want to compare all of these models.
1) What is the best way to compare various models with unique combinations
and different number
2012 Apr 01
1
NaN - trouble fixing NaN
Hi R-listers,
I am using the package plyr. I am just trying to get the hatching success
mean of each nesting event and have typed in the following and received the
below results:
> tapply(HSuccess, Aeventexhumed, mean)
A B C
0.2156265 0.1288559 NaN
What can I do about NaN? I should be able to get a result for event C
because I was able to
2012 May 26
1
Values in scatterplot??
Hi R-listers,
I am having trouble plotting the values for the R2 line, lowess smooth,
smothered conditional spread, outlier identification.
I have tried the function text and legend but was unsuccessful. Please see
below.
Your help would be appreciated. Thanks in advance.
Jean
library(car)
> SPHSHTL <- scatterplot(HSuccess ~ HTL,
+ data = data.to.analyze,
+
2012 May 15
1
Error in eval(expr, envir, enclos) : object 'Rayos' not found???
Hi R-listers,
I am trying to make a trellis boxplot with the HSuccess (y-axis) in each
Rayos (beach sections) (x-axis), for each Aeventexhumed (A, B, C) - nesting
event. I am not able to do so and keep receiving:
Error in eval(expr, envir, enclos) : object 'Rayos' not found
Please advise,
Jean
require(plyr)
resp <- read.csv("ABC Arribada R File Dec 12 Jean
2012 Jan 19
2
add1 GLM - Warning message, what does it mean?
Hi All, I am wondering if anyone can tell me what the warning message below
the model means?
J
add1(DTA.glm,~ Aeventexhumed + Veg + Berm + HTL + Estuary + Rayos)
Single term additions
Model:
cbind(MaxHatch, TotalEggs - MaxHatch) ~ Aeventexhumed + Veg +
Berm + HTL
Df Deviance AIC
<none> 488.86 4232.9
Estuary 1 454.96 4201.0
Rayos 3 258.80 4008.9
Warning
2009 Apr 22
2
function output with for loop and if statement
Hello all, turns out i'm having a bad R week. I am at my wits end with a function that I am trying to write. When I run the lines of code outside of a function, I get the desired output. When I wrap the lines of code into a function it doesn't work as expected. Not sure what is going on here. I suspected that the syntax of the if statement with the for loop was the culprit, but when I
2012 Feb 07
1
binomial vs quasibinomial
After looking at 48 glm binomial models I decided to try the quasibinomial
with the top model 25 (lowest AIC). To try to account for overdispersion
(residual deviance 2679.7/68 d.f.) After doing so the dispersion factor is
the same for the quasibinomial and less sectors of the beach were
significant by p-value. While the p-values in the binomial were more
significant for each section of the
2010 May 26
1
Fill a matrix using logical arguments?
Hello all,
I am going slightly mad trying to create a table for running
co-correspondence analysis.
What I have is seed bank and vegetation data, and my aim is to see if
the vegetation found in a site (containing several seed bank samples)
can predict the composition of a seed bank sample within that site. So
for this I need two tables with matching rows.
I have created an empty matrix,
2005 Dec 06
2
Constructing a transition matrix
Hi,
I would like to construct a transition matrix from a data frame with
annual transitions of marked plants.
plant<-c(1:6)
class<-c("seed","seed", "seed", "veg", "rep", "rep")
fate<-c("dead", "veg","veg","rep", "rep", "veg")
trans<-data.frame(plant, class, fate)
2012 Oct 08
0
Best method for comparing rectangles sections of beach
Hi R-listers,
I am trying to compare sections of the beach separated from the HTL to the
Veg (east to west), separated into indices (-5 to 30m), HTLIndex. Cross
parallel (north to south) are major beach sections (Rayos 1, 2, 3, 4 and
MNB). I am thinking to do an ANOVA for each independent rectangle of beach
(not exactly but will be treated as). The HTL from (0-5m) in Rayo 1 is one
rectangle of
2017 Nov 01
1
Creating Tag
i want to tag categories to its menuname.
i have a csv containing menu item name and in other csv i have a column
containing some strings,
i want to pick that strings from categories and look into menu items if
any menu item containing that string i want to create a new column next to
menu item name flagged as 1 otherwise 0
and the only condition is once a menu item flagged as 1 i don't need
2010 Nov 08
1
R help: reorder columns in an xyplot
Hello. First post, though I have read scores. Mostly i can solve my
problems reading the archives, but I lack the vocabulary to ask this one
right. Any social faux pais, I apologize.
I am sure this is a very basic question and I am embarrassed to ask, but I
spent several hours scouring the archives and trying various suggestions,
but I am too new to get any of them to work.
I am trying to
2008 Nov 26
1
Problem with aovlmer.fnc in languageR
Dear R list,
I have a recurring problem with the languageR package, specifically the
aovlmer.fnc function. When I try to run the following code (from R. H.
Baayen's textbook):
# Example 1:
library(languageR)
latinsquare.lmer <- lmer(RT ~ SOA + (1 | Word) + (1 | Subject),
data = latinsquare)
x <- pvals.fnc(latinsquare.lmer,
2011 Apr 21
1
Error running pvals.fnc in R version 2.13.0
Dear R-help:
I've been trying to run pvals.fnc in the newest version of R (2.13.0). The
function lmer worked fine, but when I tried to use pvals.fnc on the lmer
object, I got the following error:
"Error in pvals.fnc(elogr.subj.dys.sum.3x3.p, nsim = 10000) :
trying to get slot "coefs" from an object (class "summaryDefault") that is
not an S4 object."
How can I
2017 Dec 13
1
overlay two histograms ggplot
Hi all,
How can I overlay these two histograms?
ggplot(gg, aes(gg$Alz, fill = gg$veg)) + geom_histogram(alpha = 0.2)
ggplot(tt, aes(tt$Cont, fill = tt$veg)) + geom_histogram(alpha = 0.2)
thanks for any help!
Elahe
2008 Oct 01
1
pvals.fnc in lme4 and languageR
Hi everybody!
I was using the function pvals.fnc from package 'languageR' until April.
I do not know which version. Yesterday I updated all my packages
and tried to run my loop again. Now I get the following error message:
error in pvals.fnc(mm, nsim = 1000) :
MCMC sampling is not yet implemented in lme4_0.999375-27
for models with random correlation parameters
I guess it?s because of
2007 Jun 25
1
LanguageR pvals.fnc error message
Hi. I get an error message about not converging when I try and use the
pvals.fnc from the languageR library. The LMER analysis worked fine (See
below).
I am not an expert so I don't understand why the LMER worked but not the
pvals.fnc
Any help gratefully received.
- Mike
AIC BIC logLik MLdeviance REMLdeviance
-7324 -7254 3673 -7451 -7346
Random effects:
Groups