At 03:42 PM 26/10/2005, you wrote:
>We can double check this with yacas which is a free
>computer algebra system that supports arbitrary
>precision math:
>
>In> N(Exp(-1000))
>Out> 0.5075958898e-434
Or in Axiom, which is open source too:
numeric exp(-1000)
0.50759588975494567653E-434
Check out axis.break in the plotrix package.
On 10/26/05, Bill Hunsicker <[EMAIL PROTECTED]> wrote:
> R-Help,
>
>
>
> I am trying to do simple plots of the characteristics of cellular
> phones. I am values that fit along an axis that has many data points
> around 800 and 1800. I am not intereste
Since log base 10 is log(x)/log(10) the exponent is:
floor(log(exp(-1000))/log(10))
= floor(-1000/log(10))
which in an R session gives:
> floor(-1000/log(10))
[1] -435
and the mantissa is:
exp(-1000)/10^(-435)
= exp(-1000)/exp(-435*log(10))
Dear All,
I am a novice user of the R software package. When I try and compute,
exp(-1000) or exp(-2000), i get the answer as zero. Is there any way i can get
R to compute the answer and give me an actual number ? ( by increasing the
precision or any other method).
If I cannot get R to give me
R-Help,
I am trying to do simple plots of the characteristics of cellular
phones. I am values that fit along an axis that has many data points
around 800 and 1800. I am not interested in the "dead space' between
the two clusters of data.
I am not interested in a linear axis. I would li
The Statistical Computing Section of the American Statistical
Association announces the competition for the John M. Chambers
Statistical Software Award. In 1998 the Association for Computing
Machinery presented its Software System Award to John Chambers for the
design and development of S. Dr. Cham
> ?which
> X <- matrix(2:10, nrow=3)
> which(X == 10, arr.ind=TRUE)
row col
[1,] 3 3
HTH,
simon.
At 11:52 AM 26/10/2005, you wrote:
>Hi netters,
>
>This is probably a silly question,but I can't find the answer after
>searching the R-help archives online. ok, I have a matrix. I know
Hi netters,
This is probably a silly question,but I can't find the answer after
searching the R-help archives online. ok, I have a matrix. I know there is
a "10" somewhere in it. Now I want to
know the index of the element "10" in this matrix. That is, if X[i,j]=10, I
want to know
i and j. Is
Hi Roger and the list,
The package is working very well. What surprised me most is the speed. As I
mentioned in my previous emails, I have to find the neighbors for around
200,000 individuals. It took no more than 10 minutes for the function to
finish the searching and returned enough infor
Hi,
I am trying to use the R package 'spatstat' for generating spatial
poisson point process graphs. I can create a point pattern using the
following commands:
pp <- rpoispp(.01, win=owin(c(0,100),c(0,100)))
and also view the resulting graph by:
plot(pp)
Dear Jen,
There are lots of different ways to do what you want -- you've already had
two suggestions -- but you might consider whether you really need to do it.
In particular, R will generate its own indicator variables (and other kinds
of contrasts) in linear and other statistical models (see Sec
On Tue, 25 Oct 2005 [EMAIL PROTECTED] wrote:
> Hello everyone,
>
> I have a question regarding building the Rmysql packages from source in
> Windows.
Pretty esoteric, and the posting guide suggests you ask the maintainer.
> I am currently using: Windows XP, R-2.2.0, Mysql 3.23
MySQL 3.23 is r
Hello there,
Suppose you want to solve the following system of ODE's (a simple
Lotka-Volterra predator prey model)
dP/dt = beta*P*V - mu*P
dV/dt = r*V - beta*P*V
where P and V are the numbers of predators and prey. Now, this is
easy to do, but suppose you have a system of equations like this
I'm fairly new to R and am wondering if anybody knows of R code to
calculate confidence intervals for parameters (fixed effects and variance
components) from mixed effects models based on Sattherthwaite's method?
I'm also interested in Satterthwaite-based confidence intervals for linear
combina
x<-c(1, 2);
y<-ifelse(x==1, 1, 0);
On 10/25/05, Jennifer Dillon <[EMAIL PROTECTED]> wrote:
>
> Hello,
>
> I am almost a total novice, and I am sure there must be an easy (and
> basic) way to turn a variable of 1's and 2's into a variable of zeros
> and ones. This is in a data frame, but if I could
Hello everyone,
I have a question regarding building the Rmysql packages from source in
Windows.
I am currently using: Windows XP, R-2.2.0, Mysql 3.23
I have installed all of the tools from Installation and Administration
and after following all the instructions to the letter, I compiled the
Thank you for Prof.Ripley's suggestion. I fixed the program by adding a
lower scope, and the program ran, but I still got warning messages, and
don't know what is going on, would this affect my results?
...
Step: AIC= 12337.74
Surv(tlfup, cen) ~ MI[[j]]$trt + MI[[j]]$agem40 + MI[[j]]$agem40sq +
Le 25 Octobre 2005 15:51, Jennifer Dillon a écrit :
> Hello,
>
> I am almost a total novice, and I am sure there must be an easy (and
> basic) way to turn a variable of 1's and 2's into a variable of zeros
> and ones. This is in a data frame, but if I could do it with vectors,
> that's all I need.
I addition to all the many good options presented you can also use
trimWhiteSpace{limma} which is just a higher level call to sub()
x=scan("clipboard", what="character")
Read 6 items
x
[1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
trimWhiteSpace(x)
[1] "AIR" "ABCB" "ABXA" "ACMR" "ADCT" "A
Welcome to R. See
?merge
then
?aggregate
or
require(Hmisc)
?summarize
or
?by
You can probably find many examples in the archives, if needed.
> -Original Message-
> From: Michael Graber [mailto:[EMAIL PROTECTED]
> Sent: Tuesday, October 25, 2005 3:45 PM
> To: R-Mailingliste
Em Ter 25 Out 2005 15:24, Doran, Harold escreveu:
> There is an issue with implicit nesting in lmer. In your lme() model you
> nest block/irrigation/density/fertilizer. In lmer you need to do something
> like (I dind't include all of your variables, but I think the makes the
> point)
>
> lmer(yield
Hello,
I am almost a total novice, and I am sure there must be an easy (and
basic) way to turn a variable of 1's and 2's into a variable of zeros
and ones. This is in a data frame, but if I could do it with vectors,
that's all I need.
Can someone tell me how?
Thanks so much,
Jen
__
Dear R-List,
I am very new to R and programming itself, so my question may be easy to
answer for you.
I tried a lot and read through the manuals, but I still have the
following problem:
I have 2 data-frames:
Number<-as.numeric (Number)
Name<-as.character (Name)
TAB1<-data.frame (Name,Number)
- it
Thanks for your replies, I was unaware of these solutions.
On 10/25/05, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Does anyone know how I can set up R so that when I make a graphic, the
> graphics window remains behind the console window? It's annoying to
> have to reach for the mouse every ti
[EMAIL PROTECTED] wrote:
> Does anyone know how I can set up R so that when I make a graphic, the
> graphics window remains behind the console window? It's annoying to
> have to reach for the mouse every time I want to type another line of
> code (e.g., to add another line to the plot). Thanks.
I
On Tue, 2005-10-25 at 13:07 -0500, Marc Schwartz (via MN) wrote:
> On Tue, 2005-10-25 at 11:55 -0600, [EMAIL PROTECTED] wrote:
> > Does anyone know how I can set up R so that when I make a graphic, the
> > graphics window remains behind the console window? It's annoying to
> > have to reach for the
> Does anyone know how I can set up R so that when I make a graphic, the
> graphics window remains behind the console window? It's annoying to
> have to reach for the mouse every time I want to type another line of
> code (e.g., to add another line to the plot). Thanks.
What OS? In Windows with R
Thanks, that will certainly do the trick. I came across systemfit and sem -
packages that are useful in this respect as well.
Jukka Ruohonen
University of Helsinki
> Have you received a reply? I haven't seen one. RSiteSearch("panel
> data unit root test") produced another question on this but
On Tue, 2005-10-25 at 11:55 -0600, [EMAIL PROTECTED] wrote:
> Does anyone know how I can set up R so that when I make a graphic, the
> graphics window remains behind the console window? It's annoying to
> have to reach for the mouse every time I want to type another line of
> code (e.g., to add ano
Does anyone know how I can set up R so that when I make a graphic, the
graphics window remains behind the console window? It's annoying to
have to reach for the mouse every time I want to type another line of
code (e.g., to add another line to the plot). Thanks.
___
Dear all, someone could send me a introductory reference about Survival
Frailty Models???
Thank you very much
_
Chegou o que faltava: MSN Acesso Grátis. Instale Já!
__
R-help@stat.math.et
There is an issue with implicit nesting in lmer. In your lme() model you nest
block/irrigation/density/fertilizer. In lmer you need to do something like (I
dind't include all of your variables, but I think the makes the point)
lmer(yield~irrigation*density*fertilizer+(1|fertilizer:density)+(1|de
I'm using lm to run an ANOVA, and would like to use Ryan's Q as my
post-hoc (as recommended by Day and Quinn, 1989, Ecological
Monographs). I can't seem to find any methods in the base stats
package that implement this post-hoc. Is there a good package of
post-hoc methods out there, or has so
Em Seg 24 Out 2005 18:08, Doran, Harold escreveu:
> Ronaldo
>
> See the article on lmer pasted below for syntax. It is the only current
> source documenting the code. In lmer(), the nesting structure for the
> ranmdom effects is handled in a slightly different way. If your
> observations are nested
Sounds like you need to debug the software yourself. ?traceback ?browser
may help. See also Mark Bravington's debug package on CRAN.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
"The business of the statistician is to catalyze the scientific learning
process." - Geo
Hello,
I'm having some difficulty running Niels Waller's Maxcov Hitmax program in
the R console, and I was hoping you could provide me with some assistance.
When I attempt to run the analysis I receive the following message:
Indicators 1 3 & 2Error in if (del == 0 && to == 0) return(to) :
missi
"Guenther, Cameron" <[EMAIL PROTECTED]> writes:
> Does anyone know if there is a file size limit for inputting a SAS data
> set? The file size I am trying to import is 184 MB.
> The code is:
>
> >library("foreign")
> >
> > sashome<-"C:/Program Files/SAS Institute/SAS/V8"
> > input.data<-rea
- Weitergeleitete Nachricht von Patrick Kuss <[EMAIL PROTECTED]> -
Datum: Tue, 25 Oct 2005 17:28:33 +0200
Von: Patrick Kuss <[EMAIL PROTECTED]>
Antwort an: Patrick Kuss <[EMAIL PROTECTED]>
Betreff: running AMOVA from spreadsheet genotype data
An: [EMAIL PROTECTED]
Hi,
I p
Actually, R does not accept NA. First, the action is from lm and not from
the R base, and second, NAs are in fact deleted and not accepted.
Function lm() has a na.action argument which is by default na.omit. You
could easily make na_and_inf.omit and use it. You could also use
subset=(w > 0). T
Hi,
I plan to run amova in R from randomly generated subsets of my original dataset.
The structure of my data looks as below:
region <- rep(c("east","west"),each=8)
pop <- c(rep(1:4,each=4))
ind <- c(rep(1:4,4))
l.1 <- c(1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0)
l.2 <- sample(c(0,1),16,replace=T)
l.3 <- s
you could use:
lm( y ~ x + z + log(w), subset = w > 0)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web
Thanks, that will certainly do the trick. I came across systemfit and sem -
packages that are useful in this respect as well.
Jukka Ruohonen
University of Helsinki
> Have you received a reply? I haven't seen one. RSiteSearch("panel
> data unit root test") produced another question on this but
Hi,
Suppose I I wish to run
lm( y ~ x + z + log(w) )
where w assumes non-negative values. A problem arises when w=0, as log(0)
= -Inf, and R doesn't accept that (as it "accepts" NA). Is there a way to
tell R to do with -Inf the same it does with NA, i.e, to ignore it? (
Otherwise I have to do
What about
xx <- c("AIR ", "ABCB ", "ABXA ", "ACMR ", "ADCT ", "ADEX "," AAA")
xx=gsub("[[:blank:]]","",xx)
Stefano
On Tue, Oct 25, 2005 at 03:14:39PM +0200, Dimitris Rizopoulos wrote:
one way is to use strsplit(), i.e.,
xx <- c("AIR ", "ABCB ", "ABXA ", "ACMR ", "ADCT ", "ADEX ")
Best,
Dimi
Hi Yogesh,
try modifying the graphical parameters like
> par(cex=2,cex.axis=1.2)
for graphs in base:
or for e.g. trellis graphs in the "lattice" package the trellis parameter´s
(try ?trellis.par.get)
axis.text <- trellis.par.get("axis.text")
axis.text$cex <- 3.2
trellis.par.set("axis.text",
Dear all,
I am exporting to latex a matrix via xtable. My matrix includes e.g. 4 rows
and 4 columns. The first two rows containing real values(e.g. a laboratory
parameter's mean value), and the two rows at the bottom containing only
integers (number of cases).
>
mymat<-matrix(c(c(10.52,2.52,12
If you are interested in trimming leading and trailing whitespace (and not
interior whitespace), you can use the 'trim' function from the 'gdata'
package.
-G
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of roger bos
> Sent: Tuesday, October 25, 2005
Does anyone know if there is a file size limit for inputting a SAS data
set? The file size I am trying to import is 184 MB.
The code is:
>library("foreign")
>
> sashome<-"C:/Program Files/SAS Institute/SAS/V8"
> input.data<-read.ssd( file.path("G:/DATA/Cam/ECPATH/FIM"),"tbm_c",
+sas
Le 25.10.2005 14:59, Marc Schwartz a écrit :
>If you specifically need the plot to have a dimension measured in
>pixels, then you need to use a bitmapped format such as png and specify
>the output to be the size you require:
>
> png("test.png", width = 300, height = 300, ...)
> DoYourPlotHere()
Hello R Users,
How we can save any data which is computed in
side R . For example I am fitting a fourth
harmonic function on any time series data. We
can plot this function on each other with
timeseres data in R but suppose I want to
save this fitted function values as one
column then how to do i
Hello R Users,
I am a new user to R.
When we make a simple plot then how to
increase the font size of the data at the
each axis.
Many thanks,
Kind regards,
Yogesh
--
===
Yogesh Tiwari,
Max-Planck Institute for Biogeochemistry,
Hans-Knoell Strasse 10,
D-
roger bos wrote:
> for example:
>
>>a$tic[1:10]
>
> [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
> [7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
Not beeing an expert for regular expressions, I'd try
sub(" *$", "", sub("^ *", ""
one way is to use strsplit(), i.e.,
xx <- c("AIR ", "ABCB ", "ABXA ", "ACMR ", "ADCT ", "ADEX ")
unlist(strsplit(xx, " "))
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgiu
Try trim() in gdata package.
Wuming
On 10/25/05, roger bos <[EMAIL PROTECTED]> wrote:
> for example:
> > a$tic[1:10]
> [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
> [7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
> Thanks,
> Roger
>
Thanks to everyone for the gsub function. The white spaces were generated by
reading in the file, but its a fixed width file and I only need some of the
columns, so this was the easiest way to do it.
Thanks,
Roger
On 10/25/05, Marc Schwartz <[EMAIL PROTECTED]> wrote:
>
> On Tue, 2005-10-25 at 0
Le 25.10.2005 14:51, roger bos a écrit :
>for example:
>
>
>>a$tic[1:10]
>>
>>
>[1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
>[7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
> Thanks,
> Roger
>
>
gsub(' ','',a$tic)
--
visit the
On Tue, 2005-10-25 at 08:51 -0400, roger bos wrote:
> for example:
> > a$tic[1:10]
> [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
> [7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
> Thanks,
> Roger
Roger,
See the next to last set of
Quoting roger bos <[EMAIL PROTECTED]>:
> for example:
>> a$tic[1:10]
> [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
> [7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
> Thanks,
> Roger
>
> [[alternative HTML version deleted]]
>
> ___
On 10/25/05 8:51 AM, "roger bos" <[EMAIL PROTECTED]> wrote:
> for example:
>> a$tic[1:10]
> [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
> [7] "ABM " "AFCE " "AG " "ATG "
> Can anyone please tell me how to strip the white spaces from a$tic?
> Thanks,
> Roger
>
> [[alternative HTML version d
If you specifically need the plot to have a dimension measured in
pixels, then you need to use a bitmapped format such as png and specify
the output to be the size you require:
png("test.png", width = 300, height = 300, ...)
DoYourPlotHere()
dev.off()
Do this directly using the png() device
for example:
> a$tic[1:10]
[1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX "
[7] "ABM " "AFCE " "AG " "ATG "
Can anyone please tell me how to strip the white spaces from a$tic?
Thanks,
Roger
[[alternative HTML version deleted]]
__
R-help@sta
Dear all,
I did post a message about problems on installing a package (gstat)
which ends with a compilation error:
** libs
gcc -I/usr/lib/R/include -mieee-with-inexact -fPIC -g -O2 -c block.c
-o block.o
gcc -I/usr/lib/R/include -mieee-with-inexact -fPIC -g -O2 -c
chfactor.c -o chfact
Please send R and mgcv version numbers, as I can't replicate this
problem. best, Simon
>
>Hello, I'm a Korean researcher who have been started to learn the "R"
>package.
>
>I want to make gam model and AIC value of the model to compare several
>models.
>
>I did the GAM model,
On Mon, 2005-10-24 at 09:01 -0400, Suresh Kumar Karanam wrote:
> Hi,
> I have a dataset with four categories of data, the number of samples are not
> the same in each category. I want to find the Spearaman's Rho. Let me give
> an example.
> x=(14.22770439,26.49420624,46.7277932,19.02550707,23.373
Thanks Marc and Jim for the tips. The PDF file that I create with R looks
about the same as the one you created. However, I need to get the graphic
to be a certain size (300 pixels wide). I have been using the ImageMagick
program to do so for other graphics:
convert test.pdf -resize 300x300 out.pd
65 matches
Mail list logo