Just a remark:
> I need some help using read.ftable to read a contingency table. My columns
> are organized as follows:
> order--family--species--location--number of individuals
As the first three variables are nested, I would expect this table
to contain a lot of structural zeroes. I understand
Hello,
par(las=1) sets the orientation of the axis labels
to horizontal. That is, the tick mark labels. How
do I set the orientation of the axis label, which
annotates the variable plotted along the axis, to
horizontal?
Sorry for asking such a basic question here, but I
haven't found anything in t
the usage of MNP is described in "MNP: R Package for Fitting the
Multinomial Probit Model"
http://www.jstatsoft.org/counter.php?id=128&url=v14/i03/v14i03.pdf&ct=1
If the Dependent Variables is Unordered ,why not use Multinomial
Logistic Regression.see
http://gking.harvard.edu/zelig/docs/_TT_mlogit
I use mnp to run a multinominal probit regression model, but the summary
doesn't contain the model statistics, such as the log-likelihood and degree
of freedom, for the assessment of the goodness-of-fit of the fitted model.
Is there any way that I can generate these statistics for the fitted mod
Dear R-helpers,
While plotting a hclust tree, is it possible to mark group information on it?
Suppose I have a hclust tree "ClusTree", and I use cutree to cut the
tree into different groups as following
a=cutree(ClusTree,h=10)
How do I mark the grouping info stored in a while I plot "ClusTree"
André Anderson Carvalho wrote:
> Hi R Users...
>
> I'm trying to get the Odds Ratios for Logistic Regression Coefficients and
> their confidence intervals when there are interactions, or effect modifiers,
> in the regression model. May anyone help me how to get them?
>
> Thanks
>
> André Ander
I think it can.But if you provide more information,you will be more help.
for example,you had better give a reproducable example in you email.
2006/1/11, Naiara S. Pinto <[EMAIL PROTECTED]>:
> Hi all,
>
> I need some help using read.ftable to read a contingency table. My columns
> are organized as
The correlation between the predictions from your two model fits is
0.95. This suggests to me that the differences between the two sets of
answers have little practical importance, and anyone who disagrees may
be trying to read more from the results than can actually be supported
by
On Tue, 10 Jan 2006 20:25:23 -0500, r user <[EMAIL PROTECTED]> wrote:
> I have 2 dataframes, each with 5 columns and 20 rows.
> They are called data1 and data2.I wish to create a
> third dataframe called data3, also with 5 columns and
> 20 rows.
>
> I want data3 to contains the values in data1 whe
Hi R-help citizens,
I'm having trouble making version 0.99-6 of Matrix on FreeBSD 6.0.
The error message is:
* Installing *source* package 'Matrix' ...
** libs
gcc -I/usr/local/lib/R/include -I/usr/local/include -D__NO_MATH_INLINES -fPIC
-g -O2 -c Csparse.c -o Csparse.o
... numerous lines d
Dear list
I want to obtain the adjusted r-square given a set of coefficients (without the
intercept), and I don't know if there is a function that does it.
Exist
I know that if you make a linear regression, you enter the dataset and have in
"summary" the adjusted r-square. But
I just got 2 hits for 'RSiteSearch("significant peaks in
periodogram")', the first of which was
"http://finzi.psych.upenn.edu/R/Rhelp02a/archive/61482.html";. If you
would like more help from this group, PLEASE do read the posting guide!
"www.R-project.org/posting-guide.html". Anecd
I've done stuff like this before but
it's been a while and I'm stuck.
Suppose I have a matrix with one
column x and another column y
and both are numeric and let the
row index of the matrix be i
Starting at index i ( i would equal on the first iteration )
when the cumulative sum of x_i+1 - x_i
Take the following example:
a <- rnorm(100)
b <- trunc(3*runif(100))
g <- factor(trunc(4*runif(100)),labels=c('A','B','C','D'))
y <- rnorm(100) + a + (b+1) * (unclass(g)+2)
m <- lm(y~a+b*g)
summary(m)
Here b is discrete but not treated as a factor. I am interested in
comp
I have 2 dataframes, each with 5 columns and 20 rows.
They are called data1 and data2.I wish to create a
third dataframe called data3, also with 5 columns and
20 rows.
I want data3 to contains the values in data1 when the
value in data1 is not NA. Otherwise it should contain
the values in data2.
Hi, Matt,
See FAQ 7.22.
--sundar
Matt Sakals wrote:
> I am trying to use the levelplot function from a command file.
> Here is the code:
>
> library(sp)
> library(gstat)
> library(lattice)
>
> gatherData <- read.table("~/gather.txt", header = TRUE)
> grd = makegrid(gatherData$x, gatherData$y,
I think this is an FAQ. Anyway, levelplot() is a lattice plot and therefore
must be explicitly printed when source()ed or in a function; i.e.
print(levelplot())
is required.
-- Bert
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Matt Saka
Hi R Users...
I'm trying to get the Odds Ratios for Logistic Regression Coefficients and
their confidence intervals when there are interactions, or effect modifiers, in
the regression model. May anyone help me how to get them?
Thanks
André Anderson
Brasília, Brazil
___
I am trying to use the levelplot function from a command file.
Here is the code:
library(sp)
library(gstat)
library(lattice)
gatherData <- read.table("~/gather.txt", header = TRUE)
grd = makegrid(gatherData$x, gatherData$y, cell.size = 5)
k <- krige(z~x+y, ~x+y, data = gatherData, newdata = grd,
Hi all,
I need some help using read.ftable to read a contingency table. My columns
are organized as follows:
order--family--species--location--number of individuals
I couldn't figure out how to change the data on my text file to be
imported into R; and after you do that, is it possible to convert
Thanks to Leif Kirschenbaum, Brian Ripley and Petr Pikal,
who helped with this problem. Unfortunately it appears as if the problem
with columns beginning with NAs is deeply connected to Microsoft ODBC
DLL's automatic determination of column type based on the contents of
the first N rows of each
See 3.3.1 in the R FAQ (** not ** the Windows FAQ).
-- 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. Box
> -Original Message-
> From: [EMAIL PROTECTED]
> [mail
Leif:
Rather than trying to mimic what you might do in SAS take advantage of R's
ability to use arbitrary data structures, e.g. lists. So, one approach is:
(your.df is the data frame)
your.list<-split(your.df,your.ts[,1:3],drop=TRUE)
t(sapply(your.list,function(x)x[which.max(x$TS),]))
Cheers,
B
Hello,
Could somebody point me, is there any function in R which returns
expected values of order statistics for normal distribution? I have been
looking and couldn't find it.
Thanks!
Anna
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/m
Let's say I have a data.frame like
A B C TS other columns
1 1 1 12345
1 1 1 56789
1 2 1 23456
1 2 2 23457
2 4 7 23458
2 4 7 34567
2 4 7 45678
and I w
On 1/10/06, Antje Schüle <[EMAIL PROTECTED]> wrote:
>
> It is nearly the same example I wrote about in
> http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg54238.html. I'll
> print it out again:
>
> In the first column I have numbers from 0 to 23 (hours of a day), the second
> column contains
I reproduce from memory my exhaustive look into this issue.
RODBC uses the Microsoft ODBC DLL's developed by Microsoft.
These DLL's perform an automatic determination of column type based on the
contents of the first N rows of cells in each column, where N [0,16]. N may be
set in the Windows sy
This is a "heads up" to those using the lmer function for fitting
linear mixed models or generalized linear mixed models. In the next
few days Martin and I will upload a new version of the Matrix package
that employs a new representation of linear mixed models based on the
supernodal Cholesky fact
In
http://cran.r-project.org/contrib/extra/batchfiles/
there are Windows XP batch files that you can use that will
automatically locate R using the registry and run it so that you
don't have to change your path. Just place them anywhere
in your path. See the README for more info.
I use them fo
One other thing to try could be soundex. ITs normally used for
last names but it might work here too. Google to find the
soundex encoding rules. Reviewing the country names might
suggest minor modifications to the soundex algorithm to
improve it for your case.
On 1/10/06, Gabor Grothendieck <[E
There's a gotcha in using identical() to compare dimensions -- it also
compares names, e.g.:
> x <- array(1:14, dim=c(rows=3,cols=5))
> dim(x)
rows cols
35
> identical(dim(x)+0, c(3,5))
[1] FALSE
> identical(as.numeric(dim(x)+0), c(3,5))
[1] TRUE
>
Gabor Grothendieck wrote:
> If its
Hi,
I am trying to run a batch command on Windows XP and R CMD is not
recognized. I get the error,
"'R' is not recognized as an internal or external command, operable program
or batch file."
I have "C:\Program Files\R\rw2010\bin" in my $PATH variable and Rcmd.exe
has been installed in that fol
If its just succint you are after then this is slightly
shorter:
identical(dim(x)+0, c(3,5))
On 1/10/06, Gregory Jefferis <[EMAIL PROTECTED]> wrote:
> Thanks for suggestions. This is a simple question in principle, but there
> seem to be some wrinkles - I am always having to think quite care
You can improve it somewhat by first accepting all the largest
matches and removing the rows and columns for those and
repeatedly doing that with what is left.
On 1/10/06, Werner Wernersen <[EMAIL PROTECTED]> wrote:
> Thanks for the nice code, Gabor!
>
> Unfortunately, it seems not to work for my
Thanks to some help by Doug Bates (and the updated version of the Matrix
package), I've refined my question about fitting nested and non-nested
factors in lmer(). I can get it to work in linear regression but it
crashes in logistic regression. Here's my example:
# set up the predictors
n.age
I would throw a tolower() around s1 and s2 so that 'canada' matches with
'CANADA', and perhaps consider using a Levenshtein distance rather than
the longest common subsequence.
An algorithm for Levenshtein distance can be found here (courtesy of
Stephen Upton)
https://stat.ethz.ch/pipermail/r-help
See
http://en.wikipedia.org/wiki/Levenshtein_distance
http://thread.gmane.org/gmane.comp.lang.r.general/31499
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Werner
> Wernersen
> Sent: Tuesday, January 10, 2006 2:00 PM
> To: Gabor Grothendieck
> Cc:
To r-help:
Thanks to Mat Soukup for the help:
To draw two ablines:
xYplot(lo ~ vaR, groups = v, data = abc.fp, aspect = "xy", col = c
("red","blue"),
xlab=expression(frac(abs( bold(v) ),abs( bold(a) ))),
ylab = grid::textGrob(expression(paste(log, frac( italic(p) ( italic
(v) )
Thanks for the nice code, Gabor!
Unfortunately, it seems not to work for my purpose, confuses lots of
countries when I compare two lists of over 150 countries each.
Do you have any other suggestions?
Gabor Grothendieck <[EMAIL PROTECTED]> schrieb: If they were the same you
could
Thanks for suggestions. This is a simple question in principle, but there
seem to be some wrinkles - I am always having to think quite carefully about
how to test for equality in R. I should also have said that I would like
the check to be efficient as well safe and succinct.
One suggestion was:
Hi,
I'm having trouble with glmmPQL from the MASS package.
I'm trying to fit a model with a binary response variable, two fixed
and two random variables (nested), with a sample of about 200,000
data points.
Unfortunately, I'm getting an error message that is difficult to
understand without
For errors with:
download.file()
install.package()
update.package()
any other command that R uses to access the internet
Setting up a proxy in R 2.2.1 for windows.
I work at a Central Bank so our firewall and proxy are very particular. To
allow R to access the net I did the following:
1. In or
Amir,
if you try the following code:
require(wavelslim)
data(doppler)
dwt(doppler, wf="la6")
you get:
Error in switch(name, haar = select.haar(), d4 = select.d4(), mb4 =
select.mb4(), : Invalid selection for wave.filter
So, it seems that instead of "la8" (default) and "haar" you can choose
Hello,
I wrote a version of which.na (similar to S's) in R that returns a
vector of positions at which NA appears in the target vector v:
which.na<-function(v) {
retv <- c()
for (i in 1:length(v)) {
if (is.na(v[i])) {
retv<-append(i, retv, after = length(retv))
}
}
> (2) Is there a more simpler way to build the vector of positions (versus
> using append in a loop, as above), using functions such as match() and
> apply()?
How about
which.na <- function(x) which(is.na(x))
?
Hadley
__
R-help@stat.math.ethz.ch mailin
On 1/10/06 11:41 AM, "Seth Falcon" <[EMAIL PROTECTED]> wrote:
> On 10 Jan 2006, [EMAIL PROTECTED] wrote:
>
>> Hi everyone!!
>>
>> This is my first message to the list, so I hope not to disturb
>> anyone if the subject of my message has been already treated.
>>
>> The question is that I have
R is not thread safe, so you must not use it in a
re-entrant way.
If you want to exploit multiple processors, you can write code (e.g.,
in C) called from R (e.g., through .Call or .C) that performs
parallel/threaded computations in a thread-safe way (e.g., without
calling back into R).
Another po
> Hi everyone!!
>
> This is my first message to the list, so I hope not to disturb
> anyone if the subject of my message has been already treated.
>
> The question is that I have a tool, a GUI made with Java, connected
> to R using Rserve, and I'd like to get R and Rserve in a
> multi-processor m
You can also get access to the code chunks in vignettes
as shown here:
http://tolstoy.newcastle.edu.au/~rking/R/help/05/12/17822.html
On 1/9/06, Uwe Ligges <[EMAIL PROTECTED]> wrote:
> Mark Leeds wrote:
>
> > Sometimes I print out a package
> > and read about it and there
> > are sometimes nic
On 10 Jan 2006, [EMAIL PROTECTED] wrote:
> Hi everyone!!
>
> This is my first message to the list, so I hope not to disturb
> anyone if the subject of my message has been already treated.
>
> The question is that I have a tool, a GUI made with Java, connected
> to R using Rserve, and I'd like to g
Using R, I believe only the Haar and Daubechies wavelets are available.
wavethresh does allow for you to use different filter numbers with the
Daubechies wavelets that I think you cannot do using waveslim. You can do
cross-validation to find the best fit across different filter niumbers and
t
> "RobMcG" == McGehee, Robert <[EMAIL PROTECTED]>
> on Tue, 10 Jan 2006 10:55:47 -0500 writes:
RobMcG> Hello,
RobMcG> Why not just copy and paste the examples from the help pages? No
need to
RobMcG> type anything.
RobMcG> However, if you'd like to run the entire help
If they were the same you could use merge. To figure out
the correspondence automatically or semiautomatically, try this:
x <- c("Canada", "US", "Mexico")
y <- c("Kanada", "United States", "Mehico")
result <- outer(x, y, function(x,y) mapply(lcs2, x, y))
result[] <- sapply(result, nchar)
# try b
Andrew Gelman <[EMAIL PROTECTED]> writes:
> Error in "rownames<-"(x, value) : attempt to set rownames on object with
> no dimensions
> Am I doing something stupid here or is it actually a (minor) bug?
The latter. There's a drop=FALSE missing in the coef() method,
specifically in the line
Hi people,
I want to obtain the adjusted r-square given a set of coefficients (without
the intercept), and I don't know if there is a function that does it.
Exist
I know that if you make a linear regression, you enter the dataset and have in
"summary" the adjusted r-square. But
Hello,
Why not just copy and paste the examples from the help pages? No need to
type anything.
However, if you'd like to run the entire help section for a function,
check out the example function, i.e.:
> example(plot)
Robert
-Original Message-
From: Mark Leeds [mailto:[EMAIL PROTECTED]
[EMAIL PROTECTED] wrote:
> Dear ladies and gentlemen!
> When I use the plot funtion how can I change the size of the title for the x
> and
> y axes (xlab, ylab)and the size of the axes label ?
>
> Thank you very much.
>
> With best regards
>
> Claudia
>
> _
At 16:06 09/01/06, Gabor Grothendieck wrote:
[snip various earlier posts]
>In addition to books, the various manuals, contributed documents and
>mailing list archives, all of which one should review,
>the key thing to do if you want to really learn R is to read source code
>and lots of it. I th
I'm having another problem with lmer(), this time something simpler (I
think) involving the coef() function for a model with varying
coefficients. Here's the R code. It's a simple model with 2
observations per group and 10 groups:
# set up the predictors
n.groups <- 10
n.reps <- 2
n <- n.gro
Hi,
Before I reinvent the wheel I wanted to kindly ask you for your opinion if
there is a simple way to do it.
I want to merge a larger number of tables from different data sources in R
and the matching criterium are country names. The tables are of different size
and sometimes the c
Hi everyone!!
This is my first message to the list, so I hope not to disturb anyone if the
subject of my message has been already treated.
The question is that I have a tool, a GUI made with Java, connected to R
using Rserve, and I'd like to get R and Rserve in a multi-processor machine.
Now, wh
Dear r-helpers,
Happy New Year.
To a plot
xYplot(lo ~ vaR, groups = v, data = abc1.fp,aspect = "xy",
xlab=expression(frac(abs( bold(v) ),abs( bold(a) ))),
ylab = grid::textGrob(expression(paste(log, frac( italic(p) ( italic
(v) ), italic(p) ( italic(a) ),
abline =li
Dear ladies and gentlemen!
When I use the plot funtion how can I change the size of the title for the x and
y axes (xlab, ylab)and the size of the axes label ?
Thank you very much.
With best regards
Claudia
__
R-help@stat.math.ethz.ch mailing list
htt
> "Gregory" == Gregory Jefferis <[EMAIL PROTECTED]>
> on Tue, 10 Jan 2006 14:47:43 + writes:
Gregory> Dear R Users,
Gregory> I want to test the dimensions of an incoming
Gregory> vector, matrix or array safely
Gregory> and succinctly. Specifically I want to ch
you could use: isTRUE(all.equal(dim(obj), c(3, 5)))
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: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: ht
Dear R Users,
I want to test the dimensions of an incoming vector, matrix or array safely
and succinctly. Specifically I want to check if the unknown object has
exactly 2 dimensions with a specified number of rows and columns.
I thought that the following would work:
> obj=matrix(1,nrow=3,ncol=
Dear list,
I would like to perform an analysis on the following model:
aov(ampratio ~ Type * Place * agemF + Error(speakerid/Place) ,data=aspvotwork)
using the approach from http://www.psych.upenn.edu/~baron/rpsych/rpsych.html .
Now, I got the test results, wich indicate a significant interacti
Another year another project.
Does anyone have any experience / pointers in obtaining data from an
analog / digital converter hardware such as a National Instruments data
aquisition card using R?
I really need to obtain and process this data in a real time enviroment
and I'm hoping that R already h
We've had a small "review time" within R-core on this topic,
amd would like to state the following:
--
The R-core team welcomes proposals to develop an R-wiki.
- We would consider linking a very small number of Wikis (ideally
Dear Petr,
Thank you for your help. I have tried (and succeded) to import myfile
after conversion to text and with the help of read.table (also with the
file = 'clipboard' alternative). Both methods give correct results,
albeit dateformat turns up as a factor (minor problem). Also the
read.xls
On 09-Jan-06 Mark Leeds wrote:
> Sometimes I print out a package
> and read about it and there
> are sometimes nice examples
> that I would like to run myself.
>
> Is there a way to bring them
> into R from the package or
> are they only meant to be typed
> in manually ? If manual is the
> only w
Hi all,
I think there is something I have forgotten to emphasize. I'd like to have an
xyplot (with lattice library). A trellis graphic. For every weekday I like to
have a graphic with all hours. And inside the graphic there are the two
frequencies as lines.
The problem finally is the use of
Another year another project.
Does anyone have any experience / pointers in obtaining data from an
analog / digital converter hardware such as a National Instruments data
aquisition card using R?
I really need to obtain and process this data in a real time enviroment
and I'm hoping that R already h
Soon on CRAN a new version of package gstat will emerge, which
has a few minor changes and possible incompatibilities w.r.t. the
previous version(s).
The new gstat (0.9-23) now:
+ depends on sp, and uses internally with Spatial* classes from sp
if data are provided in the old-fashoned way (as d
On 10-Jan-06 Ted Harding wrote:
>
> On 09-Jan-06 Mark Leeds wrote:
>> Sometimes I print out a package
>> and read about it and there
>> are sometimes nice examples
>> that I would like to run myself.
>>
>> Is there a way to bring them
>> into R from the package or
>> are they only meant to be ty
On 09-Jan-06 Mark Leeds wrote:
> Sometimes I print out a package
> and read about it and there
> are sometimes nice examples
> that I would like to run myself.
>
> Is there a way to bring them
> into R from the package or
> are they only meant to be typed
> in manually ? If manual is the
> only
Jean-Christophe BOUETTE wrote:
>>-- Message transféré --
>>From: Uwe Ligges <[EMAIL PROTECTED]>
>>To: Erin Hodgess <[EMAIL PROTECTED]>
>>Date: Sun, 08 Jan 2006 15:09:06 +0100
>>Subject: Re: [R] exporting methods/classes
>>Erin Hodgess wrote:
>>
>>
>>>Dear R People:
>>>
>>>I'm still
you can use something like:
Vcov <- vcov(lmer.fit, useScale = FALSE)
betas <- fixef(lmer.fit)
se <- sqrt(diag(Vcov))
zval <- betas / se
pval <- 2 * pnorm(abs(zval), lower.tail = FALSE)
###
cbind(betas, se, zval, pval)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. S
Dear R-Helpers,
I want to compare the results of outputs from glmmPQL and lmer analyses.
I could do this if I could extract the coefficients and standard errors
from the summaries of the lmer models. This is easy to do for the glmmPQL
summaries, using
> glmm.fit <- try(glmmPQL(score ~ x*type, ran
> -- Message transféré --
> From: Uwe Ligges <[EMAIL PROTECTED]>
> To: Erin Hodgess <[EMAIL PROTECTED]>
> Date: Sun, 08 Jan 2006 15:09:06 +0100
> Subject: Re: [R] exporting methods/classes
> Erin Hodgess wrote:
>
> > Dear R People:
> >
> > I'm still struggling with sending methods a
It is nearly the same example I wrote about in
http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg54238.html. I'll print
it out again:
In the first column I have numbers from 0 to 23 (hours of a day), the second
column contains the name of a weekday (Day as factor) and the third column
c
On Tue, 10 Jan 2006, Erin Hodgess wrote:
> Dear R People:
>
> In S Plus, if you have a function which calls the plot
> function several times, you get several "pages" of graphics
> output.
Only on a graphsheet, I believe.
> Is there an eqivalent in R, please?
Pretty close, on a windows() device
Erin Hodgess wrote:
> Dear R People:
>
> In S Plus, if you have a function which calls the plot
> function several times, you get several "pages" of graphics
> output.
>
> Is there an eqivalent in R, please?
Yes, for e.g. postscript() and pdf() devices the default is to plot each
new plot on a
Dear R People:
In S Plus, if you have a function which calls the plot
function several times, you get several "pages" of graphics
output.
Is there an eqivalent in R, please?
R version 2.2.1 windows
Thanks in advance!
Sincerely,
Erin Hodgess
Associate Professor
Department of Computer and Mathem
Hi
I tried to reproduce what you have told us by copy and paste
read.delim("clipboard")
but was not successful.
Even with several blank values in each column in Excel i got correct
import to R by this process. As I do not use RODBC I do not know all
possible settings and features. If colClass
stefan semmelring wrote:
> "stefan semmelring" <[EMAIL PROTECTED]> schrieb am 09.01.06 21:31:01:
>
>>Hallo,
>>
>>
>>
>>I´m trying to find out how I can start WinEdt always directly when R is
>>started.
>>
>>Does anybody know how to do so? I tried to find it out by myself but had no
>>fortune.
>
86 matches
Mail list logo