I have been authoring a very small R package on CRAN, named
normix which implements an S3 class norMix has plot and
print methods; further, E[X] and Var[X] methods, random number
generation (r) and density evaluation.
It also provides the 16 Marron-Wand densities (known in the (1d)
density
For a more detailed answer try the American Statistician:
Title: The joy of copulas: Bivariate distributions with uniform marginals
(Com: 87V41 p248)
Author(s): Genest, Christian;MacKay, Jock;
Keywords: [Teacher's Corner];Frechet bound;Kendall $\tau$;Singular
distribution function;Kendall tau
In Jim Lindsey's repeated package,
`gausscop'Multivariate Gaussian Copula with Arbitrary
Marginals
kk
Quoting Simon Fear [EMAIL PROTECTED]:
For a more detailed answer try the American Statistician:
Title: The joy of copulas: Bivariate distributions with uniform marginals
Concerning tree graphics being trimmed, I received two
suggestions, both of which worked.
(a) plot(..., margin=0.1)
(b) set par(xpd=NA), then plot
Thanks!
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Does anyone know how I can read from a .txt file the lines that are between
two strings whose location is unknown?
My problem is that I have a .txt file with data separated by a sentence,
for example:
2.22 3.45
1.56 2.31
pattern 1
4.67 7.91
3.34 2.15
5.32 3.88
pattern 2
...
I do not know the
I'm trying to write some GMP functions but am not sure about
how and where to use pointers.
--
// alfa.c
# include R.h
# include gmp.h
void Createinteger (char *A, mpz_t *N)
{mpz_init(*N); mpz_set_str(*N,A,10);}
void
Read the lines with readLines into a character vector
Remove the lines you don't want
Use read.table on a textConnection to the character vector
If the file were very large, use an anonymous file() connection instead
(but better use OS tools such as grep to clean up the file).
On Wed, 19 Nov
With scan and match (pattern1,pattern2,..) you should be able to
build a function to perform what do you want.
A.S.
Alessandro Semeria
Models and Simulations Laboratory
Montecatini Environmental Research Center (Edison Group),
Via Ciro Menotti 48,
48023 Marina di
Fuensanta Saura Igual [EMAIL PROTECTED] Saura Igual [EMAIL PROTECTED] wrote:
Does anyone know how I can read from a .txt file the lines that
are between two strings whose location is unknown?
My problem is that I have a .txt file with data separated by a
sentence, for example:
2.22 3.45
Hi there fellow R-users,
Can anyone tell me if there exits an R package that deals with serial
correlation in the residuals of an lm model.
Perhaps, using the Cochrane Orcutt or Praise Wilson methods?
Thanks,
Wayne
Dr Wayne R. Jones
Senior Statistician / Research Analyst
KSS Limited
St
Dear expeRts,
I have lots of time series vectors and I simply want to remove certain
frequencies.
My problem is now that I can not find out how to calculate
automatic the filter coeff. for the filter() function which should
remove this certain frequencies.
Is there an elegant way to do this?
Better, arima in ts and gls in nlme can fit such models by exact maximum
likelihood.
On Wed, 19 Nov 2003, Wayne Jones wrote:
Hi there fellow R-users,
Can anyone tell me if there exits an R package that deals with serial
correlation in the residuals of an lm model.
Perhaps, using the
Dear all,
In many cases, I need a plotting region much bigger than the screen (e.g. for maps or
for graphs with many labels).
A. MS-Windows
if I try
windows(width=25, height=25, rescale=fixed)
it seems to be OK (a screen with scrollbars, exactly what I need)
but if I try then
Philippe Glaziou [EMAIL PROTECTED] writes:
Fuensanta Saura Igual [EMAIL PROTECTED] Saura Igual [EMAIL PROTECTED] wrote:
Does anyone know how I can read from a .txt file the lines that
are between two strings whose location is unknown?
My problem is that I have a .txt file with data
The Cochrane Orcutt is probably an outdated approach to deal with
autocorrelation
and it is rather easy to write code.
Why don't you use a direct likelihood-based approach?
For gaussian data see the arima() function in ts package, or the Jim
Lindsey's packages (for instance the gar() function in
Hello,
I am the president of operations at Dymaxium Inc. We are a healthcare
solutions provider for the top pharmaceutical companies all around the
world. One of our clients has asked us to recommend which statistical
packages may be best suited for some of their statistical modeling needs.
Peter Dalgaard [EMAIL PROTECTED] wrote:
My problem is that I have a .txt file with data separated by a
sentence, for example:
2.22 3.45
1.56 2.31
pattern 1
4.67 7.91
3.34 2.15
5.32 3.88
pattern 2
...
I do not know the number of lines where these separating
Philippe Glaziou [EMAIL PROTECTED] wrote:
Peter, I cannot see your point. sed can get rid of any pattern in
a text file. Fuensanta's example seemed to show that the
sentences (pattern 1, 2,...) were on separate lines from lines
containing data, thus my approach. Another one closer to your awk
Hi,
I ran a small data set from a factorial experiment through R, Minitab
and JMP... the result from R is significantly different from what
Minitab or JMP give... The data set is at the following link:
http://www.personal.psu.edu/nug107/Uploads/2x3_16repsANOVA.txt
The first 5 columns are the
My original question had the subject of making jpg-image files with the
bitmap-function.
Firstly I was told, that the problem was in my ghostscript-installation
which is called by the bitmap-function - not in R itself.
I got a lot of good suggestions (also off-list) about what problems in my
Dear all,
I have found the following (for me) incomprehensible behaviour of
ISOdate (POSIXct):
ISOdate(1900,6,16)
[1] 1900-06-15 14:00:00 Westeuropäische Sommerzeit
ISOdate(1950,6,16)
[1] 1950-06-16 14:00:00 Westeuropäische Sommerzeit
Note that in the first case I get the 15th of June back,
Dear Nirmal,
At 07:36 AM 11/19/2003 -0500, Nirmal Govind wrote:
I ran a small data set from a factorial experiment through R, Minitab and
JMP... the result from R is significantly different from what Minitab or
JMP give... The data set is at the following link:
Dear R users,
I was trying to install the package RXLisp by Duncan Temple Lang on a MDK
9.1 Linux machine running R 1.8.0 installed from a RPM.
Unfortunately I had a problem loading the shared library into R. Since
I'm a Linux newbie I was not able to solve the problem. Maybe some of
you can help
Thank you very much John and Sundar for your quick responses.
From your description, it's not possible to tell what you did or
exactly what happened -- either in R or in Minitab (and JMP).
John: in R, I fit the model using lm exactly as you did.. however, for
the ANOVA table, I used
Hello
I'm moving from using R 1.8 for Windows to using Linux (Redhat version 9)
and I cannot get any graphics. However everything else maths, models
equations is fine-- and much quicker. I built from the source file following
the commands in the install manual.
plot(1:10, 1:10)
does nothing! So
Dear Nirmal,
At 09:12 AM 11/19/2003 -0500, Nirmal Govind wrote:
Thank you very much John and Sundar for your quick responses.
From your description, it's not possible to tell what you did or
exactly what happened -- either in R or in Minitab (and JMP).
John: in R, I fit the model using lm
You probably need to set the LD_LIBRARY_PATH environment
variable to include
/usr/lib/R/library/RXLisp/libs
before you start R so that it can find the libRXLispConverters.so
file. I assume that is in that directory. If not,
use the --clean flag for R CMD INSTALL, i.e.
R CMD INSTALL
Did you build from sources or install an RPM?
If the former, you probably do not have the X11 development files
installed. Something like XFree86-devel-4.2.0-72 (from RH8.0). (My guess
is that you do have XFree86-libs.)
rm config.cache and reconfigure and rebuild after installation.
On Wed,
On Wed, 2003-11-19 at 08:47, Stephen Henderson wrote:
Hello
I'm moving from using R 1.8 for Windows to using Linux (Redhat version 9)
and I cannot get any graphics. However everything else maths, models
equations is fine-- and much quicker. I built from the source file following
the
Stephen Henderson wrote:
Hello
I'm moving from using R 1.8 for Windows to using Linux (Redhat version 9)
and I cannot get any graphics. However everything else maths, models
equations is fine-- and much quicker. I built from the source file following
the commands in the install manual.
I've
fyi, the default RH9 installation from Dell (and probably other vendors)
does not include the developmental tools. Furthermore, the default in the
RH 9.0 installation CD is not to include the developmental tools.
so, i think you are right in saying that this will become a FAQ.
On Wed, 19 Nov
Philippe Glaziou [EMAIL PROTECTED] writes:
Er, no, that wasn't the requirement. It's a job for awk or perl, e.g.
#!/usr/bin/perl -n
if (/pattern 1/){
$copy = 1;
next;
}
if (/pattern 2/){
$copy = 0;
}
print if $copy;
or
awk '/pattern
This question repeatedly comes up, e.g. see
https://www.stat.math.ethz.ch/pipermail/r-help/2003-November/040184.html
https://www.stat.math.ethz.ch/pipermail/r-help/2003-November/040348.html
from less than two weeks ago.
One thing that occurred to me is that the answer could be simplified
if
I was recently asked to perform a GLM analysis (the person comes from
the JMP world) on a repeated measures design. I have found some things
using aov but I cannot find anything with glm. In fact, multiple
regressions in general with repeated measures seems to be poorly
covered in
FWIW, there is a way to have one or more scrolled graphics devices
using the RGtk (www.omegahat.org/RGtk) and gtkDevice (from CRAN)
packages. This currently works only on Unix/Linux.
The following code creates a window with a graphics device inside a
scrolled window inside a top-level Gtk
ISOdate works, by default, in the GMT timezone. Try:
ISOdate(1900,6,16,tz=)
ISOdate(1950,6,16,tz=)
If you don't need timezones and don't want to worry about them
you can alternately use the chron library for your dates and
times.
---
Date: Wed, 19 Nov 2003 14:58:24 +0100
From: Heiko
On Wed, 19 Nov 2003, John Christie wrote:
I was recently asked to perform a GLM analysis (the person comes from
the JMP world) on a repeated measures design. I have found some things
using aov but I cannot find anything with glm. In fact, multiple
regressions in general with repeated
On Wed, 2003-11-19 at 09:29, Rafael A. Irizarry wrote:
fyi, the default RH9 installation from Dell (and probably other vendors)
does not include the developmental tools. Furthermore, the default in the
RH 9.0 installation CD is not to include the developmental tools.
so, i think you are right
Well, I live only a few seconds away from GMT and I also get
ISOdate(1900,6,16,tz=)
[1] 1900-06-15 12:00:00 GMT Daylight Time
15th, not 16th.
Was 1900 a strange leap year? I certainly haven't tested
thoroughly but note this:
ISOdate(1900,2,16,tz=)
[1] 1900-02-16 12:00:00 GMT Standard Time
#
Does this really work for you? I still get:
ISOdate(1900,6,16)
[1] 1900-06-15 14:00:00 Westeuropäische Sommerzeit
ISOdate(1900,6,16,tz=)
[1] 1900-06-15 12:00:00 Westeuropäische Sommerzeit
Obviously the time son influences the time, but it can
Not possibly account for the difference of a full
It starts on the 1st of March 1900 to go wrong I feel better now
that somebody
Else these the same effect ;-)
-Original Message-
From: Simon Fear [mailto:[EMAIL PROTECTED]
Sent: Mittwoch, 19. November 2003 17:17
To: [EMAIL PROTECTED]; [EMAIL PROTECTED];
[EMAIL PROTECTED]
Subject:
Well, one clue is that date is before the modern era, and most OSes only
go back to 1902. Some only go back to 1970! I suspect the OS does not
know that 1900 was not a leap year.
On Wed, 19 Nov 2003, Heiko Schaefer wrote:
Does this really work for you? I still get:
ISOdate(1900,6,16)
[1]
Heiko Schaefer [EMAIL PROTECTED] writes:
Does this really work for you? I still get:
ISOdate(1900,6,16)
[1] 1900-06-15 14:00:00 Westeuropäische Sommerzeit
ISOdate(1900,6,16,tz=)
[1] 1900-06-15 12:00:00 Westeuropäische Sommerzeit
Obviously the time son influences the time, but it can
You are on the right track. I suppose that linux and windows use 32 bit
For time_t and we go back beyond a valid date. Try this code:
#include stdio.h
#include time.h
int main() {
struct tm a;
time_t b;
strptime(1900-03-15 12:00:00,%Y-%m-%d %H:%M:%S,a);
printf(%s\n,asctime(a));
// now
For the record, ISOdate *is* giving the right answer, a POSIXct object.
The problem is in printing, where there was a simple coding bug: is_year
was applied to the POSIX `year' which is year-1900.
It's always worth distinguishing between the actual value and its printed
representation.
On
On Wed, 2003-11-19 at 09:26, Rick Bilonick wrote:
SNIP
I have R 1.8 installed on 3 different systems running RH 9 using the
precompiled binaries. R makes plots on all of them. One of the systems
is a laptop (Toshiba Satellite 2535cds). For some reason on this system
there is a readline
Dear r-help members
I am a newby to R and would like to estimate a simple boxcox
model, e.g.
(y^t - 1)/t ~ b0 + b1 * x1 + b2 * x2
unfortunately, R returns with an error message when I try to
perform this with the call
nls( (y^t - 1)/t ~ b0 + b1 * x1 + b2 * x2, start =
c(t=1,b0=1,b1=0,b2=0),
One other point that supports your hypothesis is that Excel 2000
on Windows 2000 has a problem in 1900 too. Enter a date into a
cell and then enter =60 and it gives the date as Feb 29, 1900 even
though 1900 is not a leap year.
Windows Excel uses 1900 as the origin for time but when
Dear all,
I am trying to compile R-1.8.0 and R-1.8.1-beta (as 19-11-2003) under
Solaris 9 using the Sun compilers. './configure' fails; the last lines of
the display are:
checking whether we can compute C Make dependencies... yes, using cc -M
checking whether cc supports -c -o FILE.lo... yes
I have a loop that increases the size of an object after each iteration. When the
Windows Task Manager shows Mem Usage about 1.8GB, the Rgui.exe process no longer
responds.
I use:
C:\Program Files\R\rw1080\bin\Rgui.exe --max-mem-size=4000M --min-vsize=10M
--max-vsize=3000M --min-nsize=500k
Based on the discussion so far, the 1950 issue is due to the
timezone and Prof. Riley has established that the 1900 issue
is a coding error in printing dates (but not in ISOdate, itself,
or in the representation of POSIXct dates).
I assume that the 1900 was found by just playing around but if
Could you compile up and try R-devel (see the FAQ)? It probably will cope
with more than 2Gb, and I've run it up to 2.5Gb.
Note that an effective limit of 1.7Gb is mentioned in the rw-FAW.
On Wed, 19 Nov 2003, Dick Beyer wrote:
I have a loop that increases the size of an object after each
Greetings all,
Every now and then, as we engage in serious discourse, it seems
appropriate to throw something into the mix that might digress and
brighten the day.
So, what is a Strouhal Number?
For an animal or insect in flight, the Strouhal number is determined by
the frequency (f) of wing
On Wed, 19 Nov 2003 10:20:16 -0800 (PST), Dick Beyer
[EMAIL PROTECTED] wrote :
I have a loop that increases the size of an object after each iteration. When the
Windows Task Manager shows Mem Usage about 1.8GB, the Rgui.exe process no longer
responds.
I use:
C:\Program
Hello.
(B
(B Is there R interface for Swarm or other MAS(Multi Agent System) tools
(Bavailable?
(B
(B Thanks.
(B
(B__
(B[EMAIL PROTECTED] mailing list
(Bhttps://www.stat.math.ethz.ch/mailman/listinfo/r-help
S-Plus has a similar feature. I once had a simulation that ran
for months. I programmed the simulator to store intermediate results
using synchronize and issue progress reports from which I could tell
how fast it was running. Every day or two, I would kill S-Plus and
restart from the
I'd like to fit a multinomial log-linear model for 4 categories of the
form
log[(P(D=i | X)/P(D=0 | X)] = alpha_i + X beta_i ; i=1,2,3
but with beta_1 constrained to zero. Is there a way to impose such a
constraint in the multinom function?
Brad
On Wed, 19 Nov 2003, Vince Forgetta wrote:
What is (or how to I change) the maximum width of a pdf image. I am
trying to make a pdf image:
pdf(file=out.pdf,width=200,height=20)
and I can't get the image to be more that 200 inches wide i.e. I get an
empty out.pdf . Is there a solution
Dear R subscribees,
First, I am new to R and I apologize if the question is naive. I am trying to
find out where arguments specifications
can be found. For example, the code below (at the very bottom of
this email) contains the argument 'digits' in the function 'lin' which
printout my output
I used dev.off(). I am certain that the problem is associated with
acrobat reader. Seeing that my collaborators use primarily acrobat
reader it looks like I'm stuck with this limit. Thanks for the help.
Vince
Marc Schwartz wrote:
On Wed, 2003-11-19 at 15:21, Prof Brian Ripley wrote:
On
On Wed, 19-Nov-2003 at 05:45PM +0100, Peter Dalgaard wrote:
| Yes, something is strange for me too (RedHat 8):
|
| ISOdate(1900,3,1)
| [1] 1900-03-01 13:00:00 CET
| ISOdate(1900,3,2)
| [1] 1900-03-01 13:00:00 CET
|
| Apparently, the one-day shift affects all dates after March 2, 1900,
| and
John Christie [EMAIL PROTECTED] writes:
OK, I found pretty much everything I needed in the short term on lme.
However, I cannot find anything on the type flag in the anova
command. When is it used and what for? In the examples I have found
type=marginal is usually entered for lme models
On Wed, 2003-11-19 at 16:34, John Christie wrote:
OK, I found pretty much everything I needed in the short term on lme.
However, I cannot find anything on the type flag in the anova
command. When is it used and what for? In the examples I have found
type=marginal is usually entered for
On Wed, 2003-11-19 at 16:54, Marc Schwartz wrote:
In your case, you are passing an nlme model object to anova, hence the
anova.nlme method is used. The details for the arguments will be listed
in the method specific help.
Quick correction, in both cases above it should be 'lme' not 'nlme'.
On Wed, 19-Nov-2003 at 05:03PM +, Prof Brian Ripley wrote:
| For the record, ISOdate *is* giving the right answer, a POSIXct object.
|
| The problem is in printing, where there was a simple coding bug: is_year
| was applied to the POSIX `year' which is year-1900.
I can't see why it doesn't
Patrick Connolly [EMAIL PROTECTED] writes:
On Wed, 19-Nov-2003 at 05:03PM +, Prof Brian Ripley wrote:
| For the record, ISOdate *is* giving the right answer, a POSIXct object.
|
| The problem is in printing, where there was a simple coding bug: is_year
| was applied to the POSIX
Peter Dalgaard [EMAIL PROTECTED] wrote:
blueberry:~/ sed -e 's/pattern 1\|pattern 2\|pattern xyz//g' tst.txt
2.22 3.45
1.56 2.31
4.67 7.91
3.34 2.15
5.32 3.88
blueberry:~/ awk '/pattern 1/{copy=1;next};/pattern 2/{copy=0};copy==1' tst.txt
4.67 7.91
3.34 2.15
5.32 3.88
blueberry:~/
Philippe Glaziou [EMAIL PROTECTED] wrote:
I see and here is my (corrected) sed solution:
cunegonde:~/tmp cat test
2.22 3.45
1.56 2.31
pattern 1
4.67 7.91
3.34 2.15
5.32 3.88
pattern 2
cunegonde:~/tmp sed -e '/pattern 1\|pattern 2/D' test
2.22 3.45
1.56 2.31
4.67 7.91
3.34 2.15
Hi R maintainers:
when I use
update.packages()
And I try to update the package maps
maps :
Version 2.0-9 in D:/rw1080/library
Version 2.0-10 on CRAN
Update (y/N)? y
The following message appears:
trying URL
`http://cran.r-project.org/bin/windows/contrib/1.8/maps_2.0-10.zip'
Error in
Kenneth -
I recall a message on the [r-pkgs] list from Ray Brownrigg
on November 1 this year that may be relevant. At that time
he was announcing the availability of maps_2.0-8.zip and
some other packages.
Therefore, 2.0-10 must be *very* recent and might still be
going through the CRAN
On Wed, 2003-11-19 at 20:27, Kenneth Cabrera wrote:
Hi R maintainers:
when I use
update.packages()
And I try to update the package maps
maps :
Version 2.0-9 in D:/rw1080/library
Version 2.0-10 on CRAN
Update (y/N)? y
The following message appears:
trying URL
Hi all,
I've ecountered a problem in the last few days with my .Rhistory file
not being able to be updated when I quit an R session because its file
attributes under Windows have been set as Hidden. Recently, I put the
following line in my global Rprofile file:-
history(max.show=Inf)
so that I
On Wed, 19 Nov 2003, Spencer Graves wrote:
Have you considered the following:
f - factor(1:2)
as.numeric(as.character(f))
[1] 1 2
See Venables and Ripley (2000) S Programming (Springer, p. 15).
Or see the FAQ, which makes the same point.
-thomas
Thomas W Blackwell [EMAIL PROTECTED] wrote:
I recall a message on the [r-pkgs] list from Ray Brownrigg
on November 1 this year that may be relevant. At that time
he was announcing the availability of maps_2.0-8.zip and
some other packages.
Therefore, 2.0-10 must be *very* recent and
Thanks to all those with helpful suggestions about my hitting the upper memory limit.
I do appreciate the help.
Thanks,
Dick
***
Richard P. Beyer, Ph.D. University of Washington
Tel.:(206) 616 7378 Env. Occ. Health
The list of arguments to a function and the default values,
if any, are given in online help:
?cov
You can also display the arguments and their defaults like this:
args(cov)
If you write your own function then there are no predefined
arguments. If you chose to use digits as an argument its
Thanks for your reply John.
this works). Applied to a linear-model object, summary() produces
coefficients, etc. (as mentioned), while anova() produces a (sequential)
ANOVA table. This seems apparent to me from the output.
What I'm having trouble with is understanding the difference between
Hi folks:
maps seems not to be available for download from
http://cran.us.r-project.org/bin/windows/contrib/1.8/. The status page
(http://cran.us.r-project.org/bin/windows/contrib/1.8/Status dated Nov 16,
2003) says though that it is OK. Any idea when it may be available?
Thanks.
Rajiv
Hi.
I could download maps package via Install Package(s) From CRAN on RGui.
maps seems not to be available for download from
http://cran.us.r-project.org/bin/windows/contrib/1.8/. The status page
(http://cran.us.r-project.org/bin/windows/contrib/1.8/Status dated Nov 16,
2003) says though
Umm... no, does not work for me yet:
local({a - CRAN.packages()
+ install.packages(select.list(a[,1],,TRUE), .libPaths()[1], available=a)})
trying URL `http://cran.r-project.org/bin/windows/contrib/1.8/PACKAGES'
Content type `text/plain; charset=iso-8859-1' length 13514 bytes
opened URL
80 matches
Mail list logo
