This was much more appropriate for R-devel, so please move any followup
there.
On Wed, 11 Apr 2007, William Constantine wrote:
I am interested in tapping into LAPACK functionality directly from R.
Using R-2.4.1 for Windows, I was able to do so ala:
dyn.load(bin/Rlapack.dll)
Hallo
[EMAIL PROTECTED] napsal dne 11.04.2007 22:45:19:
Hi
I am very new to R and have not able to find the answer to my question
in the manual or any other post. I have a dataset that has various
different sites names with data relating to each site. The data is in
one dataset. I want
Hello
I tried to load the taskPR-package in R, but it doesn't work at all.
Everytime I try /library(taskPR)/, I get the error-message:
Error in dyn.load(x, as.logical(local), as.logical(now)) :
unable to load shared library
'/home/bw135690/R-2.4.1/library/taskPR/libs/taskPR.so':
I'm doing the time-series data analysis using FFT. After fft(), I used
spectrum() to calculate the power spectrum, which is used to define the main
frequency. Now I want to get the phase information according to this frequency.
For example, input 12 time-point data,
I have a list of groups of xy positions I want to set to 0 in an
array full of 1s. When the assignments are done directly from the
list, they are incorrect, while if I use a temporary array derived
from the list the assignments are correct. The following example will
hopefully make my
It would have been helpful to have given the details the R posting guide
asks for, such as your OS.
That message probably does not mean that
/home/bw135690/R-2.4.1/library/taskPR/libs/taskPR.so does not exist.
It says 'cannot open shared object file', and this can have many causes.
One is that
Greetings list,
I have a rectangular 20 x 20 matrix with values in the range of [0 ,
1]. I'd like to plot it as an image. To that end, I have used the
image() function that seems to do what I want. Now, I just want to
tweak it to look perfect. So here is my question:
At the moment, the
try the following:
z - matrix(1, ncol = 249, nrow = 240)
zz - matrix(1, ncol = 249, nrow = 240)
for (k in seq_along(x)) {
z[rbind(x[[k]])] - 0
}
zz[xx] - 0
all.equal(z, zz)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public
Dear R-Users,
I am currently working with LME to analyse repeated measures data. I encounter
a problem when including both a random effect and a correlation structure with
different grouping levels into the LME model. The error message is:
Error in lme.formula(diameter ~ flowers*timef +
Hi
Did you by chance look at the help page of image? If you did, you could
read
x,y
locations of grid lines at which the values in z are measured. These must
be finite, non-missing and in (strictly) ascending order. By default,
equally spaced values from 0 to 1 are used. If x is a list, its
Charilaos Skiadas wrote:
A new fortune candidate perhaps?
On Apr 10, 2007, at 6:27 PM, Greg Snow wrote:
Remember, everything is better than everything else given the right
comparison.
Only if we remove the grammatical blip that turns it into an infinite
regress, i.e.
Remember, anything
Wee-Jin Goh wrote:
Greetings list,
I have a rectangular 20 x 20 matrix with values in the range of [0 ,
1]. I'd like to plot it as an image. To that end, I have used the
image() function that seems to do what I want. Now, I just want to
tweak it to look perfect. So here is my
I am hoping to obtain information regarding the reporting of
multinomial logistic regression results - either recommendations of
published studies employing R or suggestions about what output should
be reported. I am used to SPSS which routinely produces a broader
range of output than R (e.g
Hello everybody,
I'm trying to install rJava to use JRI in a Java program as JRI is a part of
rJava. There should be the option to install this package via
install.packages(rJava). Unfortunately I'm connected to the internet via a
proxy. I have no idea how to tell R to connect via the proxy.
Ron,
you're right. It's not legitimate at all. I suggest you
to take a look at the HUGE bibliography on cointegration,
as a start up.
Rogerio
Dear all R user,
Please forgive me if my question is too simple.
My question is related to Statistics rather directly to R. Suppose I have
Wow! Thanks to Sundar and to Deepayan for this selfless effort.(I
rather suspect they had a good time figuring this out.) The outcome
is that I will be able to publish just the right figure for the job.
I will be sure to aknowledge your contribtution.
Michael
On Apr 12, 2007, at 12:48 AM,
What operating system do you use?
If you use Windows, than open the Control Panel, double click on the System
icon, go to Advanced tab, press
Environment variables button, then press New to create one more new
variable, enter http_proxy as the Name of variable, and
Hello,
I have a problem with the import of a date file. I seems verry tricky.
I have a text file (end of the mail). Every file has a different number of
measurments
witch start with START OF HEIGHT DATA and ende with END OF HEIGHT DATA.
I imported the file in a matrix but the letters before the
I had that problme and ended up downloading the package and installing
locally.
there is probably a better way
Richard
Erbe, Nicole wrote:
Hello everybody,
I'm trying to install rJava to use JRI in a Java program as JRI is a part of
rJava. There should be the option to install this package
Hi,
I'm (currently forced to) use Win XP. I actually found an easy way to solve the
problem of installing the rJava package. To avoid the problems with the proxy I
just downloaded the zip file and installed the package in R from a local
source. Now I can use the JRI libraries in eclipse by
On Thu, 12 Apr 2007, Vladimir Eremeev wrote:
What operating system do you use?
If you use Windows, than open the Control Panel, double click on the System
icon, go to Advanced tab, press
Environment variables button, then press New to create one more new
variable, enter http_proxy as the
Douglas Bates writes:
One
can do data analysis by using the computer as a blunt instrument with
which to bludgeon the problem to death but one can't do elegant data
analysis like that.
One nice thing about a blunt instrument like Stata is the ability to
hold an entire dataset in memory
A re-interpretation of Zorn's lemma?
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon
Sent: Thursday, April 12, 2007 5:14 AM
To: [EMAIL PROTECTED]
Subject: Re: [R] Reasons to Use R
Charilaos Skiadas wrote:
A new fortune candidate perhaps?
Dear R-experts,
I have the following function:
userInput - function() {
ANSWER - readline(saving place of the data (example
Z:/Software/test.mat)? )
x = c(.mat)
endingTest = x %in% ANSWER
print (endingTest)
if (endingTest == ??)
cat (saving place not
At 14:39 11/04/2007, Robert Biddle wrote:
Hi:
I am doing some work that involves plotting points of interest
superimposed on photographs and maps. I can produce the plots fine
in R, but so far
I have had to do the superimposition externally, which makes it
tedious to do exploratory work.
I have
Dear all,
in the package its the print() method does not seem to correctly work in all
circumstances:
selectMethod(print, its)
Method Definition:
function (x, ...)
{
print([EMAIL PROTECTED] mailto:[EMAIL PROTECTED] , ...)
}
environment: namespace:its
Signatures:
x
Try pasting this into an R session:
Lines.raw - FILEDATE:02.02.2007
...
START OF HEIGHT DATA
S= 0 y=0.0 x=0.
S= 0 y=0.1 x=0.00055643
...
S= 9 y=4.9 x=1.67278117
S= 9 y=5.0 x=1.74873257
S=10 y=0.0 x=0.
S=10 y=0.1 x=0.00075557
...
S=99 y=5.3 x=1.94719490
END OF HEIGHT DATA
...
Without a reproducible example I can't tell what the problem is
but as a workaround you could print it as a zoo object:
library(its)
library(zoo)
# create an its object
example(its)
ii - its(mat, times)
# print it without using its' print
print(as.zoo(ii))
# or from the command line just
Albrecht, Dr. Stefan (AZ Private Equity Partner)
[EMAIL PROTECTED] writes:
Dear all,
in the package its the print() method does not seem to correctly
work in all circumstances:
My understanding is that one should define only a method for the show
generic for S4 classes and leave print
Lucke, Joseph F writes:
A re-interpretation of Zorn's lemma?
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon
Sent: Thursday, April 12, 2007 5:14 AM
To: [EMAIL PROTECTED]
Subject: Re: [R] Reasons to Use R
Charilaos Skiadas
Okay thanks, I'm going through the docs now.. and I came through this..
The named field is set and accessed by the SET_NAMED and NAMED macros, and
take values 0, 1 and 2. R has a `call by value' illusion, so an assignment
like
b - a
appears to make a copy of a and refer to it as b.
Setting S4 methods for print() works correctly, as in your examples.
Your first example does not call print(), so it is not suprising that no
method dispatch takes place.
You seem to suppose that such a method should be used for top-level
auto-printing, and it will not be as your example
Dear friends.
With file( ) to obtain a pointer of a file, every time we use scan ( ) to
read one row of it, the pointer will point to the next row of the file. In
the following example, d1 and d2 are obtained the same way but they
correspond to different rows of the same file because the pointer
Here is the contents of my testdata.txt :
-
START OF HEIGHT DATA
S= 0y=0.0 x=0.
S= 0 y=0.1 x=0.00055643
S= 9 y=4.9 x=1.67278117
S= 9 y=5.0 x=1.74873257
S=10 y=0.0 x=0.
S=10y=0.1 x=0.00075557
S=99
Thanks all of you!
But my next question is, how to improve R programming skills? I never
have time in improving it but I feel I need to.
Regards,
W
On 4/11/07, Seth Falcon [EMAIL PROTECTED] wrote:
Weiwei Shi [EMAIL PROTECTED] writes:
I forgot to add my bad solution here:
reverseList -
Schmitt, Corinna wrote:
Dear R-experts,
I have the following function:
userInput - function() {
ANSWER - readline(saving place of the data (example
Z:/Software/test.mat)? )
x = c(.mat)
endingTest = x %in% ANSWER
print (endingTest)
if (endingTest
R users:
I'd like to get some insight on an error I encounter when attempting to
work with two moderately sized sets of time series data. FYI - I'm
using the following versions of R and supporting packages on a Windows
2000 OS:
- R version 2.4.1 (2006-12-18)
- zoo version 1.2-2
- chron version
On Thu, 12 Apr 2007, Yuchen Luo wrote:
With file( ) to obtain a pointer of a file, every time we use scan ( ) to
read one row of it, the pointer will point to the next row of the file. In
the following example, d1 and d2 are obtained the same way but they
correspond to different rows of the
?seek
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Yuchen Luo
Sent: Thursday, April 12, 2007 9:33 AM
To: [EMAIL PROTECTED]
Dear Rologists,
I'm stuck with this. How would you do this efficiently:
aPGI
[[1]]
[1] 864 5576
aPGItest
[[1]]
[1] TRUE FALSE
result - [magic box involving subset)
result
[[1]]
[1] 864
Thanks for any hints,
Joh
__
[EMAIL PROTECTED] mailing
On Thu, 12 Apr 2007 [EMAIL PROTECTED] wrote:
attempted to get at this by merging the two time series (by union),
interpolating the NAs, and finally, subtracting one vector from the
other. The problem is that I can not combine the two zoo time series
objects using the merge() or cbind()
On Thu, 2007-04-12 at 18:12 +0200, Johannes Graumann wrote:
Dear Rologists,
I'm stuck with this. How would you do this efficiently:
aPGI
[[1]]
[1] 864 5576
aPGItest
[[1]]
[1] TRUE FALSE
result - [magic box involving subset)
result
[[1]]
[1] 864
Thanks for any
Johannes Graumann [EMAIL PROTECTED] writes:
Dear Rologists,
I'm stuck with this. How would you do this efficiently:
aPGI
[[1]]
[1] 864 5576
aPGItest
[[1]]
[1] TRUE FALSE
result - [magic box involving subset)
result - aPGI[aPGItest[[1]]]
result
[[1]]
[1] 864
+ seth
--
Seth
From: Marc Schwartz
On Thu, 2007-04-12 at 18:12 +0200, Johannes Graumann wrote:
Dear Rologists,
I'm stuck with this. How would you do this efficiently:
aPGI
[[1]]
[1] 864 5576
aPGItest
[[1]]
[1] TRUE FALSE
result - [magic box involving subset)
result
--- Tim Churches [EMAIL PROTECTED] wrote:
John Kane wrote:
I have 3 SAS files all in the directory F:/sas,
two
data files
and a format file :
form.ea1.sas7bdat
form.ea2.sas7bdat
sas.fmts.sas7bdat
F is a USB.
I am trying import them to R using sas.get.
I have not
Hi,
I have time series data stored in csv files (see below for an example of the
data). I understand that in order to analyze my data in R, I need to first
transform it into a ts object. Howeve, I could not find an example on how
exactly to do that. Is ts the only function I need? What are the
R plotting experts:
I have a bivariate dataset composed of 300 (x,y) continuous datapoints.
297 of these points are located within the y range of [0,10], while 2
are located at 20 and one at 55. No coding errors, real outliers.
When plotting these data with a scatterplot, I obviously have a
Cancel this! . I see that Hmisc is already changed.
When I updated today it was changed.
My thanks to Mr Dupont the maintainer.
I still have a problem with sas.get but that's another
post.
--- John Kane [EMAIL PROTECTED] wrote:
--- Tim Churches [EMAIL PROTECTED] wrote:
John Kane
On Thu, 12 Apr 2007, tom soyer wrote:
Hi,
I have time series data stored in csv files (see below for an example of the
data). I understand that in order to analyze my data in R, I need to first
transform it into a ts object. Howeve, I could not find an example on how
exactly to do that. Is
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote:
Hi,
I have time series data stored in csv files (see below for an example of the
data). I understand that in order to analyze my data in R, I need to first
transform it into a ts object. Howeve, I could not find an example on how
exactly to do
Thanks Gabor!
I think your example works, but check this out:
as.ts(z)
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1921 19.0 18.4 18.3 18.1 17.7 17.6 17.7 17.7 17.5 17.5 17.4 17.3
1922 16.9 16.9 16.7 16.7 16.7 16.7 16.8 16.6 16.6 16.7 16.8 16.9
is.ts(z)
[1] FALSE
How
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote:
Thanks Gabor!
I think your example works, but check this out:
as.ts(z)
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1921 19.0 18.4 18.3 18.1 17.7 17.6 17.7 17.7 17.5 17.5 17.4 17.3
1922 16.9 16.9 16.7 16.7 16.7 16.7 16.8
Oh yes, I forgot. Thanks!!
On 4/12/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote:
Thanks Gabor!
I think your example works, but check this out:
as.ts(z)
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1921 19.0 18.4
Gabor,
What should I do for weekly and daily data series? Are there functions
similar to yearmon() for other time intervals? I see that there is a
built-in yearqtr() function for quarterly data, but that's it.
Thanks!
On 4/12/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 4/12/07, tom
On Thu, 2007-04-12 at 13:41 -0500, Inman, Brant A. M.D. wrote:
R plotting experts:
I have a bivariate dataset composed of 300 (x,y) continuous datapoints.
297 of these points are located within the y range of [0,10], while 2
are located at 20 and one at 55. No coding errors, real outliers.
Hi all..RCounters!
I´m working with standarized time series, and i need make a counter of
consecutives positves numbers to make a cumulative experimental funtion. I
have x: the time series (0,1) and y: my counter, i have this for step. What is
wrong?.. any can help me please!
On 4/12/07, Marc Schwartz [EMAIL PROTECTED] wrote:
On Thu, 2007-04-12 at 13:41 -0500, Inman, Brant A. M.D. wrote:
R plotting experts:
I have a bivariate dataset composed of 300 (x,y) continuous datapoints.
297 of these points are located within the y range of [0,10], while 2
are located
Hi:
I've been reading Computational Methods for Multilevel Modeling by Pinheiro
and Bates,
the idea of embedding the technique in my own c-level code. The basic idea is
to rewrite
the joint density in a form to mimic a single least squares problem conditional
upon the
variance parameters.
Try something like this (modify to how you like it):
x - runif(100)
y - rnorm(100, 5, 2)
y[1:3] - c(19, 21, 50)
layout(matrix( 3:1, ncol=1 ), heights=c(2,3,4))
par(mar=c(5,4,0,2)+0.1)
plot(x,y, ylim=c(0,10), ylab='')
par(mar=c(0.5,4,0,2)+0.1)
plot(x,y, ylim=c(18,22), xlab='', xaxt='n' )
Inman, Brant A. M.D. wrote:
R plotting experts:
I have a bivariate dataset composed of 300 (x,y) continuous datapoints.
297 of these points are located within the y range of [0,10], while 2
are located at 20 and one at 55. No coding errors, real outliers.
When plotting these data with a
Hi:
I've been reading Computational Methods for Multilevel Modeling by Pinheiro
and Bates, the idea of embedding the technique in my own c-level code. The
basic idea is to rewrite the joint density in a form to mimic a single least
squares problem conditional upon the variance parameters. The
Dear Friends,
I'm trying to generate a sequence of 100 observations
with either a 1 or -1. In other words the sequence
should look something like this.
y = 1 1 -1 1 -1 -1 -1 1 1 ..
Can somebody please give me some direction on how I
can do this in R.
Thanks
Anup
Hi Anup,
(runif(100).5)*1 #would give you 0's and 1's.
sample(rep(c(-1,1),50),100) #A bit slower I think, gives you -1's and 1's
Best,
Matt
On 4/12/07, Anup Nandialath [EMAIL PROTECTED] wrote:
Dear Friends,
I'm trying to generate a sequence of 100 observations
with either a 1 or -1.
Anup Nandialath wrote:
Dear Friends,
I'm trying to generate a sequence of 100 observations
with either a 1 or -1. In other words the sequence
should look something like this.
y = 1 1 -1 1 -1 -1 -1 1 1 ..
Can somebody please give me some direction on how I
can do this in R.
There are many ways to do this.
The first that comes to my mind is sample(c(1,-1),100,TRUE). Notice
that sample also has a prob argument that may be useful for you.
Francisco
Anup Nandialath wrote:
Dear Friends,
I'm trying to generate a sequence of 100 observations
with either a 1 or
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Stephen Weigand
Sent: Thursday, April 12, 2007 2:32 PM
To: Inman, Brant A. M.D.
Cc: [EMAIL PROTECTED]
Subject: Re: [R] Putting 2 breaks on Y axis
[snip]
I was thinking plot the data without the
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Matthew Keller
Sent: Thursday, April 12, 2007 2:32 PM
To: Anup Nandialath
Cc: [EMAIL PROTECTED]
Subject: Re: [R] Random Sequence
Hi Anup,
(runif(100).5)*1 #would give you 0's and 1's.
Hi R-Users,
I have 3 replicates ('Replicate) of counts of parasites ('nor.tot.lep')
before and after an experiment ('In.Out'). I am trying to treat the
three replicates as a random effect in order to determine if the main
effect (In.Out) significantly influences my dependent variable
Hi, there:
After I upgraded my R to 2.4.1, it is my first time of trying to use
MASS and found the following error message:
install.packages(MASS)
--- Please select a CRAN mirror for use in this session ---
trying URL
'http://cran.cnr.Berkeley.edu/bin/macosx/universal/contrib/2.4/VR_7.2-33.tgz'
Hi,
I am working on a project about hospital efficiency. Due to the high
multicolinearlity of the data, I want to fit the model using ridge
regression. However, I believe that the data from large hospital(indicated
by the number of patients they treat a year) is more accurate than from
small
On 12-Apr-07 10:14:21, Jim Lemon wrote:
Charilaos Skiadas wrote:
A new fortune candidate perhaps?
On Apr 10, 2007, at 6:27 PM, Greg Snow wrote:
Remember, everything is better than everything else given the
right comparison.
Only if we remove the grammatical blip that turns it into an
Hello Brendan,
do you have any particular objection to
lmer(nor.tot.lep ~ In.Out + (1|Replicate), data=coho, family=poisson)
?
Andrew
On Thu, Apr 12, 2007 at 02:58:07PM -0700, Brendan Connors wrote:
Hi R-Users,
I have 3 replicates ('Replicate) of counts of parasites ('nor.tot.lep')
On 4/12/07, Izmirlian, Grant (NIH/NCI) [E] [EMAIL PROTECTED] wrote:
Hi:
I've been reading Computational Methods for Multilevel Modeling by Pinheiro
and Bates, the idea of embedding the technique in my own c-level code. The
basic idea is to rewrite the joint density in a form to mimic a
Hi, R-experts
I defined a new class of object, called alpha, which is basically a
data.frame. And (I think) I know how to create a method for alpha, such as
summary.alpha, plot.alpha. The problem is, when I try to access alpha
object by usual data.frame method, it won't.
For example, suppose X
Weiwei,
I have never had much success installing packages from within R.app
on MacOSX, because the location that it is supposed to save things, /
Library/Frameworks/, needs elevated priviledges, which the app
doesn't seem to try to get. So it at best ends up saving it in some
Hi Charilaos:
I think the installation has no problem as indicated by
install.packages(). Anyway, I tried your way but unfortunately it
does not help.
Thanks though,
-W
On 4/12/07, Charilaos Skiadas [EMAIL PROTECTED] wrote:
Weiwei,
I have never had much success installing packages
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote:
What should I do for weekly and daily data series? Are there functions
similar to yearmon() for other time intervals? I see that there is a
built-in yearqtr() function for quarterly data, but that's it.
ts series cannot directly represent daily
Hi gang.
So I'm actually carrying through on the threat I made a couple of months
back, and helping to start a local users group for the NYC Metropolitan
area. Here's a copy of the announcement of our introductory / planning
meeting, I also can send a flyer version of this to anyone who
I'm having trouble getting a black and white plot using the trellis device.
I have color graphics, but I'd like to use bw for a paper. For example,
trellis.device(color=F)
d - read.table(textConnection(
A,B
0,1
1,2), header=T, sep=,)
xyplot(A ~ B, data=d, type=o, lty=c(1:2), col=2, pch=19:20)
On 4/12/07, Brian Riordan [EMAIL PROTECTED] wrote:
I'm having trouble getting a black and white plot using the trellis device.
I have color graphics, but I'd like to use bw for a paper. For example,
trellis.device(color=F)
d - read.table(textConnection(
A,B
0,1
1,2), header=T, sep=,)
On Thu, 12 Apr 2007, jiho.han wrote:
Hi, R-experts
I defined a new class of object, called alpha, which is basically a
data.frame. And (I think) I know how to create a method for alpha, such as
summary.alpha, plot.alpha. The problem is, when I try to access alpha
object by usual data.frame
On Thu, 2007-04-12 at 22:30 +0200, Jose Bustos Melo wrote:
Hi all..RCounters!
Im working with standarized time series, and i need make a counter
of consecutives positves numbers to make a cumulative experimental
funtion. I have x: the time series (0,1) and y: my counter, i have
this
MASS is part of the standard R installation, as a recommended package.
I suggest you try re-installing R and not then trying to download a
different version.
I would also suggest trying 2.5.0 beta at this stage in preference to
2.4.1 if you do need to re-install R.
Beyond that, there is a
83 matches
Mail list logo