Thanks a lot.
It does work. The fitted data match the simulated ones well. Even no need the shifted
or scaled version of Chi-squared pdf. Also, I have tested the case of non-independent
RVs,generated by linear combining of independent Gaussian RVs,the result is
satisfactory too.
Regards,
J.Sun
Try ?reshape
A
-Original Message-
From: Christian Schulz [mailto:[EMAIL PROTECTED]
Sent: Wednesday, 24 September 2003 3:42 PM
To: [EMAIL PROTECTED]
Subject: [R] data.frame with duplicated id's
Hi,
is there a exstisting function (..i found nothing until now.)
what makes it possible tr
Hi,
is there a exstisting function (..i found nothing until now.)
what makes it possible transfrom a dataset:
ID AGE V.MAI V.JUNE
11 20 100 120
12 30 200 90
into
IDAGEV
1120 100
1120 120
1230200
123090
,or have i to programm th
Paul Meagher wrote:
Below is the test I ran awhile back on invoking R as a system call. It
might be faster if you had a c-extension to R but before I went that route I
would want to know 1) roughly how fast Python and Perl are in returning
results with their c-bindings/embedded stuff/dcom stuff, 2
On 24 Sep 2003 at 11:32, Patrick Connolly wrote:
read.table has an argument blank.lines.skip, which is true by
default, at least in read.table. You can try to give that to
read.delim?
Kjetil Halvorsen
> Say we have a tab delimited file called bug.txt
>
> Part Rep CageHb pupae
> 1
I am using 1.7.1 in a PC ( redhat 9, linux )
I created a subroutine in C++, mySubrutine, to be used in R. To debug this
subroutine I have a main routine in C++ that calls mySubrutine. The only thing main ()
does is to provide mySubrutine with its
arguments
(this is the easiest way for me to deb
Marc Schwartz <[EMAIL PROTECTED]> writes:
> On Tue, 2003-09-23 at 18:08, A.J. Rossini wrote:
<-- about confusion -->
Anyway, the mnemonic to use seems to be to remember that complex data
structures get simplified whenever possible. Thanks to Doug G. and
Patrick C. for that!
I'll plead an exc
Hi Zitan,
Below is the test I ran awhile back on invoking R as a system call. It
might be faster if you had a c-extension to R but before I went that route I
would want to know 1) roughly how fast Python and Perl are in returning
results with their c-bindings/embedded stuff/dcom stuff, 2) whether
On Tue, 23 Sep 2003 16:16:38 -0500, you wrote:
>Dear R People:
>
>I'm trying to install R 1.7.1 for Windows from Source.
>
>The error that I get is:
>previous declarion of 'ssize_t'
>MAKE[2]: ***[internet.o]Error 1
>MAKE[1]: ***[all]Error 1
>MAKE: *** [rmodules] Error 2
>
>Any ideas on how to proc
On Wed, 2003-09-24 at 11:04, Carmen Fridell wrote:
> I can't seem to find the command to find the least squares regression line
> for my bivariate data set. Can you please help? ~Carmen
Any time you're lost in R, you can type "help.start()". This will start
a web browser, which loads the startin
On Tue, 2003-09-23 at 18:08, A.J. Rossini wrote:
> In playing around with data.frames (and wanting a simple, cheap way to
> use the variable and case names in plots; but I've solved that with
> some hacks, yech), I noticed the following behavior with subsetting.
>
>
> testdata <- data.frame(matr
Have you investigated the drop= argument to "["? (as in the expression
testdata[,2,drop=F], which will return a dataframe).
"[.data.frame" has somewhat different behavior from "[" on matrices with
respect to the drop argument: If the result would be a dataframe with a
single column, the default
Paul Meagher <[EMAIL PROTECTED]> wrote:
>
> I am looking at developing a user modelling type app (new data points coming
> in and wanting to dynamically update regression co-efficients for each user)
> which could be viewed as a type of control problem.
>
Alan Miller's AS274 (in C, f77 or f90) doe
Say we have a tab delimited file called bug.txt
PartRep CageHb pupae
1 1 S 32
1 M 34
L 42
2 S 36
M 28
L 36
In playing around with data.frames (and wanting a simple, cheap way to
use the variable and case names in plots; but I've solved that with
some hacks, yech), I noticed the following behavior with subsetting.
testdata <- data.frame(matrix(1:20,nrow=4,ncol=5))
names(testdata) ## expect labels, ge
I can't seem to find the command to find the least squares regression line
for my bivariate data set. Can you please help? ~Carmen
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
On Tue, 2003-09-23 at 15:28, [EMAIL PROTECTED] wrote:
> Dear Peter,
>
> I don't know if this is the proper way to ask for help installing R, but if
> not, I presume you can pass this on to the appropriate place.
>
> I'm trying to install the latest binary on Redhat 9, and keep gett
Hi James,
Thanks for your response :-)
- Original Message -
> It is like anything else that you want to run as part of web services:
what
> do you want it to do? Yes, it is fast in doing computations, but what
will
> you have it do? It is probably as fast as anything else that you will
Here are all of the messages:
In file included from internet.c:858:
sock.h.27:conflicting types for 'ssize_t'
c:/mingw/include/sys/types.h:119: previous declaration of 'ssize_t'
MAKE[2]: ***[internet.o] Error 1
MAKE[1]: ***[all] Error 1
MAKE: ***[rmodules] Error 2
___
On Wed, 2003-09-24 at 03:06, Millo Giovanni wrote:
> Dear all,
>
> I was wondering whether there is a way to make loops in Sweave, i.e. for example to:
> 1) calculate a parameter, say, a=length(b)
> 2) according to that, add #a# chapters to the document, each including some
> repetitive analysis,
Dear R People:
I'm trying to install R 1.7.1 for Windows from Source.
The error that I get is:
previous declarion of 'ssize_t'
MAKE[2]: ***[internet.o]Error 1
MAKE[1]: ***[all]Error 1
MAKE: *** [rmodules] Error 2
Any ideas on how to proceed, please?
thanks in advance!
Sincerely,
Erin Hodgess
m
Dear Peter,
I don't know if this is the proper way to ask for help installing R, but if
not, I presume you can pass this on to the appropriate place.
I'm trying to install the latest binary on Redhat 9, and keep getting the same
error message no matter what I try.
[EMAIL PRO
I am posting it here because I am not sure if the behavior described below
is actually a bug.
If I do something like this:
> x1 <- 1:3
> x2 <- 5:9
> x3 <- 21:27
> ll <- list(x1, x2, x3)
> stack(x1,x2,x3)
I get the following error:
Error in rep.int(names(x), lapply(x, length)) :
invalid
On 09/23/03 15:49, K Skanes wrote:
>Hi,
>
>I have a data set with 7 years worth of data, and 2 different values (a real
>value and a model value) of interest in each year for different lengths. I
>would like a plot with the year on the y axis and an animal length along the
>x axis. For each ye
Hi,
I have a data set with 7 years worth of data, and 2 different values (a real
value and a model value) of interest in each year for different lengths. I
would like a plot with the year on the y axis and an animal length along the
x axis. For each year I would like to see two lines, one rep
Your best bet is saving as a comma-separated value file (.csv) and using
read.csv to get the data into R.
Ben
On Tue, 23 Sep 2003 [EMAIL PROTECTED] wrote:
> Hi,
>
> I would like to know how to extract the data from Excel Spreadsheet.
>
> Thank you very much.
>
> Melissa
>
>
>
Using Thomas Unternährer's handy example, one could also do:
> X <- c(1, 4.5, 2.3, 1, 7.3)
> mean(order(X, decreasing=TRUE)[1:2])
[1] 3.5
>
I think this will give the same results as Thomas Unternährer's suggested
code in almost all cases, but it is perhaps more concise and direct
(provided that
there are plenty of discretization methods, which one are you looking for?
(In Spanish
hay bastantes metodos de discretizacion cual de ellos estas busacando?
Regards (saludos)
Edgar Acuna
On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote:
> Hi R users
>
> I need to apply discretization to my conti
[EMAIL PROTECTED] wrote:
Hi,
I would like to know how to extract the data from Excel Spreadsheet.
I would like to know whether you have read the manual "R Data
Import/Export" before having posted the question. It tells about you
more than one way.
Uwe Ligges
Thank you very much.
Melissa
___
Hi,
I would like to know how to extract the data from Excel Spreadsheet.
Thank you very much.
Melissa
JLT Risk Solutions Ltd
6 Crutched Friars, London EC3N 2PH. Co Reg No 1536540
Tel: (44) (0)20 7528 4000 Fax: (44) (0)20 7528 4500
> "UweL" == Uwe Ligges <[EMAIL PROTECTED]>
> on Tue, 23 Sep 2003 11:30:02 +0200 writes:
UweL> Jan Kleinn wrote:
>> Dear all,
>>
>> I would like to make a filled contour plot without the
>> box R is generating by default around the plotting area,
>> i.e. I'm loo
On Tuesday 23 September 2003 16:38, Michael Shott wrote:
> Dear R-Help:
>
> As you can see, Prof. Friendly refers me to your site for an
> executable version of vcd. I don't mean to be obtuse, but 15
> minutes spent exploring your site failed to locate a downloadable
> version of the vcd package
This is a relatively recent article that is somewhat accessible.
Jensen, D. R., and Solomon, Herbert (1994), "Approximations to joint
distributions of definite quadratic forms", Journal of the American
Statistical Association, 89 , 480-486
It has references to previous work.
I also have an old pa
Thank you all for the reply,
the Hmisc library is exactly what i was looking for.
Johannes
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Michael Shott wrote:
Dear R-Help:
As you can see, Prof. Friendly refers me to your site for an
executable version of vcd. I don't mean to be obtuse, but 15 minutes
spent exploring your site failed to locate a downloadable version of the
vcd package to which he referred.
I know plainly what t
Dear all,
I was wondering whether there is a way to make loops in Sweave, i.e. for example to:
1) calculate a parameter, say, a=length(b)
2) according to that, add #a# chapters to the document, each including some repetitive
analysis, each time done on a particular subset of the data indexed by t
Consider:
?cut
?round
With more detail we might be able to help more ...
On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote:
> Hi R users
>
> I need to apply discretization to my continuous data.
> Is there a method in R to do this?
>
> Thanks in advance,
>
> Jaime
>
--
620B Bartram Ha
On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote:
> I need to apply discretization to my continuous data.
> Is there a method in R to do this?
See help("cut").
- tom blackwell - u michigan medical school - ann arbor -
__
[EMAIL PROTECTED] maili
Dear R-Help:
As you can see, Prof. Friendly refers me to your site for an
executable version of vcd. I don't mean to be obtuse, but 15 minutes
spent exploring your site failed to locate a downloadable version of the
vcd package to which he referred.
I know plainly what this application can do
Hi R users
I need to apply discretization to my continuous data.
Is there a method in R to do this?
Thanks in advance,
Jaime
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
On Tue, 23 Sep 2003, Federico Calboli wrote:
> Dear All,
>
> when I try to:
>
> install.packages("Design"); install.packages("Hmisc")
>
> I get the following error messages:
>
> * Installing *source* package 'Design' ...
> ** libs
> g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586
> -mcpu
On Tue, 23 Sep 2003, Jean Sun wrote:
> >From basic statistics principle,we know,given several i.i.d Gaussian
> >RVs with zero or nonzero mean,the sum of square of them is a central or
> >noncentral Chi-distributed RV.However if these Gaussian RVs have
> >different variances,what does the sum of sq
If you just want to get rid of the axes, you can do
filled.contour(x, plot.axes = { })
-roger
Uwe Ligges wrote:
Jan Kleinn wrote:
Dear all,
I would like to make a filled contour plot without the box R is
generating by default around the plotting area, i.e. I'm looking for
an option in fille
I need help computing model-predicted estimated marginal means of the dependent
variable in the cells defined by the fixed factors, including interactions. By
estimated marginal means, I just mean cell means collapsed over all other factors and
adjusted for the mean values of other covariates i
Hi,
>I would like to rank a time-series of data, extract the top ten data items from this
>series, determine the
>corresponding row numbers for each value in the sample, and take a mean of these *row
>numbers* (not the data).
>I would like to do this in R, rather than pre-process the data o
Here's one way. Suppose your "time series" is in a vector called "x".
top10 <- sort(x, decreasing=TRUE)[1:10]
mean.index <- mean(which(x %in% top10))
HTH,
Andy
> -Original Message-
> From: James Brown [mailto:[EMAIL PROTECTED]
> Sent: Tuesday, September 23, 2003 7:51 AM
> To: [EMAIL PR
> From: Murray Jorgensen [mailto:[EMAIL PROTECTED]
>
> I got a shock a few days ago when I accidentally visited
> www.r-project.com . I thought that the r-project site had been hacked
This one seems to be about some sort of city revival projects in Japan.
(The introduction starts with "Recycle
I would like to rank a time-series of data, extract the top ten data items
from this series, determine the corresponding row numbers for each value
in the sample, and take a mean of these *row numbers* (not the data).
I would like to do this in R, rather than pre-process the data on the
UNIX comm
> You need the Fortran compiler (g77) as well (as indicated by the error
> message). Is it installed? Is it in your path?
>
> Uwe Ligges
>
It was not installed. Now everything works.
Cheers,
Federico
--
=
Federico C. F. Calboli
Department of Biology
Univ
One thing you could do is to use the 'its' (irregular time-series)
package on CRAN.
e.g. using a trivial dataset
require(its)
its.format("%Y-%m-%d") #defines text format of dates in dimnames
df <- data.frame(1:3,(1:3)^2)
dimnames(df) <- list(c("2003-01-03","2003-01-06","2003-01-07"),letters[1:2]
Federico Calboli wrote:
Dear All,
when I try to:
install.packages("Design"); install.packages("Hmisc")
I get the following error messages:
* Installing *source* package 'Design' ...
** libs
g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586
-mcpu=pentiumpro -O2 -fomit-frame-pointer -p
Dear All,
when I try to:
install.packages("Design"); install.packages("Hmisc")
I get the following error messages:
* Installing *source* package 'Design' ...
** libs
g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586
-mcpu=pentiumpro -O2 -fomit-frame-pointer -pipe -march=i586
-mcpu=pen
Try this:
plot(x, y)
abline(lm(y ~ x)) #plots the fitted line
or this if you need the model results elsewhere:
mod.lm <- lm(y ~ x) # store the model
plot(x, y)# plot the data
abline(mod.lm)# plot the fitted line
see ?plot.lm and ?abline
HTH
Gav
[EMAIL PROTECTED] wrote:
H
Jan Kleinn wrote:
Dear all,
I would like to make a filled contour plot without the box R is
generating by default around the plotting area, i.e. I'm looking for an
option in filled.contour similar to 'axes=F' in 'contour' or in 'plot'.
I couldn't find any option to get rid of the box, any help
[EMAIL PROTECTED] wrote:
Hi,
I would like to enquire if by typing plot (lm(y~x)) would this show me the
plot of the fitted line? I tried this function previously but I was only
able to get the last 4 plots starting with Residuals vs fitted.
No, it shows plots for analyses of the residuals.
Try
Dear all,
I would like to make a filled contour plot without the box R is
generating by default around the plotting area, i.e. I'm looking for an
option in filled.contour similar to 'axes=F' in 'contour' or in 'plot'.
I couldn't find any option to get rid of the box, any help is welcome.
Thank
Hi,
I would like to enquire if by typing plot (lm(y~x)) would this show me the
plot of the fitted line? I tried this function previously but I was only
able to get the last 4 plots starting with Residuals vs fitted.
Thank You.
Melissa
**
Zitan Broth wrote:
Greetings All,
Been playing with R and it is very easy to get going with the UI or infile batch commands :-)
What I am wondering is how scalable and fast R is for running as part of a web service. I believe R is written in C which is a great start, but what are peoples general
"Remko Duursma" <[EMAIL PROTECTED]> writes:
> Dear R-helpers,
>
> i get some strange results using a linear mixed-effects model (lme), of the type:
>
> lme1 <- lme(y ~ x, random=~x|group, ...)
>
> For some datasets, i obtain very small standard deviations of the random effects. I
> compared th
59 matches
Mail list logo