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 explicit
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 point
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
drive
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") :
mi
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 onl
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
http://r.789695.n4.nabble.com/Re-Problem-to-resolve-a-step-for-reading-a-large-TXT-and-split-in-several-
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 t
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.
colnames(testX)<-names(coxph
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
tussock.
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 fo
> 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:
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 para
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 11
> > 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
> wrote:
> > Dear R users,
> >
> > I was wondering how I ca
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 t
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 t
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 003460306025001117101071112
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
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 func
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 an
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:
htt
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 show
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
1150305
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
Cc: R help
Sent: Thursday, May 17, 2012 2:44 PM
Subject: Re: [R] Is there R 2.4 version_
Hi Alaios,
The we
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,
s
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
1150
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 relationship
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 list
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:
http://r.789695.n4.na
sagarnikam123 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 thei
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
alternativ
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 tha
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 no
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 Nabbl
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 elimin
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, an
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 = boxplo
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.
Cfl=c(2.61,4.21,2,2.75,7.47,1.2,3.2
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.
__
R-help@r
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:
http://r.789695.n4.nabble.com/nls-a
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 o
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, Josh
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 the
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.
Cfl=c
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 tim
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
pr
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
stats-
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 c
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, ea
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, fit$
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 +
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 R
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 wrote:
> Hi al
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 bett
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 l
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, Cf
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 variou
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
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
^
T
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 Ecosyste
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 af
Your log10(HCO3 and sqrt(HCO3 seem to be missing closing brackets.
HTH,
baptiste
On 18 May 2012 11:34, Rich Shepard 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 have
> not been able t
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: mac...
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 dim
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)}
package.skeleton(
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 wrote:
> Hello,
>
> I
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 is
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
[[al
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 lis
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 us
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.
_
On Thu, May 17, 2012 at 09:43:13PM -0400, li li wrote:
> 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 int
89 matches
Mail list logo
