Re: [R] survival survfit with newdata

2012-05-17 Thread Damjan Krstajic
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-

[R] hu6800cdf

2012-05-17 Thread srod
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

Re: [R] Error on easy way for JoSAE Package

2012-05-17 Thread ana24maria
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

[R] Job opportunity in Beijing, China at Xian-Janssen Pharmaceutical Ltd

2012-05-17 Thread Fei Chen
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

[R] step function stops with subscript out of bounds

2012-05-17 Thread David A Vavra
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

[R] oldlogspline probabilities

2012-05-17 Thread Beutel, Terry S
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

[R] trouble with multi statement queries

2012-05-17 Thread tensegrity
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

Re: [R] How to sum and group data by DATE in data frame

2012-05-17 Thread Cren
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,

[R] Error in lasso regression (linear dependence issue?)

2012-05-17 Thread YaoPau
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) :

Re: [R] step function stops with subscript out of bounds

2012-05-17 Thread Prof Brian Ripley
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

[R] Complex sort problem

2012-05-17 Thread Axel Urbiz
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

Re: [R] how disable the Error massage in read.table() no lines available in input

2012-05-17 Thread Rui Barradas
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

[R] Is there R 2.4 version_

2012-05-17 Thread Alaios
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

[R] hmd.mx

2012-05-17 Thread krmartina
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

[R] MANOVA with random factor

2012-05-17 Thread David Costantini
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 __

Re: [R] survival survfit with newdata

2012-05-17 Thread David Winsemius
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.

Re: [R] oldlogspline probabilities

2012-05-17 Thread David Winsemius
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

Re: [R] Error on easy way for JoSAE Package

2012-05-17 Thread David Winsemius
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.

Re: [R] Complex sort problem

2012-05-17 Thread David Winsemius
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

Re: [R] Add column from other columns data.

2012-05-17 Thread Petr PIKAL
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

Re: [R] hu6800cdf

2012-05-17 Thread Martin Morgan
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

Re: [R] Is there R 2.4 version_

2012-05-17 Thread Sarah Goslee
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

Re: [R] Unable to install package

2012-05-17 Thread Uwe Ligges
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

Re: [R] (bug?) in survfit in R-2-15.0

2012-05-17 Thread Terry Therneau
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

Re: [R] max value

2012-05-17 Thread Petr PIKAL
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

Re: [R] replacing with NA

2012-05-17 Thread Petr PIKAL
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

[R] ctree for suvival analysis problem

2012-05-17 Thread ofraam
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

[R] [R-pkgs] new version of fda fixes bug in pca.fd and supports library(Matrix)

2012-05-17 Thread Spencer Graves
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

[R] Importing ASCII flat

2012-05-17 Thread Richard Iles
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

Re: [R] error code trying to extract second column from coeftest output

2012-05-17 Thread Petr PIKAL
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

[R] how to download source code, pdf manual for any package from internet?

2012-05-17 Thread sagarnikam123
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

Re: [R] as.function parameters

2012-05-17 Thread jackl
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

Re: [R] how to find outliers from the list of values

2012-05-17 Thread Petr PIKAL
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

[R] neighbour list in lm.LMtests

2012-05-17 Thread johnnyzn
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:

Re: [R] ctree for suvival analysis problem

2012-05-17 Thread Achim Zeileis
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

Re: [R] Importing ASCII flat

2012-05-17 Thread Uwe Ligges
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

Re: [R] Is there R 2.4 version_

2012-05-17 Thread Alaios
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]

Re: [R] Help needed for efficient way to loop through rows and columns

2012-05-17 Thread David L Carlson
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

Re: [R] Importing ASCII flat

2012-05-17 Thread David Winsemius
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

[R] nls and if statements

2012-05-17 Thread DWatts
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

Re: [R] kolmogorov-Smirnov critical values

2012-05-17 Thread aramos
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

Re: [R] Correlation Matrix

2012-05-17 Thread Cren
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:

Re: [R] how to download source code, pdf manual for any package from internet?

2012-05-17 Thread Ben Bolker
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

Re: [R] ctree for suvival analysis problem

2012-05-17 Thread ofraam
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

Re: [R] nls and if statements

2012-05-17 Thread David Winsemius
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.

Re: [R] nls and if statements

2012-05-17 Thread David Winsemius
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

Re: [R] ctree for suvival analysis problem

2012-05-17 Thread David Winsemius
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

Re: [R] ctree for suvival analysis problem

2012-05-17 Thread ofraam
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

Re: [R] how to find outliers from the list of values

2012-05-17 Thread Bert Gunter
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

[R] Covariance structure for lme

2012-05-17 Thread Charles Determan Jr
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,

Re: [R] how to find outliers from the list of values

2012-05-17 Thread Robert Baer
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 =

Re: [R] nls and if statements

2012-05-17 Thread DWatts
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.

[R] Correlation Matrix

2012-05-17 Thread mahdi
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.

Re: [R] Correlation Matrix

2012-05-17 Thread mahdi
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. __

Re: [R] nls and if statements

2012-05-17 Thread Rui Barradas
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:

Re: [R] Correlation Matrix

2012-05-17 Thread Joshua Wiley
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

Re: [R] Correlation Matrix

2012-05-17 Thread Bert Gunter
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,

[R] Windows Task Scheduler and R updates. Need basic tips

2012-05-17 Thread Paul Johnson
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

Re: [R] nls and if statements

2012-05-17 Thread Peter Ehlers
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.

Re: [R] Complex sort problem

2012-05-17 Thread Petr Savicky
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

Re: [R] hmd.mx

2012-05-17 Thread MacQueen, Don
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,

Re: [R] caret: Error when using rpart and CV != LOOCV

2012-05-17 Thread Max Kuhn
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

Re: [R] Optimization problem

2012-05-17 Thread Pacin Al
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:

Re: [R] caret: Error when using rpart and CV != LOOCV

2012-05-17 Thread Dominik Bruhn
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

Re: [R] step function stops with subscript out of bounds

2012-05-17 Thread David A Vavra
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

Re: [R] as.function parameters

2012-05-17 Thread Jason Edgecombe
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,

[R] glm convergence warning

2012-05-17 Thread Sophie Baillargeon
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,

[R] solve the equation in R

2012-05-17 Thread Pei-Ling Lin
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

Re: [R] hmd.mx

2012-05-17 Thread krmartina
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

Re: [R] solve the equation in R

2012-05-17 Thread R. Michael Weylandt
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

Re: [R] Windows Task Scheduler and R updates. Need basic tips

2012-05-17 Thread Joshua Wiley
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

Re: [R] solve the equation in R

2012-05-17 Thread Petr Savicky
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..

Re: [R] solve the equation in R

2012-05-17 Thread Ted Harding
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

Re: [R] nls and if statements

2012-05-17 Thread DWatts
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,

[R] Optimization inconsistencies

2012-05-17 Thread Nathan Stephens
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

Re: [R] Optimization inconsistencies

2012-05-17 Thread peter dalgaard
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

[R] New Eyes Needed to See Syntax Error

2012-05-17 Thread Rich Shepard
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 ^

Re: [R] New Eyes Needed to See Syntax Error

2012-05-17 Thread Rich Shepard
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

Re: [R] New Eyes Needed to See Syntax Error

2012-05-17 Thread Mercier Eloi
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

Re: [R] New Eyes Needed to See Syntax Error

2012-05-17 Thread baptiste auguie
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

Re: [R] New Eyes Needed to See Syntax Error

2012-05-17 Thread Rich Shepard
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

Re: [R] glm convergence warning

2012-05-17 Thread Duncan Mackay
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:

[R] Integration of two dimension function

2012-05-17 Thread li li
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

[R] Failure building any package

2012-05-17 Thread Noah Silverman
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)}

Re: [R] Failure building any package

2012-05-17 Thread steven mosher
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

Re: [R] Failure building any package

2012-05-17 Thread Joshua Wiley
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

[R] Financial Statements Date Subsetting

2012-05-17 Thread Brian Edmundson
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

[R] Correlation in Rattle

2012-05-17 Thread avideh yesharim
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

Re: [R] Correlation Matrix

2012-05-17 Thread mahdi
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

[R] Covariance matrix in R with non-numeric variables

2012-05-17 Thread nataraj
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

[R] Proc AutoReg (SAS like Output in R)

2012-05-17 Thread anil
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.