Displaying 20 results from an estimated 1000 matches similar to: "Difference of array and vector"
2007 Jan 02
2
R modules
Hi All,
I'd like to know what is the best way to organize R code in multiple modules=
and files. The R code we are writing is too much for a single file. Besides=
, there are a lot of reusable functions we'd like to factor out. But writing=
a package for that is quite an over-kill and might be too inflexible.
So what is the best way to organize R code into multiple files and reusable=
2006 Dec 26
4
vectorizing an iterative process.
Hi Everyone,
I am stuck with a simple problem. Suppose I have a vector x, and I want
to calculate y[i]=3Dx[i+1]-x[i], it is very easy. I just need to write
y<-x[2:length(x)]-x[1:length(x)-1].
Now if I know y, and want to know the vector x defined by
x[i]=3Dx[i-1]+y[i-1] for all i, how can I do this without a loop?
Thanks,
Geoffrey
PS. Sorry if you see a duplicate message. The previous
2009 Sep 06
5
ggplot2::qplot() -- arbitary transformations of coordinate system?
Hi,
Does anyone know how to do a coord_trans() in which the y-axis is
tranformed into (for example) -1000/y?
Thanks,
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
McCormick Road Charlottesville, VA 22903
Office: B011
2008 Feb 10
2
grep etc.
Dear R-helpers,
How do I transform
v <- c('insd-otsd', 'sppr-unsp')
into
c('insd--otsd', 'sppr--unsp')
?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
McCormick Road Charlottesville, VA 22903
2012 Mar 08
3
Packages 'effects' loads 'name' which conflicts with 'lme4'
Hi,
I would like to use the effect() function (actually a slightly modified version of it) on the output of the lmer() function in the lme4 package. But the effects package requires the nlme pacvkage, which is incompatible with lme4. Workaround?
______________________________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
for mail add: for FedEx or
2005 Aug 27
2
Defining an ex-gaussian PDF
How does one define PDFs as yet undefined in R, such as the ex-
gaussian, the sum of two RVs, one exponential, one Gaussian? The PDF
would then be the convolution of an exponential PDF, dexp(), and a
normal, dnorm().
Kindly cc me in your reply to r-help.
Thanks,
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400
2008 Apr 17
2
Conditionally swap items in a data frame
df1 <- data.frame(a = LETTERS[1:2], b = LETTERS[3:4], c = 1:2)
I am looking for an idiom that swaps the elements of df$a and df$b
when (e.g.) df$c == 2, resulting in
df2 <- data.frame(a = LETTERS[c(1, 4)], b = LETTERS[c(3, 2)], c = 1:2)
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA
2008 Jan 12
2
Lattice equivalent of par(mfrow = )
Dear r-helpers,
Does anyone have a straightforward example of putting together three
unrelated (expect for a common y-axis) xyplot() figures in what would
be in base graphics a par(mfrow = c(1, 3)) arrangement?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102
2006 Dec 19
2
Problem with glmmADMB
library(glmmADMB)
#Example for glmm.admb
data(epil2)
glmm.admb(y~Base*trt+Age
+Visit,random=~Visit,group="subject",data=epil2,family="nbinom")
Gives:
Error in glmm.admb(y ~ Base * trt + Age + Visit, random = ~Visit,
group = "subject", :
The function maximizer failed
******************
R version 2.4.1 RC (2006-12-14 r40181)
powerpc-apple-darwin8.8.0
locale:
C
2005 Feb 05
2
Problem installing Hmisc
What am I doing wrong?
> install.packages("Hmisc")
trying URL `http://cran.r-project.org/src/contrib/PACKAGES'
Content type `text/plain; charset=iso-8859-1' length 47565 bytes
opened URL
==================================================
downloaded 46Kb
trying URL `http://cran.r-project.org/src/contrib/Hmisc_3.0-1.tar.gz'
Content type `application/x-tar' length
2008 Jan 28
2
Package simex
Dear R-helpers,
It is not clear to me how you get measurement.error SD when you have a
single dataset, and it is not clear to me how sensitive SIMEX is to
errors in the estimates of measurement error.
Could someone please point me to the relevant literature?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400
2009 Jan 17
2
Concave Hull
Dear Friends,
Here is an algorithm for finding concave hulls: http://get.dsi.uminho.pt/local/
Has anyone implemented such an algorithm in R?
RSiteSearch('concave hull') didn't reveal one (I think).
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
Postal Address:
P.O.Box 400400, Charlottesville, VA 22904-4400
Express Parcels
2007 Apr 02
2
Why does lmList() fail when lm() doesn't?
Dear r-helpers,
Can anyone suggest why lm() doesn't complain here:
summary(osss.lm1 <- lm(logOdds ~ c.setSize %in% task, data = osss))
whereas in package:nlme (and in package:lme4)
osss.lmL <- lmList(logOdds ~ c.setSize %in% task | subj, data = osss)
# Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") :
# contrasts can be applied only to factors with 2 or more
2008 Sep 06
2
Sweave and/or beamer issue
Dear Friends,
I not sure whether this is an Sweave or a beamer problem.
The Rnw file:
\documentclass[compress,smaller]{beamer}
%\documentclass{article}
%\usepackage{beamerarticle}
\usepackage{Sweave}
\title{Psychophysics II}
\date{September 9, 2008}
\begin{document}
\frame{
\begin{Schunk}
\begin{Sinput}
> ro <- 0.2
> c <- seq(from = -3, to = 4, by = 0.1)
> fn <- 1 -
2006 Oct 15
4
Hide line ends behind unfilled circles?
Dear r-helpers,
xx <- c(0.000, 0.210, 0.714, 0.514, 1.000, 0.190, 0.590, 0.152)
yy <- c(0.000, 0.265, 0.256, 0.521, 0.538, 0.761, 0.821, 1.000)
aa <- c(19, 19, 19, 21, 19, 21, 21, 21)
x0 <- xx[c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 6, 6, 7)]
y0 <- yy[c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 6, 6, 7)]
x1 <- xx[c(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)]
y1 <- yy[c(2, 3, 3, 4, 6, 4, 5,
2008 Sep 26
2
ANOVA between & within variance
hi,
is there an option to calculate the 'within' & 'between' group variances
for a simple ANOVA (aov) model (2 groups, 1 trait, normally distr.) ?
or do I have to calculate them from the Sum Sq ?
thanks for your time and greetings,
gregor
--
Gregor Rolshausen
PhD Student; University of Freiburg, Germany
e-mail: gregor.rolshausen at biologie.uni-freiburg.de
tel. :
2007 Feb 08
1
Scope
Hi all,
When I write a script, I'd like to create a main() function so that I
only need to type main() t re-run it. However, I'd like all the
variables in main() to be global so that when the function terminates, I
still have access to the variables in the environment. Does anyone know
how to do that?
Best regards,
Geoffrey
_______________________________________________________=0A=
2008 Jan 02
1
Random Bernoulli sequences with given point-biserial correlation?
Dear R-listers,
Can someone suggest a method for generating a finite Bernoulli
sequence that is likely to have a given point-biserial correlation
with an existing Bernoulli sequence?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
2006 Oct 07
1
Installing Lindsey's packages
Dear r-helpers,
I downloaded http://popgen.unimaas.nl/~jlindsey/rcode/rmutil.tar (it
was originally .tgz, but got unzipped by my browser).
Can anyone give me detailed instructions on installing this and
Lindsey's other packages on R version 2.4.0 (2006-10-03)---(powerpc-
apple-darwin8.7.0, locale: C)?
_____________________________
Professor Michael Kubovy
University of Virginia
2007 Apr 26
1
xyplot() and controlling panel.polygon()
Dear R-helpers,
How do I tell panel.polygon what greoup and panel it applies to whithin
xyplot(y ~ x | c, groups = g
panel = function(x, y, groups, ...){
panel.polygon(x = xpol[c, g], y = ypol[c, g], default.units
= 'native')
panel.xYplot(x, y, groups, ...)
llines(x = c(1, 6), y = c(-24.283333, 35.941667), lwd = 2, lty =
3, col = 4)
}
x[c, g] and y[c, g] describe