Hi,
Please tell me of an R graphical routine that will create a 3d plot
something like cloud or scatterplot3d but that allows individual
labeling of each point in the plot. I've done it with SPSS. I'm sure that
there is an R routine, but I can't find it.
I'm at:
[EMAIL PROTECTED]
[EMAIL
Greetings R users,
I am an R graphics newbie trying to produce a custom trellis plot using
pairs() with R 2.7.2.
I have spatial data on which I run a geographically weighted regression
(gwr, using the -spgwr- package). I want to check the gwr coefficients
for multicollinearity and spatial
I have the book you mentioned. It basically describes the silhouette method.
I do not have it handy as I moved so it is still in some box. However I
cannot remember that book providing any other criterion to find the best
clusters number.
On the other hand I have the same problem with hierarchical
Hello,
I am having problems installing R-2.7.2 on a maching running OpenSUSE
Linux 10.2. Initially, R couldn't find a Fortran compiler, so I
downloaded and installed compat-g77-3.3.5-38.i586.rpm. I then set
FPICFLAGS=fPIC in the config.site file.
I next configured using ./configure, with no
My code is as follows:
#include iostream
#include cmath
using namespace std;
#define MATHLIB_STANDALONE 1
extern C
{
#include R_ext/Applic.h
}
typedef struct TT{
double ** tempX;
double * tempY;
int tempN;
} TT, *MM;
double fn(int N, double * beta, void * ex){
Hello, I am trying to install R on Suse linux 10.1. To get configure to
run correctly, I had to set FPICFLAGS=fPIC.
R is now configured for i686-pc-linux-gnu
Source directory: .
Installation directory:/usr/local
C compiler:gcc -std=gnu99 -g -O2
Fortran 77
Check out panel.identify.cloud in the latest version of the lattice package.
E.g.
cloud(depth ~ lat * long, quakes)
trellis.focus()
panel.identify.cloud()
2008/10/31 Charles Liard [EMAIL PROTECTED]:
Hi,
Please tell me of an R graphical routine that will create a 3d plot
something like
Indeed in the case of a lazy-loaded package it looks a bit tricky, but
it should still be possible, as I guess one can access the source code
of any function defined in a package, once loaded (am I wrong?).
When no lazy-load takes place then it is actually quite straight forward
as you said
On Thu, 30 Oct 2008, William Hudspeth wrote:
Hello, I am trying to install R on Suse linux 10.1. To get configure to
run correctly, I had to set FPICFLAGS=fPIC.
FPICFLAGS=-fPIC would have been correct (or -fpic). It is odd that
FFLAGS has been set to -g and not '-g -O2': did you change
Martin Maechler maechler at stat.math.ethz.ch writes:
Yes, indeed, PDFs *are* fine
Apologies to the original poster, I was not aware of this.
Dieter
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hi the list:
I can not create R process from Java application with my linux box,
the process never was created, or just terminated immediately after I
called Runtime.getRuntime().exec(...)
The command lines I have tried are(the paths are all correct):
/bin/sh /usr/bin/R
/bin/sh -c
tsn4867 [EMAIL PROTECTED] writes:
For the package pls, I need to understand the algorithm for simpls.fit
for Partial Least Squares. I'm not sure if simpls.fit tries to find the
weight vectors (loadings) to maximize which of the two: Cov(Xw, y) or
maximize Cov^2(Xw,y)? Are these objective
Hi ,
Dataframe Df contains more than 30 no of records.but I want only
first 30 rows of data.
How can I retrieve it.
Thanks
K.Ravichandra
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
ML == Mark Leeds [EMAIL PROTECTED]
on Wed, 29 Oct 2008 17:13:04 -0500 (CDT) writes:
ML Hi Phil: That's EXACTLY what it is. Thanks so much. It's nice to know
ML that the R Gods don't hate me. I hope it's okay that I'm going to cc
ML r-help
ML in case this thread comes up in
Df[1:30,]
--- On Fri, 31/10/08, Kurapati, Ravichandra (Ravichandra) [EMAIL PROTECTED]
wrote:
From: Kurapati, Ravichandra (Ravichandra) [EMAIL PROTECTED]
Subject: [R] how to retrieve some records from a dataframe
To: r-help@r-project.org
Date: Friday, 31 October, 2008, 2:14 PM
Hi ,
megh wrote:
Is there anything like goto loop, which exists in most computer programs?
e.g. I am looking for this kind of stuff :
if(i 6) goto step-02
Any idea?
Regards,
It doesn't exist, but it can always be replaced by if() {} else {}
constructs. (You don't usually see goto in the
Hi R,
m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5))
I want to know the methods of getting row and column positions of NA in
the above dataframe. How do I do this?
Thanks, Shubha
This e-mail may contain confidential and/or privileged i...{{dropped:13}}
Am Freitag 31 Oktober 2008 12:17:30 schrieb Shubha Vishwanath Karanth:
m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5))
? which
HTH Claudia
--
Claudia Beleites
Dipartimento dei Materiali e delle Risorse Naturali
Università degli Studi di Trieste
Via Alfonso Valerio 6/a
I-34127 Trieste
Dear friends,
there is a serious problem with RODBC on R from 2.6.2 to 2.8.0 on a postgresql
8.3 database.
Let's say we define an empty dsn called test, with user postgres and
password none. We test the connection with another application, and it
works properly.
When you call the database
Is there anything like goto loop, which exists in most computer programs?
e.g. I am looking for this kind of stuff :
if(i 6) goto step-02
Any idea?
Regards,
--
View this message in context:
http://www.nabble.com/%22Goto%22-loop-tp20263445p20263445.html
Sent from the R help mailing list
2008/10/31 megh [EMAIL PROTECTED]:
Is there anything like goto loop, which exists in most computer programs?
Really? Not since 1968, I think:
http://www.cs.utexas.edu/users/EWD/ewd02xx/EWD215.PDF
e.g. I am looking for this kind of stuff :
if(i 6) goto step-02
Any idea?
Thinking you
See head:
head(Df, n = 30)
On Fri, Oct 31, 2008 at 6:44 AM, Kurapati, Ravichandra (Ravichandra)
[EMAIL PROTECTED] wrote:
Hi ,
Dataframe Df contains more than 30 no of records.but I want only
first 30 rows of data.
How can I retrieve it.
Thanks
K.Ravichandra
Try:
which(is.na(m), arr = T)
On Fri, Oct 31, 2008 at 9:17 AM, Shubha Vishwanath Karanth
[EMAIL PROTECTED] wrote:
Hi R,
m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5))
I want to know the methods of getting row and column positions of NA in
the above dataframe. How do I do
Thanks for your rapid response. I tried your sugerences and set
C_INCLUDE_PATH=/usr/local/include. I also look in /usr/local/include
and here it appears fftw.h and also fftw3.f, fftw3.h, fftw_threads.h
and rfftw.h rfftw_threads.h
but it don't work also. I will try to update to R 2.8.0
Dear all,
I would like to create a dotplot (lattice) with a log scale as the function
http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625
mydotchart of Johannes Huesing. Is it possible?
Thanx Gianandrea
--
View this message in context:
Use the survreg function.
There are many different ways to parameterize a Weibull. The survreg function
imbeds it a general location-scale familiy, which is a different
parameterization than the rweibull function.
y - rweibull(1000, shape=2, scale=5)
survreg(Surv(y)~1, dist=weibull)
I have been using some old R scripts that were prepared for me using
R Version 1.9.0. In these, there is a call to
library(locfit)
and with this version of R, I have no problem.
But needing to upgrade to the latest version of R, I find these
scripts no longer work, and with
help(locfit)
I
My problem is that I have already a distance (similarity) matrix generated
outside R through a C++ code because the criteria to calculate the
distance between pairs of objects are none of the standard criteria
implemented in R.
If I got it right, but I might be mistaken, stepFlexclust() performs
Hi
I was wondering if anyone have done an area analysis using chi-aquare?
I have three areas (area1, area2 and area3) I need to compare them to each
other by looking at a binomial distribution factor (0=negative, 1=positive).
I have already done a logistiv regression analysis comparing these
Hi,
I'm trying to set up a prediction software, now i testing the performance
of my method, so i need to calculate a ROC curve, specially auc, cut-off,
sens and spec, i just looking at ROCH package, but it's a mass for me, i'm
not a math guy and I'm getting lost
Could any of you recommend me
On Oct 30, 2008, at 11:38 PM, FD wrote:
Hi All:
Could anyone point me to a package that can calculate stratified
kappa? My
design is like this, 4 raters, 30 types of diagnosis scores, 20
patients.
Each rater will rate each patient for each type of diagnosis score.
The
rater's value is
My guess is that you have not configured you new installation to see
(and update) your collection of packages. As a quick fix, what happens
when you go to your favorite CRAN mirror and install a new copy of
locfit? My version is 1.5-4 and runs with recent versions of R.
--
David Winsemius,
I think
subset(set, InfCode %in% x)
does it, where 'x' is your vector of ids (e.g., c(1,3,9) in your example
at the bottom)
Brar Piening wrote:
I've got a big data.frame from which I need to extract data based on a variable
number of Id's (InfCode).
Until now I've been using the following
Hi,
I'm trying to get the p-value from the 'lm' regression function as a list
object. For example, I can get r^2 from the following code by entering
summary(fm)$r.squared. Is there a way to get the p-value? If not, is there
a function where I can enter the f-value and degrees of freedom to get
Try library(ROCR)
Pau Marc Munoz Torres wrote:
Hi,
I'm trying to set up a prediction software, now i testing the performance
of my method, so i need to calculate a ROC curve, specially auc, cut-off,
sens and spec, i just looking at ROCH package, but it's a mass for me, i'm
not a math guy
eric lee wrote:
Hi,
I'm trying to get the p-value from the 'lm' regression function as a list
object. For example, I can get r^2 from the following code by entering
summary(fm)$r.squared. Is there a way to get the p-value? If not, is there
a function where I can enter the f-value and degrees
on 10/31/2008 10:59 AM eric lee wrote:
Hi,
I'm trying to get the p-value from the 'lm' regression function as a list
object. For example, I can get r^2 from the following code by entering
summary(fm)$r.squared. Is there a way to get the p-value? If not, is there
a function where I can
Hi friends,
I need suggestions/directions on how to producing a waterfall plot for present
extend of change in tumour size for a set of respondents in a study. Example
of use of waterfall plot is in the following slides presented at ASCO 2007 by
Axel Grothey. Link is
On Fri, 31 Oct 2008, Gustave Lefou wrote:
Hello,
I would like to change my R language from German to French or English.
On Windows, we presume.
I have read messages in the archives and also R Installation and
Administration paragraph 7.2 but I did not get it all.
I have to set LANGUAGE=en
Philip Twumasi-Ankrah schrieb:
Hi friends,
I need suggestions/directions on how to producing a waterfall plot for present extend of change in tumour size for a set of respondents in a study. Example of use of waterfall plot is in the following slides presented at ASCO 2007 by Axel Grothey. Link
Good morning,
I am using R to try to model the proportion of burned area in Portugal.
The dependent variable is the proportion (P). The family used is
binomial and the epsilon would be binary. They have recommended me to
use the Betareg package as it allows to work with proportions.
When I use
Quantile Regression for Longitudinal Data.
Hi,
I am trying to estimate a quantile regression using panel data. I am trying
to use the model that is described in Dr. Koenker's article. So I use the
code the that is posted in the following link:
Hi,
I am studying Kalman Filter and it seems to be difficult for me to apply the
filter on a simple ARMA.
It is easy to construct the state-space model, for instance:
dlmModARMA(ar=c(0.4,-0.2),ma=c(0.2,-0.1, sigma2=1)
but applying the dlmFilter on it, it doesn't work...
I don't know if my
2008/10/31 Peter Dalgaard [EMAIL PROTECTED]:
It doesn't exist, but it can always be replaced by if() {} else {}
constructs. (You don't usually see goto in the class of functional
programming languages to which R belongs. See also
http://en.wikipedia.org/wiki/Edsger_W._Dijkstra .)
Also see
A - matrix(-4:4, ncol=3)
A
[,1] [,2] [,3]
[1,] -4 -12
[2,] -303
[3,] -214
Apos - A; Apos[A = 0] - NA;
Apos
[,1] [,2] [,3]
[1,] NA NA2
[2,] NA NA3
[3,] NA14
Aneg - A; Aneg[A = 0] - NA
Aneg
[,1] [,2] [,3]
[1,] -4 -1 NA
Hello -
Diogo André Alagador wrote:
Dear all,
I have a matrix with positive and negative values.
From this I would like to produce 2 matrices:
1st - retaining positives and putting NA in other positions
2nd - retaining negatives and putting NA in other positions
and then apply rowMeans
The ifelse function sometimes has unwanted side effects. Your problem can
easily be done with simple subsetting. Try:
pos - neg - A
pos[ pos = 0 ] - NA
neg[ neg = 0 ] - NA
# remove the ='s if you don't want strict positive/negative
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Dear list!
I have encountered an odd problem with R-WinEdt.
The program is mostly running fine and I had no problems with the installation.
But several times a day I'm typing in R-WinEdt there appears a menu with
Chinese characters where the curser is and I'm not able to type normally. I
dont
Thanks to Mark Leeds for an implementation of exponentially weighted
moving averages as follows which solves this problem. This also resmooths
recursively.
ewma-function(x,lambda = .5, init = (1-lambda)*.raw[good.ind][1],order=1)
{
.raw - unclass(coredata(x))
good.ind - !is.na(.raw) #
On 10/31/08, glaporta [EMAIL PROTECTED] wrote:
Dear all,
I would like to create a dotplot (lattice) with a log scale as the function
http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625
mydotchart of Johannes Huesing. Is it possible?
The usual lattice syntax for this
on 10/31/2008 09:59 AM Diogo André Alagador wrote:
Dear all,
I have a matrix with positive and negative values.
From this I would like to produce 2 matrices:
1st - retaining positives and putting NA in other positions
2nd - retaining negatives and putting NA in other positions
and then
Hi,
I'm using the dotchart2-function from the Hmisc-package to plot multiple dots
on horizontal lines using the add=TRUE option.
My question is very simple: How can I change the color of the dots?
I couldn't find the right option in the help-pages of the dotchart2-function or
the
RSiteSearch(ewma,rest=func)
would have told you there already exists an ewma implementation in the qcc
package.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of [EMAIL PROTECTED]
Sent: Friday, October 31, 2008 10:04 AM
To: [EMAIL
I've been trying to set my email setting to Digest Mode for a while, but
it doesn't seem to take.
I log through to the Unsubscribe or edit options button on the Web
interface to the R-help mailing list page. I scroll down to the Set
Digest Mode box and change the radio button to On. I
Hi
I am a new user the R and I am very grateful for all your help but...
I have a problem and I can't resolve yet. I am trying to get the contrasts for
a binomial 3-way GLM (T= 4 temperature, t= 2 time and c= 2 substrate levels,
plus treatment control) in total they are 17 treatments.
If you are going to insist on doing such things you will have to learn
to read the documentation. In this case if you do a
traceback()
you will see that the error is occurring in rq.fit.slm and when you do
?rq.fit.slm
you will see that there are several storage sizes that can
Hadley et al.,
I was using the cast function to reshape some data (aggregate a melted
data frame) and I did not put in the fill and for the most part the
values that came out were fine, but there were value great than an
order of magnitude from the actual value. When I put in the fill
argument
I'm currently using R version 2.7.2. Not the newest, but not even six months
old either. Thanks for the idea.
-Original Message-
From: Erin Hodgess [mailto:[EMAIL PROTECTED]
Sent: Thursday, October 30, 2008 4:06 PM
To: Cahill, Joshua
Cc: r-help@r-project.org
Subject: Re: [R]
Saskia Sandring wrote:
Dear list!
I have encountered an odd problem with R-WinEdt.
The program is mostly running fine and I had no problems with the installation.
But several times a day I'm typing in R-WinEdt there appears a menu with
Chinese characters where the curser is and I'm not able
On Friday, 31 October 2008, 10:46 (UTC-0700), Farley, Robert wrote:
I've been trying to set my email setting to Digest Mode for a while, but
it doesn't seem to take.
I think I had such a problem at some point, and I think it was because
I didn't have cookies enabled in my browser.
m.
--
Thank you very much!
Brar
Erik Iverson schrieb:
I think
subset(set, InfCode %in% x)
does it, where 'x' is your vector of ids (e.g., c(1,3,9) in your
example at the bottom)
Brar Piening wrote:
I've got a big data.frame from which I need to extract data based on
a variable number of Id's
Hello, I have a problem with using replace() to convert a vector of dates
from -mm-dd to julian date.
For example, I type
replace(x,2004-05-14,134)
and I receive an error:
Error in `[-.data.frame`(`*tmp*`, list, value = 134) :
new columns would leave holes after
on 10/31/2008 01:07 PM [EMAIL PROTECTED] wrote:
Dear fellows,
I'm trying to extract the AIC statistic from a GLM model with quasipoisson
link.
The formula I'm referring to is
AIC = -2(maximum loglik) + 2df * phi
with phi the overdispersion parameter, as reported in:
Peng et
There are various packages that do ROC curves (others have mentioned some, you
can search for the others), but ROC curves are not that hard to do by hand in
R. Using the classic iris dataset, we can use Sepal.Length as the test/score
and species=='virginica' as the true pos/neg that we are
Perfectly... thanx Deepayan
Deepayan Sarkar wrote:
On 10/31/08, glaporta [EMAIL PROTECTED] wrote:
Dear all,
I would like to create a dotplot (lattice) with a log scale as the
function
http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625
mydotchart of Johannes
Am 30.10.2008 um 16:53 schrieb [EMAIL PROTECTED]:
Hallo! ?xml:namespace prefix = o ns =
urn:schemas-microsoft-com:office:office /o:p/o:p
[...]
[[alternative HTML version deleted]]
This is a rare case where I would have liked to see the HTML version in
order to see how
Marc Schwartz marc_schwartz at comcast.net writes:
on 10/31/2008 01:07 PM Antonio.Gasparrini at lshtm.ac.uk wrote:
I'm trying to extract the AIC statistic from a GLM model
with quasipoisson link.
The formula I'm referring to is
AIC = -2(maximum loglik) + 2df * phi
with phi
Hello,
I'd like to think all respondents for your advises and suggestions. I am sorry
for not including the ordered gene expression data frame (largedf) since it is
large.
The problem is solved now by Greg Down by including the command:
rm(.Random.seed) in the code, perhaps I used some R
Ben, I think the reference you're searching for is the one below
@ARTICLE{Lebreton1992,
author = {Lebreton, J.-D. and Burnham, K. P. and Clobert, J. and Anderson,
D. R.},
title = {Modeling survival and testing biological hypotheses using marked
animals: a unified approach with
** Sorry to repost. I forgot to include a function necessary to make
the example work **
I apologize up front for this being a little long. I hope it's
understandable. Please let me know if I need to clarify anything.
Several months ago I wrote a series of functions to help me take my R
Hi,
Similar to the R function of converting rpart into rules (see
example:
http://books.huihoo.org/data-mining-desktop-survival-guide/Convert_Tree.
html_
is there any open code to
convert ctree (from party) into rule format?
This would be very useful.
thanks
Dhruv
Dhruv
Pau Marc Munoz Torres wrote:
Hi,
I'm trying to set up a prediction software, now i testing the performance
of my method, so i need to calculate a ROC curve, specially auc, cut-off,
sens and spec, i just looking at ROCH package, but it's a mass for me, i'm
not a math guy and I'm getting lost
Maura,
No, in general it is not possible. Depending on your
goodness-of-clustering metric (and there are many besides silhouette
width), the results may demonstrate multiple peaks and be severely
non-monotonic. Perhaps more problematic, silhouette width is especially
difficult. It is
No Joy. :-(
I've tried with both IE (7.0) and Firefox (3.0.3) and the same thing
happens. Firefox has cookies enabled, and I set IE to explicitly allow
http://www.r-project.org/ cookies.
Robert Farley
Metro
www.Metro.net
-Original Message-
From: [EMAIL PROTECTED]
I want to import some target values to normalize elements of a
dataframe. I'm summarizing the existing totals using
as.matrix(xtabs(~tripid_nu+lineon, data=SurveyData)) and reading a .csv
table of target values. I can't even get them to list out as similar
data items.
My idea is to do
Dear all,
I am in need to estimate the amount of variation
explained in a variable after simulations that produce a predictor
which is in the same units as the dependent variable (numbers of
species). Since the dependent and predictor variables are the same, I would
think the most appropriate
Hi Kirsten,
The second argument of replace() is actually an index, so it's trying to
replace element number 2004-5-14 = 1985 with 134! See ?replace. I would
need more of your code to understand your error exactly.
I think you would be better off converting your dates to date format (?Date)
and
hello,
what i'm trying to achieve is shown in the comment (# ...) below, which is
surely very wrong, but probably intuitively more meaningful that i can
provide in words.
xyplot(y~x2|x1+which,data=make.groups(
a=data.a,
b=data.b,
c=data.c
),
panel=function(...){
panel.abline(a=0,b=0)
In the 'etc' directory of your installation path, e.g. C:\Program Files\R\etc
Regards,
Yihui
--
Yihui Xie [EMAIL PROTECTED]
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
School of Statistics, Room 1037, Mingde Main Building,
Renmin
On 10/31/08, erwann rogard [EMAIL PROTECTED] wrote:
hello,
what i'm trying to achieve is shown in the comment (# ...) below, which is
surely very wrong, but probably intuitively more meaningful that i can
provide in words.
xyplot(y~x2|x1+which,data=make.groups(
a=data.a,
b=data.b,
Try:
?locator
On Fri, Oct 31, 2008 at 10:06 PM, German Urrus [EMAIL PROTECTED] wrote:
Â
Hello, my name is German, I am from Argentina, and we are working atÂ
university with R. This work is about sounds and how they affect out ears,
we are using R to analice different recordings of
Hi there,
I have a data set:
a=cbind(5,2,4,7,8,3,4,11,1,20)
I want to count # of data, satistfying a[1]a[2:10].
Anyone helps me solving this case?
Thank you in advance,
Jin
--
View this message in context:
http://www.nabble.com/count-data-with-some-conditions-tp20275722p20275722.html
82 matches
Mail list logo