similar to: Lapack Routine Error

Hi there! I´m a noob when it comes to R and I´m using it to run statisc analysis. With the code for ARIMA below I´m getting this error: Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular The code is: > s.ts <- ts(x[,7], start = 2004, fre=12) > get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) + { + best.aic <- 1e8 + n <-

Hi I have a problem with this error, I have searched the archives and found previous discussion about this, can I cannot understand how the explanations apply to what I am trying to do. I am trying to do Log_rank Survival analysis, I have included tables and str command, is it a factor/integer problem? If so how do I correct this, as all my attempt to recode the data have failed. >

Dear colleagues: Despite my best efforts, I have not been able to understand/overcome an error message I received while running the fclustIndex command in package e1071, which produces validity measures for the fuzzy c-means clustering algorithm. I am relatively new to R-- and still depend heavily on Rcmdr-- and so I apologize in advance if there is an obvious answer to my question:

Hi R-users: When I run the R script (as the following), I got the error message: Error in solve.default(C[ind, ind]) : Lapack routine dgesv: system is exactly singular. Any help is appreciated. Ying library(sem) R.pw <- matrix(c( 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2137356, 0.2137356, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

Hello I am trying to evaluate an Incomplete gamma function in R. Library Zipfr gives the Igamma function. From Mathematica, I have: "Gamma[a, z] is the incomplete gamma function." In[16]: Gamma[9,11.1] Out[16]: 9000.5 Trying the same in R, I get > Igamma(9,11.1) [1] 31319.5 OR > Igamma(11.1,9) [1] 1300998 I know I have to understand the theory and the math behind it rather

#Sorry for the convoluted subject line. #I have: a=c(1,2,3) x=a #example of user supplied input #Is there any function that will tell me the "name" of the object x refers to, referring only to x itself? #i.e. the "answer" I want is "a" #I want: #fun(x) == 'a' #(I don't think this is possible, but figured I'd ask.)

Hi, I need to run a script under the variable (that comes from a csv file) that is assigned to another variable. This is very a simplified version of what I want to do: data<-read.csv('name.csv') names(data)<-("VA","VB","VC") v<-VA mn(v) Thank you in advance, Judith

Hi, everyone, I am green with R. Now I am working on a graphics projects. I have a dataset: ----------------------- a b c 1 228 83 6.69 2 274 83 7.36 3 320 83 8.86 4 366 83 7.36 5 412 83 6.81 6 228 129 5.58 --------------------- I want to draw a 3d graphic. I use the code __________________

Hi all Can I use ssh to have remote tar files from machine A to machine B? ssh from machine A to machine B tar all files in machine B to exact to machine A Thank you ____________________________________________________________________________________ Got a little couch potato? Check out fun summer activities for kids.

Hello, I took a stab at packaging RequestTracker for CentOS 5. Basically some needed FC6 packages were rebuilt as el5. The results are here: http://www.tlviewer.org/rt3/ I'm an rt3 newbie. I'm in the process of reading the DOCS and setting up the DB. The install via RPM appears to be sound. So far the login page for RT opens successfully. I've yet to setup the mail-dispatcher. --

Every time I go to the puppet website, I get the following when I try to click on any of the links, ie: http://www.reductivelabs.com/trac/puppet/wiki/SimplestPuppetInstallRecipe Thanks! ____________________________________________________________________________________ Got a little couch potato? Check out fun summer activities for kids.

Dear all I have setup with mediant 2000 with avaya now i want to install E1/PRI card with asterisk and trunk with E1 with Avaya E1 port so i want to buy E1 card for asterisk so which card is best and cast effective for my setup i want 1 port E1 card so can you suggest me which card is best for my setup and i want QSIG singaling with avaya Regards satish patel

Hi Friends, I am planning to buy "IMate PDAL" mobile phone. This contains "Wi-Fi 802.11b/g" feature. So, Is it possible to get internet using my wireless router in my office? Look forward to your response. Thank you. Regards, Chandra. --------------------------------- Got a little couch potato? Check out fun summer activities for kids. -------------- next part

Hello Asterisk-Users, I'm trying to install addons 1.2.6 on Asterisk 1.2.16 (is that OK?), but my MySQL server is installed on a different sever, so the MAKE of the addons fails with the following (truncated) error message: "app_addon_sql_mysql.c:23:19: mysql.h: No such file or directory". Is there any way to bypass/ignore the fact that MySQL is installed separately and enable the

I'm confused about something.... In Asterisk 1.4 you can collect RTP QoS metrics at the end of a call with: ${CHANNEL(rtpqos,audio,all)} Now, when your using the AMI to do a callout, like this... ACTION: Originate Async: yes Timeout: 60000 Exten: callback Channel: SIP/1000 Variable: callid=849120 Variable: destination=SIP/1001 Variable: timeout=70000 Variable: timeout_warning=60000

is it possible to run delphi apps/forms(that was build, compiled, and targeted to windows machines) on linux(fc6, suse 10.x) These delphi apps connects to MS SQL server. if it is possible, then I can easily ground all those windows xp machines, that are just use to run these delphi apps... and use linux instead Regards --------------------------------- Got a little couch potato? Check

--- Jan Engelhardt <jengelh@computergmbh.de> wrote: > Hi, > > > there is an interesting case when the FLAC encoder (using 1.2.0) is > given simple waves. Simple waves means: I have a list of {frequency, > duration, pause} tuples that define the monophonic tune. In other > words, exactly one frequency is played at a time. > > This is the original dataset from

Hi, I was not able to find a way such that R CMD build run Sweave with eval=FALSE. I have a package where I want the R code in the vignettes to run only during R CMD check but not on R CMD build. Is there any way to achieve this ? thank you. ____________________________________________________________________________________Got a little couch potato? Check out fun summer activities

Hope these help for alternatives to lm()? I show the use of a 2nd order polynomial as an example to generalize a bit. Sometimes from the subject line two separate responses can appear as reposts when in fact they are not... (though there are identical reposts too). I should probably figure a way around that. --- Stephen Tucker <brown_emu at yahoo.com> wrote: > ## Data input > input

Dear Friends, My problem is related to how to measure probabilities from a probit model by changing one independent variable keeping the others constant. A simple toy example is like this Range for my variables is defined as follows y=0 or 1, x1 = -10 to 10, x2=-40 to 100, x3 = -5 to 5 Model output <- glim(y ~ x1+x2+x3 -1, family=binomial(link="probit")) outcoef <-