Displaying 20 results from an estimated 4000 matches similar to: "Finding the effect of Box-Cox transformation using "vis.boxcoxu""
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
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
2008 Apr 30
2
arcsine transformation
I have been trying to preform both a bartlett's test and an arcsine
transformation on some average percentage data. I've tried inputting it
different ways and I keep getting the same error message:
> head(workingdata)
DYAD BEFORE AFTER
1 BG-FL 4.606772 5.787520
2 BG-LL 5.467503 7.847395
3 AD-MV 5.333735 11.107380
4 MM-FL 5.578708 12.063500
5 MM-MV 2.037605 6.415303
6 MM-RM
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
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 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
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
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
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
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 -
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. :
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,
2006 Feb 09
3
about Cox-Box transformation
Dear R-users,
I am using R version 1.8.0-1 under Suse 8.2.
I need to use the boxcox command because I want to apply a Cox-Box
transformation to a vector of rainfall values.
Within the libraries, the MASS library is present, but I don't know
whether this means that is automatically installed or not.
The command doesn't work. What do I have to do in order to make it
working? Is just a
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
2007 Apr 21
1
Changing appearance of band edges in xYplot()
Dear R-helpers,
I would like the bands in an xYplot() to be gray without dark edges.
I've tried playing around with col.bands, but that seems not to
affect the edges, although changin col.fill does change the fill as I
would expect.
I would appreciate pointers.
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box