Displaying 20 results from an estimated 6000 matches similar to: "4253H twice implementation"
2012 Aug 03
1
SEM standardized path coefficients
Hello,
I have conducted an SEM in which the resultant standardized path coefficients are much higher than would be expected from the raw correlation matrix. To explore further, I stripped the model down to a simple bivariate relationship between two variables (NDVI, and species richness), where it's my understanding that the SEM's standardized path coefficient should equal the correlation
2013 Jan 28
1
incorrect import?
Dear all,
I'm not getting what I'm doing wrong. The line below from my read.fsa.bin function throws an error when just loading my AFLP package and disappears when I load the zoo package as well.
#the line that throws the error
Index <- which(Peak == rollmax(Peak, k = 1 + 2 * floor((min(diff(SizeStandard)) * Fs - 1) / 2), fill = -Inf))
#the error
Error in UseMethod("rollmax")
2012 Nov 05
2
averaging a list of matrices element wise
Dear all,
I have a list of n matrices which all have the same dimension (r x s). What would be a fast/elegant way to calculate the element wise average? So result[1, 1] <- mean(c(raw[[1]][1, 1] , raw[[2]][1, 1], raw[[...]][1, 1], raw[[n]][1, 1]))
Here is my attempt.
#create a dummy dataset
n <- 3
r <- 5
s <- 6
raw <- lapply(seq_len(n), function(i){
matrix(rnorm(r * s), ncol =
2009 Aug 13
4
un run run...
Hi All,
I am running an Rscript with a bunch of algorithms that are UNSTABLE
under some parameter settings.
At a certain point one of them sends error massage and my whole run STOPS!
What I would like is to save the error massage in some file or variable and
carry on to the next command line without stopping this run...
Any help or ideas would be welcome, please, with a concrete example
(not
2009 Mar 11
3
Mixed models fixed effects
Dear All,
This may sound like a dumb question but I am trying to use a mixed model to
determine the predictors of bat activity along hedges within 8 sites. So my
response is continuous (bat passes) my predictors fixed effects are
continuous (height metres), width (metres) etc and the random effect is
site - can you tell me if the fixed effects can be continuous as all the
examples I have
2012 Oct 16
2
gam (mgcv) problem: Error in while (mean(ldxx/(ldxx + ldss)) > 0.4) { :, missing value where TRUE/FALSE needed
Hi All,
I'm running into a problem with GAM (in the MGCV package). When I try
to estimate the model, I get the following error message:
1> fit <-
gam(ndvi~s(rain)+s(temp)+s(rainl1)+s(rainl2)+s(rainxY)+s(rainl1xY)+s(rainl2xY)+s(tempxY),
data=dsub, weights=wvec)
Error in while (mean(ldxx/(ldxx + ldss)) > 0.4) { :
missing value where TRUE/FALSE needed
Using
2009 Sep 25
3
Problem on plotting TS using GGPLOT
Hi, I have following codes :
library(zoo); library(ggplot2); library(plyr)
dat <- rnorm(306); vv <- letters[1:6]; dat1 <- data.frame(dat, vv)
dat2 = zooreg(rnorm(51), as.yearmon(as.Date("2000-01-01")), frequency=12)
ggplot(dat1) +
geom_line(aes(y=dat, x=index(dat2), colour=vv), group=vv, size =
1.3)
However I got error while plotting them :
2009 Jul 07
6
Uncorrelated random vectors
Hello,
is it possible to create two uncorrelated random vectors for a given distribution.
In fact, I would like to have something like the function "rnorm" or "rlogis" with the extra property that they are uncorrelated.
Thanks for your help,
Luba
[[alternative HTML version deleted]]
2009 Dec 22
2
Nested For loops
Dear R experts,
Might be very simple question to ask but would be insightful. As the same story of nested "for loops". following is the code that I am using to get the autocorrelation function of the sample data. I have tried to get rid of for loops but since I am touching R after such a long time that I need to practice more but I need help to revive my skills. I know that apply() or
2009 Feb 11
2
Label bars in a faceted bar plot in ggplot2
Hi List,
I am running R 2.8.0 on a Windows XP machine, running ggplot2 version 0.8.1
I want to label the bars in a faceted grid barplot. Reproducible R
code is given below:
#### reproducible facet barplot #####
library(ggplot2)
# Dataset from which to create the barplot
ml <- rep(1:10,2)
vals <- rnorm(20,mean = 10, sd=1)
type <- c(rep("MAPE",10),rep("AIC",10))
2010 Jul 20
1
define subgroups based on position in table
Dear list,
I have a data frame with one column (group) and want to add a second column (sub) with a serial number that says to which subgroup a cell belongs. A subgroup contains the consecutive rows of the same group. The number of a subgroup is based on its position in the table. The first subgroup of A's and B's should have nr A1 and B1, the second nr A2 and B2, etc.. I hope the
2008 Aug 07
3
Bug in rnorm? (PR#12016)
Full_Name: Paul Eckermann
Version: 2.7.1
OS: Windows XP Service Pack 2
Submission from: (NULL) (129.127.183.12)
This is the first time I have submitted a bug report, so apologies if I have not
followed the correct protocol.
If I enter
length(rnorm(100*.29))
it returns 28 rather than 29.
y<-(1:100)/100
z<-sapply(y,function(x) length(rnorm(x*100)))
z!=100*y
indicates that it has
2012 Nov 14
2
aggrete data from combination
Dear R users,
A have a dataframe (matrix) with two collumns (plot, and diameter (d)). I
want all diameters values for different combination of plots.
For example I want all d values for all posible combination, 100C2 (all d
values for plot 1 with all d values in the plot 2.......with all d values
from plot 1 with all d values from plot 100, ...... with all d values from
plot 99 with all d values
2009 Jul 22
1
Link to documentation in another package
Dear all,
One of the functions that I wrote (ggsave.latex) extents the
functionality of a function (ggsave) in another package (ggplot2).
Instead of copying all the information I would like to create a link in
the helpfile of ggsave.latex to the helpfile of ggsave. I tried
\code{\link{ggsave}} and \code{\link{ggplot2::ggsave}}, but neither
worked. Both cases gave a 'missing link' warning.
2010 Apr 19
2
ecdf
Hello,
I'd like to plot an empirical cumulative distribution function, except
instead of the fraction of values < x, I'd like the fraction of values > x.
I think this can be done using the ecdf function in {Hmisc}. I installed
the package and loaded it. However, when following the example given in the
documentation, I get an error:
x <- rnorm(100)
ecdf(x,what='1-F')
2004 Mar 20
2
labels on axis(4) and adaptation of legend to size of the plot
Hi weekend R helpers,
*Is it possible to get a ylabel on the right hand side axis ==> axis(4)?
Thankx,
Jan
opar <- par(mfrow = c(2,1))
plot(ts.Origi[,1],ylab='NDII', main=name)
# Add the legend text in the right order
legend.txt <- c("NDII", "Inverse KBDI")
legend(2001.8,0.29, legend.txt,col=c(1,2), lty=1)
2010 Oct 06
1
R getting slower until it breaks...
Hello R-users,
I'm currently facing a pretty hard problem which I'm hopping you'll be able to help me with. I'm using R to create images. That alone is not the problem, the problem is that I'm using R to create 168 000 images... My code (which is given below) use different package (raster and rgdal) to import a image (size 20gig) and divide it into 168 000 pictures that are
2010 Apr 29
5
reduce size of pdf
is there a way to reduce the size of pdf files in R: ?
compression?
lower dpi ?
or some other option?
2009 Sep 08
4
barplot with lines instead of bars
Dear useRs,
I want to plot the following barplot with lines instead of bars. Is there a way?
data <- data.frame(cbind(k = 0:3, fk = c(11, 20,7,2), f0k = c(13.72, 17.64, 7.56, 1.08), fkest = c(11.85, 17.78, 8.89, 1.48)))
d <- t(data[,2:4])
barplot(d, beside=TRUE)
Regards,
Rafael.
____________________________________________________________________________________
[[elided Yahoo
2012 Apr 25
1
Strange bug in my package
Dear all,
I get a bug in the examples of my AFLP package on R-forge (https://r-forge.r-project.org/R/?group_id=1027) but only on the Linux version. The windows version compiles. The Mac version skips the examples and compiles. The strange thing is that the packages compiles on my Ubuntu 10.10 machine with R 2.15.0. Therefore I can't reproduce the error.
I have traced the problem at