Displaying 20 results from an estimated 6000 matches similar to: "Unbalanced design help"
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
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 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
2006 Mar 06
1
JGR doesn't launch after updating to R 2.2.1
Dear r-helpers,
After updating to R version 2.2.1, 2005-12-20, powerpc-apple-
darwin7.9.0 (from binary), JGR wouldn't launch. I then deleted
everything I could find related to JGR, and reinstalled it from the
latest binary. Both R and JGR are installed in /Applications/
\~LocalApps/mathStat/.
Advice on how to identify the source of the problem?
My machine:
Machine Name:
2007 Feb 18
1
Suppressing \newcolumntype{} declaration in latex.default() in Hmisc
Dear r-helpers,
When use latex() on a matrix, I set the option dcolumn = T. As a
result, in front of each tabular I get
\newcolumntype{.}{D{.}{.}{-1}}
The LaTeX compiler complains about these multiple redeclarations.
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room
2007 Mar 09
1
Problem with ci.lmer() in package:gmodels
Dear Friends,
Please note that in the following CI lower > CI higher:
> require(lmer)
> require(gmodels)
> fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject),
sleepstudy)
> ci(fm2)
Estimate CI lower CI upper Std. Error p-value
(Intercept) 251.66693 266.06895 238.630280 7.056447 0
Days 10.52773 13.63372 7.389946 1.646900
2008 Oct 06
1
Color of title in legend()
Dear R-helpers,
This produces a red title. How do I make it black, without making the
labels black?
legend(title = 'Channel Integration', 0.7, 0.3,
c('loudness','luminance'), col = c('red', 'blue'), text.col = c('red',
'blue'), lty = 1, xjust = 0.5, bty = 'n')
_____________________________
Professor Michael Kubovy
University
2009 May 05
1
Find cyclically identical binary sequences
Dear R-helpers,
I need to generate all the binary sequences of length n (here n = 8)
that start with 1 and have no fewer than two of each digit, and are
not cyclic permutations of each other. Here is what I have done:
len <- 8
df <- as.data.frame(numeric(2^(len - 1)) %o% numeric(len))
require(partitions)
for (i in 1:2^(len - 1)) df[i, ] <- binary(i, dim = len)[[1]]
df <-
2006 Apr 14
1
Dotplot x-axis
Here's a small dataset:
type <- c('hierarchical','partial','single','complete','single
+hierarchical','single+partial','partial+hierarchical','single+partial
+hierarchical')
freq <- c(1455,729,688,65,29,28,16,17)
lodds <- log(freq/(3027 - freq))
dotplot(type~lodds)
I would like to have the x-axis have ticks at nice,
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
2007 Mar 21
1
package:AlgDesign and .Random.seed
Dear r-helpers,
Could you please help me solve the following problem: When I run
require(AlgDesign)
trt <- LETTERS[1:5]
blk <- 10
trtblk <- 3
BIB <- optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk))
In response to the last command, R complains:
Error in optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk)) :
object ".Random.seed" not found
The
2007 Jan 29
1
lmer2 error under Mac OS X on PowerPC G5 but not on Dual-Core Intel Xeon
> (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy))
Error in as.double(start) : Calloc could not allocate (888475968 of
4) memory
*************************
> sessionInfo()
R version 2.4.1 (2006-12-18)
powerpc-apple-darwin8.8.0
locale:
C
attached base packages:
[1] "grid" "datasets" "stats" "graphics" "grDevices"
2007 Oct 27
1
Selectively swapping labels between factors
Dear R-helpers,
I'm trying to selectively swap labels between two factors, depending
on an indicator variable i.
Can you point me to a solution, and perhaps how I could have found it?
labels(fact1) is a character vector of r row numbers
levels(fact1) is a character vector of the n < r unique levels
How do I then get the character vector of length r of the levels of
fact1?
Once I have
2007 Apr 10
2
Positioning in xyplot
Dear R-helpers,
I have an xyplot
of the following kind:
xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4',
'3')),
strip = strip.custom(var.name = 'g', factor.levels = c(1', '5', '2',
'4', '3'),
strip.levels = T, strip.names = T, sep = ' = ', shingle.intervals =
NULL),
data =
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
2007 Oct 31
2
mosaic row and column order
Dear R-helpers,
I wasn't able to find out how to override the alphabetical ordering
of the rows and columns in a vcd::mosaic plot. I would like to have
them each ordered by numerical values in a different column of the
data frame that contains the contingency data.
I would be most grateful for a pointer toward the solution.
Thanks,
MK
_____________________________
Professor Michael
2007 Mar 17
1
Font in JavaGD() & pdf()
Dear r-helpers,
When I do an xYplot and display the result in a JavaGD() window, the
font is sans-serif (presumably Helvetica). When I send the figure to
a pdf, I get a serif font (presumably times). How do I insure that
the font in the pdf is indeed the default sans serif?
> sessionInfo()
R version 2.4.1 (2006-12-18)
i386-apple-darwin8.8.1
locale:
C
attached base packages:
[1]