G'day Spencer,
SG == Spencer Graves [EMAIL PROTECTED] writes:
SG I'm not qualified to make this suggestion since I'm incapable
SG of turning it into reality, [...]
This statement applies to me too, nevertheless I would like to point
out the following GPL library:
Hi
On 26 Jan 2006 at 8:36, stefano iacus wrote:
From: stefano iacus [EMAIL PROTECTED]
Date sent: Thu, 26 Jan 2006 08:36:49 +0100
To: Sylvain Charlat [EMAIL PROTECTED]
Copies to: r-help@stat.math.ethz.ch, r-sig-mac
Hi all,
I downloaded and installed JGR. Then when I tried to load the package in R
console, it generated the following error message:
local({pkg - select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
Loading required package: rJava
using Java
Hi all,
I successfully installed Rcmdr. And I type library(Rcmdr), nothing
happened;
or if I select menu item load package and select Rcmdr, still nothing
happened...
Why didn't Rcmdr start?
Very strangely, if I close the R console and restart R console, every first
time I load Rcmdr, it
Try using
Commander()
but, obviously, the first time.
HIH,
Ste
On Thu, Jan 26, 2006 at 01:11:23AM -0800, Michael wrote:
MichaelHi all,
Michael
MichaelI successfully installed Rcmdr. And I type library(Rcmdr), nothing
Michaelhappened;
Michael
Michaelor if I select menu item load package and
I am afraid that you must be a little bit more verbose on your
problem. You do not provide enough information to spot your problem.
Best,
Philippe Grosjean
Michael wrote:
Hi all,
I successfully installed Rcmdr. And I type library(Rcmdr), nothing
happened;
or if I select menu item load
Commander()
Error: couldn't find function Commander
On 1/26/06, Stefano Calza [EMAIL PROTECTED] wrote:
Try using
Commander()
but, obviously, the first time.
HIH,
Ste
On Thu, Jan 26, 2006 at 01:11:23AM -0800, Michael wrote:
MichaelHi all,
Michael
MichaelI successfully installed
Hi all,
Are there any R addon/pluggins with the following feature:
(1) command history? even stores the command history many days ago? Like
Matlab does?
(2) online help? for example, as I see it, the Rcmdr is a good companian for
a newbie like me who just touched R for 1.5 days. I can use Rcmdr
Ok, let's see:
1) library(Rcmdr) - RCommander should start (you say it does)
2) Close R Commander but NOT R.
3) To restart RCommander use Commander()
Is this what you are looking for?
HIH,
Ste
On Thu, Jan 26, 2006 at 01:32:01AM -0800, Michael wrote:
Michael Commander()
MichaelError: couldn't
I'm new to both R and to this list and would like to get
advice on how to build generalized additive models in R.
Based on the description of gam, which I found on the R
website, I specified the following model:
model1-gam(ST~s(MOWST1),family=binomial,data=strikes.S),
in which ST is my
Hi all,
I have a file of this kind:
entry0001:AB0032,CF32134,DF34334
entry0002:AB0033
entry0003:AB0032,CF32134,DF34334,DD343434,DD34222
entry0004:AB0032,CF32134
And I would like to read it into something like a hash, so that I can
then loop over it by keys and values.
I wonder which would be
Michael wrote:
Hi all,
Are there any R addon/pluggins with the following feature:
(1) command history? even stores the command history many days ago? Like
Matlab does?
See ?history, and also, menu entries File - Load/Save history... in
Rgui (you use R under windows, isn't it?). also, if
Which R version and which os is yours?
For me it works fine.
Lade nötiges Paket: grDevices
using Java Runtime version 1.5.0
using JAVA_HOME = C:\develope\jre
Lade nötiges Paket: rJava
Lade nötiges Paket: methods
using Java Runtime version 1.5.0
using JAVA_HOME = C:\develope\jre
R.Version()
On 1/25/06, Michael [EMAIL PROTECTED] wrote:
Hi all,
I am really new to the R language. I am a long time Matlab and C++ user
and I was forced to learn R because I am taking a statistics class.
I am seeking to reduce the learning curve to as smooth as possible.
Are there any
I upgraded and it still does not work.
I want something like this for 3d-persp() plots:
days=c(2006-01-23,2006-01-24,2006-01-25,2006-01-26,2006-01-27,2006-01-28)
sq=(1:6)^2
plot(x=as.Date(days, format=%Y-%m-%d),y=sq, type=l, main=What I
learn about R, sub=R.version.string)
sq3d=matrix(nrow = 4,
Dear Latticers,
I want to position two trellis plots of different forms
side-by-side. The plot types are slightly different,
aspect=1 required, but panels should look the same in
both plots. Current workaround uses a guessed factor.
Any way to improve this?
Dieter
library(lattice)
n1 = 20
On Thursday 26 January 2006 07:56, [EMAIL PROTECTED] wrote:
While symbolic computation is handy, I actually think a more pressing
addition to R is some kind of automatic differentiation facility,
particularly 'reverse mode' AD, which can be spectacular. There are
free tools available for it
Notice the trailing commas.
The Contains: field is space-delimited (see the example) and I surmise you
used commas.
You also seem to have a typo: 'cwhool' not 'cwhtool'.
On Thu, 26 Jan 2006, Christian Hoffmann wrote:
Hi,
Sorry to bother, but I checked around and did not succed creating a
persp(x=1:4,y=1:6,z=sq3d,theta = 30, phi = 30, expand = 0.5, col =
lightblue)
(but y=days)
of course it should be something like
persp(y=as.Date(days, format=%Y-%m-%d),x=1:4,z=sq3d, theta=30,
phi=30, expand=0.5, ticktype=detailed, col=seagreen)
But this does not what I expected (after
Hi
On 26 Jan 2006 at 11:15, Philippe Grosjean wrote:
Date sent: Thu, 26 Jan 2006 11:15:08 +0100
From: Philippe Grosjean [EMAIL PROTECTED]
Organization: SciViews UMH - EcoNum
To: Michael [EMAIL PROTECTED]
Copies to:
The following might be what you want (replace clipboard with your
filename):
mylist - strsplit(readLines(clipboard), :)
nm - sapply(mylist, [, 1)
mylist - lapply(mylist, [, -1)
names(mylist) - nm
mylist - lapply(mylist, function(s) strsplit(s, ,)[[1]])
mylist
$entry0001
[1] AB0032 CF32134
Hi,
x and y are (numeric) vectors. I wonder if one of the following is more
efficient than the other:
x%*%y
or
sum(x*y)
?
Thanks,
Dimitri Szerman
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
[EMAIL PROTECTED] wrote:
Hi,
x and y are (numeric) vectors. I wonder if one of the following is more
efficient than the other:
x%*%y
or
sum(x*y)
?
I'd try
x - rnorm(100)
y - rnorm(100)
system.time(x%*%y)
system.time(sum(x*y))
and finally (hint, hint!):
I have always been using ebitbin to set the 3GB switch in the windows
binary, but version 2.2.1 has this set as default (which I verified using
dumpbin). However, when I generate junk data to fill up my memory and read
the memory usage using gc(), it seems that I am not getting as good results
Excellent, but...
x - rnorm(100)
y - rnorm(100)
system.time(x%*%y)
[1] 0.03 0.00 0.03 NA NA
system.time(sum(x*y))
[1] 0.05 0.00 0.04 NA NA
system.time(crossprod(x, y))
[1] 0 0 0 NA NA
So, to paraphrase a well-known contributor on this mailing list:
Excellent! So,
Philippe Grosjean wrote:
Excellent, but...
x - rnorm(100)
y - rnorm(100)
system.time(x%*%y)
[1] 0.03 0.00 0.03 NA NA
system.time(sum(x*y))
[1] 0.05 0.00 0.04 NA NA
system.time(crossprod(x, y))
[1] 0 0 0 NA NA
So, to paraphrase a well-known contributor on
On 1/26/06, Dieter Menne [EMAIL PROTECTED] wrote:
Dear Latticers,
I want to position two trellis plots of different forms
side-by-side. The plot types are slightly different,
aspect=1 required, but panels should look the same in
both plots. Current workaround uses a guessed factor.
Any way
Dear Michael,
To elaborate slightly: You can't load the same package twice in an R
session. As mentioned, used Commander() to restart the Rcmdr GUI.
I hope this helps,
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
The problem is your test. You are trying to re-allocate large objects,
and memory fragmentation will take its toll, and be almost random in its
effects.
Using the CRAN binary of 2.2.1 and
for(i in 1:1000) assign(paste(r, i, sep=.), rnorm(1e6))
gc()
used (Mb) gc trigger (Mb)
Thanks, Deepayan
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
If that's the only reason, I would suggest rbind-ing them and then use
scales = list(x = free)
I probably will go for this.
If you want the first two columns to have the same x-limits, you can specify
xlim =
Dear R users
I'm having problems in putting special characters (like á, ã, ç ) in my
plots, as much in titles, as in axis names, as in legend...when I save them as
a pdf document. They don't appear...
I don't know if it is because I'm using a linux platform...
The script is the following:
Dear Alberto,
There are fisheries people also in Europe using AD Model Builder (Denmark and
England for instance),
but you are probably right that it is more widespread in North America. There
is also effort
going on where people try to make assessment models written in ADMB callable
from R.
Hi, Berton:
thanks for getting back to me.
I played around cor.rob(). Yes, I can get a robust correlation
coefficient matrix based on mcd or mve outlier detection methods.
I have two further questions:
1) How do I get a p value of the robust r?
2) What I mean by resampling is leave one out
Perhaps you guys should try to benchmark that test on an nVidia GPU?
;-)
Best regards,
Marc
On Thu, 2006-01-26 at 15:55 +0100, Uwe Ligges wrote:
Philippe Grosjean wrote:
Excellent, but...
x - rnorm(100)
y - rnorm(100)
system.time(x%*%y)
[1] 0.03 0.00 0.03 NA
Below
Hi, Berton:
thanks for getting back to me.
I played around cor.rob(). Yes, I can get a robust
correlation coefficient matrix based on mcd or mve outlier
detection methods.
I have two further questions:
1) How do I get a p value of the robust r?
A p-value for
On 1/26/06, Dieter Menne [EMAIL PROTECTED] wrote:
Thanks, Deepayan
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
If that's the only reason, I would suggest rbind-ing them and then use
scales = list(x = free)
I probably will go for this.
If you want the first two columns to
Thanks so much. The more I thought about it the more I suspected it was the
test. With your test I was able to use almost 2.9GB.
gc()
used (Mb) gc trigger (Mb) max used (Mb)
Ncells2492256.7 467875 12.5359.4
Vcells 380096322 2900.0 380460887 2902.7
Hello R-helpers,
I'm new to XML. I have been using an application for some time, but now I
wish to automate my downloads with R. When I use the web interface of their
XML application, I'm able to read the response in R with the XML package.
My problem now is to send the requests directly from R.
One more thing ...
I played around cor.rob(). Yes, I can get a robust correlation
coefficient matrix based on mcd or mve outlier detection methods.
I have two further questions:
You might call it semantics, but I prefer resistant estimation to outlier
detection methods. I recognize
You haven't told us much, for example not your R version nor your exact
platform (nor if this is a binary or self-compiled version of R) nor your
locale.
And you (or something between you and me) have wrapped the R code onto
just two lines. You talk about saving as pdf, and ran a jpeg
Hi list,
I have trained a logistic regression model on three diverse features.
Now I want to compare these values to get an idea of their influence on
the result. For most features the summary function returns p-values that
equal zero, which makes the p-values unusable.
I have the intuition it
The cor function can do spearman correlation using
method = spearman .
On 1/25/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hi, there:
As you all know, correlation is not a very robust procedure. Sometimes
correlation could be driven by a few outliers. There are a few ways to
improve the
Dear R users
I'm having problems in putting special characters (like á, ã, ç ) in my plots,
as much in titles, as in axis names, as in legend...when I save them as a pdf
OR JPEG document. They don't appear IN BOTH...
I don't know if it is because I'm using a linux platform (MANDRIVA 2006)...
Gabor:
Contrary to popular belief, rank-based procedures are **not** resistant.
Example:
x-c(1:10,100);y-c(1:10+rnorm(10,sd=.25),-100)
cor(x,y)
[1] -0.9816899 ## awful
cor(x,y,method='spearman')
[1] 0.5 ## better
require(MASS)
Loading required package: MASS
[1] TRUE
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
Documentation said that using panel.width with non-standard aspect (which
must
be 1) leads to undefined effects.
Yes, because panel.width and panel.height together define the aspect,
and will override whatever the 'aspect' argument
It's inappropriate to use the p-values associated with the coefficients for
that purpose. They tell you how sure you can be, based on the data you
have, that they are not equal to zero.
Prof. Fox's `effects' package on CRAN would be a better place to look, I
think. I believe there's an JSS
Hi group,
My case has N physicians with each seeing M patients.
One physician could have seen a group of patients, or,
a patient could have been seen by multiple number of
physicians. In order words, there are overlaps. Now,
I have the following NxM matrix
Patient#1 Patient#2
I have a data management problem which exceeds my meager R programming
skills and would greatly appreciate suggestions on how to proceed? The
data consists of a series of observation periods. Specific behaviors are
recorded for each time period in the order each is observed. Their are
8
HI R users
I have one question for using DOS command through system
I like to delete a file that is located at C:\Program
Files\DOSPROGRAM\input.dat
I can use a DOS command del on Dos prompt like this
C:\Documents and Settings del C:\Program Files\DOSPROGRAM\input.dat
to delete input.dat file.
thanks everyone, another problem is I cannot load Rcmdr at all in SCIView...
what might be the problem?
On 1/26/06, John Fox [EMAIL PROTECTED] wrote:
Dear Michael,
To elaborate slightly: You can't load the same package twice in an R
session. As mentioned, used Commander() to restart the
On Thu, 26 Jan 2006, Dale Steele wrote:
The data looks like:
--
icsrvepf
fpevrsci
ics
p
f
ic
--
I would like to convert the about to a matrix of the form:
i c s r v e p f
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 0 0 0 0 0
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1
1
Mine. What might be wrong?
R.Version()
$platform
[1] i386-pc-mingw32
$arch
[1] i386
$os
[1] mingw32
$system
[1] i386, mingw32
$status
[1]
$major
[1] 2
$minor
[1] 2.1
$year
[1] 2005
$month
[1] 12
$day
[1] 20
$svn rev
[1] 36812
$language
[1] R
On 1/26/06, Michael [EMAIL PROTECTED]
Greetings,
I have a set of bivariate data: one variable (vegetation type) which is
categorical, and one (computed annual insolation) which is continuous.
Plotting veg_type ~ insolation produces a nice overview of the patterns that
I can see in the source data. However, due to the large number
I found that read.table() does not work well on lines that is longer
than 236 bytes.
Is this know problem? Is it under fixing?
Regards,
Jincai Jiang
(Office) 212-761-3984
NOTICE: If received in error, please destroy and notify sender.
UC Davis has a statistical department, I would suggest you get consulting
help from them. Do they have a consulting service?
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P.
Jiang, Jincai (Institutional Securities Management) wrote:
I found that read.table() does not work well on lines that is longer
than 236 bytes.
Is this know problem? Is it under fixing?
Regards,
Jincai Jiang
(Office) 212-761-3984
Hi, Jincai,
What does not work well mean? 236 bytes
I can confirm that AD Model Builder is used at the Danish
Institute for Fisheries Research, and by fisheries people in
and all around the Pacific.
On a few occasions I have solved a likelihood optimization
problem in AD Model Builder, and then wrapped the binary
in an R-package, with data
Thomas Kaliwe wrote:
Hi,
I've been looking for Image Processing packages. Thresholding, Edge
Filters, Dct, Segmentation, Restoration. I'm aware, that Octave,
Matlab
etc. would be a good address but then I'm missing the statistical
power of R. Does anybody know of packages, projects etc.
Dear R Specialists,
when trying to cluster a data.frame with about 80.000 rows and 25 columns I get
the above error message. I tried hclust (using dist), agnes (entering the
data.frame directly) and pam (entering the data.frame directly). What I
actually do not want to do is generate a
Dear all,
I've a problem with bootstarp package... I want to bootstrap
correlation ratio of some data.
my code is:
f-function(d) cov(d$x,d$y)/sqrt(var(x)*var(y))
boot(d,f,R)
As a result I got an error:
Error in statistic(data, original, ...) : unused argument(s) (..
.)
d is a
You might prefer boxplot(insolation~veg_type) as a graphic. That will
give you quantiles. To get the actual numeric values you could
for (i in levels(veg_type)) {
print(i)
quantile(insolation[veg_type==i])
}
see ?quantile for more help.
Dylan Beaudette wrote:
Greetings,
I have a
Would this do?
boxplot(Sepal.Length ~ Species, iris, horizontal = TRUE)
library(Hmisc)
summary(Sepal.Length ~ Species, iris, fun = summary)
On 1/26/06, Dylan Beaudette [EMAIL PROTECTED] wrote:
Greetings,
I have a set of bivariate data: one variable (vegetation type) which is
categorical,
Why don't you use 'unlink':
unlink('c:/Program Files/DOSPROGRAM/input.dat')
If you really want to use 'del', then you have to invoke the command
processor:
system('cmd /c del c:\\Program Files\\DOSPROGRAM\\input.dat')
On 1/26/06, Taka Matzmoto [EMAIL PROTECTED] wrote:
HI R users
I have one
Here is another way:
x - scan('/tempxx.txt', what='', sep='\n', blank.lines.skip=F)
Read 7 items
x
[1] icsrvepf fpevrsci ics p f
ic
x - strsplit(x, '') # break into single characters
template - c(i=0, c=0, s=0, r=0, v=0, e=0, p=0, f=0)
mat - lapply(x, function(.l){
+
jim holtman wrote:
Why don't you use 'unlink':
unlink('c:/Program Files/DOSPROGRAM/input.dat')
If you really want to use 'del', then you have to invoke the command
processor:
system('cmd /c del c:\\Program Files\\DOSPROGRAM\\input.dat')
See also ?shell
Uwe Ligges
On 1/26/06, Taka
Hi,
I'm trying to figure out if there's an automated way to get
read.table to read in my data and *not* convert the character
columns into anything, just leave them alone. What I'm referring
to as 'character columns' are columns in the data that are quoted.
For columns of alphabetic strings
I'm having hard time understanding the computation of degrees of freedom when
runing nlme () on the following model:
formula(my data.gd)
dLt ~ Lt | ID
TasavB- function(Lt, Linf, K) (K*(Linf-Lt))
my model.nlme - nlme (dLt ~ TasavB(Lt, Linf, K),
data = my data.gd,
fixed =
On Thu, 26 Jan 2006, Douglas Grove wrote:
Hi,
I'm trying to figure out if there's an automated way to get
read.table to read in my data and *not* convert the character
columns into anything, just leave them alone. What I'm referring
to as 'character columns' are columns in the data that
Douglas Grove wrote:
Hi,
I'm trying to figure out if there's an automated way to get
read.table to read in my data and *not* convert the character
columns into anything, just leave them alone. What I'm referring
?Did you read the help page?
What about argument as.is=TRUE?
See also argument
I did read the help page, very carefully.
The colClasses argument can be used if I want
to stop and look through every data set to see
which column I need to protect. But that's what I
said that I don't want to do.
As for 'as.is', I wish it did what you suggest, but
it doesn't. If one
exportData(MU.Cost,paste(C:/RAUDSL/S,as.character(MU.Cost$Run.Id[1]),
.,as.character(MU.Cost$MU.Id[1]),.MU.PRICE.OUTPUT.txt,sep=),append
= FALSE,type=ASCII,quote=FALSE)
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
I would like to add a footnote to this graph but do not see a footnote
command in the package:lattice documentation. I would like to note the span=.8
as the footnote.
postscript(file= ”C:/Documents and Settings/dsonneborn/My
Documents/Slovak/output/pcb_tables/smooth_PCB_lines_four.ps”,
There is an `R Data Import/Export' manual to guide you on how to do data
export in R. Please consult it.
I only know vaguely what format exportData(type=ASCII) uses, and it is
probably easier to work out what you want and tailor write.table() to give
that. (I am not aware that the fine
In xyplot(), I'd like to supress the printing of the strip on certain
panels. I thought I'd do this inside the strip function, but it seems
that panel.number doesn't get passed to that. Any ideas?
Thanks,
Ben
__
R-help@stat.math.ethz.ch mailing list
After you create your xyplot use
library(grid)
panel.text(grid.locator(),label=My label)
Cheers
Francisco
PS: How is good ol' David these days?
From: Dean Sonneborn [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] footnote in postscript lattice
Date: Thu, 26 Jan 2006 15:46:28
On 1/26/06, Benjamin Tyner [EMAIL PROTECTED] wrote:
In xyplot(), I'd like to supress the printing of the strip on certain
panels. I thought I'd do this inside the strip function, but it seems
that panel.number doesn't get passed to that. Any ideas?
Well, strip has access to much more detailed
On 1/26/06, Francisco J. Zagmutt [EMAIL PROTECTED] wrote:
After you create your xyplot use
library(grid)
panel.text(grid.locator(),label=My label)
Also, look up the 'page' argument in ?xyplot.
Deepayan
Cheers
Francisco
PS: How is good ol' David these days?
From: Dean Sonneborn
Hi all,
Can someone tell me how to justify (right or left) the labels on the
branches of a dendrogram tree? I have produced a dendrogram via agnes and
plotted it with pltree. The dendrogram terminal branch labels seem to be
centre-justified by default and I was hoping to change this to left
Hi,
I have a 1-d vector x of values but it is rather noisy. I know a priori
that the density function should be monotonically decreasing, but my data
doesnt have that. Is there a way to transform the data in such a way so
that the resultant density be monotonically decreasing?
Thanks
Debayan
I see you got an error message from R. Did you have both either the
lme4 or the Matrix packages in the search path at the same time you ran
nlme to get the result you got below? If yes, please rerun with only
nlme in the search path. (This may not be necessary, but I always quite
Let's do some simple calculation: The dist object from a data set with
8 cases would have
8 * (8 - 1) / 2
elements, each takes 8 bytes to be stored in double precision. That's over
24GB if my arithmetic isn't too flaky. You'd have a devil of a time trying
to do this on a 64-bit
Dear list members,
I'm looking for a way to write nested functions similar to the
function Nest or NestList in Mathematica.
E.g.,
f-function(x) x+2*x
f(f(f(2)))
might instead be written as nest(f, 2, 3)
read as, nest function f 3 times with 2 as the initial value.
Thanks!
Manuel
Hello all R-users
_question 1_
I need to make a statistical model and respective ANOVA table
but I get distinct results for
the T-test (in summary(lm.object) function) and
the F-test (in anova(lm.object) )
shouldn't this two approach give me the same result, i.e
to indicate the same
On 1/26/2006 9:45 PM, Manuel Morales wrote:
Dear list members,
I'm looking for a way to write nested functions similar to the
function Nest or NestList in Mathematica.
E.g.,
f-function(x) x+2*x
f(f(f(2)))
might instead be written as nest(f, 2, 3)
read as, nest function f 3
From: Dave Roberts
You might prefer boxplot(insolation~veg_type) as a graphic.
That will
give you quantiles. To get the actual numeric values you could
for (i in levels(veg_type)) {
print(i)
quantile(insolation[veg_type==i])
}
see ?quantile for more help.
If you want the
From: klebyn
Hello all R-users
_question 1_
I need to make a statistical model and respective ANOVA table
but I get distinct results for
the T-test (in summary(lm.object) function) and
the F-test (in anova(lm.object) )
shouldn't this two approach give me the same result, i.e
Fellow R Enthusiasts..
I'm trying to compile R on AIX 5.2 32bit with gcc 3.3.2-5. I've
tried both the development bundle R-devel_2006-01-25.tar.gz and the
R-beta.tar.gz from about a month ago. In each instance, I'm using the
following options prior to running ./configure
Hi R users
I like to generate random numbers from the logit-normal distribution.
Are there available simple commands I can use?
I used to generate using
exp(rnorm(1000,0,1))/(1+exp(rnorm(1000,0,1)))
I am looking for a simple command such as rnorm().
Thanks in advance
TM
Code:
write.table(MU.Cost,paste(C:/RAU/R,as.character(MU.Cost$Run.Id[1]),.
,as.character(MU.Cost$MU.Id[1]),.MU.PRICE.OUTPUT.txt,sep=),append =
FALSE,quote = FALSE,sep = ,,row.names=FALSE)
Output:
Run.Id,MU.Id,Class,Av,lower,upper,Count.WO,Count.Seq,Total.WO.Value,Note
s
32,158, 1, 0.00,
Hi all,
I am wondering what's wrong with my following code?
x=read.table('c:\\.txt');
Error in file(file, r) : unable to open connection
In addition: Warning message:
cannot open file 'c:\.txt', reason 'No such file or directory'
what's wrong?
Thanks
[[alternative HTML
The error message seems pretty clear in this case. It can't
find the file. Try this instead:
x - read.table(file.choose())
and then use the resulting GUI to locate the file.
On 1/27/06, Michael [EMAIL PROTECTED] wrote:
Hi all,
I am wondering what's wrong with my following code?
Hi R users
I like to extract lower diagonal elements of a matrix in such a way like,
data[1,2], data[1,3],
, data[5,6] are extracted from a matrix called 'data'
This short script below is what I have written so far.
##
data -
try:
as.vector(as.dist(data))
Taka Matzmoto a écrit :
Hi R users
I like to extract lower diagonal elements of a matrix in such a way like,
data[1,2], data[1,3],
, data[5,6] are extracted from a matrix called 'data'
This short script below is what I have written so far.
Try
s - matrix(rnorm(36,0,1),nrow=6)
s[col(s)row(s)]
Courtesy of VR.
-Christos
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Taka Matzmoto
Sent: Friday, January 27, 2006 1:13 AM
To: r-help@stat.math.ethz.ch
Subject: [R] [Q] extracting lower diagonal
Hi R users
Is there any function or command for measureing computation time?
For example, if I like to how long it takes to generate 10 random
numbers from a normal
distribution, Is there any command I can wrap up around rnorm(10,0,1)
and returns time
in a sec unit ?
Thanks
TM
?lower.tri
hih
__
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
Selon Taka Matzmoto [EMAIL PROTECTED]:
Hi R users
I like to generate random numbers from the logit-normal distribution.
Are there available simple commands I can use?
I used to generate using
exp(rnorm(1000,0,1))/(1+exp(rnorm(1000,0,1)))
I am looking for a simple command such as
Selon Taka Matzmoto [EMAIL PROTECTED]:
Hi R users
Is there any function or command for measureing computation time?
For example, if I like to how long it takes to generate 10 random
numbers from a normal
distribution, Is there any command I can wrap up around rnorm(10,0,1)
and
On Thu, 26 Jan 2006, Taka Matzmoto wrote:
Hi R users
I like to generate random numbers from the logit-normal distribution.
which is not a standard distribution. I guess you mean a random variable
whose logit is normal (but that is not the only meaning in use).
Are there available simple
Hi all,
I am using xyplot (lattice) to generate a figure like Figure 4.18 in
MASS4, but I have the following two questions (1) how to change the font
of x(y)lab? (2) how to plot the panels for each level of the conditional
variable (a factor in my data set) in the order as occurred in the
100 matches
Mail list logo