=black,
ylab = list(expression(paste(dN/dlogD[agg], [*cm^-3*]))), xlab =
expression(paste(D[agg], [nm])),
cex.lab=1.2
)
axis(2, mgp=c(0, 0.2, -2))
dev.off()
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original
(4.5, 4.5, 1, 1) + 0.1)
###
plot(D_mean,y,type=l,yaxt=n,lty=2,lwd=2,col=black,
ylab = list(expression(paste(dN/dlogD[agg], [*cm^-3*]))),
xlab = expression(paste(D[agg], [nm])),
cex.lab=1.2
)
axis(2, mgp=c(0, 0.2, -2))
dev.off()
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561
bookseller, including the nearest university bookstore.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Donatas G.
Sent: Wednesday, July 11, 2007 9
Jessie:
How many pixels would you need to allocate for each of these 1000 parts? Is
that feasible?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED
for legend (?legend) and modify it, or use parts of
it for your own function.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Maja Schröter
Look at the definition for the transform. For example in the car package,
?box.cox
Then do the simple algebraic manipulations yourself.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL
plot( x=rnorm(25, 0.5, 0.3), y=rnorm(25, 4, 1), xlim=c(0,1), ylim=c(2,7))
# ^^
for example
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original
You might consider having R do everything: R can read the Excel sheet, do
what needs to be done, and write the results to an Excel sheet.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of LL
Sent: Saturday, May 26, 2007 1:14 PM
To: LL; r-help@stat.math.ethz.ch
Subject: [R
Perhaps format is what you are looking for:
?format
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Pietrzykowski,
Matthew (GE
on the icon, as usual.) To do that, right-click on the program's icon and
choose Run as administrator. Perhaps this will solve your RWinEdt
problems as well.
Best wishes!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original
Methods, Springer, 1999
15. Venables and Ripley, Modern Applied Statistics with S, 4th ed.,
Springer, 2002
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto
Sure.
rpois(n, lambda)
... will do it. But you should tell us something about how you want your
numbers to be distributed, since rpois() produces integers having a Poisson
distribution.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
Or you can right-click on the R icon and choose Run as administrator. That
way you won't alter the security settings and forget to re-set them. After
the packages are installed R will load in the usual way by clicking on the
icon.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax
for
implementing the statistical thinking you've acquired.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Paul Lynch
Sent: Thursday, February 15, 2007
Thank you Dan!
I dunno if I would have ever found that.
Thanks!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED
(expression(theta),cex=5),strheight(expression(
theta),
cex=5)),nrow=1)
symbols(0.5,0.5
,rectangle=dimensions,bg='white',fg='white',add=TRUE,inches=FALSE)
text(0.5,0.5,expression(theta),cex=5)
~
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax
), then install the libraries, then re-set the security levels.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Sicotte, Hugues Ph.D.
Sent
It might be helpful to those not familiar with Matlab to tell us what
function persistent does.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED
place.
Is there a way that I can use R CMD build --binary, with its zip file
output, but substitute my .First.lib function for the default?
Thanks for any help.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
--binary.
Should I just declare a namespace and use onLoad, or is there a mechanism to
use my own .First.lib?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL
, loaded the new package and there were my
new menus.
Thanks again.
Sheepishly,
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Charles Annis
.
Still struggling.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley
Sent: Thursday, November 16, 2006 4:44 AM
.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Benilton Carvalho
Sent: Thursday, November 16, 2006 12:38 AM
To: [EMAIL PROTECTED]
Cc: r-help
in the tarball it opens and shows me what I
expected.
I had hoped that I had weathered the hard part - building the package - but
I still need some help:
1) How do I get a zipped file, rather than a tarball,
2) How do I install what I've built?
Thanks in advance.
Charles Annis, P.E.
[EMAIL
_Practical Programming in Tcl
and Tk_ but would appreciate any help from fellow R users.
Many thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing
Thank you, Peter. The code provides a great example of packing the widget
nicely too.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf
of
available colors and available symbols differ, the color/symbol combinations
will be unique for the recycled usage.
Thanks for your counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
4808 0.86212 0.0066964 S
3 3 2084 0.93579 0.4407900 S
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Chris Fonnesbeck
Sent
)
tcl.tk.example(PROBLEM = FALSE)
to turn the problem on and off.
I have also attached a tcl.tk.example.txt file since these long inclusions
within an R-help note often have the line feeds missing.
Thank you for your insights.
Charles Annis, P.E.
###
tcl.tk.example - function
D'oh!
You are right of course. How very obvious these things are on the other
side of the solution. I spent the entire weekend and completely missed it.
Thanks!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
, variable=buttonValue, value=TRUE)
tkgrid(tklabel(frame.2, text=Question?yes ), answer.button1,
sticky=e)
^^^
I know this can not be as hard as I am making it.
Thanks for the help.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352
/RColorBrewer/chtml/RColorBrewer.chm'
Any help would be appreciated. Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
the restarts are not onerous.
Many thanks to you all!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Clint Bowman
Sent: Friday
?
result - repeated.measures.FACTOR.names[[1]]
for(i in 2:length(repeated.measures.FACTOR.names)) {
result - paste(result, repeated.measures.FACTOR.names[[i]], sep=,) }
result
[1] Insp1,Insp2,Insp3,Insp4,Insp5,Insp6,Insp7,Insp8,Insp9
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
] default.FACTOR.labels - c(\Probe1\, \Probe2\, \Probe3\))
^ ^ ^ ^ ^ ^
which gives me \ rather than
Is it possible to escape the character?
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
D'oh!
I've been using cat() but somehow never got the bigger picture.
Thanks!!!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]
Sent: Wednesday
.
Thanks.
Charles Annis, P.E.
PS - Yes, I have upgraded to R version 2.4.0 (2006-10-03) but do not want to
change in the middle of a project.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help
.)
You may be looking for the resulting confidence bounds on the glm fit for
which I also have code that iteratively interrogates the loglikelihood
surface without plotting it.
If any of this is interesting, please send me a note so we won't clog the
bandwidth.
Charles Annis, P.E.
[EMAIL
counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
Eternal thanks to Jim Holtman and to Gabor Grothendieck who pointed me to
this concise piece of code:
cls - function() {
require(rcom)
wsh - comCreateObject(Wscript.Shell)
comInvoke(wsh, SendKeys, \014)
invisible(wsh)
}
Perfect!
Thanks!
Charles Annis, P.E
=n)
text(0, .2, quote(H[0]:~~mu[1]==mu[2]))
text(3, .2, quote(H[1]:~~mu[1]==mu[2]+delta))
#
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL
with R_ (Springer 2002). The
writing is relaxed and succinct, not condescending as some texts might
appear to a newcomer. It's just a good book.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
How close do you think it should be, given finite resolution with digital
computing?
acos(0.5) - pi/3
[1] 2.220446e-16
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED
) exped
5.2394254 0.9294818
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Dimitrios Rizopoulos
Sent: Sunday, September 17, 2006 3
Under Windows mine is located here
C:\Program Files\R\R-2.3.1\library\base\R
The file name, however is not .Rprofile, but rather Rprofile
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
anticipate it.
My code will be used on machines with varying graphics (and memory)
capacity. Is there a way I can check the native resolution of the machine
so that I can make adjustments to my code for the possible limitations of
the machine running it?
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of bogdan romocea
Sent: Monday, August 28, 2006 10:55 AM
To: [EMAIL PROTECTED]
Cc: r-help
)))
The DOS window opened, some magic occurred in the blink of an eye, and the
DOS window closed. I haven't the foggiest idea what to do next since I can
see no evidence of having done anything.
Thanks for your patience.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
, but I don't know what.
Agin thanks for your patience.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Monday, August 28, 2006 2:21 PM
Using sub rather than gsub appears to have no effect. The DOS window opens,
then closes having had no noticeable effect.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Gabor Grothendieck
generates a
crowded jpeg on the older machine. My hope was to check the capability of
the current machine and make some modest adjustments to my R code to produce
an acceptable jpg on any machine running it.
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
Romain:
a - tempfile()
cat('htmlscript type=text/javascript document.write(screen.width) ;
/script/html', file=a)
browseURL(a)
The object a was created, but no browser opened.
ls()
[1] a
a
[1] C:\\DOCUME~1\\CHARLE~1\\LOCALS~1\\Temp\\RtmpRgWrqb\\file678418be
Charles Annis, P.E
Please visit the R site http://www.r-project.org/ and search the mailing
list for paste expression We discussed the topic recently. The trick
is that you don't paste expressions, you make an expression containing
paste.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455
some kind soul please instruct me (and perhaps subsequent searchers)
how to convert the elements of a string into numbers?
Thank you.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
.
Many thanks to all.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Charles Annis, P.E.
Sent: Saturday
This'll work.
theta - 2.1
plot(NA, xlim=c(0,1), ylim=c(0,1), xlab=bquote(theta == .(theta)),
ylab=bquote(theta == .(theta)), main=bquote(paste(Results for ,theta ==
.(theta
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
,
but closing the file with dev.off() to produce the first graphic means that
I can't add to that and produce the subsequent graphic and file.
Is there a way to do this without generating the entire plot from the
beginning for each graphic?
Thanks
Charles Annis, P.E.
[EMAIL PROTECTED]
phone
Wow! Exactly what I needed.
Thank you, Professor Ripley!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley
Sent
of a string
with \\ separators.
Can I get there from here? (I've looked through previous R-help listing of
similar problems but if it's there, I missed it.)
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
to
end before I can use the new info, which is correctly place in the
workspace, in subsequent R routines.
Is there a way I can use the updated values in the same routine that created
the widget?
Thanks for your advice and patience.
Charles Annis, P.E.
PS - I did read Prof. Ripley's post of Wed
with 2
gig RAM.
Many thanks for your counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
The trick is not to paste expressions but to make an expression of the
pastes:
composite.expression - expression(paste(f(x | , alpha, , , beta, )))
dev.off()
par(mar=c(5,5,1,1)+0.1)
plot(NA,xlim=c(0,1),
ylim=c(0,3),xlab=x,ylab=composite.expression,lwd=3,type=l,cex.lab=1.3)
Charles Annis
Why create an enormous matrix? Why not read each company's info and
immediately write it to the file using write.csv( ... append = TRUE ...)?
^
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
with the Department of Statistics at the University of South Carolina, but
as of May 2006 with the Department of Statistics at Texas A M University,
http://www.stat.sc.edu/~west/javahtml/Histogram.html and
http://www.stat.tamu.edu/~west/
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614
happy with the defaults.) And the line end, lend=2, makes the ends square.
Best wishes.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf
.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R
Many thanks to Duncan Murdoch, Thomas Lumley, Patrick Burns, and Seth
Falcon, for the illuminating advice, which will be found in the R-help
archives.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
cex=0.8 (or whatever you like) as in
legend(0.3,0.4,c(name1,name2,...),col=1:20,lty=1:20, cex=0.8)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED
Please give us an EXAMPLE of the loop you have in mind. (It's likely that
you can use simpler methods than a loop, but without an example we'd be
guessing.)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original
Vielen Dank, Uwe!
I dunno how I missed it. Looking only at the demo(plotmath), I guess.
Thanks again!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
the same but how they're constructed is rather
different.
Also, have a look at
?substitute
and
?bquote
Best Wishes,
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto
(0,1.6), xlab=, ylab=, tck=-0.015)
text(1.6, 1, expression(paste(slope = , frac(paste( f( , hat(theta),
)), hat(theta))), sep=), adj=0)
text(2.25, 1.03, \\pd, vfont=c(sans serif,plain))
text(2.3, 0.934, \\pd, vfont=c(sans serif,plain))
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561
censored, 3=interval censored.
Although unusual, the event indicator can be omitted, in
which case all subjects are assumed to have an event.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
Go to the R site and look. http://www.r-project.org/
There has been some recent traffic on this topic.
Click on search: Tobit
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL
You could substitute
lines(kernelgraf)
for your last line:
points(kernelgraf, xlab=Probability, xlim=c(0,100), ylim=c(0,.1),
col=rgb(0,0,1), main=)
See
?lines
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
website.
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Kort, Eric
Sent: Thursday, January 26, 2006 2:39 PM
To: Thomas Kaliwe; r
Please don't forget engineering! (e.g. fatigue and reliability - censored
regression and survival; quantitative nondestructive evaluation - GLM)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
You can do something similar with Microsoft's browser but it isn't quite as
easy as Foxfire:
Right-click on the frame and choose Properties. Then highlight and copy the
URL and paste into the address window and click Go.
Then save the page.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561
)
expected.cdf - ((1:n.points)-0.5)/(n.points)
qqplot(my.data, expected.cdf, las=1)
# Use the interocular trauma test for goodness-of-fit:
my.lm - lm(expected.cdf ~ sort(my.data))
abline(coef=coef(my.lm), lty=2)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
Your problem sounds like it could be modeled with logistic regression
whereby the propensity for one result or another is linked to the factors
that control it. Logistic regressions are a special case of generalized
linear models. Look at ?glm
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561
(data.df$x, data.df$y, log=xy)
This plot suggests that while Experiments B and C might have a similar
relationship between x and y, Experiment A differs.
Since I know nothing of the physical meaning of these observations I am
unqualified to comment further.
Best wishes.
Charles Annis, P.E.
[EMAIL
anyone help?
Copious Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Sundar Dorai-Raj [mailto:[EMAIL PROTECTED]
Sent: Thursday, November 10, 2005 5:00 PM
To: [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
Subject: Re
machine.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of jon butchar
Sent: Thursday, October 27, 2005 8:02 PM
To: Srinivas Iyyer
Cc: r
Works well with both IE and Firefox on my 2 year old DELL WinXP machine.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Gabor
Try
plot(Day, gene1)
lines(Day, gene2)
see
?lines
for more details.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of IAIN
= tree, logi.mod = 1, type =
dit, boxp = TRUE, rug = TRUE, las.h = 1)
which does not use the histograms but instead uses dit plots to provide a
helpful, visceral feel for the behavior of the observations.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
$coeff
(Intercept)altitude
13.43360163 -0.01220884
mod3$coef[2]
altitude
-0.01220884
as.numeric(mod3$coef[2])
[1] -0.01220884
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
Edition. Springer, 2002. ISBN 0-387-95457-0.
http://www.stats.ox.ac.uk/pub/MASS4/. This book is more demanding and
covers a broad spectrum of contemporary statistical practice.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
_Monte Carlo Statistical Methods_
by Christian P. Robert, George Casella
Springer, 2nd ed 2005
This book (I have edition 1) is a dandy. It will be rough sledding unless
you have a reasonable background in math stats but I think it is just what
you are looking for.
Charles Annis, P.E.
[EMAIL
?qqplot
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED]
Sent: Friday, July 15, 2005 9:03 PM
To: r-help
x - rep(NA, 3)
for (i in 1:length(x)){
x[i] - ...
}
will do the job, but you may be able to take advantage of R's vectorization
and do what you want with no loop at all.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
reconcile the differences.
I believe that the confidence bounds for both models should agree. After
all, both calls to survreg() produce identical parameter estimates.
So I have missed something. Would some kind soul please point out my error?
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone
included an explanation of what is happening at each iteration.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Paul Smith
Sent: Thursday, April
Vito:
Please plot your data:
y - c(1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0)
x - c(37, 35, 33, 40, 45, 41, 42, 20, 21, 25, 27, 29, 18)
plot(x, y)
You will see that ANY step function between 29 x 33 will describe these
observations perfectly.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561
for them.
Best wishes.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of
[EMAIL PROTECTED]
Sent: Thursday, January 20, 2005 10:16 AM
To: r
=agresti1998
when searching for approximate is better than exact
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Spencer Graves
Sent
pt(q, ...) returns the area to the left of q. The area to the left of 2.23
for your situation is 0.975, while the area to the left of -2.23 (which is
on the left side of zero from 2.23) is 0.025.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
Yes. Assuming you are using R in Windows, right click on the plot, and
choose Save as metafile, or more expeditiously, Copy as metafile, and then
paste directly into your WORD document.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
1 - 100 of 117 matches
Mail list logo
