This isn't really an appropriate forum, and the problem seems to be
related to how you built R (why is dmake being used?)
I've just tried RSQLite on our x86_64 Solaris 10 box (as far as I know
SunOS 5.11 is not yet released but a preview), and it worked out of the
box.
You need help from
On Wed, 31 Oct 2007, Germán Carrillo wrote:
Hello everyone!
Congratulations by this gorgeous software.
My name is Germán Carrillo, I'm writting from Colombia. I don't have much
experience with R, reason why I use the R-GUI (version 2.5.1) on Windows.
I have some problems with DBI package, I
Hi
[EMAIL PROTECTED] napsal dne 30.10.2007 21:49:37:
Excuse me for less detailed e-mail I've sent before, I guess there's a
misunderstanding here.
And still is. You even did not show us how are your data stored in R.
My problem is to use the data in the order I've typed below without
Hi
Henrique Dallazuanna wrote:
Perhaps:
strwidth(expression(bold(text)))
Or set the default font face ...
plot(1:10)
text(4, 6, testing, adj=c(0, 0))
segments(4, 6, 4 + strwidth(testing), 6, col=red)
par(font=2)
text(4, 8, testing, adj=c(0, 0))
segments(4, 8, 4 + strwidth(testing), 8,
Hi,
I get the following error message when trying to use the display
function on the ARM package:
display(model)
Error in .Internal(round(x, digits)) : no internal function round
Looks like some kind of mismatch between the ARM package and some
others? Can I somehow get around it? I have
On Wed, 31 Oct 2007, Kari Ruohonen wrote:
Hi,
I get the following error message when trying to use the display
function on the ARM package:
You seem to mean 'arm' not 'ARM'.
display(model)
Error in .Internal(round(x, digits)) : no internal function round
Looks like some kind of mismatch
Dear All,
Please respond to my personal email in case of queries -- not on the list at
the moment.
Best Marwan
The Faculty of Health Sciences of the American University of Beirut
currently seeks for its Department of Epidemiology and Population Health a
candidate with a Ph.D. degree in
thx for your help,
i checked the caret package out and the tuning works. but i can't find a way to
make a contingency table in order to see the classification result.
e.g. like:
table(outcome NaiveBayes, mydata$code)
Is there something like that?
Julia
Original-Nachricht
Jong-Hoon Kim wrote:
Dear useRs,
I am trying to draw direction fields for some differential equations.
So for I couldn't find much information on that. Could anybody give
me a hint how to draw a direction field using R?
Hi Jong-Hoon,
Here's a start on a function.
Hi
Henrique Dallazuanna wrote:
Perhaps:
strwidth(expression(bold(text)))
Or set the default font face ...
Thanks!
Things can often be so easy!
Roland
plot(1:10)
text(4, 6, testing, adj=c(0, 0))
segments(4, 6, 4 + strwidth(testing), 6, col=red)
par(font=2)
text(4, 8, testing,
Dear All,
I am trying to build a simple ggplot, but where the scale is
reversed, i.e. the largest numbers are on the bottom.
An example of the code I am using is
plotdata-data.frame(x=1:10, y=runif(10))
plot-ggplot()
plot-plot+layer(data=plotdata, mapping=aes_string(x='x',y='y'),
Hi
jiho wrote:
Hello all,
I discovered that the pdf device uses fonts to represent points
symbols (as in plot(...,type=p,...) ). Namely it uses ZapfDingbats
with symbol U+25cf. This can lead to problems when the font is not
available, or available in another version (such as points
Thank you very much for the answer. Re-installation (I did a full
reinstall of R and packages I use) helped and cured the problem. I had
somehow missed the advice to re-install packages when upgrading to 2.6.0
and had only used update.packages().
regards, Kari
On Wed, 2007-10-31 at 09:30 +,
Thank you very much for your answer, even so long after I first
posted the message.
On 2007-October-31 , at 12:00 , Paul Murrell wrote:
Hi
jiho wrote:
Hello all,
I discovered that the pdf device uses fonts to represent points
symbols (as in plot(...,type=p,...) ). Namely it uses
Hi,
I would like to ask how the paste(S1, S2, sep=) function internally
works. Are the two stings copied to a new String?
I have a program where successively strings are build up. First the
program calls an external function and depending on the result it
builds up strings to visualize the
Hi
Perhaps:
Hi,
Try this:
do.call(rbind, lapply(a, function(x)do.call(cbind, x)))
On 31/10/2007, livia [EMAIL PROTECTED] wrote:
Hello,
I would like to use the cbind() function to construct a matrix used in the
middle of a function as following
for (i in 1:1000) {
b[i] -
On 10/31/2007 8:06 AM, livia wrote:
Hello,
I would like to use the cbind() function to construct a matrix used in the
middle of a function as following
for (i in 1:1000) {
b[i] - function(cbind(a[[1]][[i]],
a[[2]][[i]],a[[3]][[i]],...a[[67]][[i]]))
}
Is there an easy way of
Dear useRs,
I am trying to draw direction fields for some differential equations. So
for I couldn't find much information on that. Could anybody give me a hint
how to draw a direction field using R?
--
734 274 2394
Jong-Hoon Kim
[[alternative HTML version deleted]]
Hi list,
I'm using barplot2 form the gplots package to plot a few numbers (I want to add
SD bars later).
However, I would like the y-axis not to start from 0 but 500. When I add the
parameters YLIM, something goes wrong. The graph is not 'cut off' at 500.
Instead the bars seems to sink trough
Mike Waters wrote:
Tony Plate wrote:
[...]
You don't say if this an R-specific problem, or if it afflicts your
computer system clock as well.
Thanks, I should have noted that my computer system clock is fine, and as
far as I can tell it (correctly) believes we are still in Daylight
Saving
Julia,
i checked the caret package out and the tuning works. but i
can't find a way to make a contingency table in order to
see the classification result.
You should read the vignettes for the package at:
http://cran.r-project.org/src/contrib/Descriptions/caret.html
these have the
I'm not sure exactly why it won't work but it could be
that conceptually a bar plot starting at 500 does not
make a lot of sense and the package writers did not
entertain such an idea?
I'd suggest having a look at ?dotchart as an
alternative why to plot the data or perhaps have a
look at
Dear Arne and Ott-Siim,
my personal opinion is that one single package is easiest both for the
useR who wants to keep track of the add-ins he might need and for the
maintainer(s) of the relevant Task Views. I'd prefer to have one single
micEcon as long as I know that I've to look there both for
Hi all,
Does anyone have the knowledge to help me identify a package capable of
forecasting a MULTIPLE regression model? i have a model with one one
dependant variable and 4 independant variables. i would like to forecast
confidence intervals for a few steps ahead...(DENSITY forecasting).
PS i
I have some problems with DBI package, I want to connect to PostgreSQL
database but I don't know the way to do it.
pg - dbDriver(PostgreSQL)
... nope. That won't work.
you need the RdbiPgSQL package from Bioconductor. Along with DBI, as I
recall.
On Wed, 2007-10-31 at 11:52 +0100, Dieter Vanderelst wrote:
Hi list,
I'm using barplot2 form the gplots package to plot a few numbers (I
want to add SD bars later).
However, I would like the y-axis not to start from 0 but 500. When I
add the parameters YLIM, something goes wrong. The
OK, I fixed it myself! Here's the code. Of course, it mostly seems simple
once one gets it working... Thanks Jim. Bryan
sample.info = read.table(input.file.name, sep=,, as.is=TRUE, nrows=3) #
get the first three lines with sample info in character format
sample.names = sample.info[1,]
On 10/31/07, Tamara Steijger [EMAIL PROTECTED] wrote:
Hi,
I would like to ask how the paste(S1, S2, sep=) function internally
works. Are the two stings copied to a new String?
I'm not 100% sure, but I'd suspect so, as this is the default
behaviour in pretty much every programming language.
Hi!
I was trying to install the RWinEdt library on a Windows Vista machine which
had R 2.6.0, WinEdt 5.4 and the RWinEdt1.7.8.zip file in the library folder
under the ...\R\2.6.0 directory. I am able to install fine form the local zip
but I can't launch winedt from within R for some reason.
Dear R-help community,
I need to find out how to make the black axis line around my plot thicker in
levelplot, I cannot find the correct command in ?levelplot or ?par - or if it
is
there I cannot get it to work - any help would be more that appreciated!
Here's an example script - I would want
I have two variables which show a typical quantile relation
I would like to fit quantile regression models based on both logarithm and
polynomial of second order functions within quantreg.
Any help appreciated...
Cheers
Duccio
Herewith the values:
--
var1var2
0.964290.00138
1
Alexander Nervedi wrote:
Hi!
I was trying to install the RWinEdt library on a Windows Vista machine which
had R 2.6.0, WinEdt 5.4 and the RWinEdt1.7.8.zip file in the library folder
under the ...\R\2.6.0 directory. I am able to install fine form the local zip
but I can't launch winedt
Has anyone been able to successfully query SAS libraries from R using RODBC?
While I believe I have configured the SAS ODBC Driver properly (I can query
SAS libraries from MS Access) when I try to query SAS libraries from R, I
only seem to get back the column names of the table I am querying.
Has anyone been able to successfully query SAS libraries from R using RODBC?
While I believe I have configured the SAS ODBC Driver properly (I can query
SAS libraries from MS Access) when I try to query SAS libraries from R, I
only seem to get back the column names of the table I am querying.
I'm not sure if i got the idea of what you want. At fist i thought that
you wanted this output:
Output:
290 380
380 440
440 443
443 468
to get it, you can use the following function:
intervals - function(Input)
{
all - c(Input$Start, Input$End)
numbers -
Dear list members,
I look for a way (or alternative) to specify initial values when estimating
linear mixed models in R, and to avoid iterative estimation.
This is a way to control specific parameter values (eg. variance parameter
values) such that the result (F-value) is based on them. This
Another approach is to use the my.symbols function in the TeachingDemos
package (in place of the points function) and define how you want your
circles represented (a polygon with enough sides is a good approximation
to a circle for most cases).
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
A simple approach is to swap x and y and rotate the theta parameter -90
degrees (play around with this until it is what you want).
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
Dear R users,
I have following code, estimating Poisson log-likelihood a number of times:
poisson.loglik - function(mu, y){
n - NROW(y)
logl - sum(y)*log(mu)-n*mu
return(-logl)
}
estimates - numeric(1e5)
for(i in seq_along(estimates)){
estimates[i] - optim(par=1,
Dear all,
Using the data set and code below, I am interested in modelling how egg
temperature (egg.temp)
is related to energy expenditure (kjday) and clutch size (treat) in
incubating birds using the
lme-function. I wish to generate the F-distribution for my model, and have
tried to do so
Helo,
please look at the log below: after loading the fSeries library, I can not
use the log function. Is this a bug or what am I doing wrong?
Because of this, I'm unable to use the garch library.
thanks a lot for any help,
Balazs Torma
log(1)
[1] 0
require(fSeries)
Loading required
Giulia Barbati wrote:
Hello,
using the Design library, and the following command (from the Harrell's book,
example at Ch.20):
dd - datadist(rx, age, wt, pf, pf.coded, heart, map, hg, sz, sg, ap, bm)
options(datadist=='dd')
That command doesn't appear in the book. The == should be a
For drawing the arrows look at the my.symbols and ms.arrows functions in
the TeachingDemos package. If you need help with determining the
directions of the arrows, then perhaps the Ryacas package or other link
to a mathematics package will help.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data
Hello all...
I am just starting to teach myself Bayesian methods, and am
interested in learning how to use UMacs. I've read the
documentation, but the single example is a bit over my head at the
level I am at right now. I was wondering if anyone has any simple
examples they'd like to
Hello,
I have a question about one of the functions from package clim.pact.
I am planning to perform canonical correlation analysis (CCA) with temperature
stations data, I want to use clim.pact package for that. As I understood I
should be able to use the function stations2field in order to
Are you using the latest version of R - with the patch?
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Balazs Torma
Sent: Wednesday, October 31, 2007 11:44 AM
To: r-help@r-project.org
Cc: [EMAIL PROTECTED]
Subject: [R] problem with package fSeries
Thanks Thomas ! I am trying to draw random sample from a household survey
which has 80,000 observations.
rural is name of the dataset, while iwt is survey weights assigned to each
observation.
the resulting error are :
z=sample(rural,5000,replace=TRUE, Prob=rural$iwt)
Error in sample(rural,
Ok, thanks to Prof Brian Ripley and to Elijah, I did the connection to
PostgeSQL database with RODBC, it's simple.
I will try with RdbiPgSQL in the next days.
Thanks a lot.
2007/10/31, [EMAIL PROTECTED] [EMAIL PROTECTED]:
I have some problems with DBI package, I want to connect to
On Oct 31, 2007, at 9:48 AM, Greg Snow wrote:
A simple approach is to swap x and y and rotate the theta parameter
-90
degrees (play around with this until it is what you want).
Thank you. That does exactly what I wanted. Oddly enough, I had
already swapped X and Y trying to get this to
Sergey Goriatchev wrote:
[snip]
I cannot run this because I always get an error message:
Error in optim(par = 1, poisson.loglik, lower = 0, method = L-BFGS-B, :
non-finite finite-difference value [1]
Could someone enlighten me what that means and why this happens?
THanks in
Dr. Lumley and Prof. Ripley,
Thank you very much for your helpful responses.
Have you found any particular distribution of Linux to work well with
64-bit R? For the cluster, I am currently considering Debian (since it
seems popular) and SUSE (since Matlab runs on it), but I remain open to
Got a little bit of a different answer:
x - Start End
+ 440 443
+ 380 443
+ 290 468
x.in - read.table(textConnection(x), header=TRUE)
# create matrix to determine queue size (overlap)
x.q - rbind(cbind(x.in$Start, 1), cbind(x.in$End, -1))
# sort
x.q - x.q[order(x.q[,1], x.q[,2]),]
This is a repeat posting from 28 Oct that generated no replies. I'm
hoping someone has some advice since I'm still stuck ...
I am trying to use tree to partition a data set. The data set has 3924
observations. Partitioning works for small subsets of the data, but when
I use the entire data set,
Hi y'all,
I wrote two functions; 1) sparsely samples a point pattern and
calculates the Clark Evans R value and Z score, 2) calls the first
function multiple times using a for loop and generates a histogram.
I would like to know if there is a more efficient way of calling the
first function
Hello,
I am using textplot function in gplots package to put some model
output inside a PDF file, but it does not seem to work properly with
PDF.
I am doing follwing:
pdf(file=C:/..., paper=a4, width=8, height=12)
.model - lm(.model.formula, data=database)
Hi all,
I've got R for windows 2.3.1
I'm trying to modify my lattice histograms.
The following does reset the y-label to blue, but I cannot turn off the label
(alpha=0), nor shift label position (I'm assuming that is what lineheight is
for). Is 'lineheight' the equivalent of 'line' in mtext?
Am
Since you do not change anything in your call to the function, you could
just use the replicate function instead of the for loop (I don't know
that it will speed things up much, but it is a bit more readable).
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
try:
trellis.par.set(axis.line=list(lwd=4))
On 10/31/07, Jenny Barnes [EMAIL PROTECTED] wrote:
Dear R-help community,
I need to find out how to make the black axis line around my plot thicker in
levelplot, I cannot find the correct command in ?levelplot or ?par - or if it
is
there I
Dear R-help,
what are the limits on xtest?
NOT_A.rf - randomForest (log10(Y[!A] ) ~ . , data = notA_desc ,
proximity=T ,xtest = A_desc)
*** caught segfault ***
address 0x9cdd000, cause 'memory not mapped'
Segmentation fault
I don't think that the matrix are large:
notA_desc is 651 obs
After installing the SciViews-R Console I am getting the following error
message when I try to run any of the commands.
Error in file.info(fn - c(...)) : invalid filename argument
See also the attached file Doc1.
-
I'm looking for a function in R similar to t.test() which was generously
pointed out to me yesterday, but which can be used for normally
distributed data.
To recap yesterday:
x - scan()
1: 62 52 68 23 34 45 27 42 83 56 40
12:
Read 11 items
alpha- .05
t.test(x)
One Sample t-test
On Wed, Oct 31, 2007 at 03:56:37PM -0400, Zembower, Kevin wrote:
I'm looking for a function in R similar to t.test() which was generously
pointed out to me yesterday, but which can be used for normally
distributed data.
...
x100-sample(x, 100, replace=TRUE)
I think that I should be able to
Given a matrix B, where B=A'A, how can I find A?
In other words, if I have a matrix B which I know is another matrix A times
its transpose, can I find matrix A?
Thanks,
Mike
__
R-help@r-project.org mailing list
Daniel, thanks, I should have remembered this, too; I've seen it and
worked with it before. Thanks.
-Kevin
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Daniel Lakeland
Sent: Wednesday, October 31, 2007 4:04 PM
To: r-help@r-project.org
Subject: Re: [R]
Dear R-helpers,
I wasn't able to find out how to override the alphabetical ordering
of the rows and columns in a vcd::mosaic plot. I would like to have
them each ordered by numerical values in a different column of the
data frame that contains the contingency data.
I would be most grateful
B is symmetric by definition; if it's also real positive-definite then A
is the upper triangular factor of the Choleski decomposition, and you can
use
chol(B)
to get A.
On Wed, 31 Oct 2007, Michael Gormley wrote:
Given a matrix B, where B=A'A, how can I find A?
In other words, if I have a
There is a z.test function in the TeachingDemos package, but it is
mainly for learning purposes to ease students into using t.test.
The rule of use the normal for n30 comes from the days when
computations were done by hand using tables and people did not want to
carry around t-tables with 100's
chol(B) doesn't give the original A, which I believe is what Mike wants.
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Katharine Mullen
Sent: Wednesday, October 31, 2007 4:08 PM
To: Michael Gormley
Cc: [EMAIL PROTECTED]
Subject: Re: [R] Find A, given B
On 1/11/2007, at 9:13 AM, Michael Gormley wrote:
Given a matrix B, where B=A'A, how can I find A?
In other words, if I have a matrix B which I know is another matrix
A times
its transpose, can I find matrix A?
You can't, because A is not unique. You can easily find ***a***
solution.
On Thu, 1 Nov 2007, Rolf Turner wrote:
On 1/11/2007, at 9:13 AM, Michael Gormley wrote:
Given a matrix B, where B=A'A, how can I find A?
In other words, if I have a matrix B which I know is another matrix
A times
its transpose, can I find matrix A?
You can't, because A is not unique.
On Wed, 31 Oct 2007, Michael Kubovy wrote:
Dear R-helpers,
I wasn't able to find out how to override the alphabetical ordering
of the rows and columns in a vcd::mosaic plot. I would like to have
them each ordered by numerical values in a different column of the
data frame that contains the
All,
I constructed a fairly simple plot using the below and scatterplot3d.
However, the axis locations in scatterplot3d seem to 'hide' quite a bit
of the points in my graphic based on the viewing angle, hence not
adhering to Tufte's rules. I tried to determine if I could rotate the
graph and
Michael Gormley wrote:
Thanks for your help, all those who submitted responses. I do not need a
specific matrix A, any solution will do. With this said, is it possible to
specify the dimensions of the A matrix in the decompostion? For example, if
A is a 2X1 matrix then A'A=B would be a
On 31/10/2007 5:50 PM, Achim Zeileis wrote:
On Wed, 31 Oct 2007, Michael Kubovy wrote:
Dear R-helpers,
I wasn't able to find out how to override the alphabetical ordering
of the rows and columns in a vcd::mosaic plot. I would like to have
them each ordered by numerical values in a
Would this be good enough:
# example using builtin BOD data frame
write.csv(cbind(Z = row.names(BOD), BOD), row.names = FALSE)
On Oct 31, 2007 7:04 PM, Earl F. Glynn [EMAIL PROTECTED] wrote:
I've always been frustrated that R never puts a label on the column name for
a column of row names.
The nice thing about R is that you can always extend it in any manner
that you want. Here is some code (that you could put in a function
wrapper) that will do what you want. It just creates the header,
writes it out, and then the dataframe itself without a column header:
out -
On 31/10/2007 5:51 PM, Bret Collier wrote:
All,
I constructed a fairly simple plot using the below and scatterplot3d.
However, the axis locations in scatterplot3d seem to 'hide' quite a bit
of the points in my graphic based on the viewing angle, hence not
adhering to Tufte's rules. I
Hi
I'm trying to use RWeka to use a NaiveBayes Classifier(the Weka
version). However it crashes whenever there is a NA in the class
Gender
Here is the.code I have with d2 as the data frame.
The first call to NB doesn't make R crash but the second call does.
NB -
Dear sir,
I am Shibu John from Thrombosis Research Institute India. It is a
multidisciplinary organisation concerned with the interrelated problems of
thrombosis and atherosclerosis.
I was searching for Cochran armitage trend test program in R. Then I had
seen your R coding for C-A trend test.
On 11/1/07, Folkes, Michael [EMAIL PROTECTED] wrote:
wohoo! A bug. I know I wasn't the first to find it, but I had a hunch this
time that it was going to be a bug and not just my poor understanding. I
wasn't really thinking as just ylab= works fine too. I think I was hoping
to squeeze
I've just updated from R-2.5.1 to R-2.6.0, under Windows XP. Is
there a simple way to reload into R-2.6.0 all the packages I had
loaded under R-2.5.1?
Thanks,
John
=
John Field Consulting Pty Ltd
10 High Street, Burnside SA 5066 Australia
Phone +61 8 8332 5294 or +61 409
Hi,
I'm new to R (and statistics) and my boss has thrown me in the deep-end with
the following task:
We want to evaluate the impact that sampling size has on our ability to create
a robust model, or evaluate how robust the model is to sample size for the
purpose of cross-validation
Shibu John wrote:
Dear sir,
I am Shibu John from Thrombosis Research Institute India. It is a
multidisciplinary organisation concerned with the interrelated problems of
thrombosis and atherosclerosis.
I was searching for Cochran armitage trend test program in R. Then I had
seen your R
83 matches
Mail list logo