Dear list,
I get random terminations of R-2.6.2 when using RODBC-1.2_3 and a remote
PostgreSQL-8.2.6 database, unixodbc-2.2.11-16,
odbc-postgresql-1:08.01.0200-2.1 on Ubuntu Linux 7.10.
The exact messages can be seen below. It seems that there is some memory
corruption going on. But it seems
Hello,
i have a set of data X_1,...,X_n of an unknown density function f.
I want to estimate the density function f with wavelets.
Can I do this with the package wavethresh ? How ?
Thanks.
--
View this message in context:
Hello list
I have a theoretical differential equation model (two coupled differential
equations) describing the acidification of dairy cultures in milk:
dX/dt = f(H, param)*X
dH/dt = g(param)*dX/dt
pH = -log10(H)
I also have actual data of the same.
I would like to estimate the parameters of
Hi,
temp1[1:6,]
_XXX BB_XXX C_XXX D_XXX E_XXX SkHep1_XXX
500.75285 0.30630 0.56330 0.45755 0.496650.74685
630.35930-0.22725 0.04725 0.53955 0.563900.51170
69 -1.57515-1.67610 -1.61695 -1.11685 -0.00085 -1.89155
910.84640
Neither of reduced_model not full_model will ever be NULL. If this
fails, they inherit from class try-error. See ?try.
So you want
inherits(reduced_model, try-error) || inherits(full_model, try-error)
On Thu, 14 Feb 2008, Juliet Hannah wrote:
Dear R Users,
I have the following glm,
Try the odesolve package or the nlmeODE package which builds on odesolve.
Best regards
Frede Aakmann Tøgersen
Scientist
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of Genetics and Biotechnology
Blichers Allé 20, P.O. BOX 50
DK-8830 Tjele
Phone: +45 8999 1900
Direct: +45
Or just read the article
Ligges, U. (2006): R Help Desk: Accessing the Sources. R News 6 (4),
43-45. http://cran.r-project.org/doc/Rnews/
which was written particularly for answering this kind of questions.
Best,
Uwe Ligges
Gabor Csardi wrote:
RSiteSearch is your friend. E.g.:
Hi.
I have a lot of plots that i need to store as components of a vector.
Here you have an example
DF-data.frame(A=rnorm(12),B=c(1:12),C=rep(c(1:4),each=3))
Q-c()
length(Q)-3
PAC-geom_point(mapping=aes(x=A,y=B),data=DF[DF$C==1,])
P-ggplot()+PAC
P #generates the plot
Q[1]-P
Warning
Ng Stanley wrote:
Hi,
temp1[1:6,]
_XXX BB_XXX C_XXX D_XXX E_XXX SkHep1_XXX
500.75285 0.30630 0.56330 0.45755 0.496650.74685
630.35930-0.22725 0.04725 0.53955 0.563900.51170
69 -1.57515-1.67610 -1.61695 -1.11685 -0.00085
[EMAIL PROTECTED] wrote:
Hi,
i'd like to plot some data that I have with the value on the x axis and
freq
on the y axis.
So, I need to calculate the freq a value is seen within my data vector
for example, say i have a vector of data
data=c(1,1,1,4,5,5,6)
I
Gallego Liberman, Matias wrote:
Hi.
I have a lot of plots that i need to store as components of a vector.
Here you have an example
DF-data.frame(A=rnorm(12),B=c(1:12),C=rep(c(1:4),each=3))
Q-c()
length(Q)-3
You need a list, hence replace the last two lines by:
Q - vector(mode=list,
Which package(s) contain something on latent class analysis?
Tom
--
++
| Tom Backer Johnsen, Psychometrics Unit, Faculty of Psychology |
| University of Bergen, Christies gt. 12, N-5015 Bergen, NORWAY |
| Tel : +47-5558-9185
You'll need to store it in a list rather than a vector.
library(ggplot2)
DF-data.frame(A=rnorm(12),B=c(1:12),C=rep(c(1:4),each=3))
Q - list()
PAC-geom_point(mapping=aes(x=A,y=B),data=DF[DF$C==1,])
P-ggplot()+PAC
Q[[length(Q) + 1]] - P
Q[[1]]
You could shorten this into
library(ggplot2)
check function lca() in package 'e1071'.
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://med.kuleuven.be/biostat/
Dear All:
Do you have any experience on calling winbugs in R enviorment? I uploaded
R2winbgus, pktools, BRugs, ect. in order to call winbugs. In the end, I try to
use bugs(). The error message still displays that bugs() can not find.
Do you have any similar issue? If so, could you give
Try this:
cbind(Row=rownames(x), stack(x)[,2:1])[cbind(Row=rownames(x),
stack(x)[,2:1])[,3]3,]
On 14/02/2008, dinesh kumar [EMAIL PROTECTED] wrote:
Hi R users
I am a new user in the field of R.
I want to subset or reshape a data.frame.
For example I have a matrix like
A B C D E
your r2 model corresponds to method = Laplace.
r4=lmer(Y~X+(1|Subject),family=binomial(link=logit),method=Laplace) is
equivalent to r2.
Bests,
Abderrahim
- Original Message -
From: Daniel Malter [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Friday, February 15, 2008 12:50 AM
Hi All,
I have a basic statistics question you might be able to help.
I have some data collected within a single population (say of chronic illness)
over a period of time (number under observation is constant, say 100 subjects).
The question is how could I test whether the frequency of an
Xin wrote:
Dear All:
Do you have any experience on calling winbugs in R enviorment? I uploaded
So you installed them, I assume?
R2winbgus, pktools, BRugs, ect. in order to call winbugs. In the end, I try
to use bugs(). The error message still displays that bugs() can not find.
What
and also check packages: flexmix, polca, depmix and a few others ...
On 15 Feb 2008, at 11:01, Dimitris Rizopoulos wrote:
check function lca() in package 'e1071'.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of
Hi,
just for the record, although I don't think it's relevant (!)
-
sessionInfo()
R version 2.6.0 (2007-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=English_United Kingdom.1252;LC_CTYPE=English_United
Kingdom.1252;LC_MONETARY=English_United
Thanks!
Ingmar Visser wrote:
and also check packages: flexmix, polca, depmix and a few others ...
On 15 Feb 2008, at 11:01, Dimitris Rizopoulos wrote:
check function lca() in package 'e1071'.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of
samsr wrote:
Hi,
I need to plot a matrix using image() such that negative values are easily
distinguishable from posittive values, while also maintaining a gradation in
color with magnitude. How can I set ranges for colors in order to achieve
this. Thanks.
Hi Sam,
This will give you red
Hi
Is there a correspondance between R commands and other software commands
(GAUSS mainly).. I'd like to translate some Gauss into R..
Thanks for your help
Naji
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Could you send us the output of sessionInfo() please so we can see
which version of the lme4 package you are using? In recent versions,
especially the development version available as
install.packages(lme4, repos = http://r-forge.r-project.org;)
the PQL algorithm is no longer used. The Laplace
Hi,
I want to add boxplots to a scatterplot:
plot(x,y, xlim=c(80,120),ylim=c(80,120))
boxplot(y,add=TRUE,at=118)
boxplot(x,add=TRUE,at=118,horizontal=TRUE)
How can I control the width of the boxes (say, I'd like them to be of
width 3 in the variables' scales). I've tried the width parameter
but
Dear lists,
Anysuggestion on how to write a loop function which has nx2 matrix as an
input, where each of the row represents an interval. the function should return
a mx2 matrix containing the no matched interval of x?? i've tried my own
function but it's only applies for 2x2 matrix.
Hello;
I've used stop() previously without any problem.
Now, I'm using a condition and getting an error when stop() is called.
Mi condition is:
if(dim(r.temp)[1] != (length(jul.names)+6)) stop(buffered data not
available)
An the R' response is:
source(my_program.R)
Error in
Try this:
x - sample(5, 100, rep=T)
y - rnorm(100)
plot(x,y, xlim=c(min(x)-.5, max(x)), ylim=c(min(y), max(y)+1))
boxplot(y, pars=list(boxwex=0.3), at=min(x)-.5, add=T)
boxplot(x, pars=list(boxwex=0.3), at=max(y)+1, add=T, horizontal=T)
On 15/02/2008, Marcin Kozak [EMAIL PROTECTED] wrote:
Hi,
Dear R admin,
I sent so many emails but none was posted.
What happened?
Thanks
Ilham
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Perhaps I am missing the obvious, but can't you simply write a small
function that, given each particle object, returns a string with
the required output format (probably a couple of paste calls with
sep and collapse set), and then use sapply and this function on the
list of particles to
Rainer M Krug r.m.krug at gmail.com writes:
Hi
I am running analysis with Kcross from the package spatstat and I am getting
a message that R can not allocate enough memory for a vector of 900MB.
R seems to be running towards the 2GB limit per process.
The dataset is not to big (ca
Dear Sir/madam,
I use lmer to extract model in your package of lme4. It seems works
well. But the problem is when I use anova/summary the extracted model,
no p-value is shown at all. In previous version(nlme), I mainly use
p-value to judge which term is significant or not, and then make a
Hi R help,
I run my data in nnet with skip layer, factor response (with 0 1
values) and explicitly put softmax=T to compare the result of the
default nnet with no softmax specification. I assume this should give
me the same result. I got the result the default one, but not the
softmax version
Hi R help,
I run my data in nnet with skip layer, factor response (with 0 1
values) and explicitly put softmax=T to compare the result of the
default nnet with no softmax specification. I assume this should give
me the same result. I got the result the default one, but not the
softmax version
Hi List,
I have a problem plotting data using the lineplot.CI command in the sciplot
package.
I want to plot the data of 2 experimental cases using different lines (traces).
Time is on the X-axis. The tricky thing is that the data collection in the
second case started later than for the first
Dear expeRt,
I would like to know whether a function similar to str_replace() in
php exists in R.
Looking forward to hearing from you.
Best,
=
Center of Excellence for Science and Innovation Studies,
Royal Institute of Technology
I don't know php very well, but perhaps you need
?sub
Gabor
On Fri, Feb 15, 2008 at 03:54:19PM +0100, Dong-hyun Oh wrote:
Dear expeRt,
I would like to know whether a function similar to str_replace() in
php exists in R.
Looking forward to hearing from you.
Best,
Thanks for your replies. My real problem is that, for my real data, I get
basically the same results from r2 and r3 (so to speak), but the coefficient
estimates and significance levels for r1 are very different from those of r2
and r3. And therefore, I do not know which of the results to trust and
Have a look at boxwex in boxplot. Something like
boxplot(y,add=TRUE,at=118, boxwex=3)
may give you what you want.
--- Marcin Kozak [EMAIL PROTECTED] wrote:
Hi,
I want to add boxplots to a scatterplot:
plot(x,y, xlim=c(80,120),ylim=c(80,120))
boxplot(y,add=TRUE,at=118)
Looks like you meant to use 'TRUE' not 'T'. The latter is probably a
variable you defined of length 1.
Double-posting at least halves your chances of getting help, so please
stop doing it.
On Fri, 15 Feb 2008, G Ilhamto wrote:
Hi R help,
I run my data in nnet with skip layer, factor
Try this:
hist(x[,1], ylim=c(0,max(x[,2])))
lines(x[,2], col=blue, lty=2)
On 15/02/2008, Xin [EMAIL PROTECTED] wrote:
Dear:
I have two series data, which are:
p1-obs p1-exp
13 13
20 23.296
22 18.139
19 14.972
10 12.576
16 10.658
8
Rather than go over all of this again, see here
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
--
View this message in context:
http://www.nabble.com/lmer-in-package-of-lme4-tp15501802p15502440.html
Sent from the R help mailing list archive at Nabble.com.
WANG ping.wang at newcastle.ac.uk writes:
I use lmer to extract model in your package of lme4. It seems works
well. But the problem is when I use anova/summary the extracted model,
no p-value is shown at all. In previous version(nlme), I mainly use
p-value to judge which term is significant
On Fri, 2008-02-15 at 15:18 +0100, Dieter Vanderelst wrote:
Hi List,
I have a problem plotting data using the lineplot.CI command in the sciplot
package.
I want to plot the data of 2 experimental cases using different lines
(traces). Time is on the X-axis. The tricky thing is that the
On Fri, Feb 15, 2008 at 8:18 AM, Dieter Vanderelst
[EMAIL PROTECTED] wrote:
Hi List,
I have a problem plotting data using the lineplot.CI command in the sciplot
package.
I want to plot the data of 2 experimental cases using different lines
(traces). Time is on the X-axis. The tricky
Have a look at
length(iris$Species)
length(prev)
--- Carla Rebelo [EMAIL PROTECTED] wrote:
Hello!
May you help me? I'm trying to do a contingency
table using this
data(iris)
library(rpart)
modelo - rpart(Species ~., iris)
prev - predict(modelo, iris)
Finally the contingency
[EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]
orld.com:
Hi All,
I have a basic statistics question you might be able to help.
I have some data collected within a single population (say of
chronic illness) over a period of time (number under observation is
constant, say 100
# Dear list,
# I am an R-beginner and
# spent the last days looking for a method to insert tables produced
# with R into a word document. I thought about SPPS: copy a table from
# an SPO-file and paste it into a word document
# (if needed do some formatting with that table).
# Annother idea was,
Hi,
first of, my apologies for the possibly off-topic message regarding a
problem when installing the XML package (I meant to post at omegahat-
help, but I get 404 at http://www.omegahat.org/mailman/listinfo)
When installing the XML package using:
install.packages(XML)
everything seems
Howdee,
*** I know that the lmList() function exists, yet I don't want to use it.
***
Would anyone be kind enough to tell how I can apply the function lm() to
each level of a given factor so to obtain the intercept and slope for each
factor level within a matrix?
For instance, suppose a
Here are 2 ways:
1. Display it on the screen and hit the PrintScreen button.
Then paste that into Word and use Word's image editor
to crop it and expand or shrink it appropriately.
2. Run this in R:
capture.output(CrossTable(...whatever...), file = clipboard)
Paste the clipboard into Word and
Hi all.
I have certain problem with the extraction of data from matrices. I know how
extract elements from a matrix, but what I would like to do is to extract a
concrete number, not an element. For example if I have this matrix X:
[,1] [,2] [,3] [,4]
[1,]3 13 23 33
[2,]9 19
Hi again,
Adding further information to my own query, this function gets to the core
of the problem, which I think lies in the behaviour of 'integrate'.
-
function (x, meanlog = 0, sdlog = 1, ...) {
require(stats)
integrand - function(t, x, meanlog,
You need a later version of libxml2 (or at least, a different one as
possibly an older one would also work). I met this on Windows -- 2.6.23
is old (more than 2 years), and I believe 2.6.28, 2.6.30 and 2.6.31 work.
The Omegahat lists have been defunct for a while, AFAIK. I'd have used
Thank you Gabor,
but with 1. or 2. I don´t get
a table, which can be changed and formatted using the WORD table functions
(like changing column width and formatting text in rows).
The result of 1. ist an bitmap and
the result of 2. an ASCII-Table
A rich formatted (*.rtf) table/object would be
That would be a good idea. I'm finishing up a much more comprehensive guide
that will become available soon.
Thanks for the suggestion,
Michael
2008/2/14 Felix Andrews [EMAIL PROTECTED]:
I just discovered a good introduction to the RGtk2 package:
RShowDoc(overview2, package=RGtk2)
I have
On 2/15/2008 11:00 AM, Marc Belisle wrote:
Howdee,
*** I know that the lmList() function exists, yet I don't want to use it.
***
Would anyone be kind enough to tell how I can apply the function lm() to
each level of a given factor so to obtain the intercept and slope for each
factor
Z is a matrix and when I run the following line, it creates a prediction
estimate using each column, how can I get it an estimate for each individual
number. I have tried changing Z to a data.frame, but this does not do it
either.
model.lm-lm(w~x)
pred.est - predict.lm(model.lm,
You will have better luck in soliciting responses if you
follow the line at the bottom of every r-help message and
post minimal reproducible code. Also read the posting
guide.
On Fri, Feb 15, 2008 at 8:29 AM, G Ilhamto [EMAIL PROTECTED] wrote:
Dear R admin,
I sent so many emails but none was
I usually use the
write.table()
function (with tab as the delimiter) to write to
a text file, then copy/paste into Word,
then use Word's convert text to table command.
Obviously, if one needs to do this for many,
many, tables, the amount of manual manipulation
is excessive. That's when
Try also:
mat[mat==29] - NA
On 15/02/2008, _Fede_ [EMAIL PROTECTED] wrote:
Hi all.
I have certain problem with the extraction of data from matrices. I know how
extract elements from a matrix, but what I would like to do is to extract a
concrete number, not an element. For example if I
Dear All:
Using the Excel formula below I am able to estimate
daily variance(Daily VarPD)for PD,but I can't figure
out how to estimate the Weekly variance for PD. Any
hints?
Daily VarPD = (PD * (1 - TD))/TD + VarTD *
((PD*(1-TD)+ PD^2*TD))/TD^3
WK TDVARTD PD Daily
On Fri, Feb 15, 2008 at 1:53 PM, Kathy Gerber [EMAIL PROTECTED] wrote:
Earlier today I sent a question to Frank Harrell as an R developer with
whom I am most familiar. He suggested also that I put my questions to
the list for additional responses. Next month I'll be giving a talk on
R as
Hi everybody!
I'm trying to use ggplot2 to return a plot from a function (so I can add
something or alter it then). Unfortunately, if I add a mapping to a
layer in the function, the variable *name* is stored in the layer,
rather than the variable's *value* - so that after the function returns
On 2/15/08, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Marlin Keith Cox said the following on 2/15/2008 11:39 AM:
Dear R Users, close to the end of this I used wireframe to create a 3D plot
from a matrix. The x and y axis tick labels (1-6) for each were created
from the matrix being a
Hi all !
I would like to simulate spatial lattice/areal data with a conditional
autoregressive (CAR) structure, for a given neighbouring matrix and for a
autocorrelation rho.
Is there any package or function in R to perform it ?
I found the function CARsimu in the hdeco library, but this is not
Dear R Users, close to the end of this I used wireframe to create a 3D plot
from a matrix. The x and y axis tick labels (1-6) for each were created
from the matrix being a 6X6 matrix. I need the axis tick labels to be the
row and column headings which you can see in the output (mat.x). I have
Marlin Keith Cox said the following on 2/15/2008 11:39 AM:
Dear R Users, close to the end of this I used wireframe to create a 3D plot
from a matrix. The x and y axis tick labels (1-6) for each were created
from the matrix being a 6X6 matrix. I need the axis tick labels to be the
row and
Look this can help you:
mtz-matrix(1:20,4,5)
mtz
[,1] [,2] [,3] [,4] [,5]
[1,]159 13 17
[2,]26 10 14 18
[3,]37 11 15 19
[4,]48 12 16 20
mtz[2,3]
[1] 10
mtz
[,1] [,2] [,3] [,4] [,5]
[1,]159 13 17
[2,]26
Dear Colleagues, please help distribute the following announcement of a
data mining competition. Thanks, -XG
=
Data Mining Competition 2008
=
Website: http://dms.stat.ucf.edu/competition08/home.htm
Department of Statistics Actuarial
On Fri, 2008-02-15 at 10:45 -0800, Tim Smith wrote:
Hi,
Is there any clustering package in R that can cluster with ordinal data?
thanks!
daisy() in recommended package 'cluster' can generate dissimilarities
for ordinal data using Gower's general (dis)similarity coefficient for
mixed data.
In that case you shouldn't be using CrossTable in the first place.
Suggest you look into the functions mentioned under See Also in
?CrossTable or look at the source of CrossTable and modify it for
your purposes.
On Fri, Feb 15, 2008 at 11:52 AM, [EMAIL PROTECTED] wrote:
Thank you Gabor,
but
?which
Does something like this help?
which(aa==29)
aa[ -which(aa==29)]
--- _Fede_ [EMAIL PROTECTED] wrote:
Hi all.
I have certain problem with the extraction of data
from matrices. I know how
extract elements from a matrix, but what I would
like to do is to extract a
concrete
On 15-Feb-08 13:29:06, G Ilhamto wrote:
Dear R admin,
I sent so many emails but none was posted.
What happened?
Thanks
Ilham
Are you subscribed to r-help? If not, you will not
receive copies of any mails sent out by the list
(though if a non-subscriber mails to the list, the
mail will be
Udo,
CrossTable() knows nothing about the proprietary formats of Word tables.
The output is designed for the R console, using a fixed width
(monospace) font. The cell boundaries are drawn using standard ASCII
characters.
There are no plans at present to modify the format of the output to
See spautolm function in spdep package
On 15/02/2008, Dae-Jin Lee [EMAIL PROTECTED] wrote:
Hi all !
I would like to simulate spatial lattice/areal data with a conditional
autoregressive (CAR) structure, for a given neighbouring matrix and for a
autocorrelation rho.
Is there any package
Dear:
I have two series data, which are:
p1-obs p1-exp
13 13
20 23.296
22 18.139
19 14.972
10 12.576
16 10.658
8 9.0794
4 7.7625
3 6.6536
5 5.7141
7 4.9146
4 4.2321
0 3.648
4 3.1471
4
Here is an example I just did. I _LOVE_ lmList(). It is such a great
convenience compared to the following. But right now it doesn't work
for family=binomial, so I had to do it differently. This is not
your example, but it should give you an idea. The [[5]] was because I
was interested only
Hi,
I am using solve.QP (from quadprog) to solve a standard quadratic
programming problem: min_w -0.5*w'Qw st ... I would like solve.QP to do two
things: 1) to start the optimization from a user-supplied initial
condition; i.e., from a vector w_0 that satisfies the constraints, and 2) to
return
For those of you who want to learn more about the history of the R
project: There will be an invited lecture by John Fox and Kurt Hornik at
this year's useR! conference in Dortmund in August (...unfortunately a bit
too late for Kathy) about The Past, Present, and Future of the R Project
see
I have R 2.6.2 installed on a Mac Probook with the lastest version of Leopoard
(including recent graphics update).
When I use either command-C or the Edit menu to copy the contents of an
active quartz graphics window, the contents of the window disappears (much like
cut).
Has anyone had
Dear Kathy,
As Achim has mentioned, I've been doing interviews with members of the R
Core team and with some other people central to the R Project. Although I
haven't entirely organized and finished reflecting on this material, the
following factors come immediately to mind:
(1) Doug has already
G'day Jorge,
On Fri, 15 Feb 2008 17:51:16 -0500
Jorge Aseff [EMAIL PROTECTED] wrote:
I am using solve.QP (from quadprog) to solve a standard quadratic
programming problem: min_w -0.5*w'Qw st ... I would like solve.QP to
do two things: 1) to start the optimization from a user-supplied
initial
84 matches
Mail list logo