similar to: Problem plotting output from tree()

?? If I haven't misunderstood, they are completely different! 1) NIR must be a matrix, or poly(NIR,...) will fail. 2) Due to the previously identified bug in poly, degree must be explicitly given as poly(NIR, degree =2,raw = TRUE). Now consider the following example: > df <-matrix(runif(60),ncol=3) > y <- runif(20) > mdl1 <-lm(y~df*I(df^2)) > mdl2

Dear All, I would like to fit some parametric survival models using left truncated, right censored data in R. However I am having problems finding a function to fit parametric survival models which can handle left truncated data. I have tested both the survreg function in package survival: fit1 <- survreg(Surv(start, stop, status) ~ X + Y + Z, data=data1) and the psm function in package

I am hoping someone can give me a couple of pointers on how to pass arguments using "..."--in this specific case in an attempt to plot an object created by ca(), a tool for correspondence analysis. #some illustrative code library(ca) set.seed(123) dat1 <- data.frame(matrix(ceiling(runif(30, 1, 100)), nc=5)) #some fake data ca1 <- ca(dat1) #the correspondence analysis plot(ca1)

Hi, Before I reinvent the wheel I wanted to kindly ask you for your opinion if there is a simple way to do it. I want to merge a larger number of tables from different data sources in R and the matching criterium are country names. The tables are of different size and sometimes the country names do differ slightly. Has anyone done this or any recommendation on what commands I

Hi all, I have two questions 1 - I have the version 1.4.1 of R, and it doesn't have the 'agrep' function in the base library. Is there a way to make this funcion avaliable in R 1.4.1? I mean, how to 'copy' it from R 1.8.1 and 'paste' it in R 1.4.1? 2 - The AGREP function doesn't give me the Levenshtein distance (edit distance). Is there a function in R that does

Hi all, I am having a problem trying to declare a slot representation to hold a file connection. With V2.8.0 I had been using: setClass("Element", representation( [other slots removed], datafileConn = setOldClass(c("file","connection")) ) ) This resulted in a warning but still ran okay. No however with V2.9.0 I am getting

After installing a new version of R, how do I make the packages that I have already installed and use on the old version available to the new version? Thank you, Roy Robertson

Michael Robertson, founder of both MP3.com and Lindows, has launched a new company to supply inexpensive SIP phones ($129 for two) and related services. See today's press release at http://www.sipphone.com/tiki-index.php?page=SIPphone%20Inc My question for the list is who will be the first to report on the compatibility and usability of the SIPphone with Asterisk? The functionality

> On Jul 13, 2017, at 7:43 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote: > > Below. > > -- Bert > Bert Gunter > > > > On Thu, Jul 13, 2017 at 3:07 AM, Luigi Biagini <luigi.biagini at gmail.com> wrote: >> I have two ideas about it. >> >> 1- >> i) Entering variables in quadratic form is done with the command I >>

Hi all, We have been fighting with telco for a entire week. Today they came here with a LITE3000 to analyze what is going on. When I configure zaptel with no external clock, E1 gets aligned/synchronized with bit rate in 2048000 bps, both me and telco. span=4,0,0,ccs,hdb3,crc4 But when I configure span4 to get clock source from telco they become unsynchronized. TElco bit rate stays in

Hi all, Prefix: I am a frustrated Java coder in R. I am coding a medium sized ecosystem modelling program in R. I have changed to using S4 objects and it has cost me an order of magnitude in execution speed over the functional model. I cannot afford this penalty and have found that it is the result of all the passing-by-value of objects. I see that you can now safely inherit from

Dear R users, I am using lmer combined with AIC model selection and averaging (in the MuMIn package) to try and assess how isotope values (which indicate diet) vary within a population of animals. I have multiple measures from individuals (variable 'Tattoo') and multiple individuals within social groups within 4 locations (A, B, C ,D) crucially I am interested if there are

I have a customer running AIX 6.1 and using pware Samba 3.0.28, and they are looking to install a SMB2/3 compatible version on their AIX server. I found this link https://wiki.samba.org/index.php/AIX_building_samba and understand that it requires xlC 13.1 and a lot of RPMs. Since this customer does not have a current AIX support contract I think it's unlikely that we'll be able to get the

Dear all, I have a basic(!) econometric question which i couldnt find the way to do it in R. Well this could be also because of my wrong interpretation of the econometric process that i am trying to implemet.so here i wanna ask if am doing a logical mistake!!! so here is the question with the explanation of the process, hope there will be someone who can help me! suppose i have a basic

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I run FreeBSD-4.5-STABLE, along with icecast 1.3.12, and I've been having two problems with ices 0.2.2. 1) When I give ices an empty playlist file and tell it to randomize, I get a segfault every time. When I give it the same playlist with Randomize set to 0, it tells me that the playlist is empty instead of crashing. 2) ices can't

A customer contacted us today concerning getting a VoIP to PSTN system with a few IP Phones setup. Asterisk should fit his needs. It is not a big job, but I think that this customer is going to need onsite work. Please contact me off list if you are an interested reseller in the Washington, DC area. Sean _______________________________________________ Sean Robertson NETXUSA p. 800-289-6389

Dear List, I would like to calculate the Goodman-Kruskal gamma for the predicted classes obtained from an ordinal regression model using lrm in the Design package. I couldn't find a way to get gamma for predicted values in Design so have found previous positings suggesting to use : Rcorr.cens(x, S outx = TRUE) in the Hmisc package My question is, will this work for predicted vs

Dear R users, Recent changes to the MuMIn package now means that the model averaging command (model.avg) no longer returns confidence intervals, but instead returns zvalues and corresponding pvalues for fixed effects included in models. Previously I have used this package for model selection/averaging following Greuber et al (2011) where it suggests that one should use confidence intervals from

Dear List, I am trying to assess the prediction accuracy of an ordinal model fit with LRM in the Design package. I used predict.lrm to predict on an independent dataset and am now attempting to assess the accuracy of these predictions. >From what I have read, the AUC is good for this because it is threshold independent. I obtained the AUC for the fit model output from the c score (c =

I'm doing a text mining project where I have to manually enter a double quote as an element inside a vector. I tried char[10]='"'#where i enclosed the double quote in a pair of single quotes. But the result is [1] "\"". Somehow a back slash is added automatically. I also tried to enclose the double quote in a pair of double quotes. That didn't work either.