For the home page, click on Search, then Searchable mail archives, and
look
for tobit There's LOTS of stuff.
I can't found what you said.
Is it assumption like OLS ? :
1. Normal
2. No autocorrelation
3. Homogeneity of variance
4. No Multicolinearity
please make sure me !
I need answer
Hello - R-experts,
Is there any way with which we can specify the number after
decimal point to take. Like I have a situation where
the values are comming 0.160325923 but I only want
4 place to decimal say 0.1603. Is there any way for that.
I am no expert in R- and this may sound simple to
have a look at ?round()
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: http://www.med.kuleuven.be/biostat/
Hallo all
Please help me. I am lost and do not know what is the problem. I have
a factor called kvartaly.
attributes(kvartaly)
$levels
[1] 1Q.04 2Q.04 3Q.04 4Q.04 1Q.05 2Q.05 3Q.05 4Q.05
$class
[1] factor
mode(kvartaly)
[1] numeric
str(kvartaly)
Factor w/ 8 levels 1Q.04,2Q.04,..: 1 1 1 1 1
I am trying to write a lesson on maximum likelihood with Sweave. I get
a surprising result with the following code, lec4.Snw:
\documentclass[a4paper,12pt]{article}
\usepackage[latin1]{inputenc}
\title{Maximum likelihood}
\author{Göran Broström}
\begin{document}
\maketitle
fig=TRUE=
##
1) The obvious test is via is.factor(), and you have not used that.
2) Your example works for me, so what versions of R and chron is this?
3) Here's my guess. split is using the C-level test isFactor. That tests
that the factor is of type integer, so please try
typeof(kvartaly)
I suspect
Dear Colleagues,
does anybody know how to add dimnames to an image. Right now I'm using
image(as.matrix(df3), col=brewer.pal(9,Blues))
where df3 is a data.frame.
dimnames(as.matrix(df3)) delivers
[[1]]
[1] RFM_A1 RFM_A2 RFM_A4 RFM_A5 RFM_A7 RFM_B3 RFM_B6 RFM_B7
RFM_B8 RFM_B9
Thank you very much.
On 13 Feb 2006 at 10:54, Prof Brian Ripley wrote:
Date sent: Mon, 13 Feb 2006 10:54:21 + (GMT)
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Petr Pikal [EMAIL PROTECTED]
Copies to: r-help@stat.math.ethz.ch
On Mon, 13 Feb 2006 11:16:14 +0100,
Göran Broström (GB) wrote:
I am trying to write a lesson on maximum likelihood with Sweave. I get
a surprising result with the following code, lec4.Snw:
\documentclass[a4paper,12pt]{article}
\usepackage[latin1]{inputenc}
\title{Maximum
On Mon, 13 Feb 2006, Markus Preisetanz wrote:
Dear Colleagues,
does anybody know how to add dimnames to an image. Right now I'm using
I think you are looking for axis labels - they are taken by default from
your missing x and y arguments. Try something like:
Hello,
Please could someone advise if it's possible to save the graphical output
from the surface3d() function? I have tried the dev.copy() function to
save as a pdf but an error message says I cannot copy from the null
device.
Laura Quinn
Institute of Atmospheric Science
School of Earth and
On Mon, 13 Feb 2006, Petr Pikal wrote:
Thank you very much.
On 13 Feb 2006 at 10:54, Prof Brian Ripley wrote:
Date sent:Mon, 13 Feb 2006 10:54:21 + (GMT)
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Petr Pikal [EMAIL PROTECTED]
Copies to:
From: Göran Broström [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date: Mon, 13 Feb 2006 11:16:14 +0100
Subject: [R] Sweave, mle and curve
I am trying to write a lesson on maximum likelihood with Sweave. I get
a surprising result with the following code, lec4.Snw:
On 2/13/2006 6:56 AM, Laura Quinn wrote:
Hello,
Please could someone advise if it's possible to save the graphical output
from the surface3d() function? I have tried the dev.copy() function to
save as a pdf but an error message says I cannot copy from the null
device.
Are you talking about
Hi!
I have a problem of finding a specific value in a column. For example, I
have a matrix with say 2 columns
X Y
1-2.0341602 9.036689e-05
2-1.4287230 1.807338e-04
3-1.1194402 2.711007e-04
4-1.0327582 3.614676e-04
5-0.8130556
Hi
I am having difficulty making apply() work as expected (and desired)
with arrays
that have zero-extent dimensions.
?apply says
If each call to 'FUN' returns a vector of length 'n', then 'apply'
returns an array of dimension 'c(n, dim(X)[MARGIN])' if 'n 1'.
If 'n' equals
maybe something like this could be of help:
mat - matrix(rnorm(20), 10, 2)
val - -0.6523
###
ind - which.min(abs(mat[, 1] - val))
mat
mat[ind, ]
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Duncan Murdoch wrote:
Are you talking about surface3d in the rgl package? It doesn't use the
standard R graphics system, so dev.copy() won't help. You need
rgl.snapshot, which can create png files. You'll need some other
utility to convert those to pdf.
Or screen-grab them somehow.
Hi,
I'm not sure how much weight you're putting on efficiency for your
algorithm, but that sounds unnecessarily complicated. If your matrices
are not too big then maybe something like this would work. Let your
matrix be M and your value be a, then
yoursearch - function(M, a) {
#
In addition to round() mentioned earlier, if you are merely looking to
*display* your results differently, you may want to check out the
digits option, e.g. in summary():
(This is the method signature for data.frame 's):
summary(object, maxsum = 7,
digits = max(3,
Hi user
I have a question about export files.
I run my program for different countries and different years.
For export files I do this in SAS:
(1) Path to data files
%let teste=C:\Documents and Settings\MEALQ\My Documents;
(2) Country code;
%let cc=pt;
(3) Survey year;
%let
Hi
On 13 Feb 2006 at 12:04, Prof Brian Ripley wrote:
Date sent: Mon, 13 Feb 2006 12:04:49 + (GMT)
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Petr Pikal [EMAIL PROTECTED]
Copies to: r-help@stat.math.ethz.ch
Subject:
for write.table in R I use the following:
# table to be stored in excel for example
sum.stat - data.frame(sum.stat)# data to be stored
write.table(sum.stat, file=D://my_table.xls, sep=\t, dec=,,
row.names=TRUE, col.names=TRUE, qmethod=double)
There is a function in the Hmisc package that will help.
Example:
require(Hmisc)
x - sort(runif(10))
x
[1] 0.1225542 0.1620869 0.2197772 0.3187375 0.5498879 0.565
0.5812717 0.7380532 0.8187384 0.9063713
whichClosest(x,.41)
[1] 4
x[4]
[1] 0.3187375
Or for your matrix (yourmat),
thanks a lot for the help. I tried some stuff and now it works.
Ita Cirovic-Donev
Risk Controlling
HYPO ALPE-ADRIA-BANK D.D.
Amruseva 6, HR-1 Zagreb
Phone: +385 1 4899 269
Fax: +385 1 6063 344
E-Mail: [EMAIL PROTECTED]
Web: www.hypo-alpe-adria.hr
I have a large dataset that I am reading into R using read.table. The first
row contains the column names, and missing observations have the standard NA
value. When I try to run principal components with the default treatment of NA
values (which is na.omit) I get the error message
Error in
iw == ivo welch [EMAIL PROTECTED] writes:
iw just a short beg for the next R 2.3 version:
iw I know it is easy to add the sd into summary() in the source bowels of
iw R---but everytime R is updated, my change disappears. :-(. I do not
iw believe that R has an easy extension mechanism for
Evidently, my R has some feature of error messages disabled.
NG == Norman Goodacre [EMAIL PROTECTED] writes:
NG the following code apprantely, for some grand old reason, induces a syntax
error:
NG if (seq[i] = A) m - trans[1,]
NG Error: syntax error in if (seq[i] =
Where Norman gets this
Hi,
I am trying to install packages in R-2.2.1 on a Redhat WS4 system. I
get the following error messages trying to install, for example, the
akima package:
R.version.string
[1] R version 2.2.1, 2005-12-20
install.packages(akima, lib=/usr/lib/R/library, repos =
http://cran.fhcrc.org/;)
.
.
.
Ana,
the following code may help:
teste=P:/test
cc=pt
yy=03
data=matrix(1:12,3,4)
dir=paste(teste,cc,yy,sep=/)
# does directory already exist? if not, then create it
if (!is.null(try(setwd(dir),silent=TRUE))) {dir.create(dir,recursive=T)}
# write data to desired location
The eta^2 you describe looks something like an R^2 (or maybe a
partial R^2), and CohensD looks like a Student's t, at least to me. The
problem with generalizing these to multi-level models is deciding which
components of variance to include where. If you can answer that, I
think
Can survreg() handle interval-censored data like the documentation
says? I ask because the command:
survreg(Surv(start, stop, event) ~ 1, data = heart)
fails with the error message
Invalid survival type
yet the documentation for Surv() states:
Presently, the only methods
Patricia J. Hawkins wrote:
Evidently, my R has some feature of error messages disabled.
NG == Norman Goodacre [EMAIL PROTECTED] writes:
NG the following code apprantely, for some grand old reason, induces a
syntax error:
NG if (seq[i] = A) m - trans[1,]
NG Error: syntax error
How have you defined event?
library(survival)
?Surv
event: The status indicator, normally 0=alive, 1=dead. Other
choices are T/F (TRUE = death) or 1/2 (2=death). For interval
censored data, the status indicator is 0=right censored, 1=
event at 'time', 2=left
Berton Gunter wrote:
Folks:
A minor (R 2.2.0) Windows package install annoyance, but I can't figure out
how to fix it. If someone could tell me how or point me to the appropriate
docs, I would appreciate it.
I set my CRAN mirror in my Rprofile.site file. After startup, I get:
Hello,
I have data in a frame:
1 2 3 4 5 6 7 8
9 10 11 12 13 14 15 16
I would like them arranged in a single column:
1
2
3
4
.
.
8
9
.
.
16
etc.
I believe this should be possible using reshape, but I can't see how
to do it.
Thanks for sparing time to help a neophyte,
Roger Mason
You might just want to try this, it is quite efficient:
as.numeric(t(as.matrix(x)))
On 2/13/06, Roger Mason [EMAIL PROTECTED] wrote:
Hello,
I have data in a frame:
1 2 3 4 5 6 7 8
9 10 11 12 13 14 15 16
I would like them arranged in a single column:
1
2
3
4
.
.
8
9
.
.
16
[BioC] ANN: Introductory R and Bioconductor Mini Course , Seattle, April
12-14
We will present a short course on using Bioconductor, primarily aimed at
analyzing microarray data, over 2 1/2 days in April. The course will be
a mix of lectures and practical labs.
Instructors:
Robert
All,
I had a question about the JGR console and whether or not I can
manipulate the location where line wrapping occurs. I have searched
'JGR' in the R listserve archives and attempted to find console
manipulation on the JGR website to no avail and could use some
direction.
TIA, Bret
As an
Thank you Tariq. That works perfectly.
Roger
¨Tariq Khan [EMAIL PROTECTED] writes:
You might just want to try this, it is quite efficient:
as.numeric(t(as.matrix(x)))
__
R-help@stat.math.ethz.ch mailing list
Hi,
Is it possible to index via a character matrix? For example, I would
like to do the following:
cont.table - table(df1$A,df1$B) # df1 a data frame with factors A and B
cont.table[cbind(df2$A,df2$B)] # df2 a smaller data frame with some
# pairings of
Hi All.
I'm new to R and trying to learn the data manipulation routines. At the
moment, this is the one that has me stumped.
Imagine n test takers provide responses to k multiple choice items. The
result is a (n+1) x k matrix of responses where the first row of the matrix
is a correct answer
Hi All.
I'm new to R and trying to learn the data manipulation routines. At the
moment, this is the one that has me stumped.
Imagine n test takers provide responses to k multiple choice items. The
result is a (n+1) x k matrix of responses where the first row of the matrix
is a correct answer
Hello there,
In the e1071 package,'lca' function only works for binary data for
Latent class analysis. Am I right?
I would like to test the agreement amongst 6 raters for nominal data
ona scale from 1-4, and conduct a latent class analysis. What package
and what function may I use?
Thank you
Assuming this test data:
iris2 - round(iris[1:10, 1:4])
Try this:
t(t(iris2) == c(iris2[1,])) + 0
On 2/13/06, Rick DeShon [EMAIL PROTECTED] wrote:
Hi All.
I'm new to R and trying to learn the data manipulation routines. At the
moment, this is the one that has me stumped.
Imagine n test
Hello,
I am new to R, so forgive my ignorance on what is probably simple.
I find package simpleboot quite useful for LOESS bootstrapping and
generation of plots.
I want to calculate the standard error for x=60 from the 100 grid points
and 50 bootstraps.
The code below gives the first fitted
On Mon, 13 Feb 2006, Waichler, Scott R wrote:
I am trying to install packages in R-2.2.1 on a Redhat WS4 system. I
What architecture is this? How did you install R?
get the following error messages trying to install, for example, the
akima package:
R.version.string
[1] R version 2.2.1,
On Mon, 13 Feb 2006, Roger Levy wrote:
Hi,
Is it possible to index via a character matrix? For example, I would
like to do the following:
cont.table - table(df1$A,df1$B) # df1 a data frame with factors A and B
cont.table[cbind(df2$A,df2$B)] # df2 a smaller data frame with some
Dear R-users,
here is an error that you might encounter when using ROracle.
library(ROracle)
Loading required package: DBI
con - dbConnect(dbDriver(Oracle), mkienzle/[EMAIL PROTECTED])
dbListTables(con)
** Internal heap ERROR 17177 addr=0x0 *
Hey All,
I am trying to create a scatter plot of x,y,z data where the
points are color coded by their z value. So far I have:
plot(BGx, BGy, pch=22, type=p, col=heat.colors(20)) is this
right?
Thank you so much for your help!
-Seth
Lisa Wang wrote:
Hello there,
In the e1071 package,'lca' function only works for binary data for
Latent class analysis. Am I right?
I would like to test the agreement amongst 6 raters for nominal data
ona scale from 1-4, and conduct a latent class analysis. What package
package irr (on
On 2/13/2006 2:49 PM, Blanchard, Seth A. wrote:
Hey All,
I am trying to create a scatter plot of x,y,z data where the
points are color coded by their z value. So far I have:
plot(BGx, BGy, pch=22, type=p, col=heat.colors(20)) is this
right?
If the z values go from 1
Hi,
This is probably documented, but I cannot find the right words or
expression for a search. My attempts failed.
I have a data frame of 3 vectors (x, y and z) and would like to
transform this so that I could use persp. Presently I have y-level
copies of each x level, and a z value for
a strategy for this that I use is just
persp(interp(x,y,z))
where interp is from the Akima package, and x,y,z are all
of the same length.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL PROTECTED]Department of Economics
vox: 217-333-4558
I'm not sure what gave you the idea that there's a limit of 99 contour
lines, but no such limit is intended. You need to read the
documentation of panel.contourplot and use it in a panel function.
Hint:
library(lattice)
levelplot(volcano,
panel = function(..., at, contour = FALSE,
Some past threads have pointed out that lists can be used as
dictionaries storing data with an associated key (environments do this
too).
I have the seemingly simpler case with integer values for the
indices--however, the integers are not necessarily contiguous low
numbers. My concern is that if
On 2/10/06, simon chou [EMAIL PROTECTED] wrote:
Hi,
I was trying to stack a topographic map readed in from esri shapefile with a
contour map and a vector map. However, the plotMap(maptools) and
contourplot(lattice) do not seem to work well on top of each other. Here is
part of my code.
Hi.
I have found stepclass {klaR}, which is similar to what I was looking
for. However, I take note of your advice en stepwise variable selection.
Greetings.
Felipe
Alexandre Santos Aguiar wrote:
Em Dom 12 Fev 2006 16:10, Felipe Martínez-Pastor escreveu:
--apparently-- not variable selection.
The obvious answer was options(width=).
Bret
Bret Collier [EMAIL PROTECTED] 2/13/2006 12:24:22 PM
All,
I had a question about the JGR console and whether or not I can
manipulate the location where line wrapping occurs. I have searched
'JGR' in the R listserve archives and attempted to find
Hi, there.
Does anyone know the R function for calculating the cdf of bivariate
normal distribution function?
Thanks.
Yulei
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
Hello,
I am new to R, so forgive my ignorance on what is probably simple.
I find package simpleboot quite useful for LOESS bootstrapping and
generation of plots.
I want to calculate the standard error for x=60 of the 100 grid points
and 50 bootstraps.
The code below gives the first fitted value.
Check out the mvtnorm package.
Cheers,
Andy
__
Andy Jaworski
518-1-01
Process Laboratory
3M Corporate Research Laboratory
-
E-mail: [EMAIL PROTECTED]
Tel: (651) 733-6092
Fax: (651) 736-3122
I'm invoking R from withing a shell script like this
R --no-save --no-restore --gui=none `hostname` 21 BYE
# various commands here
BYE
I would like to regain control from the invoking terminal at some point.
I tried source(stdin()) but got a syntax error, presumably stdin is the
little shell
The only way I can see this happening is if you installed a
ix86 RPM on a
x86_86 system. You can do that and R will run, but you will
not be able to build packages out of the box.
This is exactly what happened. Now I know to use 'uname -a' to check
machine architecture first. I was able
Prof Brian Ripley wrote:
On Mon, 13 Feb 2006, Roger Levy wrote:
Hi,
Is it possible to index via a character matrix? For example, I would
like to do the following:
cont.table - table(df1$A,df1$B) # df1 a data frame with factors A and B
cont.table[cbind(df2$A,df2$B)] # df2 a smaller
Lisa Wang wrote:
I would like to test the agreement amongst 6 raters for nominal data on
a scale from 1-4, and conduct a latent class analysis in R. How should
the data be formatted and what code should I use?
This kind of agreement test seems to call for Krippendorf's alpha. See
# Dear R list,
#
# A have a doubt about SSQ decomposition and contrasts with ANOVAs.
# So, I would like a tip from more advanced R users.
# Below my data, the basic script and my doubts:
# Data
r = paste('r', gl(3, 8), sep='')
e = paste('e', rep(gl(2, 4), 3), sep='')
tra = sort(paste('t',
Hi R users:
I got this messages when I try to use the Rcmdr library, after
I make an update with update.packages() command.
I am using R221 in windows environment.
library(Rcmdr)
Loading required package: tcltk
Loading Tcl/Tk interface ... done
Loading required package: car
Error in
I have a file named test.csv with the following 3 lines:
%y-%m-%d;VALUE
1999-01-01;100
2000-12-31;999
read.table(test.csv, header = TRUE, sep = ;)
delivers:
X.y..m..d VALUE
1 1999-01-01 100
2 2000-12-31 999
I would like to see the following ...
%y-%m-%d VALUE
1 1999-01-01
?read.table
The documentation has the parameter 'check.names'.
On 2/13/06, Diethelm Wuertz [EMAIL PROTECTED] wrote:
I have a file named test.csv with the following 3 lines:
%y-%m-%d;VALUE
1999-01-01;100
2000-12-31;999
read.table(test.csv, header = TRUE, sep = ;)
delivers:
You can do it manually by reading in the headers separately:
headers - read.table(test, header = FALSE, nrow = 1, sep = ;, as.is = TRUE)
read.table(test, header = FALSE, skip = 1, sep = ;, col.names = headers)
On 2/13/06, Diethelm Wuertz [EMAIL PROTECTED] wrote:
I have a file named test.csv
Thanks a lot that works fine!
Next problem, if I would have my own package, and the file test.csv
would be located in the data directory
How to use the function data to get
data(test)
resulting in:
test
%y-%m-%d VALUE
1 1999-01-01 100
2 2000-12-31 999
Again Thanks in advance
With repeated measures, you will want to use the correlation
argument in glmmPQL. For help with how to do that, you need to know
that glmmPQL calls lme repeatedly. Therefore, if you have any
trouble using it, I suggest you consult Pinheiro and Bates (2000)
Mixed-Effects Models in S
This is my first time making a request so please feel free to school
me on how to improve on how to make an inquiry. I am running R-1.9.1
in UNIX via a PHP script with a system command to generate a *.jpg
image from a list created by a perl script. This R script runs fine
Given a nonlinear model formula and a set of values for all the
parameters defining a point in parameter space, is there a neat way to
extract the pseudodesign matrix of the model at the point? That is the
matrix of partial derivatives of the fitted values w.r.t. the parameters
evaluated at the
Dear Kenneth,
After just running update.packages() myself, I can't duplicate this error
(using R 2.2.1 under Win XP).
Can you provide a little more information? Are any but the standard packages
loaded when you issue the library(Rcmdr) command? Have you made any
modifications to a startup file,
From the help page of format, nsmall should control the number of digits.
format(0.123456789, nsmall = 10)
[1] 0.1234567890
format(0.123456789, nsmall = 1)
[1] 0.1234568
format(0.123456789, nsmall = 2)
[1] 0.1234568
format(0.123456789, nsmall = 8)
[1] 0.12345679
It adds
The question is general for all functions, but here is a specific example -
# For the logistic regression of the following correlated variables:
C - c(457, 1371, 4113, 12339, 37017, 111051, 333153, 999459)
E - c(0.003858377, 0.014334578, 0.014092836, 0.737950754, 0.996371828,
Will your data support using lme in the 'nlme' package? If yes, I
suggest you switch. This is consistent with a response given recently
by Doug Bates to a crudely related question
(http://finzi.psych.upenn.edu/R/Rhelp02a/archive/68340.html).
You probably know that lmer
On Tue, 14 Feb 2006, Diethelm Wuertz wrote:
Thanks a lot that works fine!
Next problem, if I would have my own package, and the file test.csv
would be located in the data directory
How to use the function data to get
data(test)
resulting in:
test
%y-%m-%d VALUE
1 1999-01-01
Try this:
coef(A)
On 2/13/06, Quin Wills [EMAIL PROTECTED] wrote:
The question is general for all functions, but here is a specific example -
# For the logistic regression of the following correlated variables:
C - c(457, 1371, 4113, 12339, 37017, 111051, 333153, 999459)
E -
06.02.14whatyouwant-coef(A)
06.02.14whatyouwant
Emp C50
0.99891134.7957189 9934.6481222
2006/2/14, Quin Wills [EMAIL PROTECTED]:
The question is general for all functions, but here is a specific example -
# For the logistic regression of the following
I am trying to use the optimizing function genoud() with the snow
package on a couple of i686 machines running Redhat Linux WS4 . I don't
know anything about PVM or MPI, so I just followed the directions in
snow and rgenoud for the simplest method and started a socket cluster.
My function fn for
Dear R People:
I'm using R in a time series class. This class is being
broadcast live to 2 remote sites via closed circuit TV.
My people at the remote sites are having a terrible time
seeing the computer screen as it is broadcast(resolution issues). I have
decided to put together Power Point
Hello all,
(First of all, I'd like to thank all who replied to my previous
question. I have never encountered such a helpful community before.
Thanks for making a R so welcoming.)
To calculate a quantile for normal distributions, one simply uses
qnorm(1-a). But if I would want to do the
Dear aDVISOR,
Hope I am not disturbing you. Can you tell me how I can perform Multivariate
Garch analysis in R. Also please, it is my humble request let me know some
resource materials on Multivariate Garch analysis itself.
Sincerely yours,
--
Arun Kumar Saha, M.Sc.[C.U.]
S T A T I S T I C I A
On Mon, 13 Feb 2006, Adrian Dragulescu wrote:
From the help page of format, nsmall should control the number of digits.
That's not a quote from the help page.
format(0.123456789, nsmall = 10)
[1] 0.1234567890
format(0.123456789, nsmall = 1)
[1] 0.1234568
format(0.123456789,
On Tue, 14 Feb 2006, Jan Danielsson wrote:
Hello all,
(First of all, I'd like to thank all who replied to my previous
question. I have never encountered such a helpful community before.
Thanks for making a R so welcoming.)
To calculate a quantile for normal distributions, one simply
88 matches
Mail list logo
