I generally work in SAS but have some graphics features I would like to
run in r. I would like to do this 'automatically' from SAS.
I'm thinking of something along the lines of putting the r code in a text
file and calling system of x on it; r would expect the data in a certain place,
and my
Hello,
I feel a brief follow up to my post a few days ago is in order. Thanks
to Renaud Lancelot's truly excelent step by step guide I solved my own
difficulties in building R packages for personal use. A short while
after he sent me his document the reply was -Voilà! Je l'ai
fonctionnant
On Wed, 13 Apr 2005 09:36:21 +0100 (BST),
(Ted Harding) ((H) wrote:
On 13-Apr-05 Prof Brian Ripley wrote:
On Wed, 13 Apr 2005, BORGULYA [iso-8859-2] Gábor wrote:
Has onyone experience with pstoedit
(http://www.pstoedit.net/pstoedit)
to convert eps graphs generated by R on
Hi
I tried to cross compile R under Linux but get an
error.
i586-mingw32-gcc -isystem
/home/schouwl/unpack/mingw/include -O2 -Wall -pedantic
-I../include -I. -DHAVE_CONFIG_H -DR_DLL_BUILD -c
dynload.c -o dynload.o
dynload.c: In function `R_loadLibrary':
dynload.c:94: warning: implicit
Hi all.
I'm compiling the devel version of R on Debian GNU/Linux, and installing it
into /usr/local tree (instead of default /usr). So:
./configure --prefix=/usr/local/
make
make install
Everything works fine, but when I start R I get the following error messages
(traslated from italian,
Hi list,
I'd like to fit generalized gamma and weibull distributions to a number
of data sets. I've been searching around and found references to R and
Jim Lindsey's GNLM package, which has the gnlr and gnlr3 procedures that
can do this.
Now, I'm completely new to R, and I'm working my way
Dear R-users,
One of the main reasons I moved from GAUSS to R (as an econometrician) was
because of the existence of the library LOCFIT for local polynomial regression.
While doing some checking between my former `GAUSS code' and my new `R code', I
came to realize LOCFIT is not quite doing
I tried the versions
latest beta of R 2.1 as well
R-latest.tar.gz13-Apr-2005 17:27 11.6M
and tried out the two different mingw packages for
cross compilation
http://www.stats.ox.ac.uk/pub/Rtools/
I still get the same error.
Regards
Lars
--- Lars Schouw [EMAIL PROTECTED] wrote:
You should definitely read Loader's book. Anyway, in the meantime, you should
look an introductory paper that you will find at the Locfit web page. I think
that you can set Locfit to estimate at all the sample points, which it does
not by default, and also to use a prespecified constant
Dear Tom,
Many thanks. I think this gets me in the right direction, but
concatenates all levels into one long level. Any further thoughts?
Best regards,
Jan
Mulholland, Tom wrote:
This may not be the best way but in the past I think I have done something like
levels(x) -
Michael S wrote:
Dear ALL-R helpers,
I want to let R get vector from c ,for example :numeric array ,vector .I
saw some exmple like this :
/* useCall3.c*/
/* Getting an integer vector from C using .Call */
#include R.h
#include Rdefines.h
SEXP setInt() {
*If* you really have the header paths set correctly it does work, and has
been tested by several people.
Do read MkRules more carefully and think about how your setting differs
from the example given. You have *not* as you claim
# Set this to where the mingw32 include files are. It must be
On Thu, 14 Apr 2005, Stefano Calza wrote:
Hi all.
I'm compiling the devel version of R on Debian GNU/Linux, and installing
it into /usr/local tree (instead of default /usr). So:
The default *is* --prefix=/usr/local (no trailing space). Are you sure you
are not getting R confused with some other
From: Y Y
I generally work in SAS but have some graphics features I
would like to
run in r. I would like to do this 'automatically' from SAS.
I'm thinking of something along the lines of putting the r
code in a text
file and calling system of x on it; r would expect the data
in
Thanks to Prof. Ripley I solved it.
Actually it was my (stupid) fault. In .Renviron I actually set R_LIBS
Thanks again
Stefano
On Thu, Apr 14, 2005 at 11:07:00AM +0100, Prof Brian Ripley wrote:
ProfOn Thu, 14 Apr 2005, Stefano Calza wrote:
Prof
ProfHi all.
Prof
ProfI'm compiling the devel
Dear R-Help
I have some trouble to set the legend in a xyplot into two rows.
The code below gives me the legend in the layout I am looking for, I
just rather have it in two rows.
library(lattice)
schluessel - list(
points=list( col=red, pch=19, cex=0.5 ),
I am looking to perform a simple linear interpolation to fill a few small
gaps in a large data set.
The data set tends to be either continous with one or two gaps in it
(where I hope to perform the interpolation), or else it has large chunks
of data missing where I'd need to return an NA.
Also,
Ben Dave all,
I'm a user of ADModel (product of Otter Research)
Just a word to say that for maximisation, I always rely on Admodel.
It's really fast (amazing when you have an important number of parameters),
can be used either as a standalone application or as DLL
I do use GAUSS (Aptech), R
Hi,
I was trying like to print out data frame with POSIXct column
in html format using xtable package, but I got error message
when trying to print the table. Here is example:
aaa-data.frame(as.POSIXct(strptime('03 2005', '%d %Y')),0)
aaa.tab-xtable(aaa)
print(aaa.tab)
Error in
Hi,
I am using the following command to print to a file (I omitted the file
details):
cat( paste( paste(orientation, start, end, names,\n), paste(start, end,
exon\n), sep=))
where orientation and names are character vectors and start and
end are integer vectors.
The problem is that R coerce the
Dear all,
I've got a problem with the function atan2. For a couple of coordinates
x and y,
This function returns the angle between the vector of coordinates (x, y)
and the
abscissa axis, i.e. it is the same as atan(y/x) (as indicated on the
help page).
If we consider the vector with coordinates
On 04/13/05 21:05, Chris Bergstresser wrote:
This article is great; thanks for providing it. The authors
recommend either using ML Estimation or Multiple Imputation to fill
in the missing data. They don't talk much about which is better for
certain situations, however.
Multiple
On Thu, Apr 14, 2005 at 02:32:33PM +0300, Firas Swidan wrote:
I am using the following command to print to a file (I omitted the file
details):
cat( paste( paste(orientation, start, end, names,\n), paste(start, end,
exon\n), sep=))
where orientation and names are character vectors and
I am currently trying to get a lme analyses running to correct for the
non-independence of residuals (using e.g. corAR1, corARMA) for a larger data
set (1 obs) for an independent (lgeodisE) and dependent variable
(gendis). Previous attempts using SAS failed. In addition we were told by
SAS
Clément Calenge [EMAIL PROTECTED] writes:
Dear all,
I've got a problem with the function atan2. For a couple of
coordinates x and y,
This function returns the angle between the vector of coordinates (x,
y) and the
abscissa axis, i.e. it is the same as atan(y/x) (as indicated on the
help
On Thu, 14 Apr 2005, Miha Razinger wrote:
Hi,
I was trying like to print out data frame with POSIXct column
in html format using xtable package, but I got error message
when trying to print the table. Here is example:
aaa-data.frame(as.POSIXct(strptime('03 2005', '%d %Y')),0)
aaa.tab-xtable(aaa)
Well, you have to convert an integer to character to see it: `as is' is in
your case 64 0's and 1's.
I very much suspect that you have a double and not an integer:
10
[1] 1e+05
as.integer(10)
[1] 10
so that is one answer: actually use an `integer vector' as you claim.
A second answer
you should include the 'form' argument in corARMA(), i.e.,
corARMA(form=~1|dummy, p=1, q=1)
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:
I have noticed that after I ran a batch script multiple times I get
multiple copies of a package's name when I call search(). Is this a
problem?
search()
[1] .GlobalEnvDF DF
[4] DF DF DF
multiple copies here ...
[13] DF DF DF
Hi,
I need to compile R-2.0.1 on an opteron running suse9.1 and using
portland group compilers. Can you advise me of the environemt variables
I need to set, particulalry the FPICFLAGS.
thanks,
Julie
---
Dr Julie Harold: University of
Dear All,
I have small samples of data (between 6 and 15) for
numerious time series points. I am assuming the data
for each time point is normally distributed. The
problem is that the data arrvies sporadically and I
would like to detect the number of outliers after I
have six data points for any
I suspect you've attach()'ed `DF' multiple times in your
code (possibly inside a loop, or perhaps a function that
was called several times). Note that if it were a
`package', it would show up in search() as `package:DF'
rather than just `DF'. Also, R Core folks took care to
avoid attaching the
Building on Tom's reply, the following should work:
labels - factor(paste(This is a long label , 1:10))
labels
[1] This is a long label 1 This is a long label 2
[3] This is a long label 3 This is a long label 4
[5] This is a long label 5 This is a long label 6
[7] This is a long
On Thu, 14 Apr 2005, Liaw, Andy wrote:
I suspect you've attach()'ed `DF' multiple times in your
code (possibly inside a loop, or perhaps a function that
was called several times). Note that if it were a
`package', it would show up in search() as `package:DF'
rather than just `DF'. Also, R Core
Dear Dave,
I do not know the grubbs.test (is it a function, where can I find it?)
and probably n=6 data points are really few..
Having said that, what do you mean as outlier?
If you mean deviation from the estimated mean (of previous data), you
might have a look to the strucchange
Does anyone have nice
quality controlls for affy arrays,
Can't find any tools as are being used for
2 dye arrays.
cheers,
marinus
This e-mail and its contents are subject to the DISCLAIMER at
http://www.tno.nl/disclaimer/email.html
[[alternative HTML version deleted]]
Berton Gunter wrote:
My apologies if this is obvious:
Is there a simple way (other than simulation or bootstrapping) to obtain a
(approximate)confidence interval for the ratio of 2 variance components in a
fitted lme model? -- In particular, if there are only 2 components (1
grouping factor). I'm
Marinus
While there isn't an arrayquality package like for cDNA arrays, there
are many available tools for affy. You probably want to look at the
bioconductor packages for affy. In any case, there is a bioconductor
list that is the better list for such questions.
Sean
On Apr 14, 2005, at
On Thursday 14 April 2005 05:30, Gesmann, Markus wrote:
Dear R-Help
I have some trouble to set the legend in a xyplot into two rows.
The code below gives me the legend in the layout I am looking for, I
just rather have it in two rows.
library(lattice)
schluessel - list(
Prof Ripley has quite rightly pointed out that my previous comment could
be (quite wrongly) interpreted as denigrating the work of the R team.
This was far from my intention. I apologise to all concerned for a
terribly flippant remark that was certainly not meant as a criticism of
anyone. It
Ronaldo Reis-Jr. wrote:
Hi,
I try to make a binomial analysis using GLMM in a longitudinal data file.
Is correct to use anova(model) to access the significance of the fixed terms?
Thanks
Ronaldo
From lme4_0.95-1 on the GLMM function has been replaced by lmer with a
non-missing family argument.
Thanks Deepayan!
Your solution does excatly what I want.
Further experiments and thoughts on my side brought me also to a
solution.
If I use the option rep=FALSE, and plot the bullit with lines and
split the lines argument into two groups it gives me the same result,
as every item in the key
Hi,
thanks for the suggestions. However, for some reason the first one did not
work. Trying
cat( paste( paste(orientation, as.integer(start), as.integer(end),
names,\n), paste(as.integer(start), as.integer(end),exon\n), sep=))
resulted in the same problem.
Setting scipen in options did the
Greetings.
I've got some analysis problems I'm trying to solve, the raw data for which
are accumulated in a bunch of time-and-date-based files.
/some/path/2005-01-02-00-00-02
etc.
The best 'read all these files' method I've seen in the r-help archives comes
down to
for (df in
On Thursday 14 April 2005 10:29, Gesmann, Markus wrote:
Thanks Deepayan!
Your solution does excatly what I want.
Further experiments and thoughts on my side brought me also to a
solution.
If I use the option rep=FALSE, and plot the bullit with lines and
split the lines argument into two
On Thu, 14 Apr 2005, Firas Swidan wrote:
Hi,
thanks for the suggestions. However, for some reason the first one did not
work. Trying
cat( paste( paste(orientation, as.integer(start), as.integer(end),
names,\n), paste(as.integer(start), as.integer(end),exon\n), sep=))
resulted in the same problem.
[EMAIL PROTECTED] writes:
Greetings.
I've got some analysis problems I'm trying to solve, the raw data for which
are accumulated in a bunch of time-and-date-based files.
/some/path/2005-01-02-00-00-02
etc.
The best 'read all these files' method I've seen in the r-help archives
i have try hard to find the answer by google,but i can not find any solution.
so i wan to ask:
1,can we test the if canonical relationship is significant after using cancor?
2,if it can,how?
3,if not,is it under-developed or there is not need to do it?or there is no
good way to do it?
i hope my
On Sun, 10 Apr 2005 17:57:26 +0100 (BST)
Prof Brian Ripley [EMAIL PROTECTED] wrote:
On Sun, 10 Apr 2005, ronggui wrote:
when i use R(2.1.0) under windows,it can display Chinese well.and the
R 2.1.0 will not be out for 8 days: are you a time-traveller or careless?
(The posting guide
i have find the answer to this question.
sorry for my overlook the information about it in internet.
2,as some textbook say,when using fisher's method,we proceed by assuming
that the within-group covariance structure for our data is the same across
groups,so we need test for equailty of
__
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-project.org/posting-guide.html
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce our May-2005 2-day
S-PLUS / R: Complementing and Extending Statistical Computing for SAS
Users
www.xlsolutions-corp.com/Rsas.htm
Boston -- May 26-27
Raleigh
The Grubbs test is one of many old (1950's - '70's) and classical tests for
outliers in linear regression. Here's a link:
http://www.itl.nist.gov/div898/handbook/eda/section3/eda35h.htm
I think it fair to say that such outlier detection methods were long ago
found to be deficient and have poor
Dear John,
Thanks for the answer! In my own dataset, The
multinom() did not converge even after I had tried to
increase the maximum number of iteration (from default
100 to 1000). In this situation, there is some bigger
diffrenece in fitted probabilities under different
contrasts (e.g. 0.9687817
Dnia 2005-04-14 15:34, Uytkownik Dave Evens napisa:
Is it valid to use the grubbs.test in this way?
I'm very happy that someone is interested in my new package, but I must
worry you, that Grubbs test is probably not proper in such case. Your
data are dependent (time series) and possibly
Dear chip,
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of array chip
Sent: Thursday, April 14, 2005 1:00 PM
To: John Fox
Cc: R-help@stat.math.ethz.ch
Subject: RE: [R] multinom and contrasts
Dear John,
Thanks for the answer! In my own
Dear R-help..
I am rather new in R so i would appreciate your help in my problem..
I have 3 types of vegetation (A,B,C),50 measurements per class and 100 variables
per measurement.
I would like to perform seperability analysis between these classes meaning...
a.)create the hypercube from these
STATISTICIAN / DATA ANALYST
Centocor, an operating company of Johnson Johnson, seeks a highly
motivated
statistician/data analyst to work in its newly formed Nonclinical Statistics
group. We support innovative science in the discovery and research of
biotechnology therapies, with obligatory
On Thu, 14 Apr 2005 [EMAIL PROTECTED] wrote:
The `centers' are the means? by() can find the mean of multivariate data
by group. And dist() finds Euclidean and other distances.
However, the Jeffries-Matusita distance depends on covariance matrices,
and 50 points in 100 dims are not enough to
Hi all,
I try to develop my own R package. I have a couple of standart functions like
dim() and length() overloaded.
#Example
dim.Model - function(this) {
length(unique(this$.variables));
}
I built my package, but when I try to load it... This message appears:
Attaching package 'mudim':
Dear John,
My dataset has a response variable with 6 levels, and
12 independent variables, 10 of them are continuous
variable, one is a categorical variable with 2 levels,
the other one is also a categorical variable with 4
levels. total 206 observations. I attached my dataset
with the email
On Thu, 14 Apr 2005, [iso-8859-2] Václav Kratochvíl wrote:
I try to develop my own R package. I have a couple of standart functions
like dim() and length() overloaded.
Hmm. Not with the name dim.Model. You have defined an S3 method for
class Model, possibly not intentionally, and without
Dear array,
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of array chip
Sent: Thursday, April 14, 2005 3:35 PM
To: John Fox
Cc: R-help@stat.math.ethz.ch
Subject: RE: [R] multinom and contrasts
Dear John,
My dataset has a response variable
I've written a Smalltalk application which starts R on another machine and
then communicates with it via sockets. The functions involved are:
slaveStart - function(masterHost, masterPort){
inPort - availablePort()
sendSocket(as.character(inPort), masterHost, masterPort)
socketSlave(inPort)
Hi Andrew,
Bear with me as it's a while since I did this and I was new to R at
the time, but lme is probably what you're after. Remember that you're
actually not all that interested in _individual_ variance, because FA
is a sample-level property.
You'll want to set up something like:
Treatment
As far as I know, there are 2 libraries for latent class regression,
flexmix and mmlrc. Since I don't have experience with either one, can
someone give me some advice which library is better?
Thank you so much.
__
R-help@stat.math.ethz.ch mailing list
Hi there,
I am trying to use the source code for lowess (lowess.c in the src\appl
directory of R sources).
The problem is that some other files are included in lowess.c that I
don't know where to find them!
Any idea what I need to do to get able to use the source code?
Thanks in
[Prof Brian Ripley]
[François Pinard]
Now using this line within `~/.Renviron':
R_LIBS=/home/pinard/etc/R
my tiny package is correctly found by R. However, R does not seem to
see any library within that directory if I rather use either of:
R_LIBS=$HOME/etc/R
R_LIBS=$HOME/etc/R
If you are using only lowess, and nothing else in R, you might as well use
the code on netlib.
Andy
From: Mehrnoush Khojasteh
Hi there,
I am trying to use the source code for lowess (lowess.c in
the src\appl
directory of R sources).
The problem is that some other files are included
Dear all,
Here is a simplified version of a function I made:
plotfunc-function(x){
#x a vector
cat(please select two points,\n)
plot(x)
points-locator(2)
return(points)
}
Using R version 1.01 for mac os x (aqua GUI)
I would like to know what should I do to make
sure that
Hi there,
Is there an implementation of the Gaveh Stehfest algorithm in R
somewhere ? Or some other inversion ?
Thanks,
Tolga
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Dear all I hope you haven't received this message twice,
Here is a simplified version of a function I made:
plotfunc-function(x){
#x a vector
cat(please select two points,\n)
plot(x)
points-locator(2)
return(points)
}
Using the last R version 2.01 for mac os x (v1.01 aqua
This is just a test!
Please disregard
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Dear R help
I am having difficulty doing a biplot of the first two factors of a factor
analysis. I presume it is because the values in factor 2 for Milk and NUTS
are not displayed in the component loadings.
Loadings:
Factor1 Factor2
RedMeat0.561 -0.112
WhiteMeat 0.593 -0.432
Dear R
When I go to do the biplot
biplot(eurofood.fa$scores, eurofood$loadings)
Error in 1:p : NA/NaN argument
I think this is because the component loadings don't show values for some
variables
Loadings:
Factor1 Factor2
RedMeat0.561 -0.112
WhiteMeat 0.593 -0.432
Eggs
__
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-project.org/posting-guide.html
Hello
How do I use function 'MERGE to combine the FILE A and FILE B below to make
FILE C?
Thank you
FILE A
140151167
30.1 11.4 40
FILE B
140 167
5.7 30.3
FILE C
140
Dear Marc,
Excellent, this is exactly what I was looking for. The only thing I had
to change in your code was turning 'short.labels' into a factor.
Many thanks and best regards,
Jan
Marc Schwartz wrote:
Building on Tom's reply, the following should work:
labels - factor(paste(This is a long
Professor Ripley
I am very hourned to use R after all your hard work.
It looks as if I can't see the paths
/users/ripley/mingw even though I have set the HEADER
correct.
POINT:
I tried to include the missing header file float.h
from dynload.c directly. It can't see the file!!!
Is looks as if I
dat1 - data.frame(var1=c(140, 151, 167), var2=c(30.1, 11.4, 40))
dat2 - data.frame(var1=c(140, 167), var3=c(5.7, 30.3))
merge(dat1, dat2, all=TRUE)
var1 var2 var3
1 140 30.1 5.7
2 151 11.4 NA
3 167 40.0 30.3
Matt Austin
Statistician
Amgen
One Amgen Center Drive
M/S 24-2-C
Thousand
Briggs, Meredith M [EMAIL PROTECTED] wrote on 2005-Apr-15:
Hello
Hi,
How do I use function 'MERGE to combine the FILE A and FILE B below to make
FILE C?
Thank you
FILE A
140151167
30.1 11.4 40
FILE B
I think this might have been my code
mapply(paste,strwrap(levels(ncdata$Chapter),18,simplify = FALSE),collapse =
\n)
Tom
-Original Message-
From: Jan P. Smit [mailto:[EMAIL PROTECTED]
Sent: Thursday, 14 April 2005 5:15 PM
To: Mulholland, Tom
Cc: r-help@stat.math.ethz.ch
Subject:
Greetings:
I am trying to build R-2.0.1 from source on windows. My path is set to:
.;C:\RStools;C:\MinGW\bin;C:\perl\bin;C:\texmf\miktex\bin;C:\HTMLws\;C:\R201\R201\bin;%System
Root%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;C:\Program Files\
Common Files\Adaptec Shared\System;C:\LINGO9\
Is there a function in R that corresponds to the
function ``map'' of perl?
It could be called like:
vector.a - map( vector.b, FUN, args.for.FUN )
It should execute for each element ele.b of vector.b:
FUN( vector.b, args.for.FUN)
It should return a vector (or data.frame) of the
84 matches
Mail list logo