binomial, data = unempt )
modw <- glm( cbind( short, long ) ~ age, family = binomial, data = unempw )
library( effects )
modt.ef <- effect( 'age', modt ) # works!
modw.ef <- effect( 'age', modw ) # doesn't work!
# Error in eval(expr, envir, enclos) : object 'age
_1.2.1 ctv_0.7-4
sos_1.3-1 brew_1.0-6 Hmisc_3.9-2
survival_2.36-12
[9] MASS_7.3-17
loaded via a namespace (and not attached):
[1] cluster_1.14.2 grid_2.14.2lattice_0.20-0 tools_2.14.2
______
Professor Michael Kub
ather than
> "Depends". In this way, neither package is loaded along with effects.
>
> Best,
> John
>
> On Thu, 8 Mar 2012 20:45:37 -0500
> Michael Kubovy wrote:
>> Hi Spencer,
>>
>> In an Sweave script, which would take me too long to boil down
> packages and using namespaces, so users can record appropriately what they
> want in situations like this and not have to worry about name conflicts
> later?)
>
>
> Thanks,
> Spencer
>
>
> On 3/8/2012 4:19 AM, Michael Kubovy wrote:
>> Dear John,
>
this helps,
> John
>
>
> John Fox
> Sen. William McMaster Prof. of Social Statistics
> Department of Sociology
> McMaster University
> Hamilton, Ontario, Canada
> http://socserv.mcmaster.ca/jfox/
>
> On Wed, 7 Mar 2012 20:53:20 -
?
__
Professor Michael Kubovy
University of Virginia
Department of Psychology
for mail add: for FedEx or UPS add:
P.O.Box 400400 Gilmer Hall, Room 102
Charlottesville, VA 22904-4400 485 McCormick Road
USA
ering extended mode
\write18 enabled.
%&-line parsing enabled.
**test.tex
(./test.tex
LaTeX2e <2011/06/27>
Any thought? Please CC me when replying.
Michael Kubovy
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing li
kgs_from_ctv_or_repos(views = views, coreOnly = coreOnly, :
CRAN task view Graphics not available
Thanks,
______
Professor Michael Kubovy
University of Virginia
Department of Psychology
for mail add: for FedEx or UPS
the Zip and Unzip programs do not remedy
the problem. I am also running Windows 7.
Michael
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
On Jul 7, 2010, at 2:13 AM, Achim Zeileis wrote:
> On Tue, 6 Jul 2010, Michael Friendly wrote:
>
>> Michael Kubovy wrote:
>>> Suppose we start with
>>> data("Titanic")
>>> mosaic(Titanic, shade = TRUE)
>>> How do I combine the dashed
or bs(Pearson Resid) < 2 ?
Thanks,
Michael
__________
Professor Michael Kubovy
University of Virginia
Department of Psychology
for mail add: for FedEx or UPS add:
P.O.Box 400400
84
Done
Warning message:
closing unused connection 3 (content.xml)
When I try to open simpleOut.odt I get: "Format error discovered in the file in
sub-document content.xml at 22,124(row,col)." The file content.xml is empty.
_
On Nov 22, 2009, at 10:22 AM, Uwe Ligges wrote:
> masterinex wrote:
>> Hi guys , Im trying to do principal component analysis in R . There is 2
>> ways of doing
>> it , I believe. One is doing principal component analysis right away the
>> other way is standardizing the matrix first using s =
downloaded 1.0 Mb
The downloaded packages are in
/var/folders/TG/TG6nWu2pHYCI-BdWWiyMxU+++TI/-Tmp-//RtmpZx3zu8/downloaded_packages
_________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville,
y-data?
On Sat, Sep 5, 2009 at 8:03 PM, Michael Kubovy
wrote:
Hi,
Does anyone know how to do a coord_trans() in which the y-axis is
tranformed
into (for example) -1000/y?
Stephen Sefick
_________
Professor Michael Kubovy
University of Virginia
Departmen
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 400400Charlottesville, VA 22904-4400
Parcels
form as
mine:
data(bread)
res.uc <- smacofIndDiff(bread, itmax = 1000)
Thanks,
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McC
ale'.) In general it
seems to me that an interaction plot with lines is easier to read, and
not likely to mislead. Does anyone know if and where this has been
discussed?
_________
Professor Michael Kubovy
University of Virginia
Department of Psychology
There are quite a few others. See:
http://cran.r-project.org/web/views/ExperimentalDesign.html
On May 12, 2009, at 7:44 PM, Chris Howden wrote:
> I've spent the last week or so looking at all the experimental design
> packages I can find in R. AlgDesign, design.conf and BHH2 being the
> best o
amp;NA&100&NA\tabularnewline
5&0&5&&&&00&NA&000&NA\tabularnewline
6&0&6&&&&00&NA&011&NA\tabularnewline
7&0&7&&&&00&NA&000&NA\tabularnewline
8&0&8&&&&01&am
would like to retain the larger of the two (if considered a binary
number).
Can someone suggest an algorithm?
Thanks,
Michael
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:
15), ]
I would like to retain the larger of the two (if considered a binary
number).
Can someone suggest an algorithm?
Thanks,
Michael
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4
a summary after the three sections but tell LaTeK to
> publish it
> first? Possibly by rearranging the sections/pages?
> Currently I'm doing all the coding in a big block up front but
> that waste
> much of the advantage of weaving to me.
HTH,
MK
_
ns = 1, main = expression(pi/2))
plot(c2, stack = TRUE, bins = 1, main = expression(pi/2 + .1))
par(opar)
On Mar 29, 2009, at 8:09 AM, Duncan Murdoch wrote:
On 29/03/2009 7:39 AM, Michael Kubovy wrote:
require(circular)
c <- circular(rep(0, 20), zero = pi/2, rotation = 'clock
require(circular)
c <- circular(rep(0, 20), zero = pi/2, rotation = 'clock')
plot(c, stack = TRUE, shrink = 1.5)
Can anyone tell me why the stack is offset from 0?
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P
Take a look at the functionsnlm(), optim() in the stats package and
maxNR() in the maxLik package.
On Mar 22, 2009, at 11:15 PM, Roslina Zakaria wrote:
> Does R has a topic on newton's method?
_
Professor Michael Kubovy
University of Virginia
Depar
Thanks for the tip. I'll start with rJava, and see if that easier
approach works.
On Mar 20, 2009, at 10:50 PM, Martin Morgan wrote:
Michael Kubovy writes:
Dear R-helpers,
I run R version 2.8.1 (2008-12-22) on i386-apple-darwin8.11.1
I would like to call R from Java. I downl
rary/SJava'
Can someone tell me what I'm doing wrong? Since this is a matter of
some urgency, I would be grateful if you also cc-ed my email address.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
Postal Address:
P.O.Box 400
.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
Postal Address:
P.O.Box 400400, Charlottesville, VA 22904-4400
Express Parcels Address:
Gilmer Hall, Room 102, McCormick Road, Charlottesville, VA 22903
Office:B011;Phone
?gls
On Feb 16, 2009, at 12:28 PM, constantine wrote:
In other statistical software, such as Eviews, it is possible to
regress a model with the Least Squares method, assuming that the
residuals follow an AR(q) process.
For example the resulting regression is something like
y = 1.2154 + 0.2215
a smaller radius than the black
circle, running from 0 to 315 deg. I also thought that adding a short
horizontal line at its beginning might be helpful. I would appreciate
advice on how best to do this or anything else that would provide the
required information.
_________
Pro
Dear Friends,
require(cluster)
x <- rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)),
cbind(rnorm(15, 5, 0.5), rnorm(15, 5, 0.5)))
plot(pp <- pam(x, 2), which.plots = 1)
How can I extract the coordinates used in the plot?
_
Professor Michael Kubovy
Univers
n Jan 24, 2009, at 8:54 PM, Rolf Turner wrote:
>
> On 25/01/2009, at 2:33 PM, Michael Kubovy wrote:
>
>>
>> (1) If I type ?normal because I forgot the name dnorm() I get a long
>> list of relevant pages. Getting to right page is laborious.
>>
>> (2) If I remem
list of pages.
Advice?
It seems to me that if the output of help() listed base functions
first, it would go a long way toward improving the usefulness of this
function.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
Postal Address
7; being evaluated was:
(1, n.axes)
> coinertia(op.dudi, em3.dudi, n.axes = 2)
Error in sqrt(Dp) : Non-numeric argument to mathematical function
Can anyone tell me why the coinertia fails, and how to fix it? I would
be happy to send the two data files, if you were willing to help.
Thank
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 Virgini
"The ex-
GAUS density mu=300 ,sigma=35,nu=100")
But this does:
y<- rexGAUS(100, mu=300, nu=100, sigma=35)
m1<-gamlss(y~1, family=exGAUS)
curve(dexGAUS(y=x, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex-
GAUS density mu=300 ,sigma=35,nu=100")
__
e
efficient command).
_________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B0
c('red',
'blue'), lty = 1, xjust = 0.5, bty = 'n')
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville,
[compress,smaller]{beamer}
\documentclass{article}
\usepackage{beamerarticle}
***IT DOES NOT FAIL: *
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:
27;), each = 7), rep(c('pre', 'post'), each =
6)),
sub = c(1:7, 1:7, 8:13, 8:13))
xyplot(y ~ pp | sex, groups = sub, type = 'b', mwpp)
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 40
A cleaner approach than attach()---because you don't need to remember
to detach()---is:
my.df<-data.frame(cbind(x=runif(100), y=rnorm(100)))
head(my.df)
with(my.df, plot(y~x))
_
Professor Michael Kubovy
University of Virginia
Department of Psychol
0.249
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-squared: 0.07308,Adjusted R-squared: 0.02158
F-statistic: 1.419 on 1 and 18 DF, p-value: 0.249
Why not"g
OS details:
System Version: Mac OS X 10.5.4 (9E17)
Kernel Version: Darwin 9.4.0
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick Road
each =
3)))
xyplot(y ~ x | z, df, axis = axis.td, xlab = 'bottom', ylab = 'left',
layout = c(3, 1), aspect = 'xy', type = 'b')
Here I would like to have a different pair of top and bottom axes for
each panel (i.e., different labels at 1:3), and a differen
an appropriate key in each?
_________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011
You will find that many people are eager to help beginners, as long as
they give enough information about the problem they encountered.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-44
Perhaps
?mosaic
or
?mosaicplot
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office
groups=year,
method="filled bands" , data=s, type="l")
# this produces a bright red band at the top and a black band at the
bottom
I'll be happy to submit a bug report if this doesn't count and if you
tell me how.
_
Profess
plot2.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B019+1-434-982-475
Perhaps you need
with(workingdata, asin(BEFORE))
or
asin(workingdata$BEFORE)
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick
?boxplot
The boxplot does not show the mean but the median:
median(c(0,15,0,60,0,0,0,0,0,60,60,0,60,0,30,0))
> b <- boxplot(c(0,15,0,60,0,0,0,0,0,60,60,0,60,0,30,0))
b$stats
[,1]
[1,]0
[2,]0
[3,]0
[4,] 45
[5,] 60
_
Professor Michael
bout that.
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-43
Sorry, I meant to say: "For the moment I wonder if the solution is not
to use CIs based on the two low SEs produced by the ~ time model, and
to treat them as least-significant difference intervals."
_
Professor Michael Kubovy
University of Virginia
Dep
difference intervals.
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B01
ts on
the future direction of the debate on NHST are offered.},
url = {http://search.epnet.com/login.aspx?direct=true&db=pdh&an=met104397
}
}
_________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlo
JavaGD_0.5-1 rJava_0.5-1
The following opens a JavaGD window and freezes JGR:
super.sym <- trellis.par.get("superpose.symbol")
(from splom help page).
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Ch
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
4.1lme4_0.99875-9Matrix_0.999375-9 lattice_0.17-6
loaded via a namespace (and not attached):
[1] MASS_7.2-41 car_1.2-7 cluster_1.11.10 gdata_2.4.1
gtools_2.4.0
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS:
You might also take a look at the doBy package.
On Apr 14, 2008, at 7:30 AM, zerfetzen wrote:
> Thanks all. I will try to use both tapply and by, and have no idea
> how I
> missed the by function. Thanks again.
_
Professor Michael Kubovy
University of
g anything about search(), match() or pos).
As a former student of mine (who has long since gone into business)
once said: "I avoid learning experiences like the plague."
Thanks again,
Michael
On Mar 30, 2008, at 1:36 PM, Duncan Murdoch wrote:
> On 30/03/2008 10:06 AM, Michael
insure
that users won't get bitten by packages I had loaded but didn't do so
explicitly in the Sweave script, *and* it would solve the current
problem as well. Is there a straightforward way to do that?
On Mar 30, 2008, at 9:43 AM, Duncan Murdoch wrote:
> On 30/03/2008 8:43 AM,
yield+variety
+block+preds.sq)sq"
#
# BAD INTERACTION BETWEEN alr3 and AlgDesign
#
What now?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilme
ubset = subs, weights = wgts, na.action = naa)
Browse[1]>
Error in parse(text = x) :
unexpected symbol in "yield ~ variety + preds(yield+variety
+preds.sq)sq"
Advice?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville
= c(0.5, 3.5), panel = function(x, y){
panel.Dotplot(
x, y, col = mypal[c(1, 2, 1)],
col.segments = mypal[c(1, 2, 1)])
panel.abline(v = 50, col.line = 'red', lty = 2)}
)
_
Professor Michael Kubovy
Univers
Delightfully straightforward! Thanks.
On Feb 16, 2008, at 5:15 PM, Duncan Murdoch wrote:
> On 16/02/2008 4:51 PM, Michael Kubovy wrote:
>> Dear R-helpers,
>> > label2 <- expression(paste(italic(attraction function:), 'slope'))
>> Error: unexpected 'func
treated as a reserved word but as a string in italics?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottes
e some scholarly work on this. I can't
> remember
> where I came across it - perhaps I dreamed it up? Help, anyone?
>
> Best wishes
>
> Paul
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400
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 400400Charlott
On Feb 9, 2008, at 4:41 PM, Prof Brian Ripley wrote:
> On Sat, 9 Feb 2008, Michael Kubovy wrote:
>
>> How do I enter 'much greater than' and 'much less than' symbols in
>> an
>> expression?
>
> Those are not in the Adobe Symbol encoding used
How do I enter 'much greater than' and 'much less than' symbols in an
expression?
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 10
On Jan 28, 2008, at 11:30 AM, Liaw, Andy wrote:
> The original paper is:
>
> Cook, J. R. & Stefanski L. A. (1994) Simulation--extrapolation
> estimation in parametric measurement error models. Journal of the
> American Statistical Association, 89, 1314-1328.
>
?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B019+1-434-982-4751
Fax
hat presumably garbled the file. (Indeed, StatTransfer told me so.)
Would it be appropriate to request this as a feature of
foreign::read.dta()? (Thomas Lumley is the author.)
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville,
[1] "UTF-8" "EUC-JP"
What should he do to get to locale C?
I believe that he isn't subscribed to this list, so please include a
cc to him.
Thanks,
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 4004
Dear R-helpers,
On Jan 12, 2008, at 5:17 PM, Deepayan Sarkar wrote:
> On 1/12/08, Michael Kubovy <[EMAIL PROTECTED]> wrote:
>> Dear r-helpers,
>>
>> Does anyone have a straightforward example of putting together three
>> unrelated (expect for a common y-axis) xyp
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
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B019+1-434-982-4751
Fax:
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
el of the graphical parameters
accessible by par(), at
http://pbil.univ-lyon1.fr/R/fichestd/tdr75.pdf
On Dec 6, 2007, at 2:17 AM, Dieter Menne wrote:
> Michael Kubovy virginia.edu> writes:
>
>>
>> Dear r-helpers,
>>
>> How do I tell xYplot() and Key() t
, other = list(title = expression(italic(v)),
cex.title = 1,
labels = c(expression(italic(b)), expression(italic(c)),
expression(italic(d)
dev.off()
This works only for the title.
_________
Professor Michael Kubovy
University of Virginia
Department of Ps
ysis of cross-over trial},
publisher = {Chapman \& Hall/CRC},
year = {2003},
address = {Boca Raton, FL, USA},
edition = {2nd},
author = {Byron Jones and Michael G. Kenward}
}
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS:
e that I might adapt to this purpose?
(2) Can anyone recommend alternative ways of computing such profile
similarities?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Ro
ead? How?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesville, VA 22903
Office:B011+1-434-982-4729
Lab:B019+1-434-982-
teful for a pointer toward the solution.
Thanks,
MK
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick RoadCharlottesvill
r unique levels
How do I then get the character vector of length r of the levels of
fact1?
Once I have that, is there anything else I need to know about
swapping labels between factors in a df?
_
Professor Michael Kubovy
University of Virginia
Department of Psychol
solidate these two in (2) and make
(2) the default library (if that's adivsable) and how to make the
upgrade to future versions of R as trouble-free as possible.
_____
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P
87 matches
Mail list logo
