similar to: predict function the other way around

Hello, With this program I try to repeat analysis for different years. The results of the analysis are not printed when in the loop, except for the year sequence. What is wrong? Thanks a lot. for (i in 92:99){ cat("\n", "=============================================================\n", "YEAR =",i,"\n",

I am working on a script that takes numeric performance indicators and runs them against a series of regressors (dummy regressors, yes\no stuff via 0 and 1, e.g. Was is Christmas this week 0=no, 1=yes). The script is as follows (Written as a function): -- Begin Script -- doEnv <- function(HOUR,ENVNAME,REPORTNAME) { library(RODBC) library(forecast) library("geneplotter")

Dear list, I would like to generate random numbers from a Voigt distribution, hopefully in a way as simple as getting random numbers from a normal distribution with 'rnorm'. Is there any package to do this? Speed is an issue in this application. Or, as the Voigt distribution is a convolution of a Gaussian and a Lorentzian, can I simply combine random numbers from rnorm and rcauchy in some

Hello, Sweave seems to have trouble processing german letters in R. For example, my noweb R-input looks like this. <<>>= Oberfl?chenfehler = c(4, 11, 6, 2, 7, 9) @ If I send it through Sweave, I get the following error message. error: chunk 1 Error in parse(text = chunk) : unexpected input in "Oberfl??" extra: Warning message: In readLines(f[1]) : underfull last line in

to clean up some code I would like to make a list of arbitrary length to store?various objects for use in a loop sample code: ############ BEGIN SAMPLE ############## # You can see the need for a loop already linearModel1=lm(modelSource ~ .,mcReg) linearModel2=step(linearModel1) linearModel3=lm(modelSource ~ .-1,mcReg) linearModel4=step(linearModel3) #custom linearModel5=lm(modelSource ~ .

Dear expeRts, first of all I'd like to thank you for the quick help on my last which() problem. Here is another one I could not tackle: I have data on an absorption measurement which I want to fit with an voigt profile: fn.1 <- function(p){ for (i1 in ilong){ ff <- f[i1] ex[i1] <- exp(S*n*L*voigt(u,v,ff,p[1],p[2],p[3])[[1]]) } sum((t-ex)^2) } out <-

Hello everyone, I am trying to do within subjects repeated measures anova followed by the test of sphericity (sample dataset below). I am able to get either mixed model or linear model anova and TukeyHSD, but have no luck with Repeated-Measures Assuming Sphericity or Separate Sphericity Tests. I am trying to follow example from "car" package, but it seems that I am not getting something

Dear list, Latex/Sweave has trouble processing Sveave-output coming from the summary-command of a linear Model. >summary(lmRub) The output line causing the trouble looks in R like this Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 In my Sweaved Tex-file that line looks like this Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1? (actually in the editor the

Why are my coefficients getting appended with a 1? It borks a match I do later against the original list that doesn't have the random 1 added to the end. > linearModel[[1]] Call: lm(formula = modelSource ~ +UNITBUILD + UNITDB + ITBUILD + ITDB + UATBUILD + UATDB + HOGANCODE + RCF + ReleaseST1 + ReleaseST2 + ReleaseBLA + Small.Bank.Acquisitions + HLY.NewYear + HLY.MLK + HLY.PRES +

[corrected dataset below] Hello everyone, I am trying to do within subjects repeated measures anova followed by the test of sphericity (sample dataset below). I am able to get either mixed model or linear model anova and TukeyHSD, but have no luck with Repeated-Measures Assuming Sphericity or Separate Sphericity Tests. I am trying to follow example from "car" package, but it seems

Hi, how may I define in a callfile the CallerID presentation to be used for the requested call, eg. set it to prohibited? TIA, Bruno -------------- next part -------------- A non-text attachment was scrubbed... Name: Bruno.Voigt.vcf Type: text/x-vcard Size: 270 bytes Desc: not available Url : http://lists.digium.com/pipermail/asterisk-users/attachments/20050913/fcb5c595/Bruno.Voigt.vcf

Hi all, I'm looking for a way with any asterisk-version with TE410P (cpe EuroISDN, Q931) for sending an INFORMATION ELEMENT KeypadFacility, eg. *87, during a connected call to the PSTN switch. Are there existing functions in asterisk to generate & send such IE ? If not what existing modules would be best to derive from? TIA, Bruno -------------- next part -------------- A non-text

Dear all, I freshly installed a new smb-server that should act as file server with some smb-shares. The pdc is a Windows 2008 r2 server. My smb.conf is as follows: [global] workgroup = agroup server string = ssmbserver #passdb backend = tdbsam map to guest = Bad User usershare allow guests = No security = domain wins support = No

Hi, I am sorry for this question, but I am trying to speed up an application.... I will need to fit many x-y data sets (input from text files) to 4-parameter Pseudo-Voigt peak function. Until now I used SigmaPlot macro to do it (enclosed just in case...) peaksign(q) = if(total(q)>q[1], 1, -1) xatymin(q,r) = xatymax(q,max(r)-r) [Parameters] a = if(peaksign(y)>0, max(y), min(y))

Hi, I have a problem with the filesize of the samba 3.0.0 package. after compiling the source of the samba-3.0.0 i made the "sh makepkg.sh". the package i got was about 901 MB. does anybody know why the package is as great as mine and what can i do to get a samba package which is smaller than the package i have. i compiled under solaris 8. i made all what samba.org says that i

Hi Deepak, In econometrics there is another test very often used : the white test. The white test is based on the comparison of the estimated variances of residuals when the model is estimated by OLS under the assumption of homoscedasticity and when the model is estimated by OLS under the assumption of heteroscedastic. The White test with R install.packages("bstats") library(bstats)

Hi, Anne, assign Age and Bloodpressure in the correct order to the axes in your call to plot as in: plot(y = Age, x = BloodPressure) abline(SimpleLinearReg1) Hth -- Gerrit --------------------------------------------------------------------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eichner at math.uni-giessen.de Justus-Liebig-University Giessen Tel:

Dear Gerrit, This issue is discussed in a vignette in the car package (both for functions in the car and effects packages): vignette("embedding", package="car") . The solution suggested there is the essentially the one that you used. I hope this helps, John ----------------------------- John Fox, Professor Emeritus McMaster University Hamilton, Ontario, Canada Web:

>> Do you (or the ogg vorbis community) understand how SBR works? There is >> info on http://www.codingtechnologies.de/technology/sbr.htm >> bnut I am guessing this infop is not enough to explain the technology. Robert Voigt: >I haven't heard about SBR before. After reading that webpage I can say the >following: [ . . . ] I don't think SBR will give an

We are looking to implement a code review system that will require a review of the code before it gets to production servers. I have been looking into many different options, and would love some feedback and experiences from other users. Our puppet is in a git repo. We have 3 primary environments with a branch of puppet for each of them. All of our puppet info is in the git repo