Thanks David for prompt reply. I agree with you. However, I still fail to get
the survfit function to work with newdata. In my previous example I changed the
column names of testX matrix and I still fail.
colnames(testX)-names(coxph.model$coefficients)
sfit-
Hi,
I'm using a command in bioconductor that seems to require a package called
hu6800cdf. I've installed this properly but I still get the same error:
Could not find array definition file ' hu6800cdf.qcdef '. Simpleaffy does
not know the QC parameters for this array type.
See the package
I have attached to this message the first 20 lines of the output for
dput(amigo). Unfortunately, i can't send it the overall output.
Thank you.
On Wed, May 16, 2012 at 6:10 PM, David Winsemius [via R]
ml-node+s789695n4630253...@n4.nabble.com wrote:
On May 16, 2012, at 1:33 AM, ana24maria
The Quantitative Decision Strategies group at Janssen Research Development,
Johnson Johnson, is looking for a candidate to represent QDS in Beijing,
China in the subsidiary company of Xian-Janssen Pharmacetical Ltd. The basic
requirements for this candidate are 1) 3+ years experience in a
I've been having a problem using the step function to evaluate models. I've
simplified the code and get the same problem using the dataset Titanic. The
relevant code and output is below. The problem disappears (i.e., 'step' runs
correctly) if I rerun the code but change the 'loglm' call to
I using oldlogspline (from logspline package) to model data distributions, and
having a problem.
My data are search area sizes. They are based on circular search radii from
random points to the nearest edge of the nearest grass tussock. Search area
sizes are distributed from 0 (the random
Using both RMySQL or RODBC packages I can connect to a mySQL database and
return results from simple queries. However, I cannot return anything from
queries even with just two lines, each ending with semicolons. Using
RMySQL, the latter yields a message about an error in mysqlFetch, RS-DBI
Thank you, Michael :)
Michael Weylandt wrote
If that doesn't nail it down, I'll need you to answer the questions I
asked in my previous email.
Previously I made a mistake with *dput()*, this is the correct output:
dput(X)
new(timeSeries
, .Data = structure(c(124.3, 124.38, 124.67,
I have a rather large dataset which I just inserted hundreds of interaction
terms, and I guess R is saying there is a linear dependence when I run:
l1ce(y~., data=x, bound = .5, absolute.t = FALSE).
The error is below:
/Error in if (any(bound 0)) stop('bound'(s) must be non negative) :
On 17/05/2012 09:25, David A Vavra wrote:
I've been having a problem using the step function to evaluate models. I've
simplified the code and get the same problem using the dataset Titanic. The
relevant code and output is below. The problem disappears (i.e., 'step' runs
correctly) if I rerun the
Dear List,
Is there a way I can sort a sample based on a sort index constructed from
the data from which the sample is taken? Basically, I need to take 'many'
samples from the same source data and sort them. This can be very time
consuming for long vectors. Is there any way I can sort the data
Hello,
Apparently this is a follow-up from an earlier post.
I had answered but a misplaced comma in the subject line started another
thread, that you haven't read.
My original answer is in
Dear all,
I am trying to install the rJava package.
I am getting the following message
libjvm.so: cannot open shared object file: No such file or directory
at this web site I have found that (bottom part)
http://rwiki.sciviews.org/doku.php?id=packages:cran:rjava
that
Unix: if you encounter
I have trouble using function hmd.mx
I have all packages required.
I call function:
slovenia - hmd.mx(SVN, username, password, Slovenia) :
and I get this error:
NAs introduced by coercion
And if I then call slovenia, I get:
slovenia
Mortality data for Slovenia
Series: female male total
Dear All
I would need to perform a MANOVA with both fixed (group, sex, group*sex) and
random (brood) effects. I wonder if this is at all possible and if R does that.
At the moment, I only know that I can run a classic MANOVA with R.
Thank you
David
__
On May 17, 2012, at 2:20 AM, Damjan Krstajic wrote:
Thanks David for prompt reply. I agree with you. However, I still
fail to get the survfit function to work with newdata. In my
previous example I changed the column names of testX matrix and I
still fail.
On May 17, 2012, at 3:16 AM, Beutel, Terry S wrote:
I using oldlogspline (from logspline package) to model data
distributions, and having a problem.
My data are search area sizes. They are based on circular search
radii from random points to the nearest edge of the nearest grass
On May 17, 2012, at 2:26 AM, ana24maria wrote:
I have attached to this message the first 20 lines of the output for
dput(amigo).
If ypu attempted to send something, the mail-server scrubbed it.
Unfortunately, i can't send it the overall output.
Thank you.
I did not use the word attach.
On May 17, 2012, at 6:45 AM, Axel Urbiz wrote:
Dear List,
Is there a way I can sort a sample based on a sort index constructed
from
the data from which the sample is taken? Basically, I need to take
'many'
samples from the same source data and sort them. This can be very time
consuming
Something along the lines of
dat2 - ifelse( dat1==1 , yes, no)
Another option is in this case
dat2 - c(no, yes)[dat1+1]
Regards
Petr
should do it.
John Kane
Kingston ON Canada
-Original Message-
From: s1010...@student.hsleiden.nl
Sent: Mon, 14 May 2012
Hi Srod --
On 05/16/2012 11:04 PM, srod wrote:
Hi,
I'm using a command in bioconductor that seems to require a package called
hu6800cdf. I've installed this properly but I still get the same error:
Could not find array definition file ' hu6800cdf.qcdef '. Simpleaffy does
not know the QC
Hi Alaios,
The web page you're quoting was written in 2006. There *was* an R 2.4,
but it's no longer cutting edge, to say the least.
Earlier in the document it says R 2.4 or higher, which is what you've
got. I'd first look for newer documentation, then try the instructions
given with the current
On 17.05.2012 00:49, Rismyname wrote:
Hi,
I get the following error while installing a package. Can someone please
help?
install.packages(memisc)
Warning in install.packages :
argument 'lib' is missing: using 'C:/Users/ravi/Documents/R/R-2.15.0'
Warning in install.packages :
downloaded
The predict methods are designed to work with a dataframe. In the case
of survfit(cox model) the code has, for backwards compatability
purposes, the ability to use a vector as the newdata argument. This
feature is not documented in the help file, on purpose, and is expected
to go away one
Hi
On 2012-05-15 08:36, Melissa Rosenkranz wrote:
Here is an R problem I am struggling with:
My dataset is organized like this...
subject sessionvariable_x variable_y
01 11interger values
01 1
x[is.na(z)] - NA
This might send you a nasty bug if x and z are different lengths
though -- just a head's up.
Another option
x*!is.na(z)*z
Regards
Petr
Michael
On Wed, May 16, 2012 at 12:55 PM, Mintewab Bezabih
mintewab.beza...@economics.gu.se wrote:
Dear R users,
I was
Hi All,
I'm using the party package to grow conditional inference trees for survival
analysis.
When I used party version party_0.9-9991 everything worked well, but when I
update to party_1.0-2 (due to using 64bit R), I get an error. For simplicity
I will show the error I get for the example in
Hello, All:
fda_2.2.8 (functional data analysis) is now available on CRAN. This
revision includes the following improvements:
1. A bug in pca.fd has been fixed.
2. Many functions have a new argument returnMatrix, which if
TRUE allows the function to use
I am trying to import Indian National Sample Survey Data. It is ASCII flat
and an example is below
00146030602500111710107111201*01*00211 270104070204093
00246030602500111710107111201*02*00806104910519572 022 2600
11503055 170 4005
Hi
I want to use the standard error values in the summary that is produced
using
coeftest, but I am getting an error code- any ideas?
See what is structure of coeftest object by
str(coeftest(lmodT_WBHO))
and from this you shall deduct how to select second column.
Regards
Petr
i can get list of available packages by
r-available.packages()
No of available packages
length(r[,1])
[1] 3739
single package name
r[1,1]
[1] ACCLMA
i want download all packages starting with A / a alphabets,
with their source code, pdf manual/vignettes from internet
--
View this message
Hi..
Ok here is an example on how I wanted the tree to be implemented in R:
- the tree is, as you wrote, saved as a list of different tree levels
- each tree level is also saved as a list of different nodes in that
specific level
- and for the last part, each node is then saved as a list of
Hi
I had not see any answer yet but maybe there is nobody who wants to touch
the elusive object of outlier. Neither me, but here are some ideas how
one can proceed.
First of all its always up to you what is considered an outlier and how
will you deal with them.
I usually call an outlier any
Dear everyone! I am a new user in R. Recently, I have been struggling to work
out lagrange multiplier test in R, after read the description in R. I am
still confusing with method of assigning a neighbour list to listw.
Many thanks if anyone can give me a hint
--
View this message in context:
On Thu, 17 May 2012, ofraam wrote:
Hi All,
I'm using the party package to grow conditional inference trees for survival
analysis.
When I used party version party_0.9-9991 everything worked well, but when I
update to party_1.0-2 (due to using 64bit R), I get an error. For simplicity
I will
On 17.05.2012 14:17, Richard Iles wrote:
I am trying to import Indian National Sample Survey Data. It is ASCII flat
and an example is below
00146030602500111710107111201*01*00211 270104070204093
00246030602500111710107111201*02*00806104910519572 022 2600
115030
Thanks
Solution found
Fixed:
a. run as root the sudo R CMD javareconf
(not as simple user)
b. then install rJava from tar.gz
From: Sarah Goslee sarah.gos...@gmail.com
Cc: R help R-help@r-project.org
Sent: Thursday, May 17, 2012 2:44 PM
Subject: Re: [R]
This will accomplish what you want and should be relatively easy to modify.
# Create data.fram
names - c(S1, S2, S3, S4)
X - c(BB, AA, AB, AA)
Y - c(BB, BB, AB, AA)
Z - c(BB, BB, AA, NA)
AorB - c(A, A, B, B)
sample - data.frame(names, X, Y, Z, AorB,
stringsAsFactors=FALSE)
# Create recoded
On May 17, 2012, at 8:17 AM, Richard Iles wrote:
I am trying to import Indian National Sample Survey Data. It is
ASCII flat
and an example is below
00146030602500111710107111201*01*00211 270104070204093
00246030602500111710107111201*02*00806104910519572 022
2600
Hi All,
I have a situation where I want an 'if' variable to be parameterized. It's
entirely possible that the way I'm trying to do this is wrong, especially
given the error message I get that indicates I can't do this using an 'if'
statement.
Essentially, I have data where I think a
Thank you all for your help!! !
Alex
--
View this message in context:
http://r.789695.n4.nabble.com/kolmogorov-Smirnov-critical-values-tp4630245p4630393.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing
Hi,
unless you're dealing with heteroskedastic datas, the command *cor(x)* will
be enough, where *x* is your data matrix; in this function you can easily
select the method which has to be used: Pearson's, Kendall's or Spearman's
correlation.
--
View this message in context:
sagarnikam123 sagarnikam123 at gmail.com writes:
i can get list of available packages by
r-available.packages()
No of available packages
length(r[,1])
[1] 3739
single package name
r[1,1]
[1] ACCLMA
i want download all packages starting with A / a alphabets,
with their source
Thanks very much Achim!
I was indeed using pec which requires rms...
I see now that the old version of pec that I used didn't require rms and
therefore I had no problem.. I'm guessing I don't need the rms methods for
my use of pec, but am not sure how I can remove the dependency or
alternatively
On May 17, 2012, at 10:06 AM, DWatts wrote:
Hi All,
I have a situation where I want an 'if' variable to be
parameterized. It's
entirely possible that the way I'm trying to do this is wrong,
especially
given the error message I get that indicates I can't do this using
an 'if'
statement.
On May 17, 2012, at 11:08 AM, David Winsemius wrote:
On May 17, 2012, at 10:06 AM, DWatts wrote:
Hi All,
I have a situation where I want an 'if' variable to be
parameterized. It's
entirely possible that the way I'm trying to do this is wrong,
especially
given the error message I get
On May 17, 2012, at 11:04 AM, ofraam wrote:
Thanks very much Achim!
I was indeed using pec which requires rms...
I see now that the old version of pec that I used didn't require rms
and
therefore I had no problem.. I'm guessing I don't need the rms
methods for
my use of pec, but am not
Thanks David!!!
Sorry for the stupid question...
Seems like I still have issues with pec, but I will work on it some more...
--
View this message in context:
http://r.789695.n4.nabble.com/ctree-for-suvival-analysis-problem-tp4630362p4630408.html
Sent from the R help mailing list archive at
Petr et. al:
FWIW (probably not much).
As you know, tens of thousands of pages about outliers have been
written by statisticians. IMHO, it is another of the really terrible
ideas of our discipline and has led to much scientific abuse, as
indicated by this posting. For this reason, I have
Greetings again R users,
Some of you will likely recognize me but I hope you can help me once
more. I have tried the mixed model mailing list for this question but have
yet to find a solution. As such I hope someone will have another idea.
I have previously attempted to replicate the UN, CS,
Petr and Bert offer sound advice. At the risk of getting completely
ostracized, here's how you could find outliers using the definition of
'outlier' used by R's boxplot function and at the same time see your data.
dat = c(11489, 11008, 11873, 8000, 9558, 8645, 8024, 8371)
a =
David Winsemius wrote
and provide commented, minimal, self-contained, reproducible code.
Still untested in absence of data.
Thank you, and my apologies for not giving some sample numbers. Below is a
short subset of the data set I'm working with.
Could anyone can help me telling the way how I can build correlation matrix
in R? Thanks in advance.
--
View this message in context:
http://r.789695.n4.nabble.com/Correlation-Matrix-tp4630389.html
Sent from the R help mailing list archive at Nabble.com.
Thanks a lot. Suppose I want to use Pearson's method, then what I have to do?
--
View this message in context:
http://r.789695.n4.nabble.com/Correlation-Matrix-tp4630389p4630396.html
Sent from the R help mailing list archive at Nabble.com.
__
Hello,
I don't know if this is what you want, but your formula is equivalent to
nlstest=nls(Cfl~K*ifelse(Pw z, 1, exp(-a*Tsoil))*exp(-b*Pw)+c
,start=c(K=5.5, a=0.1, b=0.1, c=2, z=5))
Hope this helps,
Rui Barradas
--
View this message in context:
Hi Mahdi,
Look at the documentation for cor(), by typing ?cor or help(cor).
Pearson is the default and it is trivial to select the others. I
suggest you try searching google or reading R's documentation before
posting to the list. You may not understand it all, but it shows you
tried to work on
Josh:
A very nice, clear, polite, concise, and reasoned alternative to RTFM !
Probably should be templated somehow, given the volume of queries of
this sort that this list receives. (The posting guide is too involved
to serve in its stead.)
Cheers,
Bert
On Thu, May 17, 2012 at 9:31 AM,
This is a basic Windows system administrator problem, asked by a Linux
guy who is helping out in a Windows lab.
I want to keep R packages up to date on MS Windows 7 with a job in the
Task Scheduler. I have an R program that I can run (as
administrator) that updates the existing packages and then
On 2012-05-17 08:40, DWatts wrote:
David Winsemius wrote
and provide commented, minimal, self-contained, reproducible code.
Still untested in absence of data.
Thank you, and my apologies for not giving some sample numbers. Below is a
short subset of the data set I'm working with.
On Thu, May 17, 2012 at 06:45:52AM -0400, Axel Urbiz wrote:
Dear List,
Is there a way I can sort a sample based on a sort index constructed from
the data from which the sample is taken? Basically, I need to take 'many'
samples from the same source data and sort them. This can be very time
Type
str(slovenia)
or perhaps
names(slovenia)
to find out more about how the slovenia object is structured, and the
correct names of its elements.
Also
?hmd.mx
(also see,for example,
http://help.pop.psu.edu/help-by-software-package/r-project/demography-in-r/
demography-package-for-r)
Also,
Dominik,
There are a number of formulations of this statistic (see the
Kvålseth[*] reference below).
I tend to think of R^2 as the proportion of variance explained by the
model[**]. With the traditional formula, it is possible to get
negative proportions (if there are extreme outliers in the
Hi Greg,
The problem is that I also have restrictions for each variable (they must be
higher than -.07 and smaller than .2) and I'm dealing with a lot of them.
I've already tried the second approach but, as far as it seems, the function
doesn't satisfy my objective.
That's what I'm doing:
Hy Max,
thanks again for the answer.
I checked the caret implementation and you were right. If the
predictions for the model constant (or sd(pred)==0) then the
implementation returns a NA for the rSquare (in postResample). This is
mainly because the caret implementation uses `cor` (from the
Thanks.
I guess I still don't understand what's going on. It's not at all intuitive
that the table used should be in the search path. Why is it searching for
the table? Isn't the table already stored in the model? If the documentation
says this, I haven't found it. Needing to know what names I
On 05/17/2012 06:30 AM, jackl wrote:
Hi..
Ok here is an example on how I wanted the tree to be implemented in R:
- the tree is, as you wrote, saved as a list of different tree levels
- each tree level is also saved as a list of different nodes in that
specific level
- and for the last part,
Hi,
When I run the following code :
Y - c(rep(0,35),1,2,0,6,8,16,43)
cst - log(choose(42, 42:1))
beta - 42:1
tau - (beta^2)/2
fit - glm(formula = Y ~ offset(cst) + beta + tau, family = poisson)
fit
fit$converged
glm prints a warning saying that the algorithm did not converge.
However,
Hi all,
I need to solve this following equation which the unknown value (here is a)
is on the exponent place, could anybody please help me to figure it out? I
tried to use lm.sol() or solve() function but it doesn't work..
Thanks for help,
Peiling
(x, y) = (2.21, 1.01)
y = 1 + x
I used str(slovenia) and saw I was using wrong names - instead of year I
tried years and the same with age.
For hmd.mx you need package RCurl.
Thank you very much for your help!
Martina
--
View this message in context:
http://r.789695.n4.nabble.com/hmd-mx-tp4630350p4630439.html
Sent from the
You could probably find a way to do this using nls(), optim() or
uniroot() but this seems to be very much a homework problem so I won't
say more.
solve() solves matrix linear equations... not what you want here. Read
the docs.
Michael
On Thu, May 17, 2012 at 2:53 PM, Pei-Ling Lin
Hi Paul,
This is a bit OT, but here's what I would do.
1) write the R script (or if updating packages is all you want Rscript
-e update.packages(repos = 'yourrepo', ask = FALSE) would do it
without need for a script)
2) write a silly batch file (e.g., Rupdate.bat) wrapper (maybe there
are
On Thu, May 17, 2012 at 01:53:39PM -0500, Pei-Ling Lin wrote:
Hi all,
I need to solve this following equation which the unknown value (here is a)
is on the exponent place, could anybody please help me to figure it out? I
tried to use lm.sol() or solve() function but it doesn't work..
On 17-May-2012 21:50:17 Petr Savicky wrote:
On Thu, May 17, 2012 at 01:53:39PM -0500, Pei-Ling Lin wrote:
Hi all,
I need to solve this following equation which the unknown value (here is
a) is on the exponent place, could anybody please help me to figure it
out? I tried to use lm.sol() or
Peter Ehlers wrote
If the above data is all you have, I wouldn't chugg too long.
Looking at plots of Cfl~Pw and Cfl~Tsoil says to me that you
have far too much variability to fit your model.
You might find this plot instructive:
library(rgl)
plot3d(Pw, Tsoil, Cfl,
I have a very simple maximization problem where I'm solving for the vector
x:
objective function:
w'x = value to maximize
box constraints (for all elements of w):
low x high
equality constraint:
sum(x) = 1
But I get inconsistent results depending on what starting values I. I've
tried various
On May 18, 2012, at 00:14 , Nathan Stephens wrote:
I have a very simple maximization problem where I'm solving for the vector
x:
objective function:
w'x = value to maximize
box constraints (for all elements of w):
low x high
equality constraint:
sum(x) = 1
But I get
One of many scripts to produce 4 lattice plots on one page keeps throwing
an error. I've tried manipulating the file to eliminate the error, but have
not been able to do so. The error is:
source('bicarb.R')
Error in source(bicarb.R) : bicarb.R:15:1: unexpected symbol
14: 15: hco33
^
On Thu, 17 May 2012, Rich Shepard wrote:
Error in source(bicarb.R) : bicarb.R:15:1: unexpected symbol
14:
15: hco33
^
These two lines were concatenated in the message on the mail list.
Rich
--
Richard B. Shepard, Ph.D. | Integrity - Credibility - Innovation
Applied Ecosystem
On 12-05-17 04:34 PM, Rich Shepard wrote:
hco32 - qqmath(~ log10(HCO3 | factor(basin), data = surfchem.cast,
main = 'Bicarbonate (Log10)',
prepanel = prepanel.qqmathline,
panel = function(x, ...) {
panel.qqmathline(x, ...)
panel.qqmath(x, ...)
})
Missing a closing parenthesis
Your log10(HCO3 and sqrt(HCO3 seem to be missing closing brackets.
HTH,
baptiste
On 18 May 2012 11:34, Rich Shepard rshep...@appl-ecosys.com wrote:
One of many scripts to produce 4 lattice plots on one page keeps throwing
an error. I've tried manipulating the file to eliminate the error, but
On Thu, 17 May 2012, Mercier Eloi wrote:
Missing a closing parenthesis after log10.
Eloi,
A-ha! I knew new eyes would see what I kept missing.
Many thanks!
Rich
--
Richard B. Shepard, Ph.D. | Integrity - Credibility - Innovation
Applied Ecosystem Services, Inc. |Helping
Hi Sophie
It helps if you do some detective work
Try
fit1 - glm(formula = Y ~ offset(cst) + beta + tau, family =
poisson,trace = T, maxit = 200)
and compare
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home:
Dear all,
I have a function f - function(x,y,c){as.numeric(x*y c) }
I need to solve the value of c so that when I take the integral of
the function f from 0.2 to 0.8 with respect to x and from 0 to 1
with respect to y, the integral equal some prefixed value, say 0.025.
It involves two
Hello,
I'm attempting to build a package using R 2.15.0 on OS X
I am getting a generic failure when performing a cran type check on the
package. Even with a very simple test package, it still fails int he same
place.
Example:
In R:
rm(list=ls())
foo - function(x){print(x)}
Did you edit the description file and the namespace file and the Rd files?
Although my tutorial is for windows if you follow steps 8 thru 10 on mac it
should work
http://stevemosher.wordpress.com/step-8-the-sample-package/
On Thu, May 17, 2012 at 7:44 PM, Noah Silverman
Hi Noah,
Can you put your package on github or at least upload the tar ball?
Although I agree with Steven that as is, the package will not pass
cran checks, that is not the error I would expect.
Some of the experts may have seen this before or know instantly, but
for everyone else, being able to
Dear All,
I'm new at R, but I really just need a couple of things. The first thing I
need is to figure out how to get each individual financial statement
(CF,BS,IS). I need each individual one because getting them all at once
allows for formatting issues once it is a CSV. The date subsetting
Hi,
I recently installed Rattle for R 2.15.0 and all the functions work properly
except for the Correlation button. After I choose the Correlation radio button
in Explore tab, I click Execute, but nothing happens. No results at all. What
is it that I am missing?
Thank you,
Avideh
thanks a lot dear. I will keep your advice in my mind.
--
View this message in context:
http://r.789695.n4.nabble.com/Correlation-Matrix-tp4630389p4630448.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing
Dear R help forum members,
I am modeling a gaussian distribution for a computational biology application
and I am working in the statistical package R. In this regard, my problem is
that I have to construct a covariance matrix with variables (non-numeric) and
the covariance matrix is to be
Hi,
I want to find out, how can i get a SAS like output for Proc Autoreg with
AR(1) and AR(2) terms.
Thanks...
--
View this message in context:
http://r.789695.n4.nabble.com/Proc-AutoReg-SAS-like-Output-in-R-tp4630468.html
Sent from the R help mailing list archive at Nabble.com.
91 matches
Mail list logo