similar to: testing non-linear component in mgcv:gam

Dear list, I am trying to create a three-way table with percent occurrence instead of raw frequencies. However, I cannot get the results I expected: I have the following table: > ftable(table( mannerDF$agem, mannerDF$target, mannerDF$manner )) <snip> 50 bak 0 0 0 0 1 0 pak 0 0 0 0 3 0 sak

Hi, I have a function that approximates some data and indicates "segments". I'd like to plot the original data, and then the linear approximations on top of it. (Ideally, just a subset of N rows at a time, as the data set is large.) I can't figure out a clean way to do this. Suggestions? here is some sample data: ================================== Row X Seg ?.

Dear R experts, I have a Matlab code which I am translating to R in order to examine and enhance it. First of all, I need to reproduce in R the results which were already obtained in Matlab (to make sure that everything is correct). There are some matrix manipulations and '\' operation among them in the code. I have the following data frame > ABS.df Pro syn

I have a data set that I display using barplot. I don't know what you call it, but when I look at it, it looks like this: > lsu (0,0.1] (0.1,0.2] (0.2,0.3] (0.3,0.4] (0.4,0.5] (0.5,0.6] A 0.052631579 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 B 0.000000000 0.000000000 0.001007049 0.003021148 0.000000000 0.000000000 E 0.200000000 0.000000000

Hi r-sers, I have a data of relative frequencies for the interval of 0-20, 20-40,...380-400.  I would like the two data on the same graph using the same x-axis label.  My question is how to get a smooth curve using kernel density code if it possible for this data.   > cbind(rel_obs,rel_gen)           rel_obs rel_gen  [1,] 0.000000000  0.0000  [2,] 0.092534175  0.0712  [3,] 0.105152471  0.1092

I have the following matrix: > dat [,1] [,2] [,3] [,4] foo 0.7574657 0.2104075 0.02922241 0.002705617 foo 0.0000000 0.0000000 0.00000000 0.000000000 foo 0.0000000 0.0000000 0.00000000 0.000000000 foo 0.0000000 0.0000000 0.00000000 0.000000000 foo 0.0000000 0.0000000 0.00000000 0.000000000 foo 0.0000000 0.0000000 0.00000000 0.000000000 and given this:

Hi everybody, I am sorry that I am kind of spamming this forum, but I have searched for some input everywhere and cant really find a nice solution for my problem. Data looks like: price 2011-11-01 08:00:00 0.000000000 2011-11-01 08:00:00 0.000000000 2011-11-01 08:02:00 0.000000000 2011-11-01 08:03:00 -0.017033339 2011-11-01 08:13:00 0.000690001 2011-11-01

Hy Guys, My samba version: 4.0.13 ntp version: 4.2.6p3 Ubuntu 12.04.3 LTS I have a weird problem, I sat up the time service via https://wiki.samba.org/index.php/Configure_NTP On the servers I have no firewall. This port is open when i check with nmap, ntp runs. Looks like all of my Windows 7 clients works fine but w32tm /resync says Permission denied 0x8007005, so not that good :) Windows

Hi everyone, I've got a question concerning the function Box.test for testing autocorrelation in my data. My data consist of (daily) returns of several stocks over time (first row=time, all other rows=stock returns). I intend to perform a Box-Ljung test for my returns (for each stock). Since I have about 3000 stocks in my list, I'm not able to perform the test individually for each

Dear members. I am having problems to understand the kmeans- results in R. I am applying kmeans-algorithms to my big data file, and it is producing the results of the clusters. Q1) Does anybody knows how to find out in which cluster (I have fixed numberofclusters = 5 ) which data have been used? COMMAND (kmeans.results <- kmeans(mydata,centers =5, iter.max= 1000, nstart =10000)) Q2) When I

I have managed to create a data frame like this: > tsus_same_mean[1:10,] PID Grp Dist PercAln PercId 1 12638 Acidobacteria 0.000000000 1.0000000 1.0000000 2 87 Actinobacteria 0.000000000 0.9700000 0.9700000 3 92 Actinobacteria 0.008902000 1.0000000 0.9910000 4 94 Actinobacteria 0.000000000 1.0000000 1.0000000 5 189 Actinobacteria 0.005876733

Hello, I thought that I should try openblas when building a CRAN package containing lots of old (twentieth century) C-code with frequent calls to blas and lapack routines. I have the following options on my Ubuntu 20.04 machine: Selection Path Priority Status ------------------------------------------------------------ * 0

I did this: nb <- naiveBayes(users, platform) pl <- predict(nb,users) nrow(users) ==> 314781 ncol(users) ==> 109 1. naiveBayes() was quite fast (~20 seconds), while predict() was slow (tens of minutes). why? 2. the predict results were completely off the mark (quite the opposite of the expected overfitting). suffice it to show the tables: pl: android blackberry ipad

Hi, I have the binned data (observed and generated from model)  that I would like to test using the anderson-darling goodness of fit test.  But I'm not sure which package in R to use. I tried ad.test(...) but it does not recognise the test by Vito Ricci in FITTING DISTRIBUTIONS WITH R   > ad.test(hist_hume_beec[,1],hist_hume_beec[,2]) Error: could not find function "ad.test"

Hello, What would be the best set of R functions to parse and transform a file? My file looks as shown below. I would like to plot this data and I need to parse it into a single data frame that sorts of "transposes the data" with the following structure: > df <- data.frame(n=c(1,1,2,2),iter=c(1,2,1,2),step=as.factor(c('Step 1', 'Step2', 'Step 1',

Hi all, I'm new to R and have the following problem: I have a 2 factor design (a has 2 levels, b has 3 levels). I have an object kidney.aov which is an aov(y ~ a*b), and when I ask for model.tables(kidney.avo, se=T) I get the following message along with the table of effects: Design is unbalanced - use se.contrast() for se's but the design is NOT unbalanced... each fator level

Full_Name: Patrick Hodgson Version: 2.0 OS: solaris 2.9 Submission from: (NULL) (65.94.128.161) The value reported for the parameter W in the function wilcox.exact appears to be incorrect. I have checked the reference in the help file for this function (Myles & Hollander 1973, as well as 2nd ed. 1999 by same authors) and it is clear that W is the sum of the ranks of the data set with the

Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with thanks in advance. Group SI Sif SG Io I2 lol2 1 9.08e-05 9.08e+00 0.060842287 1.798556446 32.574500

I have a data frame (morgan) of hourly river flow, river levels and wind direction and speed thus: Time hour lev.morgan lev.lock2 lev.lock1 flow direction velocity 1 2009-07-06 15:00:00 15 3.266 3.274 3.240 1710.6 180.282 4.352 2 2009-07-06 16:00:00 16 3.268 3.272 3.240 1441.8 192.338 5.496 3 2009-07-06 17:00:00 17 3.268

I have a matrix which I wanted to convert to a data frame. As I could not succeed and resorted to export to csv and reimport it again. Why did I fail in the attempt and how can I achieve what I wanted without this roundabouts? The original matrix: > str(comb_model0) num [1:90, 1:4] 3.5938 0.0274 0.0342 0.0135 0.0207 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:90]