Hi,
Is there a way I can view the functions source code of a
package I installed in my PC.
For example I downloaded the great "mixtools" package.
I want to see the source code of one of its function "normalmixEM"
Is there a way to do it? Presumably from R command prompt?
I tried to take a look
> "CG" == Christophe Genolini <[EMAIL PROTECTED]>
> on Wed, 16 Jul 2008 02:40:01 +0200 writes:
CG> Hi the list, I like the R logo (grey C and blue R) very
CG> much, specialy the drawing of the letter with border and
CG> shadow. I would like to make something closed with som
Gundala Viswanath wrote:
Hi,
Is there a way I can view the functions source code of a
package I installed in my PC.
For example I downloaded the great "mixtools" package.
I want to see the source code of one of its function "normalmixEM"
Is there a way to do it? Presumably from R command prom
The R Logo was done "manually" (using software) in several
phases IIRC.
Well, that's not a surprise. Let me rephrase my question : is the logo
drawer still on this list ? Did you use gimp ? A specific filter ? or a
scriptfu ?
Christophe
__
R-hel
Dear List,
I have Multi-level Data
i= Indivitual Level
g= Group Level
var1= First Variable of interest
var2= Second Variable of interest
and I want to count the frequency of "var1" and "var2" on the group
level.
I found a way, but there must be a much simpler way.
data.ml <-
data.frame(i=c(1:
well, I'm not completely sure I understand what you want to compute
but you may try the following:
data.ml <- data.frame(i = 1:8, g = factor(c(1,1,1,2,2,3,3,3)),
var1 = c(3,3,3,4,4,4,4,4), var2 = c(8,8,8,2,2,4,4,4))
data.ml[!duplicated(data.ml$g), -1]
# or
t(sapply(split(data
Is this what you were asking about:
> aggregate(data.ml[,c('var1','var2')], list(data.ml$g), mean)
Group.1 var1 var2
1 138
2 242
3 344
On Wed, Jul 16, 2008 at 4:21 AM, Kunzler, Andreas <[EMAIL PROTECTED]> wrote:
> Dear List,
>
> I have Multi-level Data
Hi willemf,
Glad to hear that it helped. Years ago (late-90s) I Linuxed, but have since
been forced into the Windows environment (where, however, I have the great
pleasure of being able to use MiKTeX and LyX, i.e. TeX/LaTeX). I therefore
can't help you further, except to say that I have never ha
Hi willemf,
And perhaps I should have added (in case you are moving across systems) that
you should take a look at
?embedFonts
http://cran.r-project.org/doc/Rnews/Rnews_2006-2.pdf
Ghostscript should be installed on your box, so there shouldn't be a
problem. Without additional information, it's
Hmmm, I did not follow this thread closely, sorry for that,
just want to share my 2c.
If it is about quality, then I create EPS files and use the
psfrag latex package to replace the PS fonts with TeX's fonts.
This has the following advantages:
1) The figures have the same font as the text its
On Wed, Jul 16, 2008 at 04:48:28AM -0500, Gabor Csardi wrote:
[...]
> I have a little script that automates this for .fig files (this is
> based on figtex, another script that I found somewhere online and
> can't find it any more)
[...]
Ok, it is called figfrag, and it is on CTAN.
[...]
--
Cs
I solved my problem.
I had to delete stop()
It doesn't make sense to have a stop() in if. stop() is to stop the
execution of the code once the condition is met.
the else does the job of stop() in this case.
mysimbaa wrote:
>
> Hello,
> I have an error since I run my R code. This error is :
> F
On Tue, 15 Jul 2008 09:58:28 -0700 (PDT)
A Ezhil <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I didn't try automatic dependency resolution instead installed each library
> one by one. It seems that tcl is installed and the libtcl8.4.so exists.
>
> Below is the info from the system:
> > rpm -q tcl
> tcl
On 7/15/08, Angila Albaros <[EMAIL PROTECTED]> wrote:
>
> Dear Sir,
> Thanks for your reply but my data is very huge 100 x 550 (
> for x ) and 100 x 1010 ( for y). So, I think that time , I need to take one
> column of x ($x1)and do multiple regreesion with y data set.i.e x1 will be
Dear R users,
i am using the following code to produce barcharts with lattice:
Compound<-c("Glutamine", "Arginine", "Glutamate", "Glycine", "Serine",
"Glucose", "Fructose", "Raffinose",
"Glycerol", "Galacglycerol", "Threitol", "Galactinol", "Galactitol")
Class<-c("aminos","aminos","aminos","am
Mark, your suggestion results in about 75% control over the plot. This is the
best that I have managed to get it at, so thank you very much. In Linux you
create a X11() device for screen output. Specifying identical device
characteristics results in a fair degree of similarity between screen
versi
Dear R users I have a problem with the selection of the variables from a gamma
model with log link function:
step.glm.gam1.infroma<- step(glm.gam1.infroma, direction = "backward", k=2)
Errore in glm.fit(x[, jj, drop = FALSE], y, wt, offset = object$offset, :
ciclo intero 1; impossibile co
thanks that is a slim version.
Stephen
On Tue, Jul 15, 2008 at 9:03 PM, jim holtman <[EMAIL PROTECTED]> wrote:
> Typo in the last one: (resend)
>
> #GPS in Decimal Degrees in the form longitude latitude (raw data)
>
> library(maptools)
>
> # create a list of the coords
> coordList <- list(
> RM2
Hi R-listers,
>
> I would like to know how can i extract (ncomps) component no. when the
> RMSEP is lowest for PLS package?
>
> Currently, I only plot it manually, validationplot(pls) and then only feed
> the ncomp to the jack knife command. However, I would like to automate this
> step.
>
> Pleas
You could do it without a panel function using xyplot and type "h" :
df2 <- transform(df, Compound = factor(abbreviate(Compound)))
xyplot(Ratio ~ Compound | Class, df2, type = c("h", "g"), lwd = 7,
par.settings = list(grid.pars = list(lineend = 1)))
On Wed, Jul 16, 2008 at 6:48 AM, Henni
Perhaps if you follow the posting guide more closely, you might get more
(useful) replies, but without looking at your data, I doubt there's much
anyone can do for you.
The fact that the range of the outlying measures is -1 to 2 would tell
me there are no potential outliers by this measure. Pleas
Hi,
I know that if i want to see the SPLUS or R-source code of a function,i should
give the command
without the brackets.For example:
rsquared.lmRobMM
function(x)
{
str0 <- "Initial S-estimate"
str1 <- "Final M-estimate"
if(x$est == "final") {
z <- x$r.squared
attr(z, "info") <- str1
}
Thanks anyway for your answer.
That was also an option that I took into account (no potential
outliers) and I will have a look at the "value" section of ?outliers.
B.
Am 16.07.2008 um 14:11 schrieb Liaw, Andy:
Perhaps if you follow the posting guide more closely, you might get
more
(useful
On 7/16/2008 8:29 AM, arezoo bagheri wrote:
Hi,
I know that if i want to see the SPLUS or R-source code of a function,i should give the command
without the brackets.For example:
rsquared.lmRobMM
function(x)
{
str0 <- "Initial S-estimate"
str1 <- "Final M-estimate"
if(x$est == "final") {
z
Hi the list,
When we use package.skeleton, it create some file in the man
directorie. - If we use package.skeleton with namespace=FALSE, it
create a file toto-internal.Rd
- If we use package.skeleton with namespace=TRUE, it does not create
the file toto-internal.Rd
Why is that ?
Christophe
I use a different dissimlarity measure (library(analogue);Gowers Index).
I just wanted to look if there are similar values in both "tables".
I mainly try to find a way to find the best model to explain my
predefined groups (using a bunch of different variables:
factors,count,numeric, ordered
G'day Christophe,
On Wed, 16 Jul 2008 15:10:05 +0200
[EMAIL PROTECTED] wrote:
> Hi the list,
>
> When we use package.skeleton, it create some file in the man
> directorie. - If we use package.skeleton with namespace=FALSE, it
> create a file toto-internal.Rd
> - If we use package.skeleton with
Thanks, Berwin,
So the main idea is - With NAMESPACE, you do not document the
"not-for-user" because they don't have to be documented
- Witout NAMESPACE, you document the "not-for-user" with a
toto-internal.Rd that say "not for user"
That's clear.
Is it stupid to consider to use both techniq
Michael,
If you mean Manova in the car package (or manova in stats), just use
str() on the resulting object to see the slots. However, Manova()
doesn't return F values, since these depend on the type of test
statistic used for the multivariate test.
> library(car)
> ?Manova
> soils.mod <- lm
G'day Christophe,
On Wed, 16 Jul 2008 16:11:15 +0200
[EMAIL PROTECTED] wrote:
> So the main idea is - With NAMESPACE, you do not document the
> "not-for-user" because they don't have to be documented
> - Witout NAMESPACE, you document the "not-for-user" with a
> toto-internal.Rd that say "not f
Dear R-helpers,
I was curious why the output of summary (and many other functions in
R) does not add a separator between the name of a factor and the label
of a level (e.g., in the example below, 'group Trt'). If the user had
lower case labels (e.g., 'trt'), the output would be particularly
Has anybody used fuzzy time series to forecast enrollments? I have some code
that does not work so well. Thanks.
Kate
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Dear R experts,
Â
I need help to make my little program efficient which now takes 2hrs
to complete.
Â
Using arules package I developed set of rules consisted of 900 rules. Now I
want to check whether a lower rule is a subset of a higher premises rule. And
if it is a subset of higher premise
> "MC" == Mark Cowley <[EMAIL PROTECTED]>
> on Wed, 16 Jul 2008 15:32:30 +1000 writes:
MC> Dear list,
MC> I am analysing a set of quantitative proteomics data
MC> from 16 patients which has a large numbers of missing
MC> data, thus some proteins are only detected once,
I don't know, but you can change it by:
out_sum <- summary(lm.D9 <- lm(weight ~ group))
rownames(out_sum$coefficients)[2] <- "group Trt"
out_sum
On Wed, Jul 16, 2008 at 12:51 PM, Michael Kubovy <[EMAIL PROTECTED]> wrote:
> Dear R-helpers,
>
> I was curious why the output of summary (and many othe
Dear list,
I am fitting a logistic multi-level regression model and need to test the
difference between the ordinary logistic regression from a glm() fit and the
mixed effects fit from glmer(), basically I want to do a likelihood ratio test
between the two fits.
The data are like this:
My outc
I am running R 2.7.0 on a Suse 9.1 linux cluster with a job scheduler
dispatching jobs and openmpi-1.0.1. I have tried running one of the
examples at http://ace.acadiau.ca/math/ACMMaC/Rmpi/examples.html in Rmpi
and they seem to be working, except mpi.close.Rslaves() hangs. The
slaves are closed, bu
Thanks for your answer.
I guess writing
a regular expression that says "export everything that does not start
with a dot but do not export foo and bar" would be not trivial to write
(at least not for me).
The NAMESPACE created by package.skeleton contain a single line :
exportPattern("^[[:alp
Hello R Users,
I want to convert date (yr, mo, day, hr, min, sec) to decimal date,
for example:
Date (in six columns):
yr mo dy hr minsec
1993 02 13 05 52 00
Converted to Decimal date :
1993.3542
How to write a small code in R so I can convert six column date to
Hi,
I am using ns() to model the effect of time on some outcome y
[ specifically, I am using polr() in a model of the form
mod1=polr(y~x1+x2*ns(Year,df=3),...) , with x1 and x2 denoting several
covariates each ]
I understand how to use the spline basis as recorded in the model matrix in
order
Hi all,
I am trying to use rpart.predict to test rpart models built with both
numeric and character variables. The rpart trees appear to come out fine,
but then when I try to use the rpart model to predict on a different dataset
(same variables) I get error messages that claim the character variab
On 7/16/08, Henning Wildhagen <[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> i am using the following code to produce barcharts with lattice:
>
> Compound<-c("Glutamine", "Arginine", "Glutamate", "Glycine", "Serine",
> "Glucose", "Fructose", "Raffinose",
> "Glycerol", "Galacglycerol", "Threitol
Its not clear from your post what the decimal portion is supposed
to represent since Feb is not 35% of the way through the year
but in any case see R News 4/1 and the full page table at the end
of that article, in particular, for some date idioms that you can use.
On Wed, Jul 16, 2008 at 7:41 AM,
well, for computing the p-value you need to use pchisq() and dchisq()
(check ?dchisq for more info). For model fits with a logLik method you
can directly use the following simple function:
lrt <- function (obj1, obj2) {
L0 <- logLik(obj1)
L1 <- logLik(obj2)
L01 <- as.vector(- 2 *
Mike Lawrence wrote:
I'm looking to analyze a large data set: a within-Ss 2*2*1500 design
with 20 Ss. However, aov() gives me an error, reproducible as follows:
id = factor(1:20)
a = factor(1:2)
b = factor(1:2)
d = factor(1:1500)
temp = expand.grid(id=id, a=a, b=b, d=d)
temp$y = rnorm(length(te
Hey,
I am a PhD student in forestry science, and I am a brand new in R. I am
working with lidar data (cloud points with X, Y and Z value). I wish to
create a spatial map with kriging form points cloud. My problem is the Big
data-set (over 5,000,000.00 points) and I always went out of memory.
I am trying to apply the solution you mentioned to all columns in a matrix,
and output the results to another matrix and append the two using the cbind
function.
I have written something that works, using a nested For loop to go through
all the cells in the target matrix, but the trouble is that
Here's my code:
nc<-ncol(mat) #setting end point in counter to number of cols in sm
nr<-nrow(mat)
mm <- array(NA, dim=c(2, nc)) #to hold min/max ranges
sc <- array(NA, dim=c(nr, nc)) #to hold percentile scales
for (n in 1:nc) { #now calculate respect
Try this and see how fast it runs:
# x is your input matrix that will be processed by column
# y is the scaled output matrix
y <- t(apply(x, 2, function(.col){
.rang <- range(.col)
(.col - .rang[1]) / (.rang[2] - .rang[1]) * 100
}))
On Wed, Jul 16, 2008 at 3:01 PM, rcoder <[EMAIL PROTEC
?scale
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of jim holtman
Sent: Wednesday, July 16, 2008 12:14 PM
To: rcoder
Cc: r-help@r-project.org
Subject: Re: [R] Mapping data onto score
Try this and see how fast it runs:
# x is your input
On Wednesday 16 July 2008, Alessandro wrote:
> Hey,
>
>
>
> I am a PhD student in forestry science, and I am a brand new in R. I am
> working with lidar data (cloud points with X, Y and Z value). I wish to
> create a spatial map with kriging form points cloud. My problem is the Big
> data-set (over
Hi All,
I am trying to make write a table in RSQLite. And I get the error mentioned
below
mat<-as.data.frame(matrix(rnorm(n=24400),nrow=244000,ncol=1000))
> dbWriteTable(con, "array", mat)
[1] FALSE
*Warning message:
In value[[3]](cond) :
RS-DBI driver: (error in statement: too many SQL var
Hi Folks,
I'd be grateful for good suggestions about the following.
I'm plotting a family of (X,Y) curves (for different levels
of another variable, Z): say 6 curves in all but could be
more or less -- it's a rather variables situation.
I'd like to label each curve with the value of Z that it
cor
It depends on what the values for certain variables were set to when it was
compiled. See:
http://www.sqlite.org/limits.html
On Wed, Jul 16, 2008 at 4:14 PM, Vidhu Choudhary
<[EMAIL PROTECTED]> wrote:
> Hi All,
> I am trying to make write a table in RSQLite. And I get the error mentioned
> below
Hi Spencer,
I have seen your name on the web site, and perhaps you can
help me with my R problem.
I'm trying to use KernSur to put in evidence a substructure in a
bidimensional plot. My problem is that, in order to get the density
in the low density areas (in which the substructure is located) I
Dr. Harrell's Hmisc package has labcurve.
--Matt
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ted Harding
Sent: Wednesday, July 16, 2008 1:17 PM
To: [EMAIL PROTECTED]
Subject: [R] Labelling curves on graphs
Hi Folks,
I'd be grateful for good suggestio
Many thanks, Matt. for pointing that out. It certainly looks very
promising -- but also there is an awful lot to study! [Pause, while
I head for the coffee-maker]
Ted.
On 16-Jul-08 20:31:44, Austin, Matt wrote:
> Dr. Harrell's Hmisc package has labcurve.
>
> --Matt
>
> -Original Message
On Wednesday 16 July 2008, Alessandro wrote:
> Hey Dylan,
>
> Thank you. I wish to test for my PhD: TIN (ok, with Arcmap), IDW (ok, with
> Arcmap) and kriging model (and other if it is possible) to create DSM and
> DEM, and DCM (DSM-DEM). I tried with gstat in IDRISI, but my PC outs of
> memory.
>
I don't know the origin of the logo, but as far as I remember, the
current version of the logo has been manually redrawn from the
original on photoshop to increase resolution.
stefano
On 16/lug/08, at 10:50, [EMAIL PROTECTED] wrote:
The R Logo was done "manually" (using software) in sever
Hi,
I have a matrix M, quite a few (< 1/4th) of its eigen values are of
O(10^-16). Analytically I know that M is positive definite, but
numerically of course it is not. Some of the small (O(10^-16)) eigen
values (as obtained from eigen()) are negative. It is the
near-singularity that is causing th
Hello,
is there a function to check if an object exists?
Thank you,
Judith
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide co
On 16/07/2008 6:26 PM, Judith Flores wrote:
Hello,
is there a function to check if an object exists?
exists("foo") checks for foo. (There are a few other optional
parameters to tell it where to look, and what sort of thing foo might be.)
Duncan Murdoch
Thank you,
Judith
___
Please do your homework before posting: help.search("exists") --> ?exists
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Judith Flores
Sent: Wednesday, July 16, 2008 3:27 PM
To: RHelp
Subject: [R] Function to verify existence of an R obje
"Judith Flores" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> is there a function to check if an object exists?
?exists
efg
Earl F Glynn
Stowers Institute for Medical Research
__
R-help@r-project.org mailing list
https://stat.ethz.
If we knew an elegant one, we'd have an R function to do it ...
The R Logo was done "manually" (using software) in several
phases IIRC.
But yes, indeed, an 'Rfont' package on CRAN would be neat...
Martin Maechler, ETH Zurich and R Core Team.
I am not sur of what would be neat : a police wi
2008/7/16 Ted Harding <[EMAIL PROTECTED]>:
> Hi Folks,
> I'd be grateful for good suggestions about the following.
>
> I'm plotting a family of (X,Y) curves (for different levels
> of another variable, Z): say 6 curves in all but could be
> more or less -- it's a rather variables situation.
>
> I'
Dear all, I have two vector:
strata response1 11 21 32 12 22 3. .. .. .10 110 210 3. ..
.. .100 1100 2100 3
I need to get a stratified random sample of size = 100, such that:
if I select response = 1 in strata 1 then I need sample 2 or 3 in strata 2
if I select respons
Ciao Dylan,
THANKS for your help. When I arrive in this step "V <- variogram(z~1,
d.small)", appear this note:
Error in gstat(formula = object, locations = locations, data = data) :
l'argomento "data" non è specificato e non ha un valore predefinito (data
argument it's not specified and it h
I read many literatures and try to find variance estimation of Prentice
method in case cohort. Only SelfPrentice variance estimation can be
found. Does the jacknife variance using cluster(id) in coxph() only applies
to Self Prentice? I can find using [t-e,t] to estimate Beta for Prentice.
But How
How large is your matrix?
Are the very small eigenvalues well separated?
If your matrix is not very small and the lower eigenvalues are clustered, this
may be a really hard problem! You may need a special purpose algorithm and/or
higher precision arithmetic.
If your matrix is A and there exists
running R 2.5.1 on WinXP.
plot(model formula here) produces 4 useful residual diagnostic plots. Two
questions about them:
1. How can I isolate just one of the four, and have only that one be shown?
str(plot(model)) didn't give me any insight.
2. How can I make the lines (I guess they are
On 17/07/2008, at 11:45 AM, lamack lamack wrote:
Dear all, I have two vector:
strata response1 11 21 32 12 22 3. .. .. .10 110
210 3. .. .. .100 1100 2100 3
I need to get a stratified random sample of size = 100, such that:
if I select response = 1 in strata 1 then
Moshe Olshansky yahoo.com> writes:
> How large is your matrix?
Right now I am looking at sizes between 30x30 to 150x150, though it will
increase in future.
> Are the very small eigenvalues well separated?
>
> If your matrix is not very small and the lower eigenvalues are clustered,
this may
create a vector containing 1,2,3
sample one value
remove the sampled value from the initial vector
sample again
now select the observations belonging to the strata indexed by the two
sampled values
-> not difficult. Most certainly, other approaches work too, maybe even
better.
Rolf, is he ta
hi,
I am sure this is an obvious question,
I have a optimized model with coefficients from
stepAIC. I can't figure out how to find the
max value of a function with those coefficients.
Its part of a dynamic web page or else i would
have done it by hand. The code works fine, so you probably
don't nee
Hi Bjorn,
Greetings. I would like to know how can i extract component no. when the
RMSEP is lowest?
>
> Currently, I only plot it manually and then only feed the ncomp to the jack
> knife command. However, I would like to automate this step.
>
> Please let me know. Many thanks.
>
> Rgrds,
>
On 17/07/2008, at 1:11 PM, Daniel Malter wrote:
Rolf, is he taking a course of yours? :)
No, I don't actually ``teach'' these days. (If ever I did.)
It just looked like a homework problem to me.
Still does, actually.
cheers,
Your problem seems to be pretty hard (for larger dimensions)!
When speaking about the A1,A2,...,An sequence I did not mean representing A as
a product of A1*A2*...*An. What I mean is that A1,...,An is a sequence of
matrices which starts at a "good" matrix, i.e. it's eigenvalues are well
separat
How about computing the 'lhs' of the rules once and then using that
information for the subsets.
On Wed, Jul 16, 2008 at 11:56 AM, Daniel Wagner <[EMAIL PROTECTED]> wrote:
> Dear R experts,
>
> I need help to make my little program efficient which now takes 2hrs to
> complete.
>
> Using arules pa
I am using R 2.7.1 under Windows XP, installed at C:/R/R-2.7.1.
The location of the HTML SearchEngine is
file:///C:/R/R-2.7.1/doc/html/search/SearchEngine.html. Now, when I type a
phrase, say "reshape", in the search text field, the Search Results page
suggest that the location of the reshape H
Barry Rowlingson wrote:
2008/7/16 Ted Harding <[EMAIL PROTECTED]>:
Hi Folks,
I'd be grateful for good suggestions about the following.
I'm plotting a family of (X,Y) curves (for different levels
of another variable, Z): say 6 curves in all but could be
more or less -- it's a rather variables s
Hi,
could be a browser issue. What browser are you using? I see this
behavior on Firefox 3.0.1 but not Microsoft IE 7.0(.5730.11).
Looking at the source code for the generate search (of 'reshape') you
get the relative URL:
../../../library/R.utils/html/wrap.array.html
for both browsers. The
Thanks, Henrik.
I was indeed using Firefox 3.0.1, and I now see that the issue doesn't
occur with MSIE 7.0. Will use the latter for the time being.
Best regards,
Jan
"Henrik Bengtsson" <[EMAIL PROTECTED]>
Sent by: [EMAIL PROTECTED]
17/07/2008 10:45
To
"Jan Smit" <[EMAIL PROTECTED]>
cc
r-
--- On Tue, 15/7/08, haani hanni <[EMAIL PROTECTED]> wrote:
From: haani hanni <[EMAIL PROTECTED]>
Subject: how i can install Rgraphviz in R2.7.1
To: "Nabble" <[EMAIL PROTECTED]>
Cc: [EMAIL PROTECTED]
Date: Tuesday, 15 July, 2008, 1:39 PM
hello
i am a new user of R.i have window XP Proffe
84 matches
Mail list logo