similar to: Coversion from yearly to weekly data

Dear Rusers, I have met a difficult problem on explaining the differences of principal component analysis(PCA) between R,S-PLUS and SAS/STATA/SPSS, which wasn't met before. Althought they have got the same eigenvalues, their coeffiecients were different. First, I list my results from R,S-PLUS and SAS/STATA/SPSS, and then show the original dataset, hoping sb. to try and explain it.

Thank you so much Dr. Nash, I truly appreciate your kind and valuable contribution. Cheers, Paul El El s?b, 19 de ago. de 2023 a la(s) 3:35 p. m., J C Nash < profjcnash at gmail.com> escribi?: > Why bother. nlsr can find a solution from very crude start. > > Mixture <- c(17, 14, 5, 1, 11, 2, 16, 7, 19, 23, 20, 6, 13, 21, 3, 18, 15, > 26, 8, 22) > x1 <- c(69.98, 72.5,

Dear friends, Hope you are all doing well and having a great weekend. I have data that was collected on specific gravity and spectrophotometer analysis for 26 mixtures of NG (nitroglycerine), TA (triacetin), and 2 NDPA (2 - nitrodiphenylamine). In the dataset, x1 = %NG, x2 = %TA, and x3 = %2 NDPA. The response variable is the specific gravity, and the rest of the variables are the predictors.

The cautions people have given about starting values are worth heeding. That nlxb() does well in many cases is useful, but not foolproof. And John Fox has shown that the problem can be tackled very simply too. Best, JN On 2023-08-19 18:42, Paul Bernal wrote: > Thank you so much Dr. Nash, I truly appreciate your kind and valuable contribution. > > Cheers, > Paul > > El El

Hi, Since we upgraded to 2.0.9 (from 1.10 stock CentOS release), we are getting some errors with pop3. When the machines get busy, now and then it start with the following: > Mar 7 11:19:51 xxx dovecot: pop3-login: Error: net_connect_unix(pop3) failed: Resource temporarily unavailable And it generates hundreds of those before the machines dies, with the web server getting stuck as well on

I am just starting to learn Sweave (really neat tool). I am pretty early in the learning curve (I had to think a moment ago whether a # or % was the appropriate comment character). I have successfully incorporated simple graphics and outputs, but am having trouble getting a latex (xtable) table to function properly. Latex is seemingly treating the xtable code as input or verbatim text. That

Hi, I think that it would be a good thing to know more about those 2 projects (and also the future patent free format). I think that many people as me know about Lame, but not about Vorbis, and vice-versa. It would be fine that someone (perhaps the maintainer) of every project would introduce to both group of people those projects. 2 things would be interesting (to my mind): - to know about the

Hi, How can plot a frequency distribution curve for the following data.    V1      V2 1   1 160.54% 2   1 201.59% 3   1  18.45% 4   1 179.03% 5   1 274.37% 6   1   0.00% 7   1  24.52% 8   1  39.17% 9   3  43.72% 10  1  53.06% 11  1  64.97% 12  1  79.84% 13  1  98.08% 14  1 115.32% 15  1 127.96% 16  1 155.38% 17  1 157.25% 18  1 193.17% 19  1  51.53% 20 15  99.32% 21  1 106.86% 22  1 219.44%

Hi, A simple quantile question: I need to calculate the 95% and 5% quantiles (aka percentiles) for the following data: 67.12 64.51 62.06 55.45 51.41 43.78 10.74 10.14 if I use the formula: 95% quantile point= 95 (8+1)/100, I get the 8.55th point as the 95% quantile. Which does not make too much sense as I have only 8 data points. The other option is to use (95*8)/100 = 7.6th data point (which can

Hi list Maybe someone can help with the following problem (thanks in advance): I need to analisy a long list of data time \fluoscence I want to make one plot with standard deviation and the average data... Could you suggest me how to make that.. The data have this distribution: How can read and extrapolate this data? pippo.csv ps <- read.table(pippo.csv) time fluorescence 0; 0 4.013; 62.96

I'm trying to run the generalized least square approach for my 2-additive problem, unfortunately this error appeared. I have tried to figure out the error from the mailing list but couldn't find the solution. Any help is highly appreciated. This is my source code: >a1 <- c(76.18, 61.84, 60.4, 69.09) > a2 <- c(51.01, 50.39, 87.62, 52.03) > a3 <- c(80.08, 48.49, 90.86,

Dear R'ers, When I used S-plus i wrote a small program for a Mantel-Haenszel test for trend (I think it worked). Unfortunately I can't get it working in R. It appears as if my use of 'el' is the problem but I can't sort it out. Error in apply(array, c(, 2, 3), function(el) el * 1:s) : argument is missing, with no default Further down in the program I use 'el'

Hi, I've come across a strange error when using the lrm.fit function and the subsequent predict function. The model is created very quickly and can be verified by printing it on the console. Everything looks good. (In fact, the performance measures are rather nice.) Then, I want to use the model to predict some values. I get the following error: "fit was not created by a Design

Hi: I am working on a publication and I have heard about LaTEX but I haven't actually tried to learn about it until today. I've found a few examples but I can't actually make them work properly. I have a couple of questions: Does LATEX have to be installed on your computer? How does the xtable package and Sweave work together? How can I make the code below create a table as pdf file?

Hello, I have found a problem with using the glm function with a gamma family. I have a vector of data, assumed to be generated by a gamma distribution. The parameters of this gamma distribution are estimated in two ways (i) using the glm() function, (ii) "by hand", using the optim() function. I find that the -2*likelihood at the maximum found by (i) is substantially larger than that

Hi, i have that table Thesis Day A B C 1 0 83.43 90.15 22.97 1 0 85.50 94.97 16.62 1 0 83.36 95.38 20.70 1 0 84.47 92.16 23.58 1 0 83.98 95.33 19.39 1 0 82.86 93.78 24.55 1 0 83.39 92.67 19.56 1 0 85.17 95.24 17.95 1 0 81.62 93.32 28.49 1 0 82.99 92.85 19.73 1 0 81.11 95.67 27.20 1 0 83.39 94.69 16.51 1 0 79.56 89.87 30.39 1 0 80.54 93.32 21.76 1 0 82.11 92.58 22.17 1 14 85.65 94.00 19.19 1 14

Hi, I am trying to construct a logistic regression model from my data (104 patients and 25 events). I build a full model consisting of five predictors with the use of penalization by rms package (lrm, pentrace etc) because of events per variable issue. Then, I tried to approximate the full model by step-down technique predicting L from all of the componet variables using ordinary least squares

I''ve seen the problems with bug 6343667, but I haven''t seen the problem I have at the the moment. I started a scrub of a b72 system that doesn''t have any recent snapshots (none since the last scrub) and the % complete is cycling: scrub: scrub in progress, 69.08% done, 0h13m to go scrub: scrub in progress, 46.63% done, 0h28m to go scrub: scrub in progress, 6.36%

FYI for anyone using the dual ethernet ports on a Cisco 7960. I'm using a Cisco 7960 connected to an HP2524 10/100 switch, which has an asterisk box connected directly to it. No VLANs defined or in use. Measured bandwidth: PC -> HP Switch -> Asterisk : actual throughput measured at 94.1 mbps. PC -> 7960 -> HP Switch -> Asterisk : actual measured at 93.02 mbps. The

Hi All, I want to be able to change the labe size on a tree when using ape package so that the graph can look nicer, but don't know how. Playing with the "font" doesn't seem to help. (Right now the labels barely show.) Thanks for your help! Code I'm using now #------------------------- t1 = "(((1:77.37,2:77.37):9.4,4:86.77):7.56,3:94.33);" t2 =